National Library of Energy BETA

Sample records for adaptive mesh refinement

  1. CLAMR (Compute Language Adaptive Mesh Refinement)

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    CLAMR (Compute Language Adaptive Mesh Refinement) CLAMR (Compute Language Adaptive Mesh Refinement) CLAMR (Compute Language Adaptive Mesh Refinement) is being developed as a DOE...

  2. Arbitrary Lagrangian Eulerian Adaptive Mesh Refinement

    Energy Science and Technology Software Center (OSTI)

    2009-09-29

    This is a simulation code involving an ALE (arbitrary Lagrangian-Eulerian) hydrocode with AMR (adaptive mesh refinement) and pluggable physics packages for material strength, heat conduction, radiation diffusion, and laser ray tracing developed a LLNL, UCSD, and Berkeley Lab. The code is an extension of the open source SAMRAI (Structured Adaptive Mesh Refinement Application Interface) code/library. The code can be used in laser facilities such as the National Ignition Facility. The code is alsi being appliedmore » to slurry flow (landslides).« less

  3. COSMOLOGICAL ADAPTIVE MESH REFINEMENT MAGNETOHYDRODYNAMICS WITH ENZO

    SciTech Connect (OSTI)

    Collins, David C.; Xu Hao; Norman, Michael L.; Li Hui; Li Shengtai

    2010-02-01

    In this work, we present EnzoMHD, the extension of the cosmological code Enzo to include the effects of magnetic fields through the ideal magnetohydrodynamics approximation. We use a higher order Godunov method for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A second-order divergence-free reconstruction technique is used to interpolate the magnetic fields in the block-structured adaptive mesh refinement framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non-cosmological test problems to demonstrate the quality of solution resulting from this combination of solvers.

  4. Visualization Tools for Adaptive Mesh Refinement Data

    SciTech Connect (OSTI)

    Weber, Gunther H.; Beckner, Vincent E.; Childs, Hank; Ligocki,Terry J.; Miller, Mark C.; Van Straalen, Brian; Bethel, E. Wes

    2007-05-09

    Adaptive Mesh Refinement (AMR) is a highly effective method for simulations that span a large range of spatiotemporal scales, such as astrophysical simulations that must accommodate ranges from interstellar to sub-planetary. Most mainstream visualization tools still lack support for AMR as a first class data type and AMR code teams use custom built applications for AMR visualization. The Department of Energy's (DOE's) Science Discovery through Advanced Computing (SciDAC) Visualization and Analytics Center for Enabling Technologies (VACET) is currently working on extending VisIt, which is an open source visualization tool that accommodates AMR as a first-class data type. These efforts will bridge the gap between general-purpose visualization applications and highly specialized AMR visual analysis applications. Here, we give an overview of the state of the art in AMR visualization research and tools and describe how VisIt currently handles AMR data.

  5. Visualization of Scalar Adaptive Mesh Refinement Data

    SciTech Connect (OSTI)

    VACET; Weber, Gunther; Weber, Gunther H.; Beckner, Vince E.; Childs, Hank; Ligocki, Terry J.; Miller, Mark C.; Van Straalen, Brian; Bethel, E. Wes

    2007-12-06

    Adaptive Mesh Refinement (AMR) is a highly effective computation method for simulations that span a large range of spatiotemporal scales, such as astrophysical simulations, which must accommodate ranges from interstellar to sub-planetary. Most mainstream visualization tools still lack support for AMR grids as a first class data type and AMR code teams use custom built applications for AMR visualization. The Department of Energy's (DOE's) Science Discovery through Advanced Computing (SciDAC) Visualization and Analytics Center for Enabling Technologies (VACET) is currently working on extending VisIt, which is an open source visualization tool that accommodates AMR as a first-class data type. These efforts will bridge the gap between general-purpose visualization applications and highly specialized AMR visual analysis applications. Here, we give an overview of the state of the art in AMR scalar data visualization research.

  6. Structured Adaptive Mesh Refinement Application Infrastructure

    Energy Science and Technology Software Center (OSTI)

    2010-07-15

    SAMRAI is an object-oriented support library for structured adaptice mesh refinement (SAMR) simulation of computational science problems, modeled by systems of partial differential equations (PDEs). SAMRAI is developed and maintained in the Center for Applied Scientific Computing (CASC) under ASCI ITS and PSE support. SAMRAI is used in a variety of application research efforts at LLNL and in academia. These applications are developed in collaboration with SAMRAI development team members.

  7. Parallel Block Structured Adaptive Mesh Refinement on Graphics Processing Units

    SciTech Connect (OSTI)

    Beckingsale, D. A.; Gaudin, W. P.; Hornung, R. D.; Gunney, B. T.; Gamblin, T.; Herdman, J. A.; Jarvis, S. A.

    2014-11-17

    Block-structured adaptive mesh refinement is a technique that can be used when solving partial differential equations to reduce the number of zones necessary to achieve the required accuracy in areas of interest. These areas (shock fronts, material interfaces, etc.) are recursively covered with finer mesh patches that are grouped into a hierarchy of refinement levels. Despite the potential for large savings in computational requirements and memory usage without a corresponding reduction in accuracy, AMR adds overhead in managing the mesh hierarchy, adding complex communication and data movement requirements to a simulation. In this paper, we describe the design and implementation of a native GPU-based AMR library, including: the classes used to manage data on a mesh patch, the routines used for transferring data between GPUs on different nodes, and the data-parallel operators developed to coarsen and refine mesh data. We validate the performance and accuracy of our implementation using three test problems and two architectures: an eight-node cluster, and over four thousand nodes of Oak Ridge National Laboratory’s Titan supercomputer. Our GPU-based AMR hydrodynamics code performs up to 4.87× faster than the CPU-based implementation, and has been scaled to over four thousand GPUs using a combination of MPI and CUDA.

  8. Toward parallel, adaptive mesh refinement for chemically reacting flow simulations

    SciTech Connect (OSTI)

    Devine, K.D.; Shadid, J.N.; Salinger, A.G. Hutchinson, S.A.; Hennigan, G.L.

    1997-12-01

    Adaptive numerical methods offer greater efficiency than traditional numerical methods by concentrating computational effort in regions of the problem domain where the solution is difficult to obtain. In this paper, the authors describe progress toward adding mesh refinement to MPSalsa, a computer program developed at Sandia National laboratories to solve coupled three-dimensional fluid flow and detailed reaction chemistry systems for modeling chemically reacting flow on large-scale parallel computers. Data structures that support refinement and dynamic load-balancing are discussed. Results using uniform refinement with mesh sequencing to improve convergence to steady-state solutions are also presented. Three examples are presented: a lid driven cavity, a thermal convection flow, and a tilted chemical vapor deposition reactor.

  9. Advances in Patch-Based Adaptive Mesh Refinement Scalability

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Gunney, Brian T.N.; Anderson, Robert W.

    2015-12-18

    Patch-based structured adaptive mesh refinement (SAMR) is widely used for high-resolution simu- lations. Combined with modern supercomputers, it could provide simulations of unprecedented size and resolution. A persistent challenge for this com- bination has been managing dynamically adaptive meshes on more and more MPI tasks. The dis- tributed mesh management scheme in SAMRAI has made some progress SAMR scalability, but early al- gorithms still had trouble scaling past the regime of 105 MPI tasks. This work provides two critical SAMR regridding algorithms, which are integrated into that scheme to ensure efficiency of the whole. The clustering algorithm is an extensionmore » of the tile- clustering approach, making it more flexible and efficient in both clustering and parallelism. The partitioner is a new algorithm designed to prevent the network congestion experienced by its prede- cessor. We evaluated performance using weak- and strong-scaling benchmarks designed to be difficult for dynamic adaptivity. Results show good scaling on up to 1.5M cores and 2M MPI tasks. Detailed timing diagnostics suggest scaling would continue well past that.« less

  10. Advances in Patch-Based Adaptive Mesh Refinement Scalability

    SciTech Connect (OSTI)

    Gunney, Brian T.N.; Anderson, Robert W.

    2015-12-18

    Patch-based structured adaptive mesh refinement (SAMR) is widely used for high-resolution simu- lations. Combined with modern supercomputers, it could provide simulations of unprecedented size and resolution. A persistent challenge for this com- bination has been managing dynamically adaptive meshes on more and more MPI tasks. The dis- tributed mesh management scheme in SAMRAI has made some progress SAMR scalability, but early al- gorithms still had trouble scaling past the regime of 105 MPI tasks. This work provides two critical SAMR regridding algorithms, which are integrated into that scheme to ensure efficiency of the whole. The clustering algorithm is an extension of the tile- clustering approach, making it more flexible and efficient in both clustering and parallelism. The partitioner is a new algorithm designed to prevent the network congestion experienced by its prede- cessor. We evaluated performance using weak- and strong-scaling benchmarks designed to be difficult for dynamic adaptivity. Results show good scaling on up to 1.5M cores and 2M MPI tasks. Detailed timing diagnostics suggest scaling would continue well past that.

  11. ENZO: AN ADAPTIVE MESH REFINEMENT CODE FOR ASTROPHYSICS

    SciTech Connect (OSTI)

    Bryan, Greg L.; Turk, Matthew J.; Norman, Michael L.; Bordner, James; Xu, Hao; Kritsuk, Alexei G.; O'Shea, Brian W.; Smith, Britton; Abel, Tom; Wang, Peng; Skillman, Samuel W.; Wise, John H.; Reynolds, Daniel R.; Collins, David C.; Harkness, Robert P.; Kim, Ji-hoon; Kuhlen, Michael; Goldbaum, Nathan; Hummels, Cameron; Collaboration: Enzo Collaboration; and others

    2014-04-01

    This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in one, two, and three dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

  12. An adaptive mesh-moving and refinement procedure for one-dimensional conservation laws

    SciTech Connect (OSTI)

    Biswas, R.; Flaherty, J.E.; Arney, D.C. Rensselaer Polytech. Inst., Troy, NY Military Acad., West Point, NY )

    1993-01-01

    We examine the performance of an adaptive mesh-moving and /or local mesh refinement procedure for the finite difference solution of one-dimensional hyperbolic systems of conservation laws. Adaptive motion of a base mesh is designed to isolate spatially distinct phenomena, and recursive local refinement of the time step and cells of the stationary or moving base mesh is performed in regions where a refinement indicator exceeds a prescribed tolerance. These adaptive procedures are incorporated into a computer code that includes a MacCormack finite difference scheme wih Davis' artificial viscosity model and a discretization error estimate based on Richardson's extrapolation. Experiments are conducted on three problems in order to qualify the advantages of adaptive techniques relative to uniform mesh computations and the relative benefits of mesh moving and refinement. Key results indicate that local mesh refinement, with and without mesh moving, can provide reliable solutions at much lower computational cost than possible on uniform meshes; that mesh motion can be used to improve the results of uniform mesh solutions for a modest computational effort; that the cost of managing the tree data structure associated with refinement is small; and that a combination of mesh motion and refinement reliably produces solutions for the least cost per unit accuracy.

  13. RAM: a Relativistic Adaptive Mesh Refinement Hydrodynamics Code

    SciTech Connect (OSTI)

    Zhang, Wei-Qun; MacFadyen, Andrew I.; /Princeton, Inst. Advanced Study

    2005-06-06

    The authors have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. They have implemented a characteristic-wise, finite difference, weighted essentially non-oscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration they use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. They have also implemented fourth and fifth order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods and physics modules. In addition to WENO they have implemented a finite volume module with the piecewise parabolic method (PPM) for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. They examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. They show that under-resolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two and three dimensions and in Cartesian, cylindrical and spherical coordinates. they have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which they show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.

  14. A Predictive Model of Fragmentation using Adaptive Mesh Refinement and a Hierarchical Material Model

    SciTech Connect (OSTI)

    Koniges, A E; Masters, N D; Fisher, A C; Anderson, R W; Eder, D C; Benson, D; Kaiser, T B; Gunney, B T; Wang, P; Maddox, B R; Hansen, J F; Kalantar, D H; Dixit, P; Jarmakani, H; Meyers, M A

    2009-03-03

    Fragmentation is a fundamental material process that naturally spans spatial scales from microscopic to macroscopic. We developed a mathematical framework using an innovative combination of hierarchical material modeling (HMM) and adaptive mesh refinement (AMR) to connect the continuum to microstructural regimes. This framework has been implemented in a new multi-physics, multi-scale, 3D simulation code, NIF ALE-AMR. New multi-material volume fraction and interface reconstruction algorithms were developed for this new code, which is leading the world effort in hydrodynamic simulations that combine AMR with ALE (Arbitrary Lagrangian-Eulerian) techniques. The interface reconstruction algorithm is also used to produce fragments following material failure. In general, the material strength and failure models have history vector components that must be advected along with other properties of the mesh during remap stage of the ALE hydrodynamics. The fragmentation models are validated against an electromagnetically driven expanding ring experiment and dedicated laser-based fragmentation experiments conducted at the Jupiter Laser Facility. As part of the exit plan, the NIF ALE-AMR code was applied to a number of fragmentation problems of interest to the National Ignition Facility (NIF). One example shows the added benefit of multi-material ALE-AMR that relaxes the requirement that material boundaries must be along mesh boundaries.

  15. 3D Adaptive Mesh Refinement Simulations of Pellet Injection in Tokamaks

    SciTech Connect (OSTI)

    R. Samtaney; S.C. Jardin; P. Colella; D.F. Martin

    2003-10-20

    We present results of Adaptive Mesh Refinement (AMR) simulations of the pellet injection process, a proven method of refueling tokamaks. AMR is a computationally efficient way to provide the resolution required to simulate realistic pellet sizes relative to device dimensions. The mathematical model comprises of single-fluid MHD equations with source terms in the continuity equation along with a pellet ablation rate model. The numerical method developed is an explicit unsplit upwinding treatment of the 8-wave formulation, coupled with a MAC projection method to enforce the solenoidal property of the magnetic field. The Chombo framework is used for AMR. The role of the E x B drift in mass redistribution during inside and outside pellet injections is emphasized.

  16. On the Computation of Integral Curves in Adaptive Mesh Refinement Vector Fields

    SciTech Connect (OSTI)

    Deines, Eduard; Weber, Gunther H.; Garth, Christoph; Van Straalen, Brian; Borovikov, Sergey; Martin, Daniel F.; Joy, Kenneth I.

    2011-06-27

    Integral curves, such as streamlines, streaklines, pathlines, and timelines, are an essential tool in the analysis of vector field structures, offering straightforward and intuitive interpretation of visualization results. While such curves have a long-standing tradition in vector field visualization, their application to Adaptive Mesh Refinement (AMR) simulation results poses unique problems. AMR is a highly effective discretization method for a variety of physical simulation problems and has recently been applied to the study of vector fields in flow and magnetohydrodynamic applications. The cell-centered nature of AMR data and discontinuities in the vector field representation arising from AMR level boundaries complicate the application of numerical integration methods to compute integral curves. In this paper, we propose a novel approach to alleviate these problems and show its application to streamline visualization in an AMR model of the magnetic field of the solar system as well as to a simulation of two incompressible viscous vortex rings merging.

  17. Advanced numerical methods in mesh generation and mesh adaptation

    SciTech Connect (OSTI)

    Lipnikov, Konstantine; Danilov, A; Vassilevski, Y; Agonzal, A

    2010-01-01

    Numerical solution of partial differential equations requires appropriate meshes, efficient solvers and robust and reliable error estimates. Generation of high-quality meshes for complex engineering models is a non-trivial task. This task is made more difficult when the mesh has to be adapted to a problem solution. This article is focused on a synergistic approach to the mesh generation and mesh adaptation, where best properties of various mesh generation methods are combined to build efficiently simplicial meshes. First, the advancing front technique (AFT) is combined with the incremental Delaunay triangulation (DT) to build an initial mesh. Second, the metric-based mesh adaptation (MBA) method is employed to improve quality of the generated mesh and/or to adapt it to a problem solution. We demonstrate with numerical experiments that combination of all three methods is required for robust meshing of complex engineering models. The key to successful mesh generation is the high-quality of the triangles in the initial front. We use a black-box technique to improve surface meshes exported from an unattainable CAD system. The initial surface mesh is refined into a shape-regular triangulation which approximates the boundary with the same accuracy as the CAD mesh. The DT method adds robustness to the AFT. The resulting mesh is topologically correct but may contain a few slivers. The MBA uses seven local operations to modify the mesh topology. It improves significantly the mesh quality. The MBA method is also used to adapt the mesh to a problem solution to minimize computational resources required for solving the problem. The MBA has a solid theoretical background. In the first two experiments, we consider the convection-diffusion and elasticity problems. We demonstrate the optimal reduction rate of the discretization error on a sequence of adaptive strongly anisotropic meshes. The key element of the MBA method is construction of a tensor metric from hierarchical edge

  18. A Freestream-Preserving High-Order Finite-Volume Method for Mapped Grids with Adaptive-Mesh Refinement

    SciTech Connect (OSTI)

    Guzik, S; McCorquodale, P; Colella, P

    2011-12-16

    A fourth-order accurate finite-volume method is presented for solving time-dependent hyperbolic systems of conservation laws on mapped grids that are adaptively refined in space and time. Novel considerations for formulating the semi-discrete system of equations in computational space combined with detailed mechanisms for accommodating the adapting grids ensure that conservation is maintained and that the divergence of a constant vector field is always zero (freestream-preservation property). Advancement in time is achieved with a fourth-order Runge-Kutta method.

  19. Refining quadrilateral and brick element meshes

    SciTech Connect (OSTI)

    Schneiders, R.; Debye, J.

    1995-12-31

    We consider the problem of refining unstructured quadrilateral and brick element meshes. We present an algorithm which is a generalization of an algorithm developed by Cheng et. al. for structured quadrilateral element meshes. The problem is solved for the two-dimensional case. Concerning three dimensions we present a solution for some special cases and a general solution that introduces tetrahedral and pyramidal transition elements.

  20. Parallel tetrahedral mesh refinement with MOAB.

    SciTech Connect (OSTI)

    Thompson, David C.; Pebay, Philippe Pierre

    2008-12-01

    In this report, we present the novel functionality of parallel tetrahedral mesh refinement which we have implemented in MOAB. This report details work done to implement parallel, edge-based, tetrahedral refinement into MOAB. The theoretical basis for this work is contained in [PT04, PT05, TP06] while information on design, performance, and operation specific to MOAB are contained herein. As MOAB is intended mainly for use in pre-processing and simulation (as opposed to the post-processing bent of previous papers), the primary use case is different: rather than refining elements with non-linear basis functions, the goal is to increase the number of degrees of freedom in some region in order to more accurately represent the solution to some system of equations that cannot be solved analytically. Also, MOAB has a unique mesh representation which impacts the algorithm. This introduction contains a brief review of streaming edge-based tetrahedral refinement. The remainder of the report is broken into three sections: design and implementation, performance, and conclusions. Appendix A contains instructions for end users (simulation authors) on how to employ the refiner.

  1. Cubit Adaptive Meshing Algorithm Library

    Energy Science and Technology Software Center (OSTI)

    2004-09-01

    CAMAL (Cubit adaptive meshing algorithm library) is a software component library for mesh generation. CAMAL 2.0 includes components for triangle, quad and tetrahedral meshing. A simple Application Programmers Interface (API) takes a discrete boundary definition and CAMAL computes a quality interior unstructured grid. The triangle and quad algorithms may also import a geometric definition of a surface on which to define the grid. CAMAL’s triangle meshing uses a 3D space advancing front method, the quadmore » meshing algorithm is based upon Sandia’s patented paving algorithm and the tetrahedral meshing algorithm employs the GHS3D-Tetmesh component developed by INRIA, France.« less

  2. Mesh refinement for uncertainty quantification through model reduction

    SciTech Connect (OSTI)

    Li, Jing Stinis, Panos

    2015-01-01

    We present a novel way of deciding when and where to refine a mesh in probability space in order to facilitate uncertainty quantification in the presence of discontinuities in random space. A discontinuity in random space makes the application of generalized polynomial chaos expansion techniques prohibitively expensive. The reason is that for discontinuous problems, the expansion converges very slowly. An alternative to using higher terms in the expansion is to divide the random space in smaller elements where a lower degree polynomial is adequate to describe the randomness. In general, the partition of the random space is a dynamic process since some areas of the random space, particularly around the discontinuity, need more refinement than others as time evolves. In the current work we propose a way to decide when and where to refine the random space mesh based on the use of a reduced model. The idea is that a good reduced model can monitor accurately, within a random space element, the cascade of activity to higher degree terms in the chaos expansion. In turn, this facilitates the efficient allocation of computational sources to the areas of random space where they are more needed. For the Kraichnan–Orszag system, the prototypical system to study discontinuities in random space, we present theoretical results which show why the proposed method is sound and numerical results which corroborate the theory.

  3. Adaptive h -refinement for reduced-order models: ADAPTIVE h -refinement for reduced-order models

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Carlberg, Kevin T.

    2014-11-05

    Our work presents a method to adaptively refine reduced-order models a posteriori without requiring additional full-order-model solves. The technique is analogous to mesh-adaptive h-refinement: it enriches the reduced-basis space online by ‘splitting’ a given basis vector into several vectors with disjoint support. The splitting scheme is defined by a tree structure constructed offline via recursive k-means clustering of the state variables using snapshot data. This method identifies the vectors to split online using a dual-weighted-residual approach that aims to reduce error in an output quantity of interest. The resulting method generates a hierarchy of subspaces online without requiring large-scale operationsmore » or full-order-model solves. Furthermore, it enables the reduced-order model to satisfy any prescribed error tolerance regardless of its original fidelity, as a completely refined reduced-order model is mathematically equivalent to the original full-order model. Experiments on a parameterized inviscid Burgers equation highlight the ability of the method to capture phenomena (e.g., moving shocks) not contained in the span of the original reduced basis.« less

  4. Applications of automatic mesh generation and adaptive methods in computational medicine

    SciTech Connect (OSTI)

    Schmidt, J.A.; Macleod, R.S.; Johnson, C.R.; Eason, J.C.

    1995-12-31

    Important problems in Computational Medicine exist that can benefit from the implementation of adaptive mesh refinement techniques. Biological systems are so inherently complex that only efficient models running on state of the art hardware can begin to simulate reality. To tackle the complex geometries associated with medical applications we present a general purpose mesh generation scheme based upon the Delaunay tessellation algorithm and an iterative point generator. In addition, automatic, two- and three-dimensional adaptive mesh refinement methods are presented that are derived from local and global estimates of the finite element error. Mesh generation and adaptive refinement techniques are utilized to obtain accurate approximations of bioelectric fields within anatomically correct models of the heart and human thorax. Specifically, we explore the simulation of cardiac defibrillation and the general forward and inverse problems in electrocardiography (ECG). Comparisons between uniform and adaptive refinement techniques are made to highlight the computational efficiency and accuracy of adaptive methods in the solution of field problems in computational medicine.

  5. Parallel paving: An algorithm for generating distributed, adaptive, all-quadrilateral meshes on parallel computers

    SciTech Connect (OSTI)

    Lober, R.R.; Tautges, T.J.; Vaughan, C.T.

    1997-03-01

    Paving is an automated mesh generation algorithm which produces all-quadrilateral elements. It can additionally generate these elements in varying sizes such that the resulting mesh adapts to a function distribution, such as an error function. While powerful, conventional paving is a very serial algorithm in its operation. Parallel paving is the extension of serial paving into parallel environments to perform the same meshing functions as conventional paving only on distributed, discretized models. This extension allows large, adaptive, parallel finite element simulations to take advantage of paving`s meshing capabilities for h-remap remeshing. A significantly modified version of the CUBIT mesh generation code has been developed to host the parallel paving algorithm and demonstrate its capabilities on both two dimensional and three dimensional surface geometries and compare the resulting parallel produced meshes to conventionally paved meshes for mesh quality and algorithm performance. Sandia`s {open_quotes}tiling{close_quotes} dynamic load balancing code has also been extended to work with the paving algorithm to retain parallel efficiency as subdomains undergo iterative mesh refinement.

  6. Modeling, mesh generation, and adaptive numerical methods for partial differential equations

    SciTech Connect (OSTI)

    Babuska, I.; Henshaw, W.D.; Oliger, J.E.; Flaherty, J.E.; Hopcroft, J.E.; Tezduyar, T.

    1995-12-31

    Mesh generation is one of the most time consuming aspects of computational solutions of problems involving partial differential equations. It is, furthermore, no longer acceptable to compute solutions without proper verification that specified accuracy criteria are being satisfied. Mesh generation must be related to the solution through computable estimates of discretization errors. Thus, an iterative process of alternate mesh and solution generation evolves in an adaptive manner with the end result that the solution is computed to prescribed specifications in an optimal, or at least efficient, manner. While mesh generation and adaptive strategies are becoming available, major computational challenges remain. One, in particular, involves moving boundaries and interfaces, such as free-surface flows and fluid-structure interactions. A 3-week program was held from July 5 to July 23, 1993 with 173 participants and 66 keynote, invited, and contributed presentations. This volume represents written versions of 21 of these lectures. These proceedings are organized roughly in order of their presentation at the workshop. Thus, the initial papers are concerned with geometry and mesh generation and discuss the representation of physical objects and surfaces on a computer and techniques to use this data to generate, principally, unstructured meshes of tetrahedral or hexahedral elements. The remainder of the papers cover adaptive strategies, error estimation, and applications. Several submissions deal with high-order p- and hp-refinement methods where mesh refinement/coarsening (h-refinement) is combined with local variation of method order (p-refinement). Combinations of mathematically verified and physically motivated approaches to error estimation are represented. Applications center on fluid mechanics. Selected papers are indexed separately for inclusion in the Energy Science and Technology Database.

  7. Automatic Mesh Adaptivity for Hybrid Monte Carlo/Deterministic Neutronics Modeling of Fusion Energy Systems

    SciTech Connect (OSTI)

    Ibrahim, Ahmad M; Wilson, P.; Sawan, M.; Mosher, Scott W; Peplow, Douglas E.; Grove, Robert E

    2013-01-01

    Three mesh adaptivity algorithms were developed to facilitate and expedite the use of the CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques in accurate full-scale neutronics simulations of fusion energy systems with immense sizes and complicated geometries. First, a macromaterial approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as much geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm decouples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility and resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation. Additionally, because of the significant increase in the efficiency of FW-CADIS simulations, the three algorithms enabled this difficult calculation to be accurately solved on a regular computer cluster, eliminating the need for a world-class super computer.

  8. High-Resolution Numerical Simulation and Analysis of Mach Reflection Structures in Detonation Waves in Low-Pressure H 2 –O 2 –Ar Mixtures: A Summary of Results Obtained with the Adaptive Mesh Refinement Framework AMROC

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Deiterding, Ralf

    2011-01-01

    Numerical simulation can be key to the understanding of the multidimensional nature of transient detonation waves. However, the accurate approximation of realistic detonations is demanding as a wide range of scales needs to be resolved. This paper describes a successful solution strategy that utilizes logically rectangular dynamically adaptive meshes. The hydrodynamic transport scheme and the treatment of the nonequilibrium reaction terms are sketched. A ghost fluid approach is integrated into the method to allow for embedded geometrically complex boundaries. Large-scale parallel simulations of unstable detonation structures of Chapman-Jouguet detonations in low-pressure hydrogen-oxygen-argon mixtures demonstrate the efficiency of the described techniquesmore » in practice. In particular, computations of regular cellular structures in two and three space dimensions and their development under transient conditions, that is, under diffraction and for propagation through bends are presented. Some of the observed patterns are classified by shock polar analysis, and a diagram of the transition boundaries between possible Mach reflection structures is constructed.« less

  9. High-Resolution Numerical Simulation and Analysis of Mach Reflection Structures in Detonation Waves in Low-Pressure H2O2Ar Mixtures: A Summary of Results Obtained with the Adaptive Mesh Refinement Framework AMROC

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Deiterding, Ralf

    2011-01-01

    Numerical simulation can be key to the understanding of the multidimensional nature of transient detonation waves. However, the accurate approximation of realistic detonations is demanding as a wide range of scales needs to be resolved. This paper describes a successful solution strategy that utilizes logically rectangular dynamically adaptive meshes. The hydrodynamic transport scheme and the treatment of the nonequilibrium reaction terms are sketched. A ghost fluid approach is integrated into the method to allow for embedded geometrically complex boundaries. Large-scale parallel simulations of unstable detonation structures of Chapman-Jouguet detonations in low-pressure hydrogen-oxygen-argon mixtures demonstrate the efficiency of the described techniquesmorein practice. In particular, computations of regular cellular structures in two and three space dimensions and their development under transient conditions, that is, under diffraction and for propagation through bends are presented. Some of the observed patterns are classified by shock polar analysis, and a diagram of the transition boundaries between possible Mach reflection structures is constructed.less

  10. Adaptive upscaling with the dual mesh method

    SciTech Connect (OSTI)

    Guerillot, D.; Verdiere, S.

    1997-08-01

    The objective of this paper is to demonstrate that upscaling should be calculated during the flow simulation instead of trying to enhance the a priori upscaling methods. Hence, counter-examples are given to motivate our approach, the so-called Dual Mesh Method. The main steps of this numerical algorithm are recalled. Applications illustrate the necessity to consider different average relative permeability values depending on the direction in space. Moreover, these values could be different for the same average saturation. This proves that an a priori upscaling cannot be the answer even in homogeneous cases because of the {open_quotes}dynamical heterogeneity{close_quotes} created by the saturation profile. Other examples show the efficiency of the Dual Mesh Method applied to heterogeneous medium and to an actual field case in South America.

  11. Dimensional reduction as a tool for mesh refinement and trackingsingularities of PDEs

    SciTech Connect (OSTI)

    Stinis, Panagiotis

    2007-06-10

    We present a collection of algorithms which utilizedimensional reduction to perform mesh refinement and study possiblysingular solutions of time-dependent partial differential equations. Thealgorithms are inspired by constructions used in statistical mechanics toevaluate the properties of a system near a critical point. The firstalgorithm allows the accurate determination of the time of occurrence ofa possible singularity. The second algorithm is an adaptive meshrefinement scheme which can be used to approach efficiently the possiblesingularity. Finally, the third algorithm uses the second algorithm untilthe available resolution is exhausted (as we approach the possiblesingularity) and then switches to a dimensionally reduced model which,when accurate, can follow faithfully the solution beyond the time ofoccurrence of the purported singularity. An accurate dimensionallyreduced model should dissipate energy at the right rate. We construct twovariants of each algorithm. The first variant assumes that we have actualknowledge of the reduced model. The second variant assumes that we knowthe form of the reduced model, i.e., the terms appearing in the reducedmodel, but not necessarily their coefficients. In this case, we alsoprovide a way of determining the coefficients. We present numericalresults for the Burgers equation with zero and nonzero viscosity toillustrate the use of the algorithms.

  12. Automatic mesh adaptivity for CADIS and FW-CADIS neutronics modeling of difficult shielding problems

    SciTech Connect (OSTI)

    Ibrahim, A. M.; Peplow, D. E.; Mosher, S. W.; Wagner, J. C.; Evans, T. M.; Wilson, P. P.; Sawan, M. E.

    2013-07-01

    The CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques dramatically increase the efficiency of neutronics modeling, but their use in the accurate design analysis of very large and geometrically complex nuclear systems has been limited by the large number of processors and memory requirements for their preliminary deterministic calculations and final Monte Carlo calculation. Three mesh adaptivity algorithms were developed to reduce the memory requirements of CADIS and FW-CADIS without sacrificing their efficiency improvement. First, a macro-material approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as much geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm de-couples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility. Using these algorithms resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation and, additionally, increased the efficiency of the Monte Carlo simulation by a factor of at least 3.4. The three algorithms enabled this difficult calculation to be accurately solved using an FW-CADIS simulation on a regular computer cluster, obviating the need for a world-class super computer. (authors)

  13. Automatic mesh adaptivity for hybrid Monte Carlo/deterministic neutronics modeling of difficult shielding problems

    SciTech Connect (OSTI)

    Ibrahim, Ahmad M.; Wilson, Paul P.H.; Sawan, Mohamed E.; Mosher, Scott W.; Peplow, Douglas E.; Wagner, John C.; Evans, Thomas M.; Grove, Robert E.

    2015-06-30

    The CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques dramatically increase the efficiency of neutronics modeling, but their use in the accurate design analysis of very large and geometrically complex nuclear systems has been limited by the large number of processors and memory requirements for their preliminary deterministic calculations and final Monte Carlo calculation. Three mesh adaptivity algorithms were developed to reduce the memory requirements of CADIS and FW-CADIS without sacrificing their efficiency improvement. First, a macromaterial approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as much geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm decouples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility. Using these algorithms resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation and, additionally, increased the efficiency of the Monte Carlo simulation by a factor of at least 3.4. The three algorithms enabled this difficult calculation to be accurately solved using an FW-CADIS simulation on a regular computer cluster, eliminating the need for a world-class super computer.

  14. Automatic mesh adaptivity for hybrid Monte Carlo/deterministic neutronics modeling of difficult shielding problems

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Ibrahim, Ahmad M.; Wilson, Paul P.H.; Sawan, Mohamed E.; Mosher, Scott W.; Peplow, Douglas E.; Wagner, John C.; Evans, Thomas M.; Grove, Robert E.

    2015-06-30

    The CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques dramatically increase the efficiency of neutronics modeling, but their use in the accurate design analysis of very large and geometrically complex nuclear systems has been limited by the large number of processors and memory requirements for their preliminary deterministic calculations and final Monte Carlo calculation. Three mesh adaptivity algorithms were developed to reduce the memory requirements of CADIS and FW-CADIS without sacrificing their efficiency improvement. First, a macromaterial approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as muchmore » geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm decouples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility. Using these algorithms resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation and, additionally, increased the efficiency of the Monte Carlo simulation by a factor of at least 3.4. The three algorithms enabled this difficult calculation to be accurately solved using an FW-CADIS simulation on a regular computer cluster, eliminating the need for a world-class super computer.« less

  15. Local timespace mesh refinement for simulation of elastic wave propagation in multi-scale media

    SciTech Connect (OSTI)

    Kostin, Victor; Lisitsa, Vadim; Reshetova, Galina; Tcheverda, Vladimir

    2015-01-15

    This paper presents an original approach to local timespace grid refinement for the numerical simulation of wave propagation in models with localized clusters of micro-heterogeneities. The main features of the algorithm are the application of temporal and spatial refinement on two different surfaces; the use of the embedded-stencil technique for the refinement of grid step with respect to time; the use of the Fast Fourier Transform (FFT)-based interpolation to couple variables for spatial mesh refinement. The latter makes it possible to perform filtration of high spatial frequencies, which provides stability in the proposed finite-difference schemes. In the present work, the technique is implemented for the finite-difference simulation of seismic wave propagation and the interaction of such waves with fluid-filled fractures and cavities of carbonate reservoirs. However, this approach is easy to adapt and/or combine with other numerical techniques, such as finite elements, discontinuous Galerkin method, or finite volumes used for approximation of various types of linear and nonlinear hyperbolic equations.

  16. Simulating the quartic Galileon gravity model on adaptively refined meshes

    SciTech Connect (OSTI)

    Li, Baojiu; Barreira, Alexandre; Baugh, Carlton M.; Hellwing, Wojciech A.; Koyama, Kazuya; Zhao, Gong-Bo; Pascoli, Silvia E-mail: baojiu.li@durham.ac.uk E-mail: wojciech.hellwing@durham.ac.uk E-mail: silvia.pascoli@durham.ac.uk

    2013-11-01

    We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the ramses N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

  17. Enzo: An Adaptive Mesh Refinement Code for Astrophysics (Journal...

    Office of Scientific and Technical Information (OSTI)

    J. ; Reynolds, Daniel R. ; Collins, David C. ; Wang, Peng ; Skillman, Samuel W. ; Smith, Britton ; Harkness, Robert P. ; Bordner, James ; Kim, Ji-hoon ; Kuhlen, Michael ; Xu, ...

  18. ENZO: AN ADAPTIVE MESH REFINEMENT CODE FOR ASTROPHYSICS (Journal...

    Office of Scientific and Technical Information (OSTI)

    Authors: Bryan, Greg L. ; Turk, Matthew J. 1 ; Norman, Michael L. ; Bordner, James ; Xu, Hao ; Kritsuk, Alexei G. 2 ; O'Shea, Brian W. ; Smith, Britton 3 ; Abel, Tom ; Wang, ...

  19. A Predictive Model of Fragmentation using Adaptive Mesh Refinement...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    A new arena for predictive simulations in high- powered laser systems "12 NIF" Or 96 ... Numerical simulations show how to mitigate damage from debris and shrapnel Damage prior to ...

  20. Adaptive h -refinement for reduced-order models: ADAPTIVE h -refinement for reduced-order models

    SciTech Connect (OSTI)

    Carlberg, Kevin T.

    2014-11-05

    Our work presents a method to adaptively refine reduced-order models a posteriori without requiring additional full-order-model solves. The technique is analogous to mesh-adaptive h-refinement: it enriches the reduced-basis space online by ‘splitting’ a given basis vector into several vectors with disjoint support. The splitting scheme is defined by a tree structure constructed offline via recursive k-means clustering of the state variables using snapshot data. This method identifies the vectors to split online using a dual-weighted-residual approach that aims to reduce error in an output quantity of interest. The resulting method generates a hierarchy of subspaces online without requiring large-scale operations or full-order-model solves. Furthermore, it enables the reduced-order model to satisfy any prescribed error tolerance regardless of its original fidelity, as a completely refined reduced-order model is mathematically equivalent to the original full-order model. Experiments on a parameterized inviscid Burgers equation highlight the ability of the method to capture phenomena (e.g., moving shocks) not contained in the span of the original reduced basis.

  1. The Effective Combination of Mesh Adaptation and Non-linearThermo...

    Office of Scientific and Technical Information (OSTI)

    The Effective Combination of Mesh Adaptation and Non-linear Thermo-mechanical Solution Components for the Modeling of Weld Failures. Citation Details In-Document Search Title: The...

  2. The design of a parallel adaptive paving all-quadrilateral meshing algorithm

    SciTech Connect (OSTI)

    Tautges, T.J.; Lober, R.R.; Vaughan, C.

    1995-08-01

    Adaptive finite element analysis demands a great deal of computational resources, and as such is most appropriately solved in a massively parallel computer environment. This analysis will require other parallel algorithms before it can fully utilize MP computers, one of which is parallel adaptive meshing. A version of the paving algorithm is being designed which operates in parallel but which also retains the robustness and other desirable features present in the serial algorithm. Adaptive paving in a production mode is demonstrated using a Babuska-Rheinboldt error estimator on a classic linearly elastic plate problem. The design of the parallel paving algorithm is described, and is based on the decomposition of a surface into {open_quotes}virtual{close_quotes} surfaces. The topology of the virtual surface boundaries is defined using mesh entities (mesh nodes and edges) so as to allow movement of these boundaries with smoothing and other operations. This arrangement allows the use of the standard paving algorithm on subdomain interiors, after the negotiation of the boundary mesh.

  3. Solving kinetic equations with adaptive mesh in phase space for rarefied gas dynamics and plasma physics (Invited)

    SciTech Connect (OSTI)

    Kolobov, Vladimir; Arslanbekov, Robert; Frolova, Anna

    2014-12-09

    The paper describes an Adaptive Mesh in Phase Space (AMPS) technique for solving kinetic equations with deterministic mesh-based methods. The AMPS technique allows automatic generation of adaptive Cartesian mesh in both physical and velocity spaces using a Tree-of-Trees data structure. We illustrate advantages of AMPS for simulations of rarefied gas dynamics and electron kinetics on low temperature plasmas. In particular, we consider formation of the velocity distribution functions in hypersonic flows, particle kinetics near oscillating boundaries, and electron kinetics in a radio-frequency sheath. AMPS provide substantial savings in computational cost and increased efficiency of the mesh-based kinetic solvers.

  4. Minimising the error in eigenvalue calculations involving the Boltzmann transport equation using goal-based adaptivity on unstructured meshes

    SciTech Connect (OSTI)

    Goffin, Mark A.; Baker, Christopher M.J.; Buchan, Andrew G.; Pain, Christopher C.; Eaton, Matthew D.; Smith, Paul N.

    2013-06-01

    This article presents a method for goal-based anisotropic adaptive methods for the finite element method applied to the Boltzmann transport equation. The neutron multiplication factor, k{sub eff}, is used as the goal of the adaptive procedure. The anisotropic adaptive algorithm requires error measures for k{sub eff} with directional dependence. General error estimators are derived for any given functional of the flux and applied to k{sub eff} to acquire the driving force for the adaptive procedure. The error estimators require the solution of an appropriately formed dual equation. Forward and dual error indicators are calculated by weighting the Hessian of each solution with the dual and forward residual respectively. The Hessian is used as an approximation of the interpolation error in the solution which gives rise to the directional dependence. The two indicators are combined to form a single error metric that is used to adapt the finite element mesh. The residual is approximated using a novel technique arising from the sub-grid scale finite element discretisation. Two adaptive routes are demonstrated: (i) a single mesh is used to solve all energy groups, and (ii) a different mesh is used to solve each energy group. The second method aims to capture the benefit from representing the flux from each energy group on a specifically optimised mesh. The k{sub eff} goal-based adaptive method was applied to three examples which illustrate the superior accuracy in criticality problems that can be obtained.

  5. COLLABORATIVE RESEARCH: CONTINUOUS DYNAMIC GRID ADAPTATION IN A GLOBAL ATMOSPHERIC MODEL: APPLICATION AND REFINEMENT

    SciTech Connect (OSTI)

    Gutowski, William J.; Prusa, Joseph M.; Smolarkiewicz, Piotr K.

    2012-05-08

    This project had goals of advancing the performance capabilities of the numerical general circulation model EULAG and using it to produce a fully operational atmospheric global climate model (AGCM) that can employ either static or dynamic grid stretching for targeted phenomena. The resulting AGCM combined EULAG's advanced dynamics core with the "physics" of the NCAR Community Atmospheric Model (CAM). Effort discussed below shows how we improved model performance and tested both EULAG and the coupled CAM-EULAG in several ways to demonstrate the grid stretching and ability to simulate very well a wide range of scales, that is, multi-scale capability. We leveraged our effort through interaction with an international EULAG community that has collectively developed new features and applications of EULAG, which we exploited for our own work summarized here. Overall, the work contributed to over 40 peer-reviewed publications and over 70 conference/workshop/seminar presentations, many of them invited. 3a. EULAG Advances EULAG is a non-hydrostatic, parallel computational model for all-scale geophysical flows. EULAG's name derives from its two computational options: EULerian (flux form) or semi-LAGrangian (advective form). The model combines nonoscillatory forward-in-time (NFT) numerical algorithms with a robust elliptic Krylov solver. A signature feature of EULAG is that it is formulated in generalized time-dependent curvilinear coordinates. In particular, this enables grid adaptivity. In total, these features give EULAG novel advantages over many existing dynamical cores. For EULAG itself, numerical advances included refining boundary conditions and filters for optimizing model performance in polar regions. We also added flexibility to the model's underlying formulation, allowing it to work with the pseudo-compressible equation set of Durran in addition to EULAG's standard anelastic formulation. Work in collaboration with others also extended the demonstrated range of

  6. An adaptive grid refinement strategy for the simulation of negative streamers

    SciTech Connect (OSTI)

    Montijn, C. . E-mail: carolynne.montijn@cwi.nl; Hundsdorfer, W. . E-mail: willem.hundsdorfer@cwi.nl; Ebert, U. . E-mail: ute.ebert@cwi.nl

    2006-12-10

    The evolution of negative streamers during electric breakdown of a non-attaching gas can be described by a two-fluid model for electrons and positive ions. It consists of continuity equations for the charged particles including drift, diffusion and reaction in the local electric field, coupled to the Poisson equation for the electric potential. The model generates field enhancement and steep propagating ionization fronts at the tip of growing ionized filaments. An adaptive grid refinement method for the simulation of these structures is presented. It uses finite volume spatial discretizations and explicit time stepping, which allows the decoupling of the grids for the continuity equations from those for the Poisson equation. Standard refinement methods in which the refinement criterion is based on local error monitors fail due to the pulled character of the streamer front that propagates into a linearly unstable state. We present a refinement method which deals with all these features. Tests on one-dimensional streamer fronts as well as on three-dimensional streamers with cylindrical symmetry (hence effectively 2D for numerical purposes) are carried out successfully. Results on fine grids are presented, they show that such an adaptive grid method is needed to capture the streamer characteristics well. This refinement strategy enables us to adequately compute negative streamers in pure gases in the parameter regime where a physical instability appears: branching streamers.

  7. 6th International Meshing Roundtable '97

    SciTech Connect (OSTI)

    White, D.

    1997-09-01

    The goal of the 6th International Meshing Roundtable is to bring together researchers and developers from industry, academia, and government labs in a stimulating, open environment for the exchange of technical information related to the meshing process. In the pas~ the Roundtable has enjoyed significant participation born each of these groups from a wide variety of countries. The Roundtable will consist of technical presentations from contributed papers and abstracts, two invited speakers, and two invited panels of experts discussing topics related to the development and use of automatic mesh generation tools. In addition, this year we will feature a "Bring Your Best Mesh" competition and poster session to encourage discussion and participation from a wide variety of mesh generation tool users. The schedule and evening social events are designed to provide numerous opportunities for informal dialog. A proceedings will be published by Sandia National Laboratories and distributed at the Roundtable. In addition, papers of exceptionally high quaIity will be submitted to a special issue of the International Journal of Computational Geometry and Applications. Papers and one page abstracts were sought that present original results on the meshing process. Potential topics include but are got limited to: Unstructured triangular and tetrahedral mesh generation Unstructured quadrilateral and hexahedral mesh generation Automated blocking and structured mesh generation Mixed element meshing Surface mesh generation Geometry decomposition and clean-up techniques Geometry modification techniques related to meshing Adaptive mesh refinement and mesh quality control Mesh visualization Special purpose meshing algorithms for particular applications Theoretical or novel ideas with practical potential Technical presentations from industrial researchers.

  8. Interpolation Methods and the Accuracy of Lattice-Boltzmann Mesh...

    Office of Scientific and Technical Information (OSTI)

    of Lattice-Boltzmann Mesh Refinement Citation Details In-Document Search Title: Interpolation Methods and the Accuracy of Lattice-Boltzmann Mesh Refinement You are ...

  9. Development of Adaptive Model Refinement (AMoR) for Multiphysics and Multifidelity Problems

    SciTech Connect (OSTI)

    Turinsky, Paul

    2015-02-09

    This project investigated the development and utilization of Adaptive Model Refinement (AMoR) for nuclear systems simulation applications. AMoR refers to utilization of several models of physical phenomena which differ in prediction fidelity. If the highest fidelity model is judged to always provide or exceeded the desired fidelity, than if one can determine the difference in a Quantity of Interest (QoI) between the highest fidelity model and lower fidelity models, one could utilize the fidelity model that would just provide the magnitude of the QoI desired. Assuming lower fidelity models require less computational resources, in this manner computational efficiency can be realized provided the QoI value can be accurately and efficiently evaluated. This work utilized Generalized Perturbation Theory (GPT) to evaluate the QoI, by convoluting the GPT solution with the residual of the highest fidelity model determined using the solution from lower fidelity models. Specifically, a reactor core neutronics problem and thermal-hydraulics problem were studied to develop and utilize AMoR. The highest fidelity neutronics model was based upon the 3D space-time, two-group, nodal diffusion equations as solved in the NESTLE computer code. Added to the NESTLE code was the ability to determine the time-dependent GPT neutron flux. The lower fidelity neutronics model was based upon the point kinetics equations along with utilization of a prolongation operator to determine the 3D space-time, two-group flux. The highest fidelity thermal-hydraulics model was based upon the space-time equations governing fluid flow in a closed channel around a heat generating fuel rod. The Homogenous Equilibrium Mixture (HEM) model was used for the fluid and Finite Difference Method was applied to both the coolant and fuel pin energy conservation equations. The lower fidelity thermal-hydraulic model was based upon the same equations as used for the highest fidelity model but now with coarse spatial

  10. Enhancing adaptive sparse grid approximations and improving refinement strategies using adjoint-based a posteriori error estimates

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Jakeman, J. D.; Wildey, T.

    2015-01-01

    In this paper we present an algorithm for adaptive sparse grid approximations of quantities of interest computed from discretized partial differential equations. We use adjoint-based a posteriori error estimates of the interpolation error in the sparse grid to enhance the sparse grid approximation and to drive adaptivity. We show that utilizing these error estimates provides significantly more accurate functional values for random samples of the sparse grid approximation. We also demonstrate that alternative refinement strategies based upon a posteriori error estimates can lead to further increases in accuracy in the approximation over traditional hierarchical surplus based strategies. Throughout this papermore » we also provide and test a framework for balancing the physical discretization error with the stochastic interpolation error of the enhanced sparse grid approximation.« less

  11. Enhancing adaptive sparse grid approximations and improving refinement strategies using adjoint-based a posteriori error estimates

    SciTech Connect (OSTI)

    Jakeman, J.D. Wildey, T.

    2015-01-01

    In this paper we present an algorithm for adaptive sparse grid approximations of quantities of interest computed from discretized partial differential equations. We use adjoint-based a posteriori error estimates of the physical discretization error and the interpolation error in the sparse grid to enhance the sparse grid approximation and to drive adaptivity of the sparse grid. Utilizing these error estimates provides significantly more accurate functional values for random samples of the sparse grid approximation. We also demonstrate that alternative refinement strategies based upon a posteriori error estimates can lead to further increases in accuracy in the approximation over traditional hierarchical surplus based strategies. Throughout this paper we also provide and test a framework for balancing the physical discretization error with the stochastic interpolation error of the enhanced sparse grid approximation.

  12. A mesh-adaptive collocation technique for the simulation of advection-dominated single- and multiphase transport phenomena in porous media

    SciTech Connect (OSTI)

    Koch, M.

    1995-12-31

    A new mesh-adaptive 1D collocation technique has been developed to efficiently solve transient advection-dominated transport problems in porous media that are governed by a hyperbolic/parabolic (singularly perturbed) PDE. After spatial discretization a singularly perturbed ODE is obtained which is solved by a modification of the COLNEW ODE-collocation code. The latter also contains an adaptive mesh procedure that has been enhanced here to resolve linear and nonlinear transport flow problems with steep fronts where regular FD and FE methods often fail. An implicit first-order backward Euler and a third-order Taylor-Donea technique are employed for the time integration. Numerical simulations on a variety of high Peclet-number transport phenomena as they occur in realistic porous media flow situations are presented. Examples include classical linear advection-diffusion, nonlinear adsorption, two-phase Buckley-Leverett flow without and with capillary forces (Rapoport-Leas equation) and Burgers` equation for inviscid fluid flow. In most of these examples sharp fronts and/or shocks develop which are resolved in an oscillation-free manner by the present adaptive collocation method. The backward Euler method has some amount of numerical dissipation is observed when the time-steps are too large. The third-order Taylor-Donea technique is less dissipative but is more prone to numerical oscillations. The simulations show that for the efficient solution of nonlinear singularly perturbed PDE`s governing flow transport a careful balance must be struck between the optimal mesh adaptation, the nonlinear iteration method and the time-stepping procedure. More theoretical research is needed with this regard.

  13. An adaptive multiblock high-order finite-volume method for solving the shallow-water equations on the sphere

    SciTech Connect (OSTI)

    McCorquodale, Peter; Ullrich, Paul A.; Johansen, Hans; Colella, Phillip

    2015-06-16

    We present a high-order finite-volume approach for solving the shallow-water equations on the sphere, using multiblock grids on the cubed-sphere. This approach combines a Runge--Kutta time discretization with a fourth-order accurate spatial discretization, and includes adaptive mesh refinement and refinement in time. Results of tests show fourth-order convergence for the shallow-water equations as well as for advection in a highly deformational flow. Hierarchical adaptive mesh refinement allows solution error to be achieved that is comparable to that obtained with uniform resolution of the most refined level of the hierarchy, but with many fewer operations.

  14. An adaptive multiblock high-order finite-volume method for solving the shallow-water equations on the sphere

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    McCorquodale, Peter; Ullrich, Paul; Johansen, Hans; Colella, Phillip

    2015-09-04

    We present a high-order finite-volume approach for solving the shallow-water equations on the sphere, using multiblock grids on the cubed-sphere. This approach combines a Runge--Kutta time discretization with a fourth-order accurate spatial discretization, and includes adaptive mesh refinement and refinement in time. Results of tests show fourth-order convergence for the shallow-water equations as well as for advection in a highly deformational flow. Hierarchical adaptive mesh refinement allows solution error to be achieved that is comparable to that obtained with uniform resolution of the most refined level of the hierarchy, but with many fewer operations.

  15. Method of modifying a volume mesh using sheet insertion

    DOE Patents [OSTI]

    Borden, Michael J.; Shepherd, Jason F.

    2006-08-29

    A method and machine-readable medium provide a technique to modify a hexahedral finite element volume mesh using dual generation and sheet insertion. After generating a dual of a volume stack (mesh), a predetermined algorithm may be followed to modify (refine) the volume mesh of hexahedral elements. The predetermined algorithm may include the steps of locating a sheet of hexahedral mesh elements, determining a plurality of hexahedral elements within the sheet to refine, shrinking the plurality of elements, and inserting a new sheet of hexahedral elements adjacently to modify the volume mesh. Additionally, another predetermined algorithm using mesh cutting may be followed to modify a volume mesh.

  16. MeshKit

    Energy Science and Technology Software Center (OSTI)

    2010-10-05

    MeshKit is an open-source library of mesh generation functionality. MeshKit has general mesh manipulation and generation functions such as Copoy, Move, Rotate and Extrude mesh. In addition, new quad mesh and embedded boundary Cartesian mesh algorithm (EB Mesh) are included. Interfaces to several public domain meshing algorithms (TetGen, netgen, triangle, Gmsh, camal) are also offered. This library interacts with mesh data mostly through iMesh including accessing the mesh in parallel. It also can interact withmore » iGeom interface to provide geometry functionality such as importing solid model based geometries. iGeom and IMesh are implemented in the CGM and MOAB packages, respectively. For some non-existing function in iMesh such as tree-construction and ray-tracing, MeshKit also interacts with MOAB functions directly.« less

  17. Quadrilateral/hexahedral finite element mesh coarsening

    DOE Patents [OSTI]

    Staten, Matthew L; Dewey, Mark W; Scott, Michael A; Benzley, Steven E

    2012-10-16

    A technique for coarsening a finite element mesh ("FEM") is described. This technique includes identifying a coarsening region within the FEM to be coarsened. Perimeter chords running along perimeter boundaries of the coarsening region are identified. The perimeter chords are redirected to create an adaptive chord separating the coarsening region from a remainder of the FEM. The adaptive chord runs through mesh elements residing along the perimeter boundaries of the coarsening region. The adaptive chord is then extracted to coarsen the FEM.

  18. Solution-verified reliability analysis and design of bistable MEMS using error estimation and adaptivity.

    SciTech Connect (OSTI)

    Eldred, Michael Scott; Subia, Samuel Ramirez; Neckels, David; Hopkins, Matthew Morgan; Notz, Patrick K.; Adams, Brian M.; Carnes, Brian; Wittwer, Jonathan W.; Bichon, Barron J.; Copps, Kevin D.

    2006-10-01

    This report documents the results for an FY06 ASC Algorithms Level 2 milestone combining error estimation and adaptivity, uncertainty quantification, and probabilistic design capabilities applied to the analysis and design of bistable MEMS. Through the use of error estimation and adaptive mesh refinement, solution verification can be performed in an automated and parameter-adaptive manner. The resulting uncertainty analysis and probabilistic design studies are shown to be more accurate, efficient, reliable, and convenient.

  19. Mesh Morphing Pier Analysis

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Application of Mesh Morphing in STAR-CCM+ to Analysis of Scour at Cylindrical Piers Mesh morphing is a fluid structure interaction capability in STAR-CCM+ to move vertices in the computational mesh in a way that preserves mesh quality when a boundary moves. The equations being solved include terms that account for the motion of the mesh maintaining mass and property balances during the solution process. Initial work on leveraging the mesh morphing FSI capability for efficient application to

  20. A node-centered local refinement algorithm for poisson's equation in complex geometries

    SciTech Connect (OSTI)

    McCorquodale, Peter; Colella, Phillip; Grote, David P.; Vay, Jean-Luc

    2004-05-04

    This paper presents a method for solving Poisson's equation with Dirichlet boundary conditions on an irregular bounded three-dimensional region. The method uses a nodal-point discretization and adaptive mesh refinement (AMR) on Cartesian grids, and the AMR multigrid solver of Almgren. The discrete Laplacian operator at internal boundaries comes from either linear or quadratic (Shortley-Weller) extrapolation, and the two methods are compared. It is shown that either way, solution error is second order in the mesh spacing. Error in the gradient of the solution is first order with linear extrapolation, but second order with Shortley-Weller. Examples are given with comparison with the exact solution. The method is also applied to a heavy-ion fusion accelerator problem, showing the advantage of adaptivity.

  1. Rietveld Refinement

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    profile Chi-squared value 4.944 Progress so far Inverse Modeling Method 2 Fourier Method Data Model Refined Structure 40 60 80 100 120 140 160 0 1000 2000 3000 4000...

  2. Mesh Quality Improvement Toolkit

    Energy Science and Technology Software Center (OSTI)

    2002-11-15

    MESQUITE is a linkable software library to be used by simulation and mesh generation tools to improve the quality of meshes. Mesh quality is improved by node movement and/or local topological modifications. Various aspects of mesh quality such as smoothness, element shape, size, and orientation are controlled by choosing the appropriate mesh qualtiy metric, and objective function tempate, and a numerical optimization solver to optimize the quality of meshes, MESQUITE uses the TSTT mesh interfacemore » specification to provide an interoperable toolkit that can be used by applications which adopt the standard. A flexible code design makes it easy for meshing researchers to add additional mesh quality metrics, templates, and solvers to develop new quality improvement algorithms by making use of the MESQUITE infrastructure.« less

  3. Midcourse Refinements of Financing Strategies

    Broader source: Energy.gov [DOE]

    Better Buildings Neighborhood Program Financing Peer Exchange Call: Midcourse Refinements of Financing Strategies, Call Slides and Discussion Summary, March 29, 2012. During this webinar participants discussed how programs have adapted and refined their financing strategies based on initial implementation experience and learning.

  4. Adaptive methods and parallel computation for partial differential equations. Final report

    SciTech Connect (OSTI)

    Biswas, R.; Benantar, M.; Flaherty, J.E.

    1992-05-01

    Consider the adaptive solution of two-dimensional vector systems of hyperbolic and elliptic partial differential equations on shared-memory parallel computers. Hyperbolic systems are approximated by an explicit finite volume technique and solved by a recursive local mesh refinement procedure on a tree-structured grid. Local refinement of the time steps and spatial cells of a coarse base mesh is performed in regions where a refinement indicator exceeds a prescribed tolerance. Computational procedures that sequentially traverse the tree while processing solutions on each grid in parallel, that process solutions at the same tree level in parallel, and that dynamically assign processors to nodes of the tree have been developed and applied to an example. Computational results comparing a variety of heuristic processor load balancing techniques and refinement strategies are presented.

  5. AZEuS: AN ADAPTIVE ZONE EULERIAN SCHEME FOR COMPUTATIONAL MAGNETOHYDRODYNAMICS

    SciTech Connect (OSTI)

    Ramsey, Jon P.; Clarke, David A. [Institute for Computational Astrophysics, Department of Astronomy and Physics, Saint Mary's University, Halifax, Nova Scotia B3H 3C3 (Canada); Men'shchikov, Alexander B. [Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, IRFU/SAp, CEA Saclay, 91191 Gif-sur-Yvette (France)

    2012-03-01

    A new adaptive mesh refinement (AMR) version of the ZEUS-3D astrophysical magnetohydrodynamical fluid code, AZEuS, is described. The AMR module in AZEuS has been completely adapted to the staggered mesh that characterizes the ZEUS family of codes on which scalar quantities are zone-centered and vector components are face-centered. In addition, for applications using static grids, it is necessary to use higher-order interpolations for prolongation to minimize the errors caused by waves crossing from a grid of one resolution to another. Finally, solutions to test problems in one, two, and three dimensions in both Cartesian and spherical coordinates are presented.

  6. Midcourse Refinements of Financing Strategies | Department of...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    March 29, 2012. During this webinar participants discussed how programs have adapted and refined their financing strategies based on initial implementation experience and learning. ...

  7. Compact Mesh Generator

    Energy Science and Technology Software Center (OSTI)

    2007-02-02

    The CMG is a small, lightweight, structured mesh generation code. It features a simple text input parser that allows setup of various meshes via a small set of text commands. Mesh generation data can be output to text, the silo file format, or the API can be directly queried by applications. It can run serially or in parallel via MPI. The CMG includes the ability to specify varius initial conditions on a mesh via meshmore » tags.« less

  8. phdMesh

    Energy Science and Technology Software Center (OSTI)

    2008-01-01

    Parallel Heterogeneous Dynamic unstructured Mesh (phdMesh) data structure library and integration testing code that performs dynamic load balancing of the data structure and parallel geometric proximity search on a contrived test problem. The phdMesh library is intended to be module within a finite element or finite volume library or code. The integration testing code is intended to provide a compact and highly portable performance evaluation code for parallel computing systems.

  9. Refiner Crude Oil Inputs

    U.S. Energy Information Administration (EIA) Indexed Site

    Net Inputs (Refiner and Blender) of RBOB Blending Components Net Inputs (Refiner and Blender) of CBOB Blending Components Net Inputs (Refiner and Blender) of GTAB Blending ...

  10. Optimization of tetrahedral meshes

    SciTech Connect (OSTI)

    Briere De L`Isle, E.; George, P.L.

    1995-12-31

    Finite element computations are all the more exact if we start from {open_quotes}good{close_quotes} elements. We are interested in meshes where the elements are tetrahedra and we shall develop utilities allowing us to improve the quality of these meshes.

  11. Parallel partitioning strategies for the adaptive solution of conservation laws

    SciTech Connect (OSTI)

    Devine, K.D.; Flaherty, J.E.; Loy, R.M.

    1995-12-31

    We describe and examine the performance of adaptive methods for Solving hyperbolic systems of conservation laws on massively parallel computers. The differential system is approximated by a discontinuous Galerkin finite element method with a hierarchical Legendre piecewise polynomial basis for the spatial discretization. Fluxes at element boundaries are computed by solving an approximate Riemann problem; a projection limiter is applied to keep the average solution monotone; time discretization is performed by Runge-Kutta integration; and a p-refinement-based error estimate is used as an enrichment indicator. Adaptive order (p-) and mesh (h-) refinement algorithms are presented and demonstrated. Using an element-based dynamic load balancing algorithm called tiling and adaptive p-refinement, parallel efficiencies of over 60% are achieved on a 1024-processor nCUBE/2 hypercube. We also demonstrate a fast, tree-based parallel partitioning strategy for three-dimensional octree-structured meshes. This method produces partition quality comparable to recursive spectral bisection at a greatly reduced cost.

  12. Documentation for MeshKit - Reactor Geometry (&mesh) Generator

    SciTech Connect (OSTI)

    Jain, Rajeev; Mahadevan, Vijay

    2015-09-30

    This report gives documentation for using MeshKit’s Reactor Geometry (and mesh) Generator (RGG) GUI and also briefly documents other algorithms and tools available in MeshKit. RGG is a program designed to aid in modeling and meshing of complex/large hexagonal and rectilinear reactor cores. RGG uses Argonne’s SIGMA interfaces, Qt and VTK to produce an intuitive user interface. By integrating a 3D view of the reactor with the meshing tools and combining them into one user interface, RGG streamlines the task of preparing a simulation mesh and enables real-time feedback that reduces accidental scripting mistakes that could waste hours of meshing. RGG interfaces with MeshKit tools to consolidate the meshing process, meaning that going from model to mesh is as easy as a button click. This report is designed to explain RGG v 2.0 interface and provide users with the knowledge and skills to pilot RGG successfully. Brief documentation of MeshKit source code, tools and other algorithms available are also presented for developers to extend and add new algorithms to MeshKit. RGG tools work in serial and parallel and have been used to model complex reactor core models consisting of conical pins, load pads, several thousands of axially varying material properties of instrumentation pins and other interstices meshes.

  13. Geometric approaches to mesh generation

    SciTech Connect (OSTI)

    Hoffmann, C.M.

    1995-12-31

    We review three approaches to mesh generation that axe based on analyzing and accounting for the geometric structure of the domain. In the first approach, due to Armstrong, the domain is partitioned into subdomains based on the medial-axis transform, a tool for analyzing spatial structures. In the second approach, due to Cox, the design history defines a geometric structure of the domain. The design primitives of that structure are meshed separately, and mesh overlap is accounted for by coupling equations. The third approach argues that mesh generation ought to be integrated into the shape design process, by meshing design features separately and resolving overlapping meshes by standard geometric computations.

  14. Mesh2d

    Energy Science and Technology Software Center (OSTI)

    2011-12-31

    Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less

  15. SUPERIMPOSED MESH PLOTTING IN MCNP

    SciTech Connect (OSTI)

    J. HENDRICKS

    2001-02-01

    The capability to plot superimposed meshes has been added to MCNP{trademark}. MCNP4C featured a superimposed mesh weight window generator which enabled users to set up geometries without having to subdivide geometric cells for variance reduction. The variance reduction was performed with weight windows on a rectangular or cylindrical mesh superimposed over the physical geometry. Experience with the new capability was favorable but also indicated that a number of enhancements would be very beneficial, particularly a means of visualizing the mesh and its values. The mathematics for plotting the mesh and its values is described here along with a description of other upgrades.

  16. Mesh Algorithms for PDE with Sieve I: Mesh Distribution

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Knepley, Matthew G.; Karpeev, Dmitry A.

    2009-01-01

    We have developed a new programming framework, called Sieve, to support parallel numerical partial differential equation(s) (PDE) algorithms operating over distributed meshes. We have also developed a reference implementation of Sieve in C++ as a library of generic algorithms operating on distributed containers conforming to the Sieve interface. Sieve makes instances of the incidence relation, or arrows, the conceptual first-class objects represented in the containers. Further, generic algorithms acting on this arrow container are systematically used to provide natural geometric operations on the topology and also, through duality, on the data. Finally, coverings and duality are used to encode notmore » only individual meshes, but all types of hierarchies underlying PDE data structures, including multigrid and mesh partitions. In order to demonstrate the usefulness of the framework, we show how the mesh partition data can be represented and manipulated using the same fundamental mechanisms used to represent meshes. We present the complete description of an algorithm to encode a mesh partition and then distribute a mesh, which is independent of the mesh dimension, element shape, or embedding. Moreover, data associated with the mesh can be similarly distributed with exactly the same algorithm. The use of a high level of abstraction within the Sieve leads to several benefits in terms of code reuse, simplicity, and extensibility. We discuss these benefits and compare our approach to other existing mesh libraries.« less

  17. Mesh Oriented datABase

    Energy Science and Technology Software Center (OSTI)

    2004-04-01

    MOAB is a component for representing and evaluating mesh data. MOAB can store stuctured and unstructured mesh, consisting of elements in the finite element "zoo". The functional interface to MOAB is simple yet powerful, allowing the representation of many types of metadata commonly found on the mesh. MOAB is optimized for efficiency in space and time, based on access to mesh in chunks rather than through individual entities, while also versatile enough to support individualmore » entity access. The MOAB data model consists of a mesh interface instance, mesh entities (vertices and elements), sets, and tags. Entities are addressed through handles rather than pointers, to allow the underlying representation of an entity to change without changing the handle to that entity. Sets are arbitrary groupings of mesh entities and other sets. Sets also support parent/child relationships as a relation distinct from sets containing other sets. The directed-graph provided by set parent/child relationships is useful for modeling topological relations from a geometric model or other metadata. Tags are named data which can be assigned to the mesh as a whole, individual entities, or sets. Tags are a mechanism for attaching data to individual entities and sets are a mechanism for describing relations between entities; the combination of these two mechanisms isa powerful yet simple interface for representing metadata or application-specific data. For example, sets and tags can be used together to describe geometric topology, boundary condition, and inter-processor interface groupings in a mesh. MOAB is used in several ways in various applications. MOAB serves as the underlying mesh data representation in the VERDE mesh verification code. MOAB can also be used as a mesh input mechanism, using mesh readers induded with MOAB, or as a t’anslator between mesh formats, using readers and writers included with MOAB.« less

  18. ITP Petroleum Refining: Energy Bandwidth for Petroleum Refining...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Bandwidth for Petroleum Refining Processes ITP Petroleum Refining: Energy Bandwidth for Petroleum Refining Processes bandwidth.pdf (1.16 MB) More Documents & Publications ITP ...

  19. ITP Petroleum Refining: Profile of the Petroleum Refining Industry...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Profile of the Petroleum Refining Industry in California: California Industries of the Future Program ITP Petroleum Refining: Profile of the Petroleum Refining Industry in ...

  20. Creating Interoperable Meshing and Discretization Software: The Terascale Simulation Tools and Technology Center

    SciTech Connect (OSTI)

    Brown, D.; Freitag, L.; Glimm, J.

    2002-03-28

    We present an overview of the technical objectives of the Terascale Simulation Tools and Technologies center. The primary goal of this multi-institution collaboration is to develop technologies that enable application scientists to easily use multiple mesh and discretization strategies within a single simulation on terascale computers. The discussion focuses on our efforts to create interoperable mesh generation tools, high-order discretization techniques, and adaptive meshing strategies.

  1. CREATING INTEROPERABLE MESHING AND DISCRETIZATION SOFTWARE: THE TERASCALE SIMULATION TOOLS AND TECHNOLOGY CENTER.

    SciTech Connect (OSTI)

    BROWN,D.; FREITAG,L.; GLIMM,J.

    2002-06-02

    We present an overview of the technical objectives of the Terascale Simulation Tools and Technologies center. The primary goal of this multi-institution collaboration is to develop technologies that enable application scientists to easily use multiple mesh and discretization strategies within a single simulation on terascale computers. The discussion focuses on our efforts to create interoperable mesh generation tools, high-order discretization techniques, and adaptive meshing strategies.

  2. Step 3: Project Refinement

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    3: Project Refinement 2 1 Potential 3 Refinement 4 Implementation 5 Operations & Maintenance 2 Options 3 Refinement 1/28/2016 2 3 FUNDING AND FINANCING OPTIONS Project Ownership Financing structure is highly dependent on size of the project and the capital available for a given project: * Tribe owns the project (cash purchase or debt) * Tribe hosts the project and buys the electricity (power purchase agreement) * Tribe partners with private sector and co-owns the project (uncertainties about

  3. Quadrilateral finite element mesh coarsening

    DOE Patents [OSTI]

    Staten, Matthew L; Dewey, Mark W; Benzley, Steven E

    2012-10-16

    Techniques for coarsening a quadrilateral mesh are described. These techniques include identifying a coarsening region within the quadrilateral mesh to be coarsened. Quadrilateral elements along a path through the coarsening region are removed. Node pairs along opposite sides of the path are identified. The node pairs along the path are then merged to collapse the path.

  4. Spherical geodesic mesh generation (Conference) | SciTech Connect

    Office of Scientific and Technical Information (OSTI)

    In ALE simulations with moving meshes, mesh topology has a direct influence on feature ... This mesh topology is fundamentally different from traditional mesh topologies and ...

  5. Collaboration Topics - Meshing | National Nuclear Security Administration |

    National Nuclear Security Administration (NNSA)

    (NNSA) Meshing This work focuses on automatic algorithms for 3-D mesh generation and their extension to scalable parallel environments. The objective is to explore new theoretical advancements and develop new algorithms in mesh generation that will lead to practical advances in automatic meshing with an emphasis on all-hexahedral mesh generation technology. The development of technologies to facilitate model preparation for analysis such as CAD support for meshing and geometry simplification

  6. Refiners trade hydroprocessing experience

    SciTech Connect (OSTI)

    Not Available

    1984-04-16

    Hydrogen treating and conversion processes less severe than hydrocracking abound in most refineries and therefore were a subject of high interest at the National Petroleum Refiners Association question and answer session on refining technology. The present paper, which is the second abstract of the transcript of the most recent meeting, covers hydroprocessing and some of its mechanical, process, and catalytic aspects.

  7. Collaboration Topics - Meshing | National Nuclear Security Administrat...

    National Nuclear Security Administration (NNSA)

    The objective is to explore new theoretical advancements and develop new algorithms in mesh generation that will lead to practical advances in automatic meshing with an emphasis on ...

  8. Nanowire mesh solar fuels generator

    DOE Patents [OSTI]

    Yang, Peidong; Chan, Candace; Sun, Jianwei; Liu, Bin

    2016-05-24

    This disclosure provides systems, methods, and apparatus related to a nanowire mesh solar fuels generator. In one aspect, a nanowire mesh solar fuels generator includes (1) a photoanode configured to perform water oxidation and (2) a photocathode configured to perform water reduction. The photocathode is in electrical contact with the photoanode. The photoanode may include a high surface area network of photoanode nanowires. The photocathode may include a high surface area network of photocathode nanowires. In some embodiments, the nanowire mesh solar fuels generator may include an ion conductive polymer infiltrating the photoanode and the photocathode in the region where the photocathode is in electrical contact with the photoanode.

  9. Toward Interoperable Mesh, Geometry and Field Components for PDE Simulation Development

    SciTech Connect (OSTI)

    Chand, K K; Diachin, L F; Li, X; Ollivier-Gooch, C; Seol, E S; Shephard, M; Tautges, T; Trease, H

    2005-07-11

    Mesh-based PDE simulation codes are becoming increasingly sophisticated and rely on advanced meshing and discretization tools. Unfortunately, it is still difficult to interchange or interoperate tools developed by different communities to experiment with various technologies or to develop new capabilities. To address these difficulties, we have developed component interfaces designed to support the information flow of mesh-based PDE simulations. We describe this information flow and discuss typical roles and services provided by the geometry, mesh, and field components of the simulation. Based on this delineation for the roles of each component, we give a high-level description of the abstract data model and set of interfaces developed by the Department of Energy's Interoperable Tools for Advanced Petascale Simulation (ITAPS) center. These common interfaces are critical to our interoperability goal, and we give examples of several services based upon these interfaces including mesh adaptation and mesh improvement.

  10. Multilevel adaptive solution procedure for material nonlinear problems in visual programming environment

    SciTech Connect (OSTI)

    Kim, D.; Ghanem, R.

    1994-12-31

    Multigrid solution technique to solve a material nonlinear problem in a visual programming environment using the finite element method is discussed. The nonlinear equation of equilibrium is linearized to incremental form using Newton-Rapson technique, then multigrid solution technique is used to solve linear equations at each Newton-Rapson step. In the process, adaptive mesh refinement, which is based on the bisection of a pair of triangles, is used to form grid hierarchy for multigrid iteration. The solution process is implemented in a visual programming environment with distributed computing capability, which enables more intuitive understanding of solution process, and more effective use of resources.

  11. NAFTA opportunities: Petroleum refining

    SciTech Connect (OSTI)

    Not Available

    1993-01-01

    The North American Free Trade Agreement (NAFTA) creates a more transparent environment for the sale of refined petroleum products to Mexico, and locks in access to Canada's relatively open market for these products. Canada and Mexico are sizable United States export markets for refined petroleum products, with exports of $556 million and $864 million, respectively, in 1992. These markets represent approximately 24 percent of total U.S. exports of these goods.

  12. User Manual for the PROTEUS Mesh Tools

    SciTech Connect (OSTI)

    Smith, Micheal A.; Shemon, Emily R.

    2015-06-01

    This report describes the various mesh tools that are provided with the PROTEUS code giving both descriptions of the input and output. In many cases the examples are provided with a regression test of the mesh tools. The most important mesh tools for any user to consider using are the MT_MeshToMesh.x and the MT_RadialLattice.x codes. The former allows the conversion between most mesh types handled by PROTEUS while the second allows the merging of multiple (assembly) meshes into a radial structured grid. Note that the mesh generation process is recursive in nature and that each input specific for a given mesh tool (such as .axial or .merge) can be used as “mesh” input for any of the mesh tools discussed in this manual.

  13. SAMRAI V3.1

    Energy Science and Technology Software Center (OSTI)

    001825MLTPL01 Structured Adaptive Mesh Refinement Application Infrastructure https://computation.llnl.gov/casc/samrai

  14. A Discontinuous Petrov-Galerkin Methodology for Adaptive Solutions to the Incompressible Navier-Stokes Equations

    SciTech Connect (OSTI)

    Roberts, Nathan V.; Demkowiz, Leszek; Moser, Robert

    2015-11-15

    The discontinuous Petrov-Galerkin methodology with optimal test functions (DPG) of Demkowicz and Gopalakrishnan [18, 20] guarantees the optimality of the solution in an energy norm, and provides several features facilitating adaptive schemes. Whereas Bubnov-Galerkin methods use identical trial and test spaces, Petrov-Galerkin methods allow these function spaces to differ. In DPG, test functions are computed on the fly and are chosen to realize the supremum in the inf-sup condition; the method is equivalent to a minimum residual method. For well-posed problems with sufficiently regular solutions, DPG can be shown to converge at optimal rates—the inf-sup constants governing the convergence are mesh-independent, and of the same order as those governing the continuous problem [48]. DPG also provides an accurate mechanism for measuring the error, and this can be used to drive adaptive mesh refinements. We employ DPG to solve the steady incompressible Navier-Stokes equations in two dimensions, building on previous work on the Stokes equations, and focusing particularly on the usefulness of the approach for automatic adaptivity starting from a coarse mesh. We apply our approach to a manufactured solution due to Kovasznay as well as the lid-driven cavity flow, backward-facing step, and flow past a cylinder problems.

  15. Self-Organizing Mesh Generation

    Energy Science and Technology Software Center (OSTI)

    1991-11-01

    A set of five programs which make up a self organizing mesh generation package. QMESH generates meshes having quadrilateral elements on arbitrarily shaped two-dimensional (planar or axisymmetric) bodies. It is designed for use with two-dimensional finite element analysis applications. A flexible hierarchal input scheme is used to describe bodies to QMESH as collections of regions. A mesh for each region is developed independently, with the final assembly and bandwidth minimization performed by the independent program,more » RENUM or RENUM8. RENUM is applied when four-node elements are desired. Eight node elements (with mid side nodes) may be obtained with RENUM8. QPLOT and QPLOT8 are plot programs for meshes generated by the QMESH/RENUM and QMESH/RENUM8 program pairs respectively. QPLOT and QPLOT8 automatically section the mesh into appropriately-sized sections for legible display of node and element numbers, An overall plot showing the position of the selected plot areas is produced.« less

  16. MeshMaker: Configurable Meshing Framework for Eco-Hydrology Models

    DOE Data Explorer [Office of Scientific and Technical Information (OSTI)]

    Nathan Collier; Jitendra Kumar

    2016-02-09

    MeshMaker is a Python-based framework for generation of high quality structured and unstructured grid computational meshes for Eco-Hydrological models.

  17. Spherical geodesic mesh generation (Conference) | SciTech Connect

    Office of Scientific and Technical Information (OSTI)

    This mesh topology is fundamentally different from traditional mesh topologies and displays superior qualities such as topological symmetry. This work describes the geodesic mesh ...

  18. Simulation of single mode Richtmyer-Meshkov instability using the adaptive free Lagrange method

    SciTech Connect (OSTI)

    Burton, D.E.; Harrison, A.K.

    1991-06-01

    Because of extreme deformation, instability problems have been traditionally modeled using Eulerian techniques or Arbitrary Lagrange Euler (ALE) techniques. The free-Lagrange (FL) method first introduced by Crowley offers a third basic alternative. The specific FL method described in this paper consists of: (1) a 2D staggered-grid hydrodynamics (SGH) differencing scheme appropriate to an unstructured mesh consisting of triangular and quadrilateral zones; (2) a set of primitive mesh optimization algorithms; and (3) a set of rules for adaptive mesh refinement which we refer to as the adaptive free Lagrange (AFL) method. In a previous study, Burton showed that the spurious vorticity commonly associated with SGH differencing results from a failure to explicity conserve angular momentum in the difference equations. This led to a significantly improved technique for removing the spurious vorticity, called spurious vorticity damping (SVD). SVD was shown to produce superior results in a test problem which had no intrinsic vorticity. In the present study, SVD is shown to also produce superior results in a problem with intrinsic vorticity. The AFL and SVD methods are applied to a low Mach number shock tube problem and compared with other methods. AFL seems to perform as well as Eulerian methods and is somewhat better at resolving fine structures in the calculation. The SVD method produced results far superior to other methods tried. 21 refs., 10 figs.

  19. A Parallel Ghosting Algorithm for The Flexible Distributed Mesh Database

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Mubarak, Misbah; Seol, Seegyoung; Lu, Qiukai; Shephard, Mark S.

    2013-01-01

    Critical to the scalability of parallel adaptive simulations are parallel control functions including load balancing, reduced inter-process communication and optimal data decomposition. In distributed meshes, many mesh-based applications frequently access neighborhood information for computational purposes which must be transmitted efficiently to avoid parallel performance degradation when the neighbors are on different processors. This article presents a parallel algorithm of creating and deleting data copies, referred to as ghost copies, which localize neighborhood data for computation purposes while minimizing inter-process communication. The key characteristics of the algorithm are: (1) It can create ghost copies of any permissible topological order inmore » a 1D, 2D or 3D mesh based on selected adjacencies. (2) It exploits neighborhood communication patterns during the ghost creation process thus eliminating all-to-all communication. (3) For applications that need neighbors of neighbors, the algorithm can create n number of ghost layers up to a point where the whole partitioned mesh can be ghosted. Strong and weak scaling results are presented for the IBM BG/P and Cray XE6 architectures up to a core count of 32,768 processors. The algorithm also leads to scalable results when used in a parallel super-convergent patch recovery error estimator, an application that frequently accesses neighborhood data to carry out computation.« less

  20. SHARP/PRONGHORN Interoperability: Mesh Generation

    SciTech Connect (OSTI)

    Avery Bingham; Javier Ortensi

    2012-09-01

    Progress toward collaboration between the SHARP and MOOSE computational frameworks has been demonstrated through sharing of mesh generation and ensuring mesh compatibility of both tools with MeshKit. MeshKit was used to build a three-dimensional, full-core very high temperature reactor (VHTR) reactor geometry with 120-degree symmetry, which was used to solve a neutron diffusion critical eigenvalue problem in PRONGHORN. PRONGHORN is an application of MOOSE that is capable of solving coupled neutron diffusion, heat conduction, and homogenized flow problems. The results were compared to a solution found on a 120-degree, reflected, three-dimensional VHTR mesh geometry generated by PRONGHORN. The ability to exchange compatible mesh geometries between the two codes is instrumental for future collaboration and interoperability. The results were found to be in good agreement between the two meshes, thus demonstrating the compatibility of the SHARP and MOOSE frameworks. This outcome makes future collaboration possible.

  1. Local, Optimization-based Simplicial Mesh Smoothing

    Energy Science and Technology Software Center (OSTI)

    1999-12-09

    OPT-MS is a C software package for the improvement and untangling of simplicial meshes (triangles in 2D, tetrahedra in 3D). Overall mesh quality is improved by iterating over the mesh vertices and adjusting their position to optimize some measure of mesh quality, such as element angle or aspect ratio. Several solution techniques (including Laplacian smoothing, "Smart" Laplacian smoothing, optimization-based smoothing and several combinations thereof) and objective functions (for example, element angle, sin (angle), and aspectmore » ratio) are available to the user for both two and three-dimensional meshes. If the mesh contains invalid elements (those with negative area) a different optimization algorithm for mesh untangling is provided.« less

  2. Minimally refined biomass fuel

    DOE Patents [OSTI]

    Pearson, Richard K.; Hirschfeld, Tomas B.

    1984-01-01

    A minimally refined fluid composition, suitable as a fuel mixture and derived from biomass material, is comprised of one or more water-soluble carbohydrates such as sucrose, one or more alcohols having less than four carbons, and water. The carbohydrate provides the fuel source; water solubilizes the carbohydrates; and the alcohol aids in the combustion of the carbohydrate and reduces the vicosity of the carbohydrate/water solution. Because less energy is required to obtain the carbohydrate from the raw biomass than alcohol, an overall energy savings is realized compared to fuels employing alcohol as the primary fuel.

  3. Introducing Enabling Computational Tools to the Climate Sciences: Multi-Resolution Climate Modeling with Adaptive Cubed-Sphere Grids

    SciTech Connect (OSTI)

    Jablonowski, Christiane

    2015-07-14

    The research investigates and advances strategies how to bridge the scale discrepancies between local, regional and global phenomena in climate models without the prohibitive computational costs of global cloud-resolving simulations. In particular, the research explores new frontiers in computational geoscience by introducing high-order Adaptive Mesh Refinement (AMR) techniques into climate research. AMR and statically-adapted variable-resolution approaches represent an emerging trend for atmospheric models and are likely to become the new norm in future-generation weather and climate models. The research advances the understanding of multi-scale interactions in the climate system and showcases a pathway how to model these interactions effectively with advanced computational tools, like the Chombo AMR library developed at the Lawrence Berkeley National Laboratory. The research is interdisciplinary and combines applied mathematics, scientific computing and the atmospheric sciences. In this research project, a hierarchy of high-order atmospheric models on cubed-sphere computational grids have been developed that serve as an algorithmic prototype for the finite-volume solution-adaptive Chombo-AMR approach. The foci of the investigations have lied on the characteristics of both static mesh adaptations and dynamically-adaptive grids that can capture flow fields of interest like tropical cyclones. Six research themes have been chosen. These are (1) the introduction of adaptive mesh refinement techniques into the climate sciences, (2) advanced algorithms for nonhydrostatic atmospheric dynamical cores, (3) an assessment of the interplay between resolved-scale dynamical motions and subgrid-scale physical parameterizations, (4) evaluation techniques for atmospheric model hierarchies, (5) the comparison of AMR refinement strategies and (6) tropical cyclone studies with a focus on multi-scale interactions and variable-resolution modeling. The results of this research project

  4. Message Passing for Linux Clusters with Gigabit Ethernet Mesh Connections

    SciTech Connect (OSTI)

    Jie Chen; Chip Watson; Robert Edwards; Weizhen Mao

    2005-04-01

    Multiple copper-based commodity Gigabit Ethernet (GigE) interconnects (adapters) on a single host can lead to Linux clusters with mesh/torus connections without using expensive switches and high speed network interconnects (NICs). However traditional message passing systems based on TCP for GigE will not perform well for this type of clusters because of the overhead of TCP for multiple GigE links. In this paper, we present two os-bypass message passing systems that are based on a modified M-VIA (an implementation of VIA specification) for two production GigE mesh clusters: one is constructed as a 4 x 8 x 8 (256 nodes) torus and has been in production use for a year; the other is constructed as a 6 x 8 x 8 (384 nodes) torus and was deployed recently. One of the message passing systems targets to a specific application domain and is called QMP and the other is an implementation of MPI specification 1.1. The GigE mesh clusters using these two message passing systems achieve about 18.5 {micro}s half-way round trip latency and 400MB/s total bandwidth, which compare reasonably well to systems using specialized high speed adapters in a switched architecture at much lower costs.

  5. TAP - Tools for Adaptive Partitioning v. 0.99 Beta

    Energy Science and Technology Software Center (OSTI)

    2008-11-19

    TAP is a set of tools which are essential for conducting research on adaptive partitioners. The basic premise is that a single partitioner may not be a good choice for adaptive mesh simulations; rather one must match a partitioner (obtained from a partitioning package like Zoltan, ParMetis etc) with the mesh being partitioned. TAP provides the tools that can judge the suitability of a partitioning algorithm to a given mesh.

  6. Refines Efficiency Improvement

    SciTech Connect (OSTI)

    WRI

    2002-05-15

    Refinery processes that convert heavy oils to lighter distillate fuels require heating for distillation, hydrogen addition or carbon rejection (coking). Efficiency is limited by the formation of insoluble carbon-rich coke deposits. Heat exchangers and other refinery units must be shut down for mechanical coke removal, resulting in a significant loss of output and revenue. When a residuum is heated above the temperature at which pyrolysis occurs (340 C, 650 F), there is typically an induction period before coke formation begins (Magaril and Aksenova 1968, Wiehe 1993). To avoid fouling, refiners often stop heating a residuum before coke formation begins, using arbitrary criteria. In many cases, this heating is stopped sooner than need be, resulting in less than maximum product yield. Western Research Institute (WRI) has developed innovative Coking Index concepts (patent pending) which can be used for process control by refiners to heat residua to the threshold, but not beyond the point at which coke formation begins when petroleum residua materials are heated at pyrolysis temperatures (Schabron et al. 2001). The development of this universal predictor solves a long standing problem in petroleum refining. These Coking Indexes have great potential value in improving the efficiency of distillation processes. The Coking Indexes were found to apply to residua in a universal manner, and the theoretical basis for the indexes has been established (Schabron et al. 2001a, 2001b, 2001c). For the first time, a few simple measurements indicates how close undesired coke formation is on the coke formation induction time line. The Coking Indexes can lead to new process controls that can improve refinery distillation efficiency by several percentage points. Petroleum residua consist of an ordered continuum of solvated polar materials usually referred to as asphaltenes dispersed in a lower polarity solvent phase held together by intermediate polarity materials usually referred to as

  7. ITP Petroleum Refining: Energy Bandwidth for Petroleum Refining Processes

    Broader source: Energy.gov [DOE]

    This energy bandwidth analysis provides a realistic estimate of the energy that may be saved in petroleum refining processes by quantifying measures of energy consumption.

  8. RGG: Reactor geometry (and mesh) generator

    SciTech Connect (OSTI)

    Jain, R.; Tautges, T.

    2012-07-01

    The reactor geometry (and mesh) generator RGG takes advantage of information about repeated structures in both assembly and core lattices to simplify the creation of geometry and mesh. It is released as open source software as a part of the MeshKit mesh generation library. The methodology operates in three stages. First, assembly geometry models of various types are generated by a tool called AssyGen. Next, the assembly model or models are meshed by using MeshKit tools or the CUBIT mesh generation tool-kit, optionally based on a journal file output by AssyGen. After one or more assembly model meshes have been constructed, a tool called CoreGen uses a copy/move/merge process to arrange the model meshes into a core model. In this paper, we present the current state of tools and new features in RGG. We also discuss the parallel-enabled CoreGen, which in several cases achieves super-linear speedups since the problems fit in available RAM at higher processor counts. Several RGG applications - 1/6 VHTR model, 1/4 PWR reactor core, and a full-core model for Monju - are reported. (authors)

  9. Mesh infrastructure for coupled multiprocess geophysical simulations

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Garimella, Rao V.; Perkins, William A.; Buksas, Mike W.; Berndt, Markus; Lipnikov, Konstantin; Coon, Ethan; Moulton, John D.; Painter, Scott L.

    2014-01-01

    We have developed a sophisticated mesh infrastructure capability to support large scale multiphysics simulations such as subsurface flow and reactive contaminant transport at storage sites as well as the analysis of the effects of a warming climate on the terrestrial arctic. These simulations involve a wide range of coupled processes including overland flow, subsurface flow, freezing and thawing of ice rich soil, accumulation, redistribution and melting of snow, biogeochemical processes involving plant matter and finally, microtopography evolution due to melting and degradation of ice wedges below the surface. In addition to supporting the usual topological and geometric queries about themore » mesh, the mesh infrastructure adds capabilities such as identifying columnar structures in the mesh, enabling deforming of the mesh subject to constraints and enabling the simultaneous use of meshes of different dimensionality for subsurface and surface processes. The generic mesh interface is capable of using three different open source mesh frameworks (MSTK, MOAB and STKmesh) under the hood allowing the developers to directly compare them and choose one that is best suited for the application's needs. We demonstrate the results of some simulations using these capabilities as well as present a comparison of the performance of the different mesh frameworks.« less

  10. Mesh infrastructure for coupled multiprocess geophysical simulations

    SciTech Connect (OSTI)

    Garimella, Rao V.; Perkins, William A.; Buksas, Mike W.; Berndt, Markus; Lipnikov, Konstantin; Coon, Ethan; Moulton, John D.; Painter, Scott L.

    2014-01-01

    We have developed a sophisticated mesh infrastructure capability to support large scale multiphysics simulations such as subsurface flow and reactive contaminant transport at storage sites as well as the analysis of the effects of a warming climate on the terrestrial arctic. These simulations involve a wide range of coupled processes including overland flow, subsurface flow, freezing and thawing of ice rich soil, accumulation, redistribution and melting of snow, biogeochemical processes involving plant matter and finally, microtopography evolution due to melting and degradation of ice wedges below the surface. In addition to supporting the usual topological and geometric queries about the mesh, the mesh infrastructure adds capabilities such as identifying columnar structures in the mesh, enabling deforming of the mesh subject to constraints and enabling the simultaneous use of meshes of different dimensionality for subsurface and surface processes. The generic mesh interface is capable of using three different open source mesh frameworks (MSTK, MOAB and STKmesh) under the hood allowing the developers to directly compare them and choose one that is best suited for the application's needs. We demonstrate the results of some simulations using these capabilities as well as present a comparison of the performance of the different mesh frameworks.

  11. Method of generating a surface mesh

    DOE Patents [OSTI]

    Shepherd, Jason F.; Benzley, Steven; Grover, Benjamin T.

    2008-03-04

    A method and machine-readable medium provide a technique to generate and modify a quadrilateral finite element surface mesh using dual creation and modification. After generating a dual of a surface (mesh), a predetermined algorithm may be followed to generate and modify a surface mesh of quadrilateral elements. The predetermined algorithm may include the steps of generating two-dimensional cell regions in dual space, determining existing nodes in primal space, generating new nodes in the dual space, and connecting nodes to form the quadrilateral elements (faces) for the generated and modifiable surface mesh.

  12. MOAB : a mesh-oriented database.

    SciTech Connect (OSTI)

    Tautges, Timothy James; Ernst, Corey; Stimpson, Clint; Meyers, Ray J.; Merkley, Karl

    2004-04-01

    A finite element mesh is used to decompose a continuous domain into a discretized representation. The finite element method solves PDEs on this mesh by modeling complex functions as a set of simple basis functions with coefficients at mesh vertices and prescribed continuity between elements. The mesh is one of the fundamental types of data linking the various tools in the FEA process (mesh generation, analysis, visualization, etc.). Thus, the representation of mesh data and operations on those data play a very important role in FEA-based simulations. MOAB is a component for representing and evaluating mesh data. MOAB can store structured and unstructured mesh, consisting of elements in the finite element 'zoo'. The functional interface to MOAB is simple yet powerful, allowing the representation of many types of metadata commonly found on the mesh. MOAB is optimized for efficiency in space and time, based on access to mesh in chunks rather than through individual entities, while also versatile enough to support individual entity access. The MOAB data model consists of a mesh interface instance, mesh entities (vertices and elements), sets, and tags. Entities are addressed through handles rather than pointers, to allow the underlying representation of an entity to change without changing the handle to that entity. Sets are arbitrary groupings of mesh entities and other sets. Sets also support parent/child relationships as a relation distinct from sets containing other sets. The directed-graph provided by set parent/child relationships is useful for modeling topological relations from a geometric model or other metadata. Tags are named data which can be assigned to the mesh as a whole, individual entities, or sets. Tags are a mechanism for attaching data to individual entities and sets are a mechanism for describing relations between entities; the combination of these two mechanisms is a powerful yet simple interface for representing metadata or application

  13. Petroleum Refining | Department of Energy

    Broader source: Energy.gov (indexed) [DOE]

    Energy Efficiency Roadmap for Petroleum Refineries in California (2004) ... Photo credit: iStockphoto.com15906809 Petroleum Refining The United States is the largest, most ...

  14. Massively parallel mesh generation for physics codes

    SciTech Connect (OSTI)

    Hardin, D.D.

    1996-06-01

    Massively parallel processors (MPPs) will soon enable realistic 3-D physical modeling of complex objects and systems. Work is planned or presently underway to port many of LLNL`s physical modeling codes to MPPs. LLNL`s DSI3D electromagnetics code already can solve 40+ million zone problems on the 256 processor Meiko. However, the author lacks the software necessary to generate and manipulate the large meshes needed to model many complicated 3-D geometries. State-of-the-art commercial mesh generators run on workstations and have a practical limit of several hundred thousand elements. In the foreseeable future MPPs will solve problems with a billion mesh elements. The objective of the Parallel Mesh Generation (PMESH) Project is to develop a unique mesh generation system that can construct large 3-D meshes (up to a billion elements) on MPPs. Such a capability will remove a critical roadblock to unleashing the power of MPPs for physical analysis and will put LLNL at the forefront of mesh generation technology. PMESH will ``front-end`` a variety of LLNL 3-D physics codes, including those in the areas of electromagnetics, structural mechanics, thermal analysis, and hydrodynamics. The DSI3D and DYNA3D codes are already running on MPPs. The primary goal of the PMESH project is to provide the robust generation of large meshes for complicated 3-D geometries through the appropriate distribution of the generation task between the user`s workstation and the MPP. Secondary goals are to support the unique features of LLNL physics codes (e.g., unusual elements) and to minimize the user effort required to generate different meshes for the same geometry. PMESH`s capabilities are essential because mesh generation is presently a major limiting factor in simulating larger and more complex 3-D geometries. PMESH will significantly enhance LLNL`s capabilities in physical simulation by advancing the state-of-the-art in large mesh generation by 2 to 3 orders of magnitude.

  15. 7th International Meshing Roundtable '98

    SciTech Connect (OSTI)

    Eldred, T.J.

    1998-10-01

    The goal of the 7th International Meshing Roundtable is to bring together researchers and developers from industry, academia, and government labs in a stimulating, open environment for the exchange of technical information related to the meshing process. In the past, the Roundtable has enjoyed significant participation from each of these groups from a wide variety of countries.

  16. Crystal structure refinement with SHELXL

    SciTech Connect (OSTI)

    Sheldrick, George M., E-mail: gsheldr@shelx.uni-ac.gwdg.de [Department of Structural Chemistry, Georg-August Universitt Gttingen, Tammannstrae 4, Gttingen 37077 (Germany)

    2015-01-01

    New features added to the refinement program SHELXL since 2008 are described and explained. The improvements in the crystal structure refinement program SHELXL have been closely coupled with the development and increasing importance of the CIF (Crystallographic Information Framework) format for validating and archiving crystal structures. An important simplification is that now only one file in CIF format (for convenience, referred to simply as a CIF) containing embedded reflection data and SHELXL instructions is needed for a complete structure archive; the program SHREDCIF can be used to extract the .hkl and .ins files required for further refinement with SHELXL. Recent developments in SHELXL facilitate refinement against neutron diffraction data, the treatment of H atoms, the determination of absolute structure, the input of partial structure factors and the refinement of twinned and disordered structures. SHELXL is available free to academics for the Windows, Linux and Mac OS X operating systems, and is particularly suitable for multiple-core processors.

  17. Feature recognition applications in mesh generation

    SciTech Connect (OSTI)

    Tautges, T.J.; Liu, S.S.; Lu, Y.; Kraftcheck, J.; Gadh, R.

    1997-06-01

    The use of feature recognition as part of an overall decomposition-based hexahedral meshing approach is described in this paper. The meshing approach consists of feature recognition, using a c-loop or hybrid c-loop method, and the use of cutting surfaces to decompose the solid model. These steps are part of an iterative process, which proceeds either until no more features can be recognized or until the model has been completely decomposed into meshable sub-volumes. This method can greatly reduce the time required to generate an all-hexahedral mesh, either through the use of more efficient meshing algorithms on more of the geometry or by reducing the amount of manual decomposition required to mesh a volume.

  18. Update on Development of Mesh Generation Algorithms in MeshKit

    SciTech Connect (OSTI)

    Jain, Rajeev; Vanderzee, Evan; Mahadevan, Vijay

    2015-09-30

    MeshKit uses a graph-based design for coding all its meshing algorithms, which includes the Reactor Geometry (and mesh) Generation (RGG) algorithms. This report highlights the developmental updates of all the algorithms, results and future work. Parallel versions of algorithms, documentation and performance results are reported. RGG GUI design was updated to incorporate new features requested by the users; boundary layer generation and parallel RGG support were added to the GUI. Key contributions to the release, upgrade and maintenance of other SIGMA1 libraries (CGM and MOAB) were made. Several fundamental meshing algorithms for creating a robust parallel meshing pipeline in MeshKit are under development. Results and current status of automated, open-source and high quality nuclear reactor assembly mesh generation algorithms such as trimesher, quadmesher, interval matching and multi-sweeper are reported.

  19. Method of modifying a volume mesh using sheet extraction

    DOE Patents [OSTI]

    Borden, Michael J.; Shepherd, Jason F.

    2007-02-20

    A method and machine-readable medium provide a technique to modify a hexahedral finite element volume mesh using dual generation and sheet extraction. After generating a dual of a volume stack (mesh), a predetermined algorithm may be followed to modify the volume mesh of hexahedral elements. The predetermined algorithm may include the steps of determining a sheet of hexahedral mesh elements, generating nodes for merging, and merging the nodes to delete the sheet of hexahedral mesh elements and modify the volume mesh.

  20. Introducing the target-matrix paradigm for mesh optimization...

    Office of Scientific and Technical Information (OSTI)

    are automatically constructed to enable definitions of application-specific mesh quality. ... that address mesh quality issues such as those shown in the examples of the present paper. ...

  1. Label-invariant Mesh Quality Metrics. (Conference) | SciTech...

    Office of Scientific and Technical Information (OSTI)

    Label-invariant Mesh Quality Metrics. Citation Details In-Document Search Title: Label-invariant Mesh Quality Metrics. Abstract not provided. Authors: Knupp, Patrick Publication ...

  2. Electronic structure of superconductivity refined

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Electronic structure of superconductivity refined Electronic structure of superconductivity refined A team of physicists propose a new model that expands on a little understood aspect of the electronic structure in high-temperature superconductors. July 10, 2008 Los Alamos National Laboratory sits on top of a once-remote mesa in northern New Mexico with the Jemez mountains as a backdrop to research and innovation covering multi-disciplines from bioscience, sustainable energy sources, to plasma

  3. ,"U.S. Refiner Petroleum Product Prices"

    U.S. Energy Information Administration (EIA) Indexed Site

    ...EMAEPPRLPTGNUSDPG","EMAEPPRHPTGNUSDPG" "Date","U.S. Total Gasoline Retail Sales by Refiners (Dollars per Gallon)","U.S. Aviation Gasoline Retail Sales by Refiners (Dollars ...

  4. Renewable Energy Financial Instruments Guidance Tool (REFINe...

    Open Energy Info (EERE)

    Guidance Tool (REFINe) Jump to: navigation, search Tool Summary LAUNCH TOOL Name: Renewable Energy Financial Instruments Guidance Tool (REFINe) Focus Area: Renewable Energy...

  5. MECS 2006 - Petroleum Refining | Department of Energy

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Manufacturing Energy and Carbon Footprint Petroleum Refining (123.98 KB) More Documents & Publications Petroleum Refining (2010 MECS) MECS 2006 - Alumina and Aluminum MECS 2006 - ...

  6. ,"U.S. Refiner Petroleum Product Prices"

    U.S. Energy Information Administration (EIA) Indexed Site

    ...AEPPRLPTGNUSDPG","EMAEPPRHPTGNUSDPG" "Date","U.S. Total Gasoline Retail Sales by Refiners (Dollars per Gallon)","U.S. Aviation Gasoline Retail Sales by Refiners (Dollars ...

  7. Mesh-free data transfer algorithms for partitioned multiphysics problems: Conservation, accuracy, and parallelism

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Slattery, Stuart R.

    2015-12-02

    In this study we analyze and extend mesh-free algorithms for three-dimensional data transfer problems in partitioned multiphysics simulations. We first provide a direct comparison between a mesh-based weighted residual method using the common-refinement scheme and two mesh-free algorithms leveraging compactly supported radial basis functions: one using a spline interpolation and one using a moving least square reconstruction. Through the comparison we assess both the conservation and accuracy of the data transfer obtained from each of the methods. We do so for a varying set of geometries with and without curvature and sharp features and for functions with and without smoothnessmore » and with varying gradients. Our results show that the mesh-based and mesh-free algorithms are complementary with cases where each was demonstrated to perform better than the other. We then focus on the mesh-free methods by developing a set of algorithms to parallelize them based on sparse linear algebra techniques. This includes a discussion of fast parallel radius searching in point clouds and restructuring the interpolation algorithms to leverage data structures and linear algebra services designed for large distributed computing environments. The scalability of our new algorithms is demonstrated on a leadership class computing facility using a set of basic scaling studies. Finally, these scaling studies show that for problems with reasonable load balance, our new algorithms for both spline interpolation and moving least square reconstruction demonstrate both strong and weak scalability using more than 100,000 MPI processes with billions of degrees of freedom in the data transfer operation.« less

  8. Mesh-free data transfer algorithms for partitioned multiphysics problems: Conservation, accuracy, and parallelism

    SciTech Connect (OSTI)

    Slattery, Stuart R.

    2015-12-02

    In this study we analyze and extend mesh-free algorithms for three-dimensional data transfer problems in partitioned multiphysics simulations. We first provide a direct comparison between a mesh-based weighted residual method using the common-refinement scheme and two mesh-free algorithms leveraging compactly supported radial basis functions: one using a spline interpolation and one using a moving least square reconstruction. Through the comparison we assess both the conservation and accuracy of the data transfer obtained from each of the methods. We do so for a varying set of geometries with and without curvature and sharp features and for functions with and without smoothness and with varying gradients. Our results show that the mesh-based and mesh-free algorithms are complementary with cases where each was demonstrated to perform better than the other. We then focus on the mesh-free methods by developing a set of algorithms to parallelize them based on sparse linear algebra techniques. This includes a discussion of fast parallel radius searching in point clouds and restructuring the interpolation algorithms to leverage data structures and linear algebra services designed for large distributed computing environments. The scalability of our new algorithms is demonstrated on a leadership class computing facility using a set of basic scaling studies. Finally, these scaling studies show that for problems with reasonable load balance, our new algorithms for both spline interpolation and moving least square reconstruction demonstrate both strong and weak scalability using more than 100,000 MPI processes with billions of degrees of freedom in the data transfer operation.

  9. NERSC Selects 20 NESAP Code Teams

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    particle-mesh algorithms, adaptive mesh refinement techniques, molecular dynamics, eigenvalue problems for complex molecular systems and ab initio computational chemistry methods. ...

  10. Parallel octree-based hexahedral mesh generation for eulerian to lagrangian conversion.

    SciTech Connect (OSTI)

    Staten, Matthew L.; Owen, Steven James

    2010-09-01

    Computational simulation must often be performed on domains where materials are represented as scalar quantities or volume fractions at cell centers of an octree-based grid. Common examples include bio-medical, geotechnical or shock physics calculations where interface boundaries are represented only as discrete statistical approximations. In this work, we introduce new methods for generating Lagrangian computational meshes from Eulerian-based data. We focus specifically on shock physics problems that are relevant to ASC codes such as CTH and Alegra. New procedures for generating all-hexahedral finite element meshes from volume fraction data are introduced. A new primal-contouring approach is introduced for defining a geometric domain. New methods for refinement, node smoothing, resolving non-manifold conditions and defining geometry are also introduced as well as an extension of the algorithm to handle tetrahedral meshes. We also describe new scalable MPI-based implementations of these procedures. We describe a new software module, Sculptor, which has been developed for use as an embedded component of CTH. We also describe its interface and its use within the mesh generation code, CUBIT. Several examples are shown to illustrate the capabilities of Sculptor.

  11. Detecting Translation Errors in CAD Surfaces and Preparing Geometries for Mesh Generation

    SciTech Connect (OSTI)

    Petersson, N Anders; Chand, K K

    2001-08-27

    The authors have developed tools for the efficient preparation of CAD geometries for mesh generation. Geometries are read from IGES files and then maintained in a boundary-representation consisting of a patchwork of trimmed and untrimmed surfaces. Gross errors in the geometry can be identified and removed automatically while a user interface is provided for manipulating the geometry (such as correcting invalid trimming curves or removing unwanted details). Modifying the geometry by adding or deleting surfaces and/or sectioning it by arbitrary planes (e.g. symmetry planes) is also supported. These tools are used for robust and accurate geometry models for initial mesh generation and will be applied to in situ mesh generation requirements of moving and adaptive grid simulations.

  12. Mesh Generation for SHARP: Unprecedented Complexity | Department of Energy

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Mesh Generation for SHARP: Unprecedented Complexity Mesh Generation for SHARP: Unprecedented Complexity January 29, 2013 - 1:36pm Addthis SHARP Supporting Elements During this quarter, the framework team was involved in two primary efforts, mesh generation and implementation of a MOAB-based coupled multi-physics simulation. For mesh generation, finishing touches were put on three major, high-complexity hexahedral meshes, and support was provided for their use in various simulations: MATiS-H, an

  13. GRAIN REFINEMENT OF URANIUM BILLETS

    DOE Patents [OSTI]

    Lewis, L.

    1964-02-25

    A method of refining the grain structure of massive uranium billets without resort to forging is described. The method consists in the steps of beta- quenching the billets, annealing the quenched billets in the upper alpha temperature range, and extrusion upset of the billets to an extent sufficient to increase the cross sectional area by at least 5 per cent. (AEC)

  14. Method for refining contaminated iridium

    DOE Patents [OSTI]

    Heshmatpour, B.; Heestand, R.L.

    1982-08-31

    Contaminated iridium is refined by alloying it with an alloying agent selected from the group consisting of manganese and an alloy of manganese and copper, and then dissolving the alloying agent from the formed alloy to provide a purified iridium powder.

  15. ITP Petroleum Refining: Energy Bandwidth for Petroleum Refining Processes

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Bandwidth for Petroleum Refining Processes Prepared by Energetics Incorporated for the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Industrial Technologies Program October 2006 Foreword The Industrial Technologies Program (ITP) is a research and development (R&D) program within the U.S. Department of Energy, Office of Energy Efficiency and Renewable Energy (EERE). This program works in collaboration with U.S. industry to improve industrial energy efficiency and

  16. On combining Laplacian and optimization-based mesh smoothing techniques

    SciTech Connect (OSTI)

    Freitag, L.A.

    1997-07-01

    Local mesh smoothing algorithms have been shown to be effective in repairing distorted elements in automatically generated meshes. The simplest such algorithm is Laplacian smoothing, which moves grid points to the geometric center of incident vertices. Unfortunately, this method operates heuristically and can create invalid meshes or elements of worse quality than those contained in the original mesh. In contrast, optimization-based methods are designed to maximize some measure of mesh quality and are very effective at eliminating extremal angles in the mesh. These improvements come at a higher computational cost, however. In this article the author proposes three smoothing techniques that combine a smart variant of Laplacian smoothing with an optimization-based approach. Several numerical experiments are performed that compare the mesh quality and computational cost for each of the methods in two and three dimensions. The author finds that the combined approaches are very cost effective and yield high-quality meshes.

  17. Reformulated Gasoline Market Affected Refiners Differently, 1995

    Reports and Publications (EIA)

    1996-01-01

    This article focuses on the costs of producing reformulated gasoline (RFG) as experienced by different types of refiners and on how these refiners fared this past summer, given the prices for RFG at the refinery gate.

  18. The Uranium Processing Facility (UPF) Finite Element Meshing Discussion

    Office of Energy Efficiency and Renewable Energy (EERE)

    The Uranium Processing Facility (UPF) Finite Element Meshing Discussion Loring Wyllie Arne Halterman Degenkolb Engineers, San Francisco

  19. Cubit Mesh Generation Toolkit V11.1

    Energy Science and Technology Software Center (OSTI)

    2009-03-25

    CUBIT prepares models to be used in computer-based simulation of real-world events. CUBIT is a full-featured software toolkit for robust generation of two- and three-dimensional finite element meshes (grids) and geometry preparation. Its main goal is to reduce the time to generate meshes, particularly large hex meshes of complicated, interlocking assemblies.

  20. Firing of pulverized solvent refined coal

    DOE Patents [OSTI]

    Derbidge, T. Craig; Mulholland, James A.; Foster, Edward P.

    1986-01-01

    An air-purged burner for the firing of pulverized solvent refined coal is constructed and operated such that the solvent refined coal can be fired without the coking thereof on the burner components. The air-purged burner is designed for the firing of pulverized solvent refined coal in a tangentially fired boiler.

  1. Automatic meshing of curved three-dimensional domains: Curving finite elements and curvature-based mesh control

    SciTech Connect (OSTI)

    Shephard, M.S.; Dey, S.; Georges, M.K.

    1995-12-31

    Specific issues associated with the automatic generation of finite element meshes for curved geometric domains axe considered. A review of the definition of when a triangulation is a valid mesh, a geometric triangulation, for curved geometric domains is given. Consideration is then given to the additional operations necessary to maintain the validity of a mesh when curved finite elements are employed. A procedure to control the mesh gradations based on the curvature of the geometric model faces is also given.

  2. INGRID; 3-D Mesh Generation Nonlinear Systems

    SciTech Connect (OSTI)

    Stillman, D.W.; Rainsberger, R.

    1985-07-01

    INGRID is a general-purpose, three-dimensional mesh generator developed for use with finite element, nonlinear, structural dynamics codes. INGRID generates the large and complex input data files for DYNA3D, NIKE3D, FACET, and TOPAZ3D. One of the greatest advantages of INGRID is that virtually any shape can be described without resorting to wedge elements, tetrahedrons, triangular elements or highly distorted quadrilateral or hexahedral elements. Other capabilities available are in the areas of geometry and graphics. Exact surface equations and surface intersections considerably improve the ability to deal with accurate models, and a hidden line graphics algorithm is included which is efficient on the most complicated meshes. The primary new capability is associated with the boundary conditions, loads, and material properties required by nonlinear mechanics programs. Commands have been designed for each case to minimize user effort. This is particularly important since special processing is almost always required for each load or boundary condition.

  3. INGRID. 3-D Mesh Generation Nonlinear Systems

    SciTech Connect (OSTI)

    Christon, M.A.; Dovey, D.; Stillman, D.W.; Hallquist, J.O.; Rainsberger, R.B.

    1992-09-01

    INGRID is a general-purpose, three-dimensional mesh generator developed for use with finite element, nonlinear, structural dynamics codes. INGRID generates the large and complex input data files for DYNA3D, NIKE3D, FACET, and TOPAZ3D. One of the greatest advantages of INGRID is that virtually any shape can be described without resorting to wedge elements, tetrahedrons, triangular elements or highly distorted quadrilateral or hexahedral elements. Other capabilities available are in the areas of geometry and graphics. Exact surface equations and surface intersections considerably improve the ability to deal with accurate models, and a hidden line graphics algorithm is included which is efficient on the most complicated meshes. The primary new capability is associated with the boundary conditions, loads, and material properties required by nonlinear mechanics programs. Commands have been designed for each case to minimize user effort. This is particularly important since special processing is almost always required for each load or boundary condition.

  4. INGRID. 3-D Mesh Generation Nonlinear Systems

    SciTech Connect (OSTI)

    Stillman, D.W.; Rainsberger, R.

    1985-07-01

    INGRID is a general-purpose, three-dimensional mesh generator developed for use with finite element, nonlinear, structural dynamics codes. INGRID generates the large and complex input data files for DYNA3D, NIKE3D, FACET, and TOPAZ3D (ESTSC. One of the greatest advantages of INGRID is that virtually any shape can be described without resorting to wedge elements, tetrahedrons, triangular elements or highly distorted quadrilateral or hexahedral elements. Other capabilities available are in the areas of geometry and graphics. Exact surface equations and surface intersections considerably improve the ability to deal with accurate models, and a hidden line graphics algorithm is included which is efficient on the most complicated meshes. The primary new capability is associated with the boundary conditions, loads, and material properties required by nonlinear mechanics programs. Commands have been designed for each case to minimize user effort. This is particularly important since special processing is almost always required for each load or boundary condition.

  5. INGRID. 3-D Mesh Generation Nonlinear Systems

    SciTech Connect (OSTI)

    Stillman, D.W.; Rainsberger, R.

    1985-07-01

    INGRID is a general-purpose, three-dimensional mesh generator developed for use with finite element, nonlinear, structural dynamics codes. INGRID generates the large and complex input data files for DYNA3D, NIKE3D, FACET, and TOPAZ3D. One of the greatest advantages of INGRID is that virtually any shape can be described without resorting to wedge elements, tetrahedrons, triangular elements or highly distorted quadrilateral or hexahedral elements. Other capabilities available are in the areas of geometry and graphics. Exact surface equations and surface intersections considerably improve the ability to deal with accurate models, and a hidden line graphics algorithm is included which is efficient on the most complicated meshes. The primary new capability is associated with the boundary conditions, loads, and material properties required by nonlinear mechanics programs. Commands have been designed for each case to minimize user effort. This is particularly important since special processing is almost always required for each load or boundary condition.

  6. 3-D Mesh Generation Nonlinear Systems

    SciTech Connect (OSTI)

    1994-04-07

    INGRID is a general-purpose, three-dimensional mesh generator developed for use with finite element, nonlinear, structural dynamics codes. INGRID generates the large and complex input data files for DYNA3D, NIKE3D, FACET, and TOPAZ3D. One of the greatest advantages of INGRID is that virtually any shape can be described without resorting to wedge elements, tetrahedrons, triangular elements or highly distorted quadrilateral or hexahedral elements. Other capabilities available are in the areas of geometry and graphics. Exact surface equations and surface intersections considerably improve the ability to deal with accurate models, and a hidden line graphics algorithm is included which is efficient on the most complicated meshes. The primary new capability is associated with the boundary conditions, loads, and material properties required by nonlinear mechanics programs. Commands have been designed for each case to minimize user effort. This is particularly important since special processing is almost always required for each load or boundary condition.

  7. Vertical-Axis Wind Turbine Mesh Generator

    Energy Science and Technology Software Center (OSTI)

    2014-01-24

    VAWTGen is a mesh generator for creating a finite element beam mesh of arbitrary vertical-axis wind turbines (VAWT). The software accepts input files specifying tower and blade structural and aerodynamic descriptions and constructs a VAWT using a minimal set of inputs. VAWTs with an arbitrary number of blades can be constructed with or without a central tower. Strut connections between the tower and blades can be specified in an arbitrary manner. The software also facilitatesmore » specifying arbitrary joints between structural components and concentrated structural tenns (mass and stiffness). The output files which describe the VAWT configuration are intended to be used with the Offshore Wind ENergy Simulation (OWENS) Toolkit software for structural dynamics analysis of VAWTs. Furthermore, VAWTGen is useful for visualizing output from the OWENS analysis software.« less

  8. Vertical-Axis Wind Turbine Mesh Generator

    SciTech Connect (OSTI)

    2014-01-24

    VAWTGen is a mesh generator for creating a finite element beam mesh of arbitrary vertical-axis wind turbines (VAWT). The software accepts input files specifying tower and blade structural and aerodynamic descriptions and constructs a VAWT using a minimal set of inputs. VAWTs with an arbitrary number of blades can be constructed with or without a central tower. Strut connections between the tower and blades can be specified in an arbitrary manner. The software also facilitates specifying arbitrary joints between structural components and concentrated structural tenns (mass and stiffness). The output files which describe the VAWT configuration are intended to be used with the Offshore Wind ENergy Simulation (OWENS) Toolkit software for structural dynamics analysis of VAWTs. Furthermore, VAWTGen is useful for visualizing output from the OWENS analysis software.

  9. Diffusive mesh relaxation in ALE finite element numerical simulations

    SciTech Connect (OSTI)

    Dube, E.I.

    1996-06-01

    The theory for a diffusive mesh relaxation algorithm is developed for use in three-dimensional Arbitary Lagrange/Eulerian (ALE) finite element simulation techniques. This mesh relaxer is derived by a variational principle for an unstructured 3D grid using finite elements, and incorporates hourglass controls in the numerical implementation. The diffusive coefficients are based on the geometric properties of the existing mesh, and are chosen so as to allow for a smooth grid that retains the general shape of the original mesh. The diffusive mesh relaxation algorithm is then applied to an ALE code system, and results from several test cases are discussed.

  10. Wireless experiments on a Motorola mesh testbed.

    SciTech Connect (OSTI)

    Riblett, Loren E., Jr.; Wiseman, James M.; Witzke, Edward L.

    2010-06-01

    Motomesh is a Motorola product that performs mesh networking at both the client and access point levels and allows broadband mobile data connections with or between clients moving at vehicular speeds. Sandia National aboratories has extensive experience with this product and its predecessors in infrastructure-less mobile environments. This report documents experiments, which characterize certain aspects of how the Motomesh network performs when obile units are added to a fixed network infrastructure.

  11. ITP Petroleum Refining: Energy and Environmental Profile of the...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    and Environmental Profile of the U.S. Petroleum Refining Industry (November 2007) ITP Petroleum Refining: Energy and Environmental Profile of the U.S. Petroleum Refining Industry ...

  12. An h-adaptive finite element method for turbulent heat transfer

    SciTech Connect (OSTI)

    Carriington, David B [Los Alamos National Laboratory

    2009-01-01

    A two-equation turbulence closure model (k-{omega}) using an h-adaptive grid technique and finite element method (FEM) has been developed to simulate low Mach flow and heat transfer. These flows are applicable to many flows in engineering and environmental sciences. Of particular interest in the engineering modeling areas are: combustion, solidification, and heat exchanger design. Flows for indoor air quality modeling and atmospheric pollution transport are typical types of environmental flows modeled with this method. The numerical method is based on a hybrid finite element model using an equal-order projection process. The model includes thermal and species transport, localized mesh refinement (h-adaptive) and Petrov-Galerkin weighting for the stabilizing the advection. This work develops the continuum model of a two-equation turbulence closure method. The fractional step solution method is stated along with the h-adaptive grid method (Carrington and Pepper, 2002). Solutions are presented for 2d flow over a backward-facing step.

  13. Data-Parallel Mesh Connected Components Labeling and Analysis

    SciTech Connect (OSTI)

    Harrison, Cyrus; Childs, Hank; Gaither, Kelly

    2011-04-10

    We present a data-parallel algorithm for identifying and labeling the connected sub-meshes within a domain-decomposed 3D mesh. The identification task is challenging in a distributed-memory parallel setting because connectivity is transitive and the cells composing each sub-mesh may span many or all processors. Our algorithm employs a multi-stage application of the Union-find algorithm and a spatial partitioning scheme to efficiently merge information across processors and produce a global labeling of connected sub-meshes. Marking each vertex with its corresponding sub-mesh label allows us to isolate mesh features based on topology, enabling new analysis capabilities. We briefly discuss two specific applications of the algorithm and present results from a weak scaling study. We demonstrate the algorithm at concurrency levels up to 2197 cores and analyze meshes containing up to 68 billion cells.

  14. Improving the accuracy of macromolecular structure refinement

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    is challenging at low resolution. We compared refinement methods using synchrotron diffraction data of photosystem I at 7.4 resolution, starting from different initial models...

  15. Refiners Increasingly Employing Catalyst Regeneration as Alternative...

    Open Energy Info (EERE)

    million by the end of 2019. Refiners Benefit from Catalyst Regeneration Technology via Price Reductions and Lower Maintenance Costs The catalyst regeneration technology is the...

  16. Right rock: Finding/refining customer expectations

    SciTech Connect (OSTI)

    Ashby, R.

    1997-11-01

    This report contains viewgraphs on methods of finding customer expectations of software and then refining the documentation and software to fit these expectations.

  17. MeSh ToolKit v1.2

    Energy Science and Technology Software Center (OSTI)

    2004-05-15

    MSTK or Mesh Toolkit is a mesh framework that allows users to represent, manipulate and query unstructured 3D arbitrary topology meshes in a general manner without the need to code their own data structures. MSTK is a flexible framework in that is allows (or will eventually allow) a wide variety of underlying representations for the mesh while maintaining a common interface. It will allow users to choose from different mesh representations either at initialization ormore » during the program execution so that the optimal data structures are used for the particular algorithm. The interaction of users and applications with MSTK is through a functional interface that acts as through the mesh always contains vertices, edges, faces and regions and maintains connectivity between all these entities.« less

  18. 3-D Mesh Generation Nonlinear Systems

    Energy Science and Technology Software Center (OSTI)

    1994-04-07

    INGRID is a general-purpose, three-dimensional mesh generator developed for use with finite element, nonlinear, structural dynamics codes. INGRID generates the large and complex input data files for DYNA3D, NIKE3D, FACET, and TOPAZ3D. One of the greatest advantages of INGRID is that virtually any shape can be described without resorting to wedge elements, tetrahedrons, triangular elements or highly distorted quadrilateral or hexahedral elements. Other capabilities available are in the areas of geometry and graphics. Exact surfacemore » equations and surface intersections considerably improve the ability to deal with accurate models, and a hidden line graphics algorithm is included which is efficient on the most complicated meshes. The primary new capability is associated with the boundary conditions, loads, and material properties required by nonlinear mechanics programs. Commands have been designed for each case to minimize user effort. This is particularly important since special processing is almost always required for each load or boundary condition.« less

  19. Introducing the Target-Matrix Paradigm for Mesh Optimization...

    Office of Scientific and Technical Information (OSTI)

    Title: Introducing the Target-Matrix Paradigm for Mesh Optimization via Node-movement. Abstract not provided. Authors: Knupp, Patrick Michael Publication Date: 2010-10-01 OSTI ...

  20. Incorporating the min-max mesh optimization method within the...

    Office of Scientific and Technical Information (OSTI)

    Incorporating the min-max mesh optimization method within the Target-Matrix Paradigm. ... Country of Publication: United States Language: English Word Cloud More Like This Full ...

  1. Incorporating the min-max mesh optimization method within the...

    Office of Scientific and Technical Information (OSTI)

    Title: Incorporating the min-max mesh optimization method within the Target-Matrix ... Country of Publication: United States Language: English Word Cloud More Like This Full ...

  2. Hexahedral finite element mesh coarsening using pillowing technique

    DOE Patents [OSTI]

    Staten, Matthew L.; Woodbury, Adam C.; Benzley, Steven E.; Shepherd, Jason F.

    2012-06-05

    A techniques for coarsening a hexahedral mesh is described. The technique includes identifying a coarsening region within a hexahedral mesh to be coarsened. A boundary sheet of hexahedral elements is inserted into the hexahedral mesh around the coarsening region. A column of hexahedral elements is identified within the boundary sheet. The column of hexahedral elements is collapsed to create an extraction sheet of hexahedral elements contained within the coarsening region. Then, the extraction sheet of hexahedral elements is extracted to coarsen the hexahedral mesh.

  3. A Lagrangian-Eulerian method with zoomable hidden fine-mesh approach to solving advection-dispersion equations

    SciTech Connect (OSTI)

    Yeh, G.T. )

    1990-06-01

    A Lagrangian-Eulerian method with zoomable hidden fine-mesh approach (LEZOOM), that can be adapted with either finite element or finite difference methods, is used to solve the advection-dispersion equation. The approach is based on automatic adaptation of zooming a hidden fine mesh in regions where the sharp front is located. Application of LEZOOM to four bench mark problems indicates that it can handle the advection-dispersion/diffusion problems with mesh Peclet numbers ranged from 0 to {infinity} and with mesh Courant numbers well in excess of 1. Difficulties that can be resolved with LEZOOM include numerical dispersion, oscillations, the clipping of peaks, and the effect of grid orientation. Nonuniform grid as well as spatial temporally variable flow pose no problems with LEZOOM. Both initial and boundary value problems can be solved accurately with LEZOOM. It is shown that although the mixed Lagrangian-Eulerian (LE) approach (LEZOOM without zooming) also produces excessive numerical dispersion as the upstream finite element (UFE) method, the LE approach is superior to the UFE method.

  4. Pneumatic conveying of pulverized solvent refined coal

    DOE Patents [OSTI]

    Lennon, Dennis R.

    1984-11-06

    A method for pneumatically conveying solvent refined coal to a burner under conditions of dilute phase pneumatic flow so as to prevent saltation of the solvent refined coal in the transport line by maintaining the transport fluid velocity above approximately 95 ft/sec.

  5. MeshVoro: A Three-Dimensional Voronoi Mesh Building Tool for the TOUGH Family of Codes

    SciTech Connect (OSTI)

    Freeman, C. M.; Boyle, K. L.; Reagan, M.; Johnson, J.; Rycroft, C.; Moridis, G. J.

    2013-09-30

    Few tools exist for creating and visualizing complex three-dimensional simulation meshes, and these have limitations that restrict their application to particular geometries and circumstances. Mesh generation needs to trend toward ever more general applications. To that end, we have developed MeshVoro, a tool that is based on the Voro (Rycroft 2009) library and is capable of generating complex threedimensional Voronoi tessellation-based (unstructured) meshes for the solution of problems of flow and transport in subsurface geologic media that are addressed by the TOUGH (Pruess et al. 1999) family of codes. MeshVoro, which includes built-in data visualization routines, is a particularly useful tool because it extends the applicability of the TOUGH family of codes by enabling the scientifically robust and relatively easy discretization of systems with challenging 3D geometries. We describe several applications of MeshVoro. We illustrate the ability of the tool to straightforwardly transform a complex geological grid into a simulation mesh that conforms to the specifications of the TOUGH family of codes. We demonstrate how MeshVoro can describe complex system geometries with a relatively small number of grid blocks, and we construct meshes for geometries that would have been practically intractable with a standard Cartesian grid approach. We also discuss the limitations and appropriate applications of this new technology.

  6. North Dakota Refining Capacity Study

    SciTech Connect (OSTI)

    Dennis Hill; Kurt Swenson; Carl Tuura; Jim Simon; Robert Vermette; Gilberto Marcha; Steve Kelly; David Wells; Ed Palmer; Kuo Yu; Tram Nguyen; Juliam Migliavacca

    2011-01-05

    According to a 2008 report issued by the United States Geological Survey, North Dakota and Montana have an estimated 3.0 to 4.3 billion barrels of undiscovered, technically recoverable oil in an area known as the Bakken Formation. With the size and remoteness of the discovery, the question became 'can a business case be made for increasing refining capacity in North Dakota?' And, if so what is the impact to existing players in the region. To answer the question, a study committee comprised of leaders in the region's petroleum industry were brought together to define the scope of the study, hire a consulting firm and oversee the study. The study committee met frequently to provide input on the findings and modify the course of the study, as needed. The study concluded that the Petroleum Area Defense District II (PADD II) has an oversupply of gasoline. With that in mind, a niche market, naphtha, was identified. Naphtha is used as a diluent used for pipelining the bitumen (heavy crude) from Canada to crude markets. The study predicted there will continue to be an increase in the demand for naphtha through 2030. The study estimated the optimal configuration for the refinery at 34,000 barrels per day (BPD) producing 15,000 BPD of naphtha and a 52 percent refinery charge for jet and diesel yield. The financial modeling assumed the sponsor of a refinery would invest its own capital to pay for construction costs. With this assumption, the internal rate of return is 9.2 percent which is not sufficient to attract traditional investment given the risk factor of the project. With that in mind, those interested in pursuing this niche market will need to identify incentives to improve the rate of return.

  7. Volume Decomposition and Feature Recognition for Hexahedral Mesh Generation

    SciTech Connect (OSTI)

    GADH,RAJIT; LU,YONG; TAUTGES,TIMOTHY J.

    1999-09-27

    Considerable progress has been made on automatic hexahedral mesh generation in recent years. Several automatic meshing algorithms have proven to be very reliable on certain classes of geometry. While it is always worth pursuing general algorithms viable on more general geometry, a combination of the well-established algorithms is ready to take on classes of complicated geometry. By partitioning the entire geometry into meshable pieces matched with appropriate meshing algorithm the original geometry becomes meshable and may achieve better mesh quality. Each meshable portion is recognized as a meshing feature. This paper, which is a part of the feature based meshing methodology, presents the work on shape recognition and volume decomposition to automatically decompose a CAD model into meshable volumes. There are four phases in this approach: (1) Feature Determination to extinct decomposition features, (2) Cutting Surfaces Generation to form the ''tailored'' cutting surfaces, (3) Body Decomposition to get the imprinted volumes; and (4) Meshing Algorithm Assignment to match volumes decomposed with appropriate meshing algorithms. The feature determination procedure is based on the CLoop feature recognition algorithm that is extended to be more general. Results are demonstrated over several parts with complicated topology and geometry.

  8. Firing of pulverized solvent refined coal

    DOE Patents [OSTI]

    Lennon, Dennis R.; Snedden, Richard B.; Foster, Edward P.; Bellas, George T.

    1990-05-15

    A burner for the firing of pulverized solvent refined coal is constructed and operated such that the solvent refined coal can be fired successfully without any performance limitations and without the coking of the solvent refined coal on the burner components. The burner is provided with a tangential inlet of primary air and pulverized fuel, a vaned diffusion swirler for the mixture of primary air and fuel, a center water-cooled conical diffuser shielding the incoming fuel from the heat radiation from the flame and deflecting the primary air and fuel steam into the secondary air, and a watercooled annulus located between the primary air and secondary air flows.

  9. U.S. Refining Capacity Utilization

    Reports and Publications (EIA)

    1995-01-01

    This article briefly reviews recent trends in domestic refining capacity utilization and examines in detail the differences in reported crude oil distillation capacities and utilization rates among different classes of refineries.

  10. Hydrogen Piping Experience in Chevron Refining

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Piping Experience in Chevron Refining Ned Niccolls Materials Engineer Chevron Energy Technology Company Hydrogen Pipeline Working Group Workshop August 30-31, 2005 Outline 2 Overall perspectives from long term use of hydrogen piping in refining. Piping specifications and practices. The (few) problem areas. Related industry work: American Petroleum Institute corrosion and materials work on high temperature hydrogen attack. Overall Perspectives 3 Few problems with hydrogen piping operating at

  11. Refining Bio-Oil alongside Petroleum | Department of Energy

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Refining Bio-Oil alongside Petroleum Refining Bio-Oil alongside Petroleum April 9, 2013 - 12:00am Addthis W.R. Grace, a leading provider of refining technologies, and Pacific ...

  12. Method and apparatus for connecting finite element meshes and performing simulations therewith

    DOE Patents [OSTI]

    Dohrmann, Clark R.; Key, Samuel W.; Heinstein, Martin W.

    2003-05-06

    The present invention provides a method of connecting dissimilar finite element meshes. A first mesh, designated the master mesh, and a second mesh, designated the slave mesh, each have interface surfaces proximal the other. Each interface surface has a corresponding interface mesh comprising a plurality of interface nodes. Each slave interface node is assigned new coordinates locating the interface node on the interface surface of the master mesh. The slave interface surface is further redefined to be the projection of the slave interface mesh onto the master interface surface.

  13. Source: Energy Information Administration, Form EIA-782A, "Refiners...

    U.S. Energy Information Administration (EIA) Indexed Site

    Motor Gasoline No. 2 Distillate Residual Fuel Oil 5. U.S. Refiner Wholesale Petroleum Product Volumes Figure Percentages of Refiner Wholesale Volumes 1995 Annual Averages Motor...

  14. Source: Energy Information Administration, Form EIA-782A, "Refiners...

    U.S. Energy Information Administration (EIA) Indexed Site

    Day Motor Gasoline No. 2 Distillate Residual Fuel Oil 3. U.S. Refiner Retail Petroleum Product Volumes Figure Percentages of Refiner Retail Volumes 1997 Annual Averages Motor...

  15. Source: Energy Information Administration, Form EIA-782A, "Refiners...

    U.S. Energy Information Administration (EIA) Indexed Site

    Motor Gasoline No. 2 Distillate Residual Fuel Oil 5. U.S. Refiner Wholesale Petroleum Product Volumes Figure Percentages of Refiner Wholesale Volumes 1997 Annual Averages Motor...

  16. Source: Energy Information Administration, Form EIA-782A, "Refiners...

    U.S. Energy Information Administration (EIA) Indexed Site

    Day Motor Gasoline No. 2 Distillate Residual Fuel Oil 3. U.S. Refiner Retail Petroleum Product Volumes Figure Percentages of Refiner Retail Volumes 1995 Annual Averages Motor...

  17. Source: Energy Information Administration, Form EIA-782A, "Refiners...

    U.S. Energy Information Administration (EIA) Indexed Site

    Day Motor Gasoline No. 2 Distillate Residual Fuel Oil 3. U.S. Refiner Retail Petroleum Product Volumes Figure Percentages of Refiner Retail Volumes 1996 Annual Averages Motor...

  18. Source: Energy Information Administration, Form EIA-782A, "Refiners...

    U.S. Energy Information Administration (EIA) Indexed Site

    Motor Gasoline No. 2 Distillate Residual Fuel Oil 5. U.S. Refiner Wholesale Petroleum Product Volumes Figure Percentages of Refiner Wholesale Volumes 1996 Annual Averages Motor...

  19. ,"U.S. Motor Gasoline Refiner Sales Volumes"

    U.S. Energy Information Administration (EIA) Indexed Site

    Data for" ,"Data 1","U.S. Motor Gasoline Refiner Sales ... AM" "Back to Contents","Data 1: U.S. Motor Gasoline Refiner Sales Volumes" ...

  20. ITP Petroleum Refining: Impacts of Condition Assessment on Energy...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Impacts of Condition Assessment on Energy Use: Selected Applications in Chemicals Processing and Petroleum Refining ITP Petroleum Refining: Impacts of Condition Assessment on ...

  1. Analyte separation utilizing temperature programmed desorption of a preconcentrator mesh

    DOE Patents [OSTI]

    Linker, Kevin L.; Bouchier, Frank A.; Theisen, Lisa; Arakaki, Lester H.

    2007-11-27

    A method and system for controllably releasing contaminants from a contaminated porous metallic mesh by thermally desorbing and releasing a selected subset of contaminants from a contaminated mesh by rapidly raising the mesh to a pre-determined temperature step or plateau that has been chosen beforehand to preferentially desorb a particular chemical specie of interest, but not others. By providing a sufficiently long delay or dwell period in-between heating pulses, and by selecting the optimum plateau temperatures, then different contaminant species can be controllably released in well-defined batches at different times to a chemical detector in gaseous communication with the mesh. For some detectors, such as an Ion Mobility Spectrometer (IMS), separating different species in time before they enter the IMS allows the detector to have an enhanced selectivity.

  2. Finite element meshing approached as a global minimization process

    SciTech Connect (OSTI)

    WITKOWSKI,WALTER R.; JUNG,JOSEPH; DOHRMANN,CLARK R.; LEUNG,VITUS J.

    2000-03-01

    The ability to generate a suitable finite element mesh in an automatic fashion is becoming the key to being able to automate the entire engineering analysis process. However, placing an all-hexahedron mesh in a general three-dimensional body continues to be an elusive goal. The approach investigated in this research is fundamentally different from any other that is known of by the authors. A physical analogy viewpoint is used to formulate the actual meshing problem which constructs a global mathematical description of the problem. The analogy used was that of minimizing the electrical potential of a system charged particles within a charged domain. The particles in the presented analogy represent duals to mesh elements (i.e., quads or hexes). Particle movement is governed by a mathematical functional which accounts for inter-particles repulsive, attractive and alignment forces. This functional is minimized to find the optimal location and orientation of each particle. After the particles are connected a mesh can be easily resolved. The mathematical description for this problem is as easy to formulate in three-dimensions as it is in two- or one-dimensions. The meshing algorithm was developed within CoMeT. It can solve the two-dimensional meshing problem for convex and concave geometries in a purely automated fashion. Investigation of the robustness of the technique has shown a success rate of approximately 99% for the two-dimensional geometries tested. Run times to mesh a 100 element complex geometry were typically in the 10 minute range. Efficiency of the technique is still an issue that needs to be addressed. Performance is an issue that is critical for most engineers generating meshes. It was not for this project. The primary focus of this work was to investigate and evaluate a meshing algorithm/philosophy with efficiency issues being secondary. The algorithm was also extended to mesh three-dimensional geometries. Unfortunately, only simple geometries were tested

  3. Metal decontamination for waste minimization using liquid metal refining technology

    SciTech Connect (OSTI)

    Joyce, E.L. Jr.; Lally, B.; Ozturk, B.; Fruehan, R.J.

    1993-09-01

    The current Department of Energy Mixed Waste Treatment Project flowsheet indicates that no conventional technology, other than surface decontamination, exists for metal processing. Current Department of Energy guidelines require retrievable storage of all metallic wastes containing transuranic elements above a certain concentration. This project is in support of the National Mixed Low Level Waste Treatment Program. Because of the high cost of disposal, it is important to develop an effective decontamination and volume reduction method for low-level contaminated metals. It is important to be able to decontaminate complex shapes where surfaces are hidden or inaccessible to surface decontamination processes and destruction of organic contamination. These goals can be achieved by adapting commercial metal refining processes to handle radioactive and organic contaminated metal. The radioactive components are concentrated in the slag, which is subsequently vitrified; hazardous organics are destroyed by the intense heat of the bath. The metal, after having been melted and purified, could be recycled for use within the DOE complex. In this project, we evaluated current state-of-the-art technologies for metal refining, with special reference to the removal of radioactive contaminants and the destruction of hazardous organics. This evaluation was based on literature reports, industrial experience, plant visits, thermodynamic calculations, and engineering aspects of the various processes. The key issues addressed included radioactive partitioning between the metal and slag phases, minimization of secondary wastes, operability of the process subject to widely varying feed chemistry, and the ability to seal the candidate process to prevent the release of hazardous species.

  4. Optimization-based mesh correction with volume and convexity constraints

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    D'Elia, Marta; Ridzal, Denis; Peterson, Kara J.; Bochev, Pavel; Shashkov, Mikhail

    2016-02-24

    Here, we consider the problem of finding a mesh such that 1) it is the closest, with respect to a suitable metric, to a given source mesh having the same connectivity, and 2) the volumes of its cells match a set of prescribed positive values that are not necessarily equal to the cell volumes in the source mesh. Also, this volume correction problem arises in important simulation contexts, such as satisfying a discrete geometric conservation law and solving transport equations by incremental remapping or similar semi-Lagrangian transport schemes. In this paper we formulate volume correction as a constrained optimization problemmore » in which the distance to the source mesh defines an optimization objective, while the prescribed cell volumes, mesh validity and/or cell convexity specify the constraints. We solve this problem numerically using a sequential quadratic programming (SQP) method whose performance scales with the mesh size. To achieve scalable performance we develop a specialized multigrid-based preconditioner for optimality systems that arise in the application of the SQP method to the volume correction problem. Numerical examples illustrate the importance of volume correction, and showcase the accuracy, robustness and scalability of our approach.« less

  5. Trends in petroleum refining process technology

    SciTech Connect (OSTI)

    Kowalczyk, D.

    1995-12-31

    In the 1990`s, the shift toward reformulated fuels and the unrelenting economic pressures on the petroleum refining industry have led to the ongoing development of a series of technological advances to improve fuels quality and industry operating efficiency. In this paper, ten of the most innovative and high impact recent developments in petroleum refining process technology will be highlighted. Process improvements and innovations have occurred in all facets of petroleum refining operations including fluid catalytic cracking, ether production, desulfurization, hydrocracking, gas processing, environmental control and heavy oil processing. Discussed will be the technical and economic impact of each of these new technologies on the petroleum refinery of the late 20th and early 21st century.

  6. JAPAN: Refining options and liberalization plans

    SciTech Connect (OSTI)

    Totto, L.; Isaak, D.T.

    1988-01-01

    The reformulation of Japan's petroleum import policy is important to Japan and to the international energy community. The technical and economic factors involved in opening the Japanese market is unknown. Assuming that foreign refiners have the technical capability to meet Japanese product specifications, the prices and volumes that will prove economical to both parties are also little known. Japanese refiners have been investigating the possibility of exporting gasoline and processing crude for China and Malaysia. The study demonstrates that under certain trade policies, Japan could become a significant product exporter. This study's purpose is to investigate the economic efficiency of alternative petroleum supply options, i.e., combining Japan's refining and product imports. An evaluation of the current import program and the technically feasible alternatives will play a major role in a new import policy. For this evaluation, a linear programming (LP) model of Japan's domestic refinery capacity and the petroleum import system was constructed. 3 figs., 23 tabs.

  7. A wavelet-MRA-based adaptive semi-Lagrangian method for the relativistic Vlasov-Maxwell system

    SciTech Connect (OSTI)

    Besse, Nicolas Latu, Guillaume Ghizzo, Alain Sonnendruecker, Eric Bertrand, Pierre

    2008-08-10

    In this paper we present a new method for the numerical solution of the relativistic Vlasov-Maxwell system on a phase-space grid using an adaptive semi-Lagrangian method. The adaptivity is performed through a wavelet multiresolution analysis, which gives a powerful and natural refinement criterion based on the local measurement of the approximation error and regularity of the distribution function. Therefore, the multiscale expansion of the distribution function allows to get a sparse representation of the data and thus save memory space and CPU time. We apply this numerical scheme to reduced Vlasov-Maxwell systems arising in laser-plasma physics. Interaction of relativistically strong laser pulses with overdense plasma slabs is investigated. These Vlasov simulations revealed a rich variety of phenomena associated with the fast particle dynamics induced by electromagnetic waves as electron trapping, particle acceleration, and electron plasma wavebreaking. However, the wavelet based adaptive method that we developed here, does not yield significant improvements compared to Vlasov solvers on a uniform mesh due to the substantial overhead that the method introduces. Nonetheless they might be a first step towards more efficient adaptive solvers based on different ideas for the grid refinement or on a more efficient implementation. Here the Vlasov simulations are performed in a two-dimensional phase-space where the development of thin filaments, strongly amplified by relativistic effects requires an important increase of the total number of points of the phase-space grid as they get finer as time goes on. The adaptive method could be more useful in cases where these thin filaments that need to be resolved are a very small fraction of the hyper-volume, which arises in higher dimensions because of the surface-to-volume scaling and the essentially one-dimensional structure of the filaments. Moreover, the main way to improve the efficiency of the adaptive method is to

  8. Using supercritical fluids to refine hydrocarbons

    DOE Patents [OSTI]

    Yarbro, Stephen Lee

    2015-06-09

    A system and method for reactively refining hydrocarbons, such as heavy oils with API gravities of less than 20 degrees and bitumen-like hydrocarbons with viscosities greater than 1000 cp at standard temperature and pressure, using a selected fluid at supercritical conditions. A reaction portion of the system and method delivers lightweight, volatile hydrocarbons to an associated contacting unit which operates in mixed subcritical/supercritical or supercritical modes. Using thermal diffusion, multiphase contact, or a momentum generating pressure gradient, the contacting unit separates the reaction products into portions that are viable for use or sale without further conventional refining and hydro-processing techniques.

  9. The strategic outlook for petroleum refiners

    SciTech Connect (OSTI)

    Cobb, C.

    1994-12-31

    After several years of acceptable and relatively reliable profitability, refiners were plunged into uncertainty once again following passage of the Clean Air Act Amendments (CAAA) in 1990. While many of the original uncertainties that attended these regulations have since been resolved, the industry now must prepare for the new manufacturing, distributions, and transportation challenges that undoubtedly will accompany the first stages of reformulated gasoline (RFG) compliance in early 1995. This impending challenge introduces several fundamental questions: (1) How has the industry changed and adjusted to respond to anticipated future needs? (2) What strategies are refiners employing today? (3) What are industry performance expectations over the next several years?

  10. Unconstrained plastering : all-hexahedral mesh generation via advancing front geometry decomposition (2004-2008).

    SciTech Connect (OSTI)

    Blacker, Teddy Dean; Staten, Matthew L.; Kerr, Robert A.; Owen, Steven James

    2010-03-01

    The generation of all-hexahedral finite element meshes has been an area of ongoing research for the past two decades and remains an open problem. Unconstrained plastering is a new method for generating all-hexahedral finite element meshes on arbitrary volumetric geometries. Starting from an unmeshed volume boundary, unconstrained plastering generates the interior mesh topology without the constraints of a pre-defined boundary mesh. Using advancing fronts, unconstrained plastering forms partially defined hexahedral dual sheets by decomposing the geometry into simple shapes, each of which can be meshed with simple meshing primitives. By breaking from the tradition of previous advancing-front algorithms, which start from pre-meshed boundary surfaces, unconstrained plastering demonstrates that for the tested geometries, high quality, boundary aligned, orientation insensitive, all-hexahedral meshes can be generated automatically without pre-meshing the boundary. Examples are given for meshes from both solid mechanics and geotechnical applications.

  11. Focus on Venezuelan heavy crude: refining margins

    SciTech Connect (OSTI)

    Not Available

    1984-01-25

    Of six crudes refined in the US Gulf Coast, heavy Venezuelan crude Lagunillas (15/sup 0/ API) provides the best margin per barrel. Data for end of December 1983 and the first three weeks of January show that margins on all crudes are on the rise in this market, due to a turnaround in product prices. The lighter crudes are showing the greatest increase in Gross Product Worth. This is having a modest shrinking effect on the margin differential between light and heavy crudes in this market. The domestic crude West Texas Intermediate, at 40/sup 0/ API, provides the highest GPW in this crude slate sample, over US $31 per barrel, compared to GPW of under US $28 per barrel for Lagunillas. Still, as Lagunillas cost about US $8 less than does WTI, refiners with sufficient residue conversion capacity can be earning about US $3.50 more in margin per barrel than they can with WTI. Although few refiners would be using a 15/sup 0/ API crude exclusively for any length of time, heavier oil's inclusion in modern refiners' diets is enhancing their competitive position more than any other single factor. This issue of Energy Detente presents the fuel price/tax series and industrial fuel prices for January 1984 for countries of the Western Hemisphere.

  12. Energy Bandwidth for Petroleum Refining Processes

    SciTech Connect (OSTI)

    none,

    2006-10-01

    The petroleum refining energy bandwidth report analyzes the most energy-intensive unit operations used in U.S. refineries: crude oil distillation, fluid catalytic cracking, catalytic hydrotreating, catalytic reforming, and alkylation. The "bandwidth" provides a snapshot of the energy losses that can potentially be recovered through best practices and technology R&D.

  13. First-Principles Simulations of High-Speed Combustion and Detonation...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Physical model includes multi-species equation of state, viscosity, diffusion, and 8 ... a uniform grid as well as static and dynamic adaptive mesh refinement capabilities

  14. BoxLib

    Energy Science and Technology Software Center (OSTI)

    2000-10-01

    BoxLib is a C++ foundation library used to aid in constructing parallel grid based Partial Differential Equation (PDE) solvers, particularly for Adaptive Mesh Refinement (AMR) applications.

  15. Using supercritical fluids to refine hydrocarbons

    DOE Patents [OSTI]

    Yarbro, Stephen Lee

    2014-11-25

    This is a method to reactively refine hydrocarbons, such as heavy oils with API gravities of less than 20.degree. and bitumen-like hydrocarbons with viscosities greater than 1000 cp at standard temperature and pressure using a selected fluid at supercritical conditions. The reaction portion of the method delivers lighter weight, more volatile hydrocarbons to an attached contacting device that operates in mixed subcritical or supercritical modes. This separates the reaction products into portions that are viable for use or sale without further conventional refining and hydro-processing techniques. This method produces valuable products with fewer processing steps, lower costs, increased worker safety due to less processing and handling, allow greater opportunity for new oil field development and subsequent positive economic impact, reduce related carbon dioxide, and wastes typical with conventional refineries.

  16. Refining New-Physics Searches in B→Dτν with Lattice QCD | Argonne

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Data Series: Refiner Crude Oil Inputs Refiner Gross Inputs Refiner Operable Capacity (Calendar Day) Refiner Percent Operable Utilization Net Inputs (Refiner and Blender) of Motor Gasoline Blending Comp Net Inputs (Refiner and Blender) of RBOB Blending Components Net Inputs (Refiner and Blender) of CBOB Blending Components Net Inputs (Refiner and Blender) of GTAB Blending Components Net Inputs (Refiner and Blender) of All Other Blending Component Net Inputs (Refiner and Blender) of Fuel Ethanol

  17. Turbine component cooling channel mesh with intersection chambers

    DOE Patents [OSTI]

    Lee, Ching-Pang; Marra, John J

    2014-05-06

    A mesh (35) of cooling channels (35A, 35B) with an array of cooling channel intersections (42) in a wall (21, 22) of a turbine component. A mixing chamber (42A-C) at each intersection is wider (W1, W2)) than a width (W) of each of the cooling channels connected to the mixing chamber. The mixing chamber promotes swirl, and slows the coolant for more efficient and uniform cooling. A series of cooling meshes (M1, M2) may be separated by mixing manifolds (44), which may have film cooling holes (46) and/or coolant refresher holes (48).

  18. Secretary Bodman Tours Refinery and Calls for More Domestic Refining...

    Office of Environmental Management (EM)

    Bodman today renewed the call for expanded oil refining capacity in the United States and ... and one way to do that is to strengthen and expand our domestic oil refining capacity. ...

  19. Diameter-Refined Metallic Carbon Nanotubes as Optically Tunable Transparent

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Conductors | ANSER Center | Argonne-Northwestern National Laboratory Diameter-Refined Metallic Carbon Nanotubes as Optically Tunable Transparent Conductors Home > Research > ANSER Research Highlights > Diameter-Refined Metallic Carbon Nanotubes as Optically Tunable Transparent Conductors

  20. ,"Residual Fuel Oil Sales to End Users Refiner Sales Volumes...

    U.S. Energy Information Administration (EIA) Indexed Site

    Data for" ,"Data 1","Residual Fuel Oil Sales to End Users Refiner Sales ... "Back to Contents","Data 1: Residual Fuel Oil Sales to End Users Refiner Sales Volumes" ...

  1. ,"U.S. Residual Fuel Oil Refiner Sales Volumes"

    U.S. Energy Information Administration (EIA) Indexed Site

    Data for" ,"Data 1","U.S. Residual Fuel Oil Refiner Sales Volumes",2,"Monthly","5... "Back to Contents","Data 1: U.S. Residual Fuel Oil Refiner Sales Volumes" ...

  2. Grain Refinement of Permanent Mold Cast Copper Base Alloys

    SciTech Connect (OSTI)

    M.Sadayappan; J.P.Thomson; M.Elboujdaini; G.Ping Gu; M. Sahoo

    2005-04-01

    Grain refinement is a well established process for many cast and wrought alloys. The mechanical properties of various alloys could be enhanced by reducing the grain size. Refinement is also known to improve casting characteristics such as fluidity and hot tearing. Grain refinement of copper-base alloys is not widely used, especially in sand casting process. However, in permanent mold casting of copper alloys it is now common to use grain refinement to counteract the problem of severe hot tearing which also improves the pressure tightness of plumbing components. The mechanism of grain refinement in copper-base alloys is not well understood. The issues to be studied include the effect of minor alloy additions on the microstructure, their interaction with the grain refiner, effect of cooling rate, and loss of grain refinement (fading). In this investigation, efforts were made to explore and understand grain refinement of copper alloys, especially in permanent mold casting conditions.

  3. Bandwidth Study U.S. Petroleum Refining | Department of Energy

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Petroleum Refining Bandwidth Study U.S. Petroleum Refining Bandwidth Study U.S. Petroleum Refining Energy bandwidth studies of U.S. manufacturing sectors can serve as foundational references in framing the range (or bandwidth) of potential energy savings opportunities. This bandwidth study examines energy consumption and potential energy savings opportunities in U.S. petroleum refining. The study relies on multiple sources to estimate the energy used in nine individual process areas,

  4. Parallel performance optimizations on unstructured mesh-based simulations

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Sarje, Abhinav; Song, Sukhyun; Jacobsen, Douglas; Huck, Kevin; Hollingsworth, Jeffrey; Malony, Allen; Williams, Samuel; Oliker, Leonid

    2015-06-01

    This paper addresses two key parallelization challenges the unstructured mesh-based ocean modeling code, MPAS-Ocean, which uses a mesh based on Voronoi tessellations: (1) load imbalance across processes, and (2) unstructured data access patterns, that inhibit intra- and inter-node performance. Our work analyzes the load imbalance due to naive partitioning of the mesh, and develops methods to generate mesh partitioning with better load balance and reduced communication. Furthermore, we present methods that minimize both inter- and intranode data movement and maximize data reuse. Our techniques include predictive ordering of data elements for higher cache efficiency, as well as communication reduction approaches.more » We present detailed performance data when running on thousands of cores using the Cray XC30 supercomputer and show that our optimization strategies can exceed the original performance by over 2×. Additionally, many of these solutions can be broadly applied to a wide variety of unstructured grid-based computations.« less

  5. Solves Poisson's Equation in Axizymmetric Geometry on a Rectangular Mesh

    Energy Science and Technology Software Center (OSTI)

    1996-09-10

    DATHETA4.0 computes the magnetostatic field produced by multiple point current sources in the presence of perfect conductors in axisymmetric geometry. DATHETA4.0 has an interactive user interface and solves Poisson''s equation using the ADI method on a rectangular finite-difference mesh. DATHETA4.0 uncludes models specific to applied-B ion diodes.

  6. Full Core Multiphysics Simulation with Offline Mesh Deformation

    SciTech Connect (OSTI)

    Merzari, E.; Shemon, E. R.; Yu, Y.; Thomas, J. W.; Obabko, A.; Jain, Rajeev; Mahadevan, Vijay; Solberg, Jerome; Ferencz, R.; Whitesides, R.

    2015-12-21

    In this report, building on previous reports issued in FY13 we describe our continued efforts to integrate thermal/hydraulics, neutronics, and structural mechanics modeling codes to perform coupled analysis of a representative fast sodium-cooled reactor core. The focus of the present report is a full core simulation with off-line mesh deformation.

  7. Generating unstructured nuclear reactor core meshes in parallel

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Jain, Rajeev; Tautges, Timothy J.

    2014-10-24

    Recent advances in supercomputers and parallel solver techniques have enabled users to run large simulations problems using millions of processors. Techniques for multiphysics nuclear reactor core simulations are under active development in several countries. Most of these techniques require large unstructured meshes that can be hard to generate in a standalone desktop computers because of high memory requirements, limited processing power, and other complexities. We have previously reported on a hierarchical lattice-based approach for generating reactor core meshes. Here, we describe efforts to exploit coarse-grained parallelism during reactor assembly and reactor core mesh generation processes. We highlight several reactor coremore » examples including a very high temperature reactor, a full-core model of the Korean MONJU reactor, a ¼ pressurized water reactor core, the fast reactor Experimental Breeder Reactor-II core with a XX09 assembly, and an advanced breeder test reactor core. The times required to generate large mesh models, along with speedups obtained from running these problems in parallel, are reported. A graphical user interface to the tools described here has also been developed.« less

  8. Generating unstructured nuclear reactor core meshes in parallel

    SciTech Connect (OSTI)

    Jain, Rajeev; Tautges, Timothy J.

    2014-10-24

    Recent advances in supercomputers and parallel solver techniques have enabled users to run large simulations problems using millions of processors. Techniques for multiphysics nuclear reactor core simulations are under active development in several countries. Most of these techniques require large unstructured meshes that can be hard to generate in a standalone desktop computers because of high memory requirements, limited processing power, and other complexities. We have previously reported on a hierarchical lattice-based approach for generating reactor core meshes. Here, we describe efforts to exploit coarse-grained parallelism during reactor assembly and reactor core mesh generation processes. We highlight several reactor core examples including a very high temperature reactor, a full-core model of the Korean MONJU reactor, a ¼ pressurized water reactor core, the fast reactor Experimental Breeder Reactor-II core with a XX09 assembly, and an advanced breeder test reactor core. The times required to generate large mesh models, along with speedups obtained from running these problems in parallel, are reported. A graphical user interface to the tools described here has also been developed.

  9. Genealogy of major U.S. refiners - Energy Information Administration

    U.S. Energy Information Administration (EIA) Indexed Site

    See full Genealogy of Major U.S. Refiners Previous Release Genealogy of Major U.S. Refiners Genealogy of Major U.S. Refiners Release date: September 18, 2013 figre 1. World energy consumption, 1990-2040. The structure of the U.S. petroleum refining industry has changed substantially over the past several years. In the diagram the companies shown on the right side are presently active in U.S. refining. The transactions over the past 25 years that created these companies also are shown. The

  10. Deformable elastic network refinement for low-resolution macromolecular crystallography

    SciTech Connect (OSTI)

    Schrder, Gunnar F.; Levitt, Michael; Brunger, Axel T.

    2014-09-01

    An overview of applications of the deformable elastic network (DEN) refinement method is presented together with recommendations for its optimal usage. Crystals of membrane proteins and protein complexes often diffract to low resolution owing to their intrinsic molecular flexibility, heterogeneity or the mosaic spread of micro-domains. At low resolution, the building and refinement of atomic models is a more challenging task. The deformable elastic network (DEN) refinement method developed previously has been instrumental in the determinion of several structures at low resolution. Here, DEN refinement is reviewed, recommendations for its optimal usage are provided and its limitations are discussed. Representative examples of the application of DEN refinement to challenging cases of refinement at low resolution are presented. These cases include soluble as well as membrane proteins determined at limiting resolutions ranging from 3 to 7 . Potential extensions of the DEN refinement technique and future perspectives for the interpretation of low-resolution crystal structures are also discussed.

  11. Pollution prevention in the petroleum refining industry

    SciTech Connect (OSTI)

    Fromm, C.H.; White, S.L.

    1995-09-01

    Pollution prevention (P2) as applied to petroleum refining should seek opportunities to reduce waste by preventing oil/hydrocarbon loss, by decreasing consumption of auxiliary input materials, and by improving conversion of incoming impurities into useful products. This chapter will focus on P2 techniques that have found or could find applications in minimizing or eliminating reducible solid waste and wastewater from the petroleum refining process. Air emissions are not covered. The following commonly encountered wastes are considered in this chapter: oily sludges; spent caustics; spent catalysts; miscellaneous process wastes; wastewater; maintenance and materials handling wastes. Following a brief description of waste components and sources, specific P2 techniques are presented in tabular form for each of these wastestreams. None of the P2 techniques presented is discussed here in any detail--the intent is to give the reader a menu of potentially effective P2 options to consider, along with the references where a more detailed discussion may be found. Some of the options presented were advanced in the original references merely as suggestions or plans for improvement. No effort was made to verify their efficacy or applicability in this compilation.

  12. A Novel Coarsening Method for Scalable and Efficient Mesh Generation

    SciTech Connect (OSTI)

    Yoo, A; Hysom, D; Gunney, B

    2010-12-02

    In this paper, we propose a novel mesh coarsening method called brick coarsening method. The proposed method can be used in conjunction with any graph partitioners and scales to very large meshes. This method reduces problem space by decomposing the original mesh into fixed-size blocks of nodes called bricks, layered in a similar way to conventional brick laying, and then assigning each node of the original mesh to appropriate brick. Our experiments indicate that the proposed method scales to very large meshes while allowing simple RCB partitioner to produce higher-quality partitions with significantly less edge cuts. Our results further indicate that the proposed brick-coarsening method allows more complicated partitioners like PT-Scotch to scale to very large problem size while still maintaining good partitioning performance with relatively good edge-cut metric. Graph partitioning is an important problem that has many scientific and engineering applications in such areas as VLSI design, scientific computing, and resource management. Given a graph G = (V,E), where V is the set of vertices and E is the set of edges, (k-way) graph partitioning problem is to partition the vertices of the graph (V) into k disjoint groups such that each group contains roughly equal number of vertices and the number of edges connecting vertices in different groups is minimized. Graph partitioning plays a key role in large scientific computing, especially in mesh-based computations, as it is used as a tool to minimize the volume of communication and to ensure well-balanced load across computing nodes. The impact of graph partitioning on the reduction of communication can be easily seen, for example, in different iterative methods to solve a sparse system of linear equation. Here, a graph partitioning technique is applied to the matrix, which is basically a graph in which each edge is a non-zero entry in the matrix, to allocate groups of vertices to processors in such a way that many of

  13. Extended Virtual Spring Mesh (EVSM): The Distributed Self-Organizing Mobile Ad Hoc Network for Area Exploration

    SciTech Connect (OSTI)

    Kurt Derr

    2011-12-01

    Mobile Ad hoc NETworks (MANETs) are distributed self-organizing networks that can change locations and configure themselves on the fly. This paper focuses on an algorithmic approach for the deployment of a MANET within an enclosed area, such as a building in a disaster scenario, which can provide a robust communication infrastructure for search and rescue operations. While a virtual spring mesh (VSM) algorithm provides scalable, self-organizing, and fault-tolerant capabilities required by aMANET, the VSM lacks the MANET's capabilities of deployment mechanisms for blanket coverage of an area and does not provide an obstacle avoidance mechanism. This paper presents a new technique, an extended VSM (EVSM) algorithm that provides the following novelties: (1) new control laws for exploration and expansion to provide blanket coverage, (2) virtual adaptive springs enabling the mesh to expand as necessary, (3) adapts to communications disturbances by varying the density and movement of mobile nodes, and (4) new metrics to assess the performance of the EVSM algorithm. Simulation results show that EVSM provides up to 16% more coverage and is 3.5 times faster than VSM in environments with eight obstacles.

  14. Surface micromachined counter-meshing gears discrimination device

    DOE Patents [OSTI]

    Polosky, Marc A.; Garcia, Ernest J.; Allen, James J.

    2000-12-12

    A surface micromachined Counter-Meshing Gears (CMG) discrimination device which functions as a mechanically coded lock. Each of two CMG has a first portion of its perimeter devoted to continuous driving teeth that mesh with respective pinion gears. Each EMG also has a second portion of its perimeter devoted to regularly spaced discrimination gear teeth that extend outwardly on at least one of three levels of the CMG. The discrimination gear teeth are designed so as to pass each other without interference only if the correct sequence of partial rotations of the CMG occurs in response to a coded series of rotations from the pinion gears. A 24 bit code is normally input to unlock the device. Once unlocked, the device provides a path for an energy or information signal to pass through the device. The device is designed to immediately lock up if any portion of the 24 bit code is incorrect.

  15. GMINC: a mesh generator for flow simulations in fractured reservoirs

    SciTech Connect (OSTI)

    Pruess, K.

    1983-03-01

    GMINC is a pre-processor computer program for generating geometrical meshes to be used in modeling fluid and heat flow in fractured porous media. It is based on the method of multiple interacting continua (MINC) as developed by Pruess and Narasimhan. The meshes generated by GMINC are in integral finite difference form, and are compatible with the simulators SHAFT79 and MULKOM. Applications with other integral finite difference simulators are possible, and require slight modifications in input/output formats. This report describes methodology and application of GMINC, including preparation of input decks and sample problems. A rather comprehensive overview of the MINC-method is also provided to make the presentation self-contained as a guide for modeling of flow in naturally fractured media.

  16. MESH2D GRID GENERATOR DESIGN AND USE

    SciTech Connect (OSTI)

    Flach, G.; Smith, F.

    2012-01-20

    Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j{sub 0}) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations. The overall mesh is constructed from grid zones that are typically then subdivided into a collection of smaller grid cells. The grid zones usually correspond to distinct materials or larger-scale geometric shapes. The structured grid zones are identified through uppercase indices (I,J). Subdivision of zonal regions into grid cells can be done uniformly, or nonuniformly using either a polynomial or geometric skewing algorithm. Grid cells may be concentrated backward, forward, or toward both ends. Figure 1 illustrates the above concepts in the context of a simple four zone grid.

  17. A novel three-dimensional mesh deformation method based on sphere relaxation

    SciTech Connect (OSTI)

    Zhou, Xuan; Li, Shuixiang

    2015-10-01

    In our previous work (2013) [19], we developed a disk relaxation based mesh deformation method for two-dimensional mesh deformation. In this paper, the idea of the disk relaxation is extended to the sphere relaxation for three-dimensional meshes with large deformations. We develop a node based pre-displacement procedure to apply initial movements on nodes according to their layer indices. Afterwards, the nodes are moved locally by the improved sphere relaxation algorithm to transfer boundary deformations and increase the mesh quality. A three-dimensional mesh smoothing method is also adopted to prevent the occurrence of the negative volume of elements, and further improve the mesh quality. Numerical applications in three-dimension including the wing rotation, bending beam and morphing aircraft are carried out. The results demonstrate that the sphere relaxation based approach generates the deformed mesh with high quality, especially regarding complex boundaries and large deformations.

  18. Adaptive Thresholds

    Energy Science and Technology Software Center (OSTI)

    2014-08-26

    ADAPT is a topological analysis code that allow to compute local threshold, in particular relevance based thresholds for features defined in scalar fields. The initial target application is vortex detection but the software is more generally applicable to all threshold based feature definitions.

  19. ITP Petroleum Refining: Energy and Environmental Profile of the U.S.

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Petroleum Refining Industry (November 2007) | Department of Energy and Environmental Profile of the U.S. Petroleum Refining Industry (November 2007) ITP Petroleum Refining: Energy and Environmental Profile of the U.S. Petroleum Refining Industry (November 2007) profile.pdf (1.62 MB) More Documents & Publications ITP Petroleum Refining: Energy Bandwidth for Petroleum Refining Processes ITP Petroleum Refining: Energy and Environmental Profile of the U.S. Petroleum Refining Industry

  20. ITP Petroleum Refining: Petroleum Technology Vision 2020 | Department...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Petroleum Technology Vision 2020 ITP Petroleum Refining: Petroleum Technology Vision 2020 techvision.pdf (684.96 KB) More Documents & Publications Manufacturing Energy and Carbon ...

  1. The Space-Based Telescopes for Actionable Refinement of Ephemeris...

    Office of Scientific and Technical Information (OSTI)

    Title: The Space-Based Telescopes for Actionable Refinement of Ephemeris Pathfinder Mission Authors: Simms, L ; De Vries, W ; RIot, V ; Olivier, S ; Pertica, A ; Bauman, B ; ...

  2. ,"U.S. Total Refiner Petroleum Product Prices"

    U.S. Energy Information Administration (EIA) Indexed Site

    NUSDPG","EMAEPPRPTGNUSDPG","EMAEPPRLPTGNUSDPG","EMAEPPRHPTGNUSDPG" "Date","U.S. Total Gasoline Retail Sales by Refiners (Dollars per Gallon)","U.S. Aviation Gasoline...

  3. The Space-Based Telescopes for Actionable Refinement of Ephemeris...

    Office of Scientific and Technical Information (OSTI)

    Journal Article: The Space-Based Telescopes for Actionable Refinement of Ephemeris Pathfinder Mission Citation Details In-Document Search Title: The Space-Based Telescopes for ...

  4. ,"U.S. Reformulated, Average Refiner Gasoline Prices"

    U.S. Energy Information Administration (EIA) Indexed Site

    ...www.eia.govdnavpetpetprirefmg2cnusepm0rdpgalm.htm" ,"Source:","Energy Information ... Reformulated Gasoline Retail Sales by Refiners (Dollars per ...

  5. ,"U.S. Conventional, Average Refiner Gasoline Prices"

    U.S. Energy Information Administration (EIA) Indexed Site

    ...www.eia.govdnavpetpetprirefmg2cnusepm0udpgalm.htm" ,"Source:","Energy Information ... Conventional Gasoline Retail Sales by Refiners (Dollars per ...

  6. Refiners react to changes in the pipeline infrastructure

    SciTech Connect (OSTI)

    Giles, K.A.

    1997-06-01

    Petroleum pipelines have long been a critical component in the distribution of crude and refined products in the U.S. Pipelines are typically the most cost efficient mode of transportation for reasonably consistent flow rates. For obvious reasons, inland refineries and consumers are much more dependent on petroleum pipelines to provide supplies of crude and refined products than refineries and consumers located on the coasts. Significant changes in U.S. distribution patterns for crude and refined products are reshaping the pipeline infrastructure and presenting challenges and opportunities for domestic refiners. These changes are discussed.

  7. EERE Success Story—Refining Bio-Oil alongside Petroleum

    Office of Energy Efficiency and Renewable Energy (EERE)

    W.R. Grace and Pacific Northwest National Laboratory are working to establish a bio-oil refining process that users existing petroleum refinery infrastructure.

  8. New Process for Grain Refinement of Aluminum. Final Report

    SciTech Connect (OSTI)

    Dr. Joseph A. Megy

    2000-09-22

    A new method of grain refining aluminum involving in-situ formation of boride nuclei in molten aluminum just prior to casting has been developed in the subject DOE program over the last thirty months by a team consisting of JDC, Inc., Alcoa Technical Center, GRAS, Inc., Touchstone Labs, and GKS Engineering Services. The Manufacturing process to make boron trichloride for grain refining is much simpler than preparing conventional grain refiners, with attendant environmental, capital, and energy savings. The manufacture of boride grain refining nuclei using the fy-Gem process avoids clusters, salt and oxide inclusions that cause quality problems in aluminum today.

  9. The Space-based Telescopes for Actionable Refinement of Ephemeris...

    Office of Scientific and Technical Information (OSTI)

    Title: The Space-based Telescopes for Actionable Refinement of Ephemeris (STARE) mission Authors: Riot, V ; Devries, W ; Bauman, B ; Simms, L ; Carter, D ; Phillion, D ; Olivier, S ...

  10. Table 35. Refiner Motor Gasoline Prices by Grade, Sales Type...

    U.S. Energy Information Administration (EIA) Indexed Site

    Energy Information Administration Petroleum Marketing Annual 1995 Table 35. Refiner Motor Gasoline Prices by Grade, Sales Type, PAD District, and State (Cents per Gallon...

  11. Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales...

    U.S. Energy Information Administration (EIA) Indexed Site

    250 Energy Information AdministrationPetroleum Marketing Annual 1999 Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales Type, PAD District, and State (Thousand Gallons...

  12. Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales...

    U.S. Energy Information Administration (EIA) Indexed Site

    Energy Information Administration Petroleum Marketing Annual 1995 Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales Type, PAD District, and State (Thousand Gallons...

  13. Table 35. Refiner Motor Gasoline Prices by Grade, Sales Type...

    U.S. Energy Information Administration (EIA) Indexed Site

    134 Energy Information AdministrationPetroleum Marketing Annual 1998 Table 35. Refiner Motor Gasoline Prices by Grade, Sales Type, PAD District, and State (Cents per Gallon...

  14. Table 35. Refiner Motor Gasoline Prices by Grade, Sales Type...

    U.S. Energy Information Administration (EIA) Indexed Site

    134 Energy Information AdministrationPetroleum Marketing Annual 1999 Table 35. Refiner Motor Gasoline Prices by Grade, Sales Type, PAD District, and State (Cents per Gallon...

  15. Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type...

    U.S. Energy Information Administration (EIA) Indexed Site

    220 Energy Information AdministrationPetroleum Marketing Annual 1998 Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State (Thousand Gallons per...

  16. Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales...

    U.S. Energy Information Administration (EIA) Indexed Site

    - - - - W W - - - - - - See footnotes at end of table. 44. Refiner Motor Gasoline Volumes by Formulation, Sales Type, PAD District, and State 292 Energy...

  17. Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales...

    U.S. Energy Information Administration (EIA) Indexed Site

    250 Energy Information AdministrationPetroleum Marketing Annual 1998 Table 44. Refiner Motor Gasoline Volumes by Formulation, Sales Type, PAD District, and State (Thousand Gallons...

  18. Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type...

    U.S. Energy Information Administration (EIA) Indexed Site

    220 Energy Information AdministrationPetroleum Marketing Annual 1999 Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State (Thousand Gallons per...

  19. Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type...

    U.S. Energy Information Administration (EIA) Indexed Site

    Energy Information Administration Petroleum Marketing Annual 1995 Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State (Thousand Gallons per...

  20. Petroleum Products Table 43. Refiner Motor Gasoline Volumes...

    U.S. Energy Information Administration (EIA) Indexed Site

    of table. 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State 262 Energy Information Administration Petroleum Marketing Annual 1997 Table 43....

  1. Petroleum Products Table 43. Refiner Motor Gasoline Volumes...

    U.S. Energy Information Administration (EIA) Indexed Site

    1995 Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State (Thousand Gallons per Day) - Continued Geographic Area Month Premium All Grades Sales...

  2. Petroleum Products Table 43. Refiner Motor Gasoline Volumes...

    U.S. Energy Information Administration (EIA) Indexed Site

    2000 Table 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State (Thousand Gallons per Day) - Continued Geographic Area Month Premium All Grades Sales...

  3. Petroleum Products Table 43. Refiner Motor Gasoline Volumes...

    U.S. Energy Information Administration (EIA) Indexed Site

    of table. 43. Refiner Motor Gasoline Volumes by Grade, Sales Type, PAD District, and State 262 Energy Information Administration Petroleum Marketing Annual 1996 Table 43....

  4. ITP Petroleum Refining: Energy Efficiency Roadmap for Petroleum...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Energy Efficiency Roadmap for Petroleum Refineries in California ITP Petroleum Refining: Energy Efficiency Roadmap for Petroleum Refineries in California refiningroadmap.pdf (1.34 ...

  5. Adaptive reconnection-based arbitrary Lagrangian Eulerian method

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Bo, Wurigen; Shashkov, Mikhail

    2015-07-21

    We present a new adaptive Arbitrary Lagrangian Eulerian (ALE) method. This method is based on the reconnection-based ALE (ReALE) methodology of Refs. [35], [34] and [6]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh in which the solution and positions of grid nodes are updated; a rezoning phase in which a new grid is defined by changing the connectivity (using Voronoi tessellation) but not the number of cells; and a remapping phase in which the Lagrangian solution is transferred onto the new grid. Furthermore, in the standard ReALEmore » method, the rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way.« less

  6. Adaptive sampler

    DOE Patents [OSTI]

    Watson, Bobby L.; Aeby, Ian

    1982-01-01

    An adaptive data compression device for compressing data having variable frequency content, including a plurality of digital filters for analyzing the content of the data over a plurality of frequency regions, a memory, and a control logic circuit for generating a variable rate memory clock corresponding to the analyzed frequency content of the data in the frequency region and for clocking the data into the memory in response to the variable rate memory clock.

  7. Adaptive sampler

    DOE Patents [OSTI]

    Watson, B.L.; Aeby, I.

    1980-08-26

    An adaptive data compression device for compressing data is described. The device has a frequency content, including a plurality of digital filters for analyzing the content of the data over a plurality of frequency regions, a memory, and a control logic circuit for generating a variable rate memory clock corresponding to the analyzed frequency content of the data in the frequency region and for clocking the data into the memory in response to the variable rate memory clock.

  8. Steel refining with an electrochemical cell

    DOE Patents [OSTI]

    Blander, Milton; Cook, Glenn M.

    1988-01-01

    Apparatus for processing a metallic fluid containing iron oxide, container for a molten metal including an electrically conductive refractory disposed for contact with the molten metal which contains iron oxide, an electrolyte in the form of a basic slag on top of the molten metal, an electrode in the container in contact with the slag electrically separated from the refractory, and means for establishing a voltage across the refractory and the electrode to reduce iron oxide to iron at the surface of the refractory in contact with the iron oxide containing fluid. A process is disclosed for refining an iron product containing not more than about 10% by weight oxygen and not more than about 10% by weight sulfur, comprising providing an electrolyte of a slag containing one or more of calcium oxide, magnesium oxide, silica or alumina, providing a cathode of the iron product in contact with the electrolyte, providing an anode in contact with the electrolyte electrically separated from the cathode, and operating an electrochemical cell formed by the anode, the cathode and the electrolyte to separate oxygen or sulfur present in the iron product therefrom.

  9. Steel refining with an electrochemical cell

    DOE Patents [OSTI]

    Blander, M.; Cook, G.M.

    1988-05-17

    Apparatus is described for processing a metallic fluid containing iron oxide, container for a molten metal including an electrically conductive refractory disposed for contact with the molten metal which contains iron oxide, an electrolyte in the form of a basic slag on top of the molten metal, an electrode in the container in contact with the slag electrically separated from the refractory, and means for establishing a voltage across the refractory and the electrode to reduce iron oxide to iron at the surface of the refractory in contact with the iron oxide containing fluid. A process is disclosed for refining an iron product containing not more than about 10% by weight oxygen and not more than about 10% by weight sulfur, comprising providing an electrolyte of a slag containing one or more of calcium oxide, magnesium oxide, silica or alumina, providing a cathode of the iron product in contact with the electrolyte, providing an anode in contact with the electrolyte electrically separated from the cathode, and operating an electrochemical cell formed by the anode, the cathode and the electrolyte to separate oxygen or sulfur present in the iron product therefrom. 2 figs.

  10. Steel refining with an electrochemical cell

    DOE Patents [OSTI]

    Blander, M.; Cook, G.M.

    1985-05-21

    Disclosed is an apparatus for processing a metallic fluid containing iron oxide, container for a molten metal including an electrically conductive refractory disposed for contact with the molten metal which contains iron oxide, an electrolyte in the form of a basic slag on top of the molten metal, an electrode in the container in contact with the slag electrically separated from the refractory, and means for establishing a voltage across the refractory and the electrode to reduce iron oxide to iron at the surface of the refractory in contact with the iron oxide containing fluid. A process is disclosed for refining an iron product containing not more than about 10% by weight sulfur, comprising providing an electrolyte of a slag containing one or more of calcium oxide, magnesium oxide, silica or alumina, providing a cathode of the iron product in contact with the electrolyte, providing an anode in contact with the electrolyte electrically separated from the cathode, and operating an electrochemical cell formed by the anode, the cathode and the electrolyte to separate oxygen or sulfur present in the iron product therefrom.