Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational geometry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

TU-FG-201-04: Computer Vision in Autonomous Quality Assurance of Linear Accelerators

Purpose: Routine quality assurance (QA) of linear accelerators represents a critical and costly element of a radiation oncology center. Recently, a system was developed to autonomously perform routine quality assurance on linear accelerators. The purpose of this work is to extend this system and contribute computer vision techniques for obtaining quantitative measurements for a monthly multi-leaf collimator (MLC) QA test specified by TG-142, namely leaf position accuracy, and demonstrate extensibility for additional routines. Methods: Grayscale images of a picket fence delivery on a radioluminescent phosphor coated phantom are captured using a CMOS camera. Collected images are processed to correct for camera distortions, rotation and alignment, reduce noise, and enhance contrast. The location of each MLC leaf is determined through logistic fitting and a priori modeling based on knowledge of the delivered beams. Using the data collected and the criteria from TG-142, a decision is made on whether or not the leaf position accuracy of the MLC passes or fails. Results: The locations of all MLC leaf edges are found for three different picket fence images in a picket fence routine to 0.1mm/1pixel precision. The program to correct for image alignment and determination of leaf positions requires a runtime of 21– 25 seconds for a single picket, and 44 – 46 seconds for a group of three pickets on a standard workstation CPU, 2.2 GHz Intel Core i7. Conclusion: MLC leaf edges were successfully found using techniques in computer vision. With the addition of computer vision techniques to the previously described autonomous QA system, the system is able to quickly perform complete QA routines with minimal human contribution.

60 APPLIED LIFE SCIENCES↗

WE-AB-207A-09: Optimization of the Design of a Moving Blocker for Cone-Beam CT Scatter Correction: Experimental Evaluation

Purpose: A moving blocker based strategy has shown promising results for scatter correction in cone-beam computed tomography (CBCT). Different geometry designs and moving speeds of the blocker affect its performance in image reconstruction accuracy. The goal of this work is to optimize the geometric design and moving speed of the moving blocker system through experimental evaluations. Methods: An Elekta Synergy XVI system and an anthropomorphic pelvis phantom CIRS 801-P were used for our experiment. A blocker consisting of lead strips was inserted between the x-ray source and the phantom moving back and forth along rotation axis to measure the scatter signal. Accoriding to our Monte Carlo simulation results, three blockers were used, which have the same lead strip width 3.2mm and different gap between neighboring lead strips, 3.2, 6.4 and 9.6mm. For each blocker, three moving speeds were evaluated, 10, 20 and 30 pixels per projection (on the detector plane). Scatter signal in the unblocked region was estimated by cubic B-spline based interpolation from the blocked region. CBCT image was reconstructed by a total variation (TV) based algebraic iterative reconstruction (ART) algorithm from the partially blocked projection data. Reconstruction accuracy in each condition is quantified as CT number error of region of interest (ROI) by comparing to a CBCT reconstructed image from analytically simulated unblocked and scatter free projection data. Results: Highest reconstruction accuracy is achieved when the blocker width is 3.2 mm, the gap between neighboring lead strips is 9.6 mm and the moving speed is 20 pixels per projection. RMSE of the CT number of ROIs can be reduced from 436 to 27. Conclusions: Image reconstruction accuracy is greatly affected by the geometry design of the blocker. The moving speed does not have a very strong effect on reconstruction result if it is over 20 pixels per projection.

60 APPLIED LIFE SCIENCES↗

DIF3D-VARIANT 11.0: A Decade of Updates

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 30 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems. VARIANT performs nodal neutron transport calculations using PN or SPN theory in Cartesian and hexagonal two- and threedimensional geometries. Because of the computing capabilities at that time, VARIANT was restricted to 33 group P3 flux approximations with P1 scattering. Clearly computer memory capabilities have increased since then and thus large space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original manual and details more than a decade worth of changes made to DIF3D to make DIF3D 11.0.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Final report for Texas A&M University Group Contribution to DE-FG02-09ER25949/DE-SC0002505: Topology for Statistical Modeling of Petascale Data (and ASCR-funded collaboration between Sandia National Labs, Texas A&M University and University of Utah)

We summarize the contributions of the Texas A\&M University Group to the project (DE-FG02-09ER25949/DE-SC0002505: Topology for Statistical Modeling of Petascale Data - an ASCR-funded collaboration between Sandia National Labs, Texas A\&M U, and U Utah) during 6/9/2011 -- 2/27/2013.

97 MATHEMATICS AND COMPUTING↗

Partnership for Edge Physics Simulation (EPSI)

We propose to develop advanced simulation codes, based upon an extreme parallelism, first principles kinetic approach, to address the challenges associated with the edge region of magnetically confined plasmas. This work is relevant to both existing magnetic fusion facilities and essential for next-generation burning plasma experiments, such as ITER where success is critically dependent upon H-mode operation achieving an edge pedestal of sufficient height for good core plasma performance without producing deleterious large scale edge localized instabilities. The plasma edge presents a well-known set of multi-physics, multi-scale problems involving complex 3D magnetic geometry. Perhaps the greatest computational challenge is the lack of scale separation – temporal scales for drift waves, Alfven waves, ELM dynamics for example have strong overlap. Similar overlap occurs on the spatial scales for the ion poloidal gyro-radius, drift wave and pedestal width. The traditional approach of separating fusion problems into weakly interacting spatial or temporal domains clearly breaks down in the edge. A full kinetic model (full-f model) must be solved to understand and predict the edge physics including non-equilibrium thermodynamic issues arising from the magnetic topology (the open field lines producing a spatially sensitive velocity hole), plasma wall interactions, neutral and atomic physics. The plan here is to model these phenomena within a comprehensive first principles set of equations without the need for the insurmountable multiple-codes coupling issues by building on the XGC1 code developed under the SciDAC Proto-FSP Center for Plasma Edge Simulation (CPES). This proposal includes the critical participants in the XGC1 development. We propose enhancing the capability of XGC1 by including all the important turbulence physics contained in kinetic ion and electron electromagnetic dynamics, by extending the PIC technology to incorporate several positive features found in the Eulerian grid technology, combining the coarse grained XGC0 code into it with a multiscale projective time integration technique to be developed in the project for experimental time scale simulation of the edge multi-physics, and by implementing modern computational technologies to move forward with the extreme scale heterogeneous hardware/software platform development. Ultra-fast edge localized instabilities, involving compressional Alven modes, will be studied using a proper MHD or two-fluid codes employing kinetic closure information from XGC. The proposed research integrates the most important physical processes on overlapping temporal and spatial scales, to study (i) transitions and thresholds for enhanced confinement regimes, (ii) the predictive understanding of the edge pedestal formation, structure and dynamics, (iii) the effect of the edge plasma on the core confinement, (iv) the physics of Edge Localized Modes and their suppression or mitigation via external control techniques, and (v) the heat load on the material wall. Our approach will implement physical models that are valid in their relevant collisionality regimes for current experiments and future burning plasma devices. Comparisons with experimental data will be made at multiple levels including characterization of turbulence dynamics, flux-gradient responses in the energy, particle and momentum channels, and observations of global phenomena such as L-H transitions, momentum generation, ELM dynamics, ELM-free regimes and the tokamak density limit. The proposal team has established strong connections, through joint membership, with other organizations that will be critical for success. These include FASTMath, SUPER and QUEST Institutes respectively for meshing tools and solvers; code performance tuning and uncertainty quantification as well as a proposed Data Management and Visualization Institute.. We will also be teaming with ongoing fusion SciDAC’s: CSWPI for RF-edge plasmas interactions; CEMM for the M3D MHD simulations and with PEPSC for incorporation of energetic particle physics. Finally, the team includes representatives from each of the major U.S. experimental facilities who will coordinate the validation program and the development of synthetic diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Technical Report for "Nuclear Technologies for Near Term Fusion Devices"

Over approximately 18 years, this project evolved to focus on a number of related topics, all tied to the nuclear analysis of fusion energy systems. For the earliest years, the University of Wisconsin (UW)’s effort was in support of the Advanced Power Extraction (APEX) study to investigate high power density first wall and blanket systems. A variety of design concepts were studied before this study gave way to a design effort for a US Test Blanket Module (TBM) to be installed in ITER. Simultaneous to this TBM project, nuclear analysis supported the conceptual design of a number of fusion nuclear science facilities that might fill a role in the path to fusion energy. Beginning in approximately 2005, this project added a component focused on the development of novel radiation transport software capability in support of the above nuclear analysis needs. Specifically, a clear need was identified to support neutron and photon transport on the complex geometries associated with Computer-Aided Design (CAD). Following the initial development of the Direct Accelerated Geoemtry Monte Carlo (DAGMC) capability, additional features were added, including unstructured mesh tallies and multi-physics analysis such as the Rigorous 2-Step (R2S) methodology for Shutdown Dose Rate (SDR) prediction. Throughout the project, there were also smaller tasks in support of the fusion materials community and for the testing of changes to the nuclear data that is fundamental to this kind of nuclear analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LLNL Mercury Project Trinity Open Science Final Report

The Mercury Monte Carlo particle transport code is used to simulate the transport of radiation through urban environments. These challenging calculations include complicated geometries and require significant computational resources to complete. In the proposed Trinity Open Science calculations, I will investigate computer science aspects of the code which are relevant to convergence of the simulation quantities with increasing Monte Carlo particle counts.

97 MATHEMATICS AND COMPUTING↗

Adhesive-bonded scarf and stepped-lap joints

Continuum mechanics solutions are derived for the static load-carrying capacity of scarf and stepped-lap adhesive-bonded joints. The analyses account for adhesive plasticity and adherend stiffness imbalance and thermal mismatch. The scarf joint solutions include a simple algebraic formula which serves as a close lower bound, within a small fraction of a per cent of the true answer for most practical geometries and materials. Digital computer programs were developed and, for the stepped-lap joints, the critical adherend and adhesive stresses are computed for each step. The scarf joint solutions exhibit grossly different behavior from that for double-lap joints for long overlaps inasmuch as that the potential bond shear strength continues to increase with indefinitely long overlaps on the scarf joints. The stepped-lap joint solutions exhibit some characteristics of both the scarf and double-lap joints. The stepped-lap computer program handles arbitrary (different) step lengths and thickness and the solutions obtained have clarified potentially weak design details and the remedies. The program has been used effectively to optimize the joint proportions.

Hart-Smith, L. J.↗

Fracture of composite orthotropic plates for materials type 2

The fracture problem of laminated plates which consist of orthotropic layers is considered. Symmetrical cracks are located normal to the bimaterial interfaces. The external loads are applied away from the crack region. Three cases are considered: (a) the case of internal cracks; (b) the case of broken laminates; and (c) case of a crack crossing the interface. A general formulation of the problem is given for plane strain and generalized plane stress cases. The singular behavior of stresses at the crack tips and at the interfaces is studied. In each case the stress intensity factors are computed for various crack geometries.

Delale, F.↗

A computer-aided design system geared toward conceptual design in a research environment

A computer-aided design system has recently been developed specifically for the small research group environment. The system is implemented on a Prime 400 minicomputer linked with a CDC 6600 computer. The goal was to assign the minicomputer specific tasks, such as data input and graphics, thereby reserving the large mainframe computer for time-consuming analysis codes. The basic structure of the design system consists of GEMPAK, a computer code that generates detailed configuration geometry from a minimum of input; interface programs that reformat GEMPAK geometry for input to the analysis codes; and utility programs that simplify computer access and data interpretation. The working system has had a large positive impact on the quantity and quality of research performed by the originating group. This paper describes the system, the major factors that contributed to its particular form, and presents examples of its application.

STACK S. H.↗

Use of interactive graphics to analyze QUICK-geometry

The QUICK InterActive Graphics Analysis (QUIAGA) program and its advantages for displaying aircraft QUICK-geometry to aid in detection and analysis of errors are described. The QUICK-geometry system generates a completely analytical aircraft geometry description for use by finite-difference flow codes. The QUIAGA program was developed to exercise the QUICK-geometry subroutines to examine the analytic definition of a configuration by plotting cross sections and body lines on a graphics terminal. A number of options are available, including multiple cross-section views, hidden-line removal, and display of control point locations. Use of these options for the detection and analysis of errors in the QUICK-geometry definition can be of great assistance in speedily arriving at a correct analytical geometry description for flow-field computation. The QUIAGA program has been used in developing a QUICK-geometry model of the NASA Space Shuttle Orbiter, and examples from this experience are given to show some of the program's features. Details of program usage and an example session are given in the appendixes.

Townsend, J. C.↗

Final mission design for IRAS

The Infrared Astronomical Satellite (IRAS) was operated on orbit from 26 January to 22 November 1983 (GMT). Its primary purpose, successfully accomplished, was to conduct an all-sky survey in the infrared from 8 to 120 microns. The mission design for this project featured orbit selection; numerous exercises in the geometry of the sphere; computer simulation of mission, sky, and telescope; and an active interface with the IRAS Science Team. In addition to mission design, the subject of mission planning after launch is addressed. The paper makes extensive use of reference to other works on the topic and seeks to derive some general conclusions on the subject based upon the IRAS experience.

Mclaughlin, W. I.↗

Predictions and measurements of isothermal flowfields in axisymmetric combustor geometries

Numerical predictions, flow visualization experiments and time-mean velocity measurements were obtained for six basic nonreacting flowfields (with inlet swirl vane angles of 0 (swirler removed), 45 and 70 degrees and sidewall expansion angles of 90 and 45 degrees) in an idealized axisymmetric combustor geometry. A flowfield prediction computer program was developed which solves appropriate finite difference equations including a conventional two equation k-epsilon eddy viscosity turbulence model. The wall functions employed were derived from previous swirling flow measurements, and the stairstep approximation was employed to represent the sloping wall at the inlet to the test chamber. Recirculation region boundaries have been sketched from the entire flow visualization photograph collection. Tufts, smoke, and neutrally buoyant helium filled soap bubbles were employed as flow tracers. A five hole pitot probe was utilized to measure the axial, radial, and swirl time mean velocity components.

Rhodes, D. L.↗

A 3-D chimera grid embedding technique

A three-dimensional (3-D) chimera grid-embedding technique is described. The technique simplifies the construction of computational grids about complex geometries. The method subdivides the physical domain into regions which can accommodate easily generated grids. Communication among the grids is accomplished by interpolation of the dependent variables at grid boundaries. The procedures for constructing the composite mesh and the associated data structures are described. The method is demonstrated by solution of the Euler equations for the transonic flow about a wing/body, wing/body/tail, and a configuration of three ellipsoidal bodies.

Benek, J. A.↗

A 3-component laser-Doppler velocimeter data acquisition and reduction system

A laser doppler velocimeter capable of measuring all three components of velocity simultaneously in low-speed flows is described. All the mean velocities, Reynolds stresses, and higher-order products can be evaluated. The approach followed is to split one of the two colors used in a 2-D system, thus creating a third set of beams which is then focused in the flow from an off-axis direction. The third velocity component is computed from the known geometry of the system. The laser optical hardware and the data acquisition electronics are described in detail. In addition, full operating procedures and listings of the software (written in BASIC and ASSEMBLY languages) are also included. Some typical measurements obtained with this system in a vortex/mixing layer interaction are presented and compared directly to those obtained with a cross-wire system.

Rodman, L. C.↗

A 3-component laser-Doppler velocimeter data acquisition and reduction system

This report describes a laser Doppler velocimeter capable of measuring all three components of velocity simultaneously in low-speed flows. All the mean velocities, Reynolds stresses, and higher-order products can then be evaluated. The approach followed is to split one of the colors used in a 2-D system, thus creating a third set of beams which is then focused in the flow from an off-axis direction. The third velocity component is computed from the known geometry of the system. In this report, the laser optical hardware and the data acquisition electronics are described in detail. In addition, full operating procedures and listings of the software (written in BASIC and assembly languages) are also included. Some typical measurements obtained with this system in a vortex/mixing layer interaction are presented and compared directly to those obtained with a cross-wire system.

Rodman, L. C.↗

Aspects of the design and performance of a 3-component LDV system

A conventional two color, two component Laser Doppler Velocimeter (LDV) has been converted into one capable of measuring all three components of velocity simultaneously in low-speed flows. All the Reynolds stresses and higher-order products can then be evaluated. The approach followed is to split one of the two colors and to focus a third pair of beams from an off axis direction. The third velocity component is computed from the known geometry of the system. This paper contains details of the system design with emphasis on the additional problems that are encountered through the introduction of the third component. Some preliminary measurements obtained in a vortex/mixing layer interaction are presented and compared directly with those obtained with a cross-wire system.

Bell, J. H.↗

Supersonic second order analysis and optimization program user's manual

Approximate nonlinear inviscid theoretical techniques for predicting aerodynamic characteristics and surface pressures for relatively slender vehicles at supersonic and moderate hypersonic speeds were developed. Emphasis was placed on approaches that would be responsive to conceptual configuration design level of effort. Second order small disturbance theory was utilized to meet this objective. Numerical codes were developed for analysis and design of relatively general three dimensional geometries. Results from the computations indicate good agreement with experimental results for a variety of wing, body, and wing-body shapes. Case computational time of one minute on a CDC 176 are typical for practical aircraft arrangement.

Clever, W. C.↗