Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal sizing method”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-I-19: Optimization of Low Contrast Detectability Across Two CT Manufacturers

Purpose: Our goal was to evaluate low contrast detectability (LCD) for abdominal CT protocols across two CT scanner manufacturers, while producing a similar noise texture and CTDIvol for acquired images. Methods: A CIRS tissue equivalent LCD phantom containing three columns of 7 spherical targets, ranging from 10 mm to 2.4 mm, that are 5, 10, and 20 HU below the background matrix (HUBB) was scanned using two a GE HD750 64 slice scanner and a Siemens Somatom Definition AS 64 slice scanner. Protocols were designed to deliver a CTDIvol of 12.26 mGy and images were reconstructed with FBP, ASIR and Sapphire. Comparisons were made with those algorithms that had matching noise power spectrum peaks (NPS). NPS information was extracted from a previously published article that matched NPS peak frequencies across manufacturers by calculating the NPS from uniform phantom images reconstructed with several IR algorithms. Results: The minimum detectable lesion size in the 20 HUBB and 10 HUBB column was 6.3 mm, and 10 mm in the 5 HUBB column for the GE HD 750 scanner. The minimum detectable lesion size in the 20 HUBB column was 4.8 mm, in the 10 HUBB column, 9.5 mm, and the 5 HUBB column, 10 mm for the Siemens Somatom Definition AS. Conclusion: Reducing radiation dose while improving or maintaining LCD is possible with application of IR. However, there are several different IR algorithms, with each generating a different resolution and noise texture. In multi-manufacturer settings, matching only the CTDIvol between manufacturers may Result in a loss of clinically relevant information.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

A data-driven framework for error estimation and mesh-model optimization in system-level thermal-hydraulic simulation

We report over the past decades, several computer codes have been developed for simulation and analysis of thermal-hydraulics and system response in nuclear reactors under operating, abnormal transient, and accident conditions. However, simulation errors and uncertainties still inevitably exist even while these codes have been extensively assessed and used. In this work, a data-driven framework (Optimal Mesh/Model Information System, OMIS) is formulated and demonstrated to estimate simulation error and suggest optimal selection of computational mesh size (i.e., nodalization) and constitutive correlations (e.g., wall functions and turbulence models) for low-fidelity, coarse-mesh thermal-hydraulic simulation, in order to achieve accuracy comparable to that of high-fidelity simulation. Using results from high-fidelity simulations and experimental data with many fast-running low-fidelity simulations, an error database is built and used to train a machine learning model that can determine the relationship between local simulation error and local physical features. This machine learning model is then used to generate insight and help correct low-fidelity simulations for similar physical conditions. The OMIS framework is designed as a modularized six-step procedure and accomplished with state-of-the-art methods and algorithms. A mixed-convection case study was performed to illustrate the entire framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Cross-scale efficient tensor contractions for coupled cluster computations through multiple programming model backends

Coupled-cluster methods provide highly accurate models of molecular structure through explicit numerical calculation of tensors representing the correlation between electrons. These calculations are dominated by a sequence of tensor contractions, motivating the development of numerical libraries for such operations. While based on matrix–matrix multiplication, these libraries are specialized to exploit symmetries in the molecular structure and in electronic interactions, and thus reduce the size of the tensor representation and the complexity of contractions. The resulting algorithms are irregular and their parallelization has been previously achieved via the use of dynamic scheduling or specialized data decompositions. We introduce our efforts to extend the Libtensor framework to work in the distributed memory environment in a scalable and energy-efficient manner. We achieve up to 240× speedup compared with the optimized shared memory implementation of Libtensor. We attain scalability to hundreds of thousands of compute cores on three distributed-memory architectures (Cray XC30 and XC40, and IBM Blue Gene/Q), and on a heterogeneous GPU-CPU system (Cray XK7). As the bottlenecks shift from being compute-bound DGEMM's to communication-bound collectives as the size of the molecular system scales, we adopt two radically different parallelization approaches for handling load-imbalance, tasking and bulk synchronous models. Nevertheless, we preserve a unified interface to both programming models to maintain the productivity of computational quantum chemists.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Dislocation half-loop control for optimal V-defect density in GaN-based light emitting diodes

V-defects are morphological defects that typically form on threading dislocations during epitaxial growth of $(0001)$-oriented GaN layers. A V-defect is a hexagonal pyramid-shaped depression with six {$10\bar{1}1$}-oriented sidewalls. These semipolar sidewalls have a lower polarization barrier than the polarization barriers present between the polar c-plane quantum wells and quantum barriers and can laterally inject carriers directly into quantum wells in GaN-based light emitting diodes (LEDs). This is especially important, as the high polarization field in c-plane GaN is a significant factor in the high forward voltage of GaN LEDs. The optimal V-defect density for efficient lateral carrier injection in a GaN LED (∼10 9 cm −2 ) is typically an order of magnitude higher than the threading dislocation density of GaN grown on patterned sapphire substrates (∼10 8 cm −2 ). Pure-edge dislocation loops have been known to exist in GaN, and their formation into large V-defects via low-temperature growth with high Si-doping has recently been studied. Here, in this work, we develop a method for pure-edge threading dislocation half-loop formation and density control via disilane flow, growth temperature, and thickness of the half-loop generation layer. We also develop a method of forming the threading dislocation half-loops into V-defects of comparable size to those originating from substrate threading dislocations.

Crystallographic defects↗
DOE OSTI2024

Integrated floating wind farm layout design and mooring system optimization to increase annual energy production

As we cluster wind turbines in wind farms to gain energy from sites with high wind speeds, wake losses occur within the wind farm. Wake loss is a term used to describe the lower energy production of a downwind turbine that is totally or partially in the wake of an upwind turbine. To decrease wake losses inside the wind farm, the wind farm’s layout is optimized. However, a variety of factors constrain the wind farm layout optimization, such as the size of the lease area relative to the number of turbines to be placed, or the shape of the lease area. Therefore, many wind farms end up with a regular grid layout, such as the Horns Rev 1 wind farm in the North Sea. The ability of a floating offshore wind turbine (FOWT) to change its position based on the wind direction and its mooring system design presents an opportunity to further decrease wake losses in floating wind farms. In this work, we integrate the design of the FOWT mooring systems with the floating wind farm layout design with the goal of increasing the farm’s annual energy production. We use the Horns Rev 1 wind farm as a case study to demonstrate our method. The results show that allowing the FOWT to relocate can decrease wake losses up to 18%. Moreover, the newly developed mooring systems are less stiff and therefore allow larger motion of the FOWT; hence, the material cost of the mooring system decreases by an estimated 17%.

17 WIND ENERGY↗
DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗
DOE OSTI2018

UAVs for Radiological Monitoring: Determining Ideal Operational Parameters From GEANT4 Simulations - 18025

Owing in no small part to the enormous advances in radiation detector technology (size, sensitivity and spectral-resolution), there has been an increased drive and requirement to accurately map both the spatial extent and magnitude of radiological contamination at the greatest possible spatial resolutions. To undertake this, a range of different methods have been devised and include the use of cars and human, but also planes (and more recently) unmanned aerial vehicles (UAVs). However, critical to the remote sensing using such UAVs in contrast to ground-based sensing, is (i) the selection of the ideal operational parameters to attain the best on-ground resolution in addition to, (ii) understanding the correct normalisation procedures and other influences on results that surveys at these reduced altitudes might present. Resulting from earlier works, through employing combined modelling and experimental data, the optimal height for a UAV survey using the more realistic 'area' radioactive source term (opposed to an infinitely-small 'point-source') is determined to be approximately 5 m - suitably maximizing cross detector counts. At greater altitudes than this, the influence of variations in neighboring topology / topography and nearby structures are observed to be negligible. However, at this and further-reduced altitudes (where optimal count-rates exists), small height-related variations are observed to result, and impact considerably on the complexity of the data normalisation. Whilst not occurring as a tool to entirely replace ground-based surveys, UAV systems are able to perform in challenging or dangerous environments, inadvisable or problematic to send humans. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

A new method for in situ structural investigations of nano-sized amorphous and crystalline materials using mixed-flow reactors

Structural investigations of amorphous and nanocrystalline phases forming in solution are historically challenging. Few methods are capable of in situ atomic structural analysis and rigorous control of the system. A mixed-flow reactor (MFR) is used for total X-ray scattering experiments to examine the short- and long-range structure of phases in situ with pair distribution function (PDF) analysis. The adaptable experimental setup enables data collection for a range of different system chemistries, initial supersaturations and residence times. The age of the sample during analysis is controlled by adjusting the flow rate. Faster rates allow for younger samples to be examined, but if flow is too fast not enough data are acquired to average out excess signal noise. Slower flow rates form older samples, but at very slow speeds particles settle and block flow, clogging the system. Proper background collection and subtraction is critical for data optimization. Overall, this MFR method is an ideal scheme for analyzing the in situ structures of phases that form during crystal growth in solution. As a proof of concept, high-resolution total X-ray scattering data of amorphous and crystalline calcium phosphates and amorphous calcium carbonate were collected for PDF analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Corrosion Behavior of Zn/Graphene Oxide Composite Coatings Fabricated by Direct Current Electrodeposition

The zinc (Zn) and Zn/graphene oxide (Zn/GO) composite coatings on the mild steel were fabricated by direct current electrodeposition from an acidic chloride bath without GO and that with four different concentrations (50, 100, 200 and 300 mg/L) of GO, respectively. The GO first is synthesized by using the typical Hummer method. Physical characterizations including the surface and cross-sectional morphology, chemical composition and crystal structure of the Zn and Zn/GO composite coatings were done by scanning electron microscope, energy-dispersive x-ray and x-ray diffraction. The crystalline size, texture coefficient, water contact angle and hardness of the Zn and Zn/GO composite coatings were measured. The corrosion behavior of the Zn and Zn/GO composite coatings was studied by the electrochemical corrosion tests including the potentiodynamic polarization test and the electrochemical impedance spectroscopy test and the immersion test. The result reveals that the corrosion resistance of the Zn coating is remarkably increased by the incorporation of GO. Furthermore, it does not show a linear increase with the increase in the amount of incorporated GO and an optimal corrosion resistance is provided by the Zn/GO composite coating attained from the electrolyte containing 100 mg/L GO. The influence of the incorporation of GO on the corrosion resistance of the Zn coating is explained in combination with their physical characterizations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Identifying the Atomic-Level Effects of Metal Composition on the Structure and Catalytic Activity of Peptide-Templated Materials

Bioinspired approaches for the formation of metallic nanomaterials have been extensively employed for a diverse range of applications including diagnostics and catalysis. These materials can often be used under sustainable conditions; however, it is challenging to control the material size, morphology, and composition simultaneously. Here we have employed the R5 peptide, which forms a 3D scaffold to direct the size and linear shape of bimetallic PdAu nanomaterials for catalysis. The materials were prepared at varying Pd:Au ratios to probe optimal compositions to achieve maximal catalytic efficiency. These materials were extensively characterized at the atomic level using transmission electron microscopy, extended X-ray absorption fine structure spectroscopy, and atomic pair distribution function analysis derived from high-energy X-ray diffraction patterns to provide highly resolved structural information. The results confirmed PdAu alloy formation, but also demonstrated that significant surface structural disorder was present. The catalytic activity of the materials was studied for olefin hydrogenation, which demonstrated enhanced reactivity from the bimetallic structures.These results present a pathway to the bioinspired production of multimetallic materials with enhanced properties, which can be assessed via a suite of characterization methods to fully ascertain structure/function relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Cross-scale Efficient Tensor Contractions for Coupled Cluster Computations Through Multiple Programming Model Backends

Coupled-cluster methods provide highly accurate models of molecular structure by explicit numerical calculation of tensors representing the correlation between electrons. These calculations are dominated by a sequence of tensor contractions, motivating the development of numerical libraries for such operations. While based on matrix-matrix multiplication, these libraries are specialized to exploit symmetries in the molecular structure and in electronic interactions, and thus reduce the size of the tensor representation and the complexity of contractions. The resulting algorithms are irregular and their parallelization has been previously achieved via the use of dynamic scheduling or specialized data decompositions. We introduce our efforts to extend the Libtensor framework to work in the distributed memory environment in a scalable and energy efficient manner. We achieve up to 240 speedup compared with the best optimized shared memory implementation. We attain scalability to hundreds of thousands of compute cores on three distributed-memory architectures, (Cray XC30&XC40, BlueGene/Q), and on a heterogeneous GPU-CPU system (Cray XK7). As the bottlenecks shift from being compute-bound DGEMM's to communication-bound collectives as the size of the molecular system scales, we adopt two radically different parallelization approaches for handling load-imbalance. Nevertheless, we preserve a uni ed interface to both programming models to maintain the productivity of computational quantum chemists.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Image-Based Fracture Surface Defect Characterization Methods for Additively Manufactured Ti-6Al-4V Tested in Fatigue

Abstract Fatigue initiation in additively manufactured samples/parts often occurs at processed-induced defects such as lack-of-fusion (LoF), keyhole, or other morphological/microstructural defects that have unique characteristics and measurable qualities. Attempts at identifying and minimizing such defects have utilized optimized processing conditions along with in situ and ex situ characterization that includes metallography and/or X-ray computed tomography (XCT). This paper highlights the benefits of using fracture surface analyses to detect and quantify defects that may not be detected by metallography/XCT due to sectioning and resolution limits. In addition to using manual quantification of fatigue initiating LoF and keyhole defects on fracture surfaces, image-based machine learning using convolutional neural networks such as U-Net were also used to automate the process. Statistical analyses were used to identify the extreme cases of defects that initiated and accelerated fatigue and to model the distribution of defect size and shape characteristics to distinguish the type of defect. Initial results show agreement between trained machine learning models and ground truth data in defect segmentation, and the distributions of defect characteristics are distinguishable to particular process-induced defect types.

Materials Science↗
DOE OSTI2015

MO-FG-204-06: A New Algorithm for Gold Nano-Particle Concentration Identification in Dual Energy CT

Purpose: Gold nano-particle (GNP) has recently attracted a lot of attentions due to its potential as an imaging contrast agent and radiotherapy sensitiser. Imaging the GNP at its low contraction is a challenging problem. We propose a new algorithm to improve the identification of GNP based on dual energy CT (DECT). Methods: We consider three base materials: water, bone, and gold. Determining three density images from two images in DECT is an under-determined problem. We propose to solve this problem by exploring image domain sparsity via an optimization approach. The objective function contains four terms. A data-fidelity term ensures the fidelity between the identified material densities and the DECT images, while the other three terms enforces the sparsity in the gradient domain of the three images corresponding to the density of the base materials by using total variation (TV) regularization. A primal-dual algorithm is applied to solve the proposed optimization problem. We have performed simulation studies to test this model. Results: Our digital phantom in the tests contains water, bone regions and gold inserts of different sizes and densities. The gold inserts contain mixed material consisting of water with 1g/cm3 and gold at a certain density. At a low gold density of 0.0008 g/cm3, the insert is hardly visible in DECT images, especially for those with small sizes. Our algorithm is able to decompose the DECT into three density images. Those gold inserts at a low density can be clearly visualized in the density image. Conclusion: We have developed a new algorithm to decompose DECT images into three different material density images, in particular, to retrieve density of gold. Numerical studies showed promising results.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

U.S. Solar Photovoltaic System and Energy Storage Cost Benchmark: Q1 2021 [Slides]

Our benchmarking method includes bottom-up accounting for all necessary system and projectdevelopment costs incurred when installing residential, commercial, and utility-scale systems, and it models the Q1 2021 costs for such systems, excluding any previous supply agreements or contracts. In general, we attempt to model the typical installation techniques and business operations from an installed-cost perspective, and our benchmarks are national averages. The residential PV-only benchmark and the commercial rooftop PV-only benchmark average costs by inverter type (string inverters, string inverters with direct current [DC] optimizers, and microinverters), weighted by inverter market share. The residential PV-only benchmark is further averaged across small installer and national integrator business models, weighted by market share. All benchmarks include variations—accounting for the differences in size, equipment, and operational use (particularly for storage)—that are currently available in the marketplace. All benchmarks assume nonunionized construction labor; residential and commercial PV systems predominantly use nonunionized labor, and the type of labor required for utility-scale PV systems depends heavily on the development process. All benchmarks assume the use of monofacial monocrystalline silicon PV modules. Benchmarks using cadmium telluride or bifacial modules could result in significantly different results. The data in this annual benchmark report inform the formulation of and track progress toward the U.S. Department of Energy (DOE) Solar Energy Technologies Office’s Government Performance and Reporting Act cost targets.

14 SOLAR ENERGY↗
DOE OSTI2019

Multidimensional sum-up rounding for integer programming in optimal experimental design

Here we present a numerical method for approximating the solution of convex integer programs stemming from optimal experimental design. The statistical setup consists of a Bayesian framework for linear inverse problems for which the direct relationship is described by a discretized integral equation. Specifically, we aim to find the optimal sensor placement from a set of candidate locations where data are collected with measurement error. The convex objective function is a measure of the uncertainty, described here by the trace or log-determinant of the posterior covariance matrix, for the discretized linear inverse problem solution. The resulting convex integer program is relaxed, producing a lower bound. An upper bound is obtained by extending the sum-up rounding approach to multiple dimensions. For this extension, we analyze its accuracy as a function of the discretization mesh size for a rectangular domain. We show asymptotic optimality of the integer solution defining the upper bound for different experimental design criteria (A- and D-optimal), by proving the convergence to zero of the gap between the upper and lower bounds as the mesh size goes to zero. The technique is illustrated on a two-dimensional gravity surveying problem for both A-optimal and D-optimal sensor placement where our designs yield better results compared with a thresholding rounding approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Preface to the Special Issue on TOUGH Symposium 2015

The TOUGH Symposium 2015 was held in Berkeley, California, September 28-30, 2015. The TOUGH family of codes, developed at the Energy Geosciences Division of Lawrence Berkeley National Laboratory (LBNL), is a suite of computer programs for the simulation of multiphase and multicomponent fluid and heat flows in porous and fractured media with applications in many geosciences fields, such as geothermal reservoir engineering, nuclear waste disposal, geological carbon sequestration, oil and gas reservoirs, gas hydrate research, vadose zone hydrology and environmental remediation. Since the first release in the 1980s, many modifications and enhancements have been continuously made to TOUGH and its various descendants (iTOUGH2, TOUGH+, TOUGH-MP, TOUGHREACT, TOUGH+HYDRATE, TMVOC...), at LBNL and elsewhere. Today, these codes are used worldwide in academia, government organizations and private companies in problems involving coupled hydrological, thermal, biogeochemical and geomechanical processes. The Symposia, organized every 2-3 years, bring together developers and users for an open exchange on recent code enhancements and applications. In 2015, the Symposium was attended by one hundred participants, representing thirty-four nationalities. This Special Issue in Computers & Geosciences gathers extended versions of selected Symposium proceedings related to (i) recent enhancements to the TOUGH family of codes and (ii) coupled flow and geomechanics processes modeling. In this special issue, a paper by Jung et al. provides a new base version of TOUGH, TOUGH3. This latest version addresses the computational robustness and efficiency needed for the increasing complexity of the simulated processes and the growing size of model domains. Moreover, it consolidates the serial and parallel versions into a single code, hence, inconsistencies between different versions are less likely to occur, and code maintenance and development become more efficient. Finsterle et al. present a review of the capabilities of iTOUGH2 and illustrate the workflow of an integrated simulation-optimization framework. Liu et al. present an approach to reduce the computational cost of performing Bayesian inverse modeling and demonstrate their approach using a synthetic experiment simulated with TOUGH2. Regarding the mesh generation process, the paper by Doughty presents a generalization of the grid-generation algorithm embedded in TOUGH2 (MESHMAKER) to create one-column grids with fractal dimension, useful for representing flow through fractures or highly heterogeneous media with fractal geometry, and also for inverse methods. Sentís and Gable present a module to interface the open source code Lagrit with TOUGH2, that assists in the processes of mesh generation and model setup, and that can be particularly useful for complex geometries and large models. Bonduà et al. present a pre-processing tool (VORO2MESH) that features a full 3D Voronoi approach for structured and unstructured grids, and an enhanced post-processor (TOUGH2Viewer) for the visualization of 3D Voronoi grids. The literature review shows that coupled fluid flow and geomechanics has gained increasing attention in the last years in many geosciences fields. In this context, several codes of the TOUGH family have been successfully coupled with geomechanics software. The paper by Rutqvist presents a review of some of these TOUGH-based geomechanics models, of which the TOUGH-FLAC simulator – presented for the first time in the TOUGH Symposium 2003 – is well-known worldwide and has been used in a large variety of applications, encompassing the application domains of the TOUGH codes. Blanco-Martín et al. present an extension of TOUGH-FLAC to the finite strain framework, motivated by the existence of large strains associated with rock salt deformation in the context of heat-generating nuclear waste disposal. Following a paper published in the Special Issue of the TOUGH Symposium 2012 in this journal, Kim et al. present the current status of the TOUGH-RSBN simulator for coupled hydro-mechanical modeling of hydraulic fracturing using a discrete fracture network approach, and present some validation and demonstration cases. Rinaldi and Nespoli present TOUGH-SEED, an injection-induced seismicity modeling approach using TOUGH2 and a geomechanics stochastic model, and apply it to the Basel EGS project. Finally, Rinaldi et al. use the PEST protocol implemented in iTOUGH2, combined with TOUGH-FLAC as the forward model, to perform inverse modeling for the In Salah CO2 demonstration site in Algeria. Using pore pressures and InSAR displacements as observations, this approach allows for a multi-physics joint inversion and takes advantage of the inverse capabilities of iTOUGH2 (parameter estimation, sensitivity analysis, uncertainty estimation …), applied to coupled fluid flow and geomechanics. On behalf of the members of the TOUGH Symposium 2015 Organizing Committee, I hope that this collection of papers will provide the readers of Computers & Geosciences with a recent view of the capabilities of the TOUGH family of codes, as well as with current modeling research domains in the geosciences. Last but not least, I would like to express my gratitude to all the reviewers who kindly contributed to this Special Issue, as well as the Journal staff who have made this publication possible.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2015

Numerical modeling of the SNS H{sup −} ion source

Ion source rf antennas that produce H- ions can fail when plasma heating causes ablation of the insulating coating due to small structural defects such as cracks. Reducing antenna failures that reduce the operating capabilities of the Spallation Neutron Source (SNS) accelerator is one of the top priorities of the SNS H- Source Program at ORNL. Numerical modeling of ion sources can provide techniques for optimizing design in order to reduce antenna failures. There are a number of difficulties in developing accurate models of rf inductive plasmas. First, a large range of spatial and temporal scales must be resolved in order to accurately capture the physics of plasma motion, including the Debye length, rf frequencies on the order of tens of MHz, simulation time scales of many hundreds of rf periods, large device sizes on tens of cm, and ion motions that are thousands of times slower than electrons. This results in large simulation domains with many computational cells for solving plasma and electromagnetic equations, short time steps, and long-duration simulations. In order to reduce the computational requirements, one can develop implicit models for both fields and particle motions (e.g. divergence-preserving ADI methods), various electrostatic models, or magnetohydrodynamic models. We have performed simulations using all three of these methods and have found that fluid models have the greatest potential for giving accurate solutions while still being fast enough to perform long timescale simulations in a reasonable amount of time. We have implemented a number of fluid models with electromagnetics using the simulation tool USim and applied them to modeling the SNS H- ion source. We found that a reduced, single-fluid MHD model with an imposed magnetic field due to the rf antenna current and the confining multi-cusp field generated increased bulk plasma velocities of > 200 m/s in the region of the antenna where ablation is often observed in the SNS source. We report here on comparisons of simulated plasma parameters and code performance using more accurate physical models, such as two-temperature extended MHD models, for both a related benchmark system describing a inductively coupled plasma reactor, and for the SNS ion source. We also present results from scaling studies for mesh generation and solvers in the USim simulation code.

07 ISOTOPE AND RADIATION SOURCES↗