Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “permit”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Geologic carbon sequestration injection wells in overpressured storage reservoirs: estimating area of review

The Area of Review (AoR) under the US Environmental Protection Agency's (EPA) Class VI CO 2 injection permit is defined as the region surrounding the geologic carbon sequestration (GCS) project where underground sources of drinking water (USDWs) may be endangered. Estimation of the AoR is based on the calculated reservoir pressurization due to CO 2 injection and the associated potential to lift saline water into potable groundwater aquifers through open flow paths (e.g. wells) assuming the system is hydrostatic. In cases where the storage reservoirs are not initially hydrostatic, and in particular where they are overpressured, AoR estimation methods need to be altered. Here in this paper, we present and apply an approach to evaluating potential endangerment of USDW based on comparing brine leakage through a hypothetical open flow path in a no-injection scenario and brine leakage in a CO 2 -injection scenario. We present six possible ways to normalize injection-related leakage relative to no-injection leakage. We calculate leakage using semi-analytical solutions for single-phase flow and model reservoir pressurization and flow up (single) leaky wells located progressively farther from the injection well. For an example case of relative overpressure and using an injection-rate-based approach, results show 50-60% larger open-well-leakage rates for wells located at 2 km and 10% increase for wells located at 10 km from the injection well relative to the no-injection case. If total brine leakage is considered, the results depend strongly on the assumed pre-injection to post-injection time frames and on the methods of normalization used to calculate incremental leakage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mapping existing wellbore locations to compare technical risks between onshore and offshore CCS activities in Texas

Carbon dioxide capture and geologic storage (CCS; geologic sequestration) is a promising technology for reducing anthropogenic greenhouse gas emissions to the atmosphere from industrial point sources. Aspects of CCS have been investigated for over two decades, and many large- and small-scale geologic storage field demonstration projects are now underway globally. Interest in offshore CCS has been increasing in recent years (e.g., European Union, Australia, Japan, and the United States). Deep geologic storage in offshore settings is analogous to onshore CCS activities in many respects (i.e., geologic and geotechnical aspects), but is distinct from previously explored seabed sediment CO 2 storage) or deep marine dissolution). Given the large subsurface geologic storage volumes available in offshore settings, much discussion of offshore CCS is focused on the benefits and risks of such activity compared to onshore settings. Similar to onshore settings, existing (legacy) wells likely present the most direct migration pathway and largest risk of noncontainment in offshore settings. As part of current studies to evaluate geologic storage options in offshore settings along the Texas coast and greater Gulf of Mexico (GoM), mapping of the geographic distribution and ages of wells in a region containing coastal counties and extending 30 miles offshore Texas indicates that both well spatial density and well age decrease moving from onshore to offshore. Results suggest reduced risk of leakage owing to more rigorous and documented well completion and abandonment practices for these generally younger wells (although many are decades old). A result of decreased well density is that larger areas are available for leasing for CCS projects that avoid legacy wells altogether (> 1 mile from any existing well). The one-mile designation is used as an arbitrary convention, and while it is recognized that this is smaller than a typical area of review (AoR) for permitting, each site will have a different AoR radius for consideration. Finally, the combination of large subsurface storage volumes under control of a single landowner and reduced risks from legacy wells makes offshore CCS attractive in the GoM.

58 GEOSCIENCES↗
DOE OSTI2012

Stereoselection in the diels alderase ribozyme: A molecular dynamics study

The Diels-Alderase ribozyme is an in vitro -evolved ribonucleic acid enzyme that catalyzes a [4 + 2] cycloaddition reaction between an anthracene diene and a maleimide dienophile. The ribozyme can in principle be used to selectively synthesize only one product enantiomer, depending on which of the two entrances to the catalytic pocket, front or back , the substrate is permitted to use. Here, we investigate stereoselection and substrate recognition in the ribozyme by means of multiple molecular dynamics simulations, performed on each of the two substrates individually in the pocket, on the reactant state, and on the product state. The results are consistent with a binding mechanism in which the maleimide likely binds first followed by the anthracene, which enters preferentially through the front door. The free energy profiles for anthracene binding indicate that the pre-(R,R)-enantiomer conformation is slightly preferred, in agreement with the experimentally observed small enantiomeric excess of the (R,R)-enantiomer of the product. In conclusion, the reactant state is stabilized by the simultaneous presence of both substrates bound to their binding sites in the hydrophobic pocket as well as by stacking interactions between them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Retrospective evaluation of dosimetric quality for prostate carcinomas treated with 3D conformal, intensity modulated and volumetric modulated arc radiotherapy

This study examines and compares the dosimetric quality of radiotherapy treatment plans for prostate carcinoma across a cohort of 163 patients treated across five centres: 83 treated with three-dimensional conformal radiotherapy (3DCRT), 33 treated with intensity modulated radiotherapy (IMRT) and 47 treated with volumetric modulated arc therapy (VMAT). Treatment plan quality was evaluated in terms of target dose homogeneity and organs at risk (OAR), through the use of a set of dose metrics. These included the mean, maximum and minimum doses; the homogeneity and conformity indices for the target volumes; and a selection of dose coverage values that were relevant to each OAR. Statistical significance was evaluated using two-tailed Welch's T-tests. The Monte Carlo DICOM ToolKit software was adapted to permit the evaluation of dose metrics from DICOM data exported from a commercial radiotherapy treatment planning system. The 3DCRT treatment plans offered greater planning target volume dose homogeneity than the other two treatment modalities. The IMRT and VMAT plans offered greater dose reduction in the OAR: with increased compliance with recommended OAR dose constraints, compared to conventional 3DCRT treatments. When compared to each other, IMRT and VMAT did not provide significantly different treatment plan quality for like-sized tumour volumes. This study indicates that IMRT and VMAT have provided similar dosimetric quality, which is superior to the dosimetric quality achieved with 3DCRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Supercell precipitation contribution to the United States hydroclimate

Abstract This research seeks to understand simulated supercell precipitation characteristics across the conterminous United States (CONUS) using high‐resolution, convection‐permitting, dynamically downscaled simulations for three 15‐year epochs. Epochs include a historical end‐of‐20th‐century period (1990–2005) and two end‐of‐21st‐century (2085–2100) scenarios for intermediate and pessimistic greenhouse gas concentration trajectories. Simulated updraft helicity, which measures the corkscrew flow within a storm's updraft, is used as a proxy for supercells. An algorithm tracks and catalogues updraft helicity swaths that, when buffered, are used to acquire simulated precipitation from supercells. The historical epoch provides a baseline climatology of supercell precipitation for a contemporary climate, which is then compared against the two future epochs to assess how supercell precipitation may change during the 21st century. Despite their relatively small size, supercells provide critical precipitation to the Wheat and Corn Belts, large expanses of CONUS pasture and rangeland, regional aquifers and several large river basins. Many areas in the central CONUS receive upwards of 3%–6% of their annual and 5%–8% of their warm‐season precipitation from these storms. Results suggest that precipitation contribution from supercells will decrease in the future across most of the High Plains and Central and Northern Great Plains with robust increases likely across the south‐central and Southeast regions. Supercell precipitation rates are expected to increase for large portions of the CONUS by the end‐of‐the‐21st‐century, suggesting a growing threat for flash floods from these storms as they become more efficient precipitation producers. This research provides an initial perspective on the magnitude of supercell precipitation and potential changes to this important hydrologic input to assist water‐sensitive industries, private and public insurance markets, agriculture entities, as well as inform plans to mitigate and build resilience to rapid environmental and societal change.

Meteorology & Atmospheric Sciences↗
DOE OSTI2019

Spatially tracking carbon through the root-rhizosphere-soil system using laser ablation-IRMS

The intimate relationships between plant roots, rhizosphere, and soil are fostered by the release of organic compounds from the plant into soil through various forms of rhizodeposition and the simultaneous harvestingof inorganic nutrients from the soil to the plant. Here we present a method to spatially track and map the migration of plant-derived carbon (C) through roots into the rhizosphere and surrounding soil using laser ablation-isotope ratio mass spectrometry (LA-IRMS). Weused switchgrass microcosms containing soil from field plots at the Kellogg Biological Station (Hickory Corners, Michigan, USA), which have been cropped with switchgrass since 2008. We used a 13 CO 2 tracerto isotopically label switchgrass plants for two diel cycles and tracked subsequent movement of labeled C using the spatially specific (<100 µm resolution) δ 13 C analysis enabled by LA-IRMS.This approach permitted assessment of variable C flow through different roots and enabled mapping of spatial variability of C allocation to the rhizosphere. Highly 13 C-enriched C (consistent with production during the 13 CO 2 application period) extended ~0.5 – 1 mm from the root into the soil, suggesting that the majority of recent plant-derived C was within this distance of the root after 48 hours. Tracking the physical extent of root exudation into the rhizosphere can help evaluate the localization of plant-microbe interactions in highly variable subsurface environments,and the use of the isotopic label can differentiate freshly fixed C (presumably from root exudates) from other types of subsurface C (e.g., plant necromass and microbial turnover). The LA-IRMS technique mayalso serve as a valuable screening technique to identify areas of high activity for additional microbial or geochemicalassays.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Laser-induced annealing of aged PuO 2

Plutonium dioxide (PuO 2 ) is an important compound used in nuclear fuel, irradiation targets, and heat sources. As such, improved understanding of its structural and spectroscopic properties has numerous applications. Alpha particle-induced damage of a PuO 2 crystal lattice modifies several properties of its Raman spectrum, including band intensities, positions, and widths. The decay also induces growth of new bands and creates electron-trapped defects with luminescent properties. Herein, we show for the first time that laser-induced heating can reverse damage to an aged and damaged PuO 2 lattice. Using automated instrumentation to heat a single 6–10 μm spot to temperatures above 1,300°C, we show that laser-induced annealing of aged PuO 2 results in restoration of T2g band intensity with a concomitant decrease in defect band intensity, a result that indicates laser annealing can be used to reverse age-damage in PuO 2 . This laser annealing approach permitted in situ observation of temperature-dependent Raman spectral changes, thereby providing insight into the thermodynamics of structural alterations in a radiolytically damaged PuO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Advancing wildlife connectivity in land use planning: a case study with four-toed salamanders

Stable habitat connections that wildlife can safely traverse are essential to biodiversity conservation and healthy ecosystems. We developed high-resolution landscape connectivity models to predict resistance to movement by a threatened wetland-obligate amphibian, the four-toed salamander ( Hemidactylium scutatum ), and identified priority management areas on the 13,000-ha Department of Energy Oak Ridge Reservation (ORR) from 2019 to 2022. We developed a resistance surface based on aerial light detection and ranging data (LiDAR), >30 years of field-based mapping of forest, hydrologic, and geologic features, and contemporary population surveys, alongside derived predictors at <1-m resolution. We then modeled predicted movement corridors using a circuit theory-based modeling approach. We worked closely with land management and natural resources personnel to integrate ecological modeling with broader land use priorities, monetary costs, and feasibility. We identified important terrestrial and aquatic areas on ORR and simulated management scenarios to promote stable connections for four-toed salamanders. This approach allowed us to narrow down a list of 438 potential habitat manipulation sites to 10 sites where open-bottomed culverts and buffers could be implemented. This smaller-scale restoration approach produced a similar increase in landscape connectivity while costing <20% of a larger-scale approach based on barrier removal. We successfully identified feasible, cost-effective management strategies that integrated knowledge from a variety of sources. In conclusion, we offer a strategy that permitted integration of wildlife management goals into infrastructure upgrades wherein wildlife was not an initial consideration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Accelerating technology development to monitor and minimize effects from land‐based wind energy on birds and bats

While wind energy is a key sector of domestic energy production for the United States, operation of wind turbines directly and indirectly adversely affects certain species of birds and bats. The cumulative effect of wind turbine strikes can have both biological and regulatory consequences, and, in some cases, delay permitting and construction or affect ongoing operations. Technology can help quantify and minimize these effects, but the pace of development, acceptance, and adoption of technological solutions is slow. Although adopting cost‐effective technologies may reduce negative effects on wildlife and help achieve both energy production and conservation goals, consensus is lacking among developers, regulators, and the conservation community regarding how to define technology effectiveness and acceptance and how to develop a standardized process for doing so. Removing barriers to technology advancement requires deviating from the status quo. Changes include 1) creating incentives to mitigate impacts, 2) establishing options for research as mitigation, 3) rethinking how research is funded, 4) increasing stakeholder coordination, and 5) increasing the efficiency of research and development. We recommend the creation of a national framework to establish clear criteria and protocols for technology evaluation and adoption.

17 WIND ENERGY↗
DOE OSTI2017

Hydrothermal Conditioning of Physical Hydrogels Prepared from a Midblock-Sulfonated Multiblock Copolymer

Since nanostructured amphiphilic macromolecules capable of affording high ion and water transport are becoming increasingly important in a wide range of contemporary energy and environmental technologies, the swelling kinetics and temperature dependence of water uptake are investigated in a series of midblock-sulfonated thermoplastic elastomers. Upon self-assembly, these materials maintain a stable hydrogel network in the presence of a polar liquid. In this study, real-time water-sorption kinetics in copolymer films prepared by different casting solvents are elucidated by synchrotron small-angle X-ray scattering and gravimetric measurements, which directly correlate nanostructural changes with macroscopic swelling to establish fundamental structure-property behavior. By monitoring the equilibrium swelling capacity of these materials over a range of temperatures, an unexpected transition in the vicinity of 50 degrees C has been discovered. Furthermore, depending on copolymer morphology and degree of sulfonation, hydrothermal conditioning of specimens to temperatures above this transition permits retention of superabsorbent swelling at ambient temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Quantitative impact of changes in marrow cellularity, skeletal size, and bone mineral density on active marrow dosimetry based upon a reference model

Purpose The hematopoietically active tissues of skeletal bone marrow are a prime target for computational dosimetry given potential risks of leukemia and, at higher dose levels, acute marrow toxicity. The complex three‐dimensional geometry of trabecular spongiosa, however, complicates schema for dose assessment in such a way that only a few reference skeletal models have been developed to date, and which are based upon microimaging of a limited number of cadaveric bone spongiosa cores. The question then arises as to what degree of accuracy is achievable from reference skeletal dose models when applied to individual patients or specific exposed populations? Methods Patient variability in marrow dosimetry were quantified for three skeletal sites — the ribs, lumbar vertebrae, and cranium — for the beta‐emitters 45 Ca, 153 Sm, and 90 Y, and the alpha‐particle emitters 223 Ra, 219 Rn, and 215 Po, the latter two being the immediate progeny of the former. For each radionuclide and bone site, three patient parameters were altered from their values in the reference model: (1) bone size as a surrogate for patient stature, (2) marrow cellularity as a surrogate for age‐ or disease‐related changes in marrow adiposity, and (3) the trabecular bone volume fraction as a surrogate for bone mineral density. Marrow dose variability is expressed as percent differences in the radionuclide S value given by the reference model and the patient‐parameterized model. The impact of radionuclide biokinetics on marrow dosimetry was not considered. Results Variations in overall bone size play a very minor role in active marrow dose variability. Marrow cellularity is a significant factor in dose variability for active marrow self‐irradiation, but it plays no role for radionuclides localized to the trabecular bone matrix. Variations in trabecular bone volume fractions impact the active marrow dose variability for short‐range particle emitters 45 Ca, 223 Ra, 219 Rn, and 215 Po in the vertebrae and ribs, skeletal sites with small spongiosa proportions of trabecular bone. In the cranium, with its relative high proportion of trabecular bone, significant differences in marrow dosimetry from the reference model were noted for all radionuclides. Conclusions Skeletal models of active marrow dosimetry should be more fully parameterized to permit closer matching to patient bone density and marrow cellularity, particularly when considering short‐range particle emitters localized to either the bone trabeculae or active marrow, respectively.

Geyer, Amy M.↗
DOE OSTI2019

Parallel beamlet dose calculation via beamlet contexts in a distributed multi‐GPU framework

Purpose Dose calculation is one of the most computationally intensive, yet essential tasks in the treatment planning process. With the recent interest in automatic beam orientation and arc trajectory optimization techniques, there is a great need for more efficient model‐based dose calculation algorithms that can accommodate hundreds to thousands of beam candidates at once. Foundational work has shown the translation of dose calculation algorithms to graphical processing units ( GPU s), lending to remarkable gains in processing efficiency. But these methods provide parallelization of dose for only a single beamlet, serializing the calculation of multiple beamlets and under‐utilizing the potential of modern GPU s. In this paper, the authors propose a framework enabling parallel computation of many beamlet doses using a novel beamlet context transformation and further embed this approach in a scalable network of multi‐ GPU computational nodes. Methods The proposed context‐based transformation separates beamlet‐local density and TERMA into distinct beamlet contexts that independently provide sufficient data for beamlet dose calculation. Beamlet contexts are arranged in a composite context array with dosimetric isolation, and the context array is subjected to a GPU collapsed‐cone convolution superposition procedure, producing the set of beamlet‐specific dose distributions in a single pass. Dose from each context is converted to a sparse representation for efficient storage and retrieval during treatment plan optimization. The context radius is a new parameter permitting flexibility between the speed and fidelity of the dose calculation process. A distributed manager‐worker architecture is constructed around the context‐based GPU dose calculation approach supporting an arbitrary number of worker nodes and resident GPU s. Phantom experiments were executed to verify the accuracy of the context‐based approach compared to Monte Carlo and a reference CPU ‐ CCCS implementation for single beamlets and broad beams composed by addition of beamlets. Dose for representative 4 π beam sets was calculated in lung and prostate cases to compare its efficiency with that of an existing beamlet‐sequential GPU ‐ CCCS implementation. Code profiling was also performed to evaluate the scalability of the framework across many networked GPU s. Results The dosimetric accuracy of the context‐based method displays <1.35% and 2.35% average error from the existing serialized CPU ‐ CCCS algorithm and Monte Carlo simulation for beamlet‐specific PDD s in water and slab phantoms, respectively. The context‐based method demonstrates substantial speedup of up to two orders of magnitude over the beamlet‐sequential GPU ‐ CCCS method in the tested configurations. The context‐based framework demonstrates near linear scaling in the number of distributed compute nodes and GPU s employed, indicating that it is flexible enough to meet the performance requirements of most users by simply increasing the hardware utilization. Conclusions The context‐based approach demonstrates a new expectation of performance for beamlet‐based dose calculation methods. This approach has been successful in accelerating the dose calculation process for very large‐scale treatment planning problems ‐ such as automatic 4 π IMRT beam orientation and VMAT arc trajectory selection, with hundreds of thousands of beamlets ‐ in clinically feasible timeframes. The flexibility of this framework makes it as a strong candidate for use in a variety of other very large‐scale treatment planning tasks and clinical workflows.

Neph, Ryan↗
DOE OSTI2018

Passage Route and Upstream Migration Success: A Case Study of Snake River Salmonids Ascending Lower Granite Dam

Abstract Lower Granite Dam is the last dam that federally protected Snake River salmonids Oncorhynchus spp. must ascend during their spawning migration. The dam has an adult fish ladder equipped with a trapping system to facilitate fisheries research and hatchery broodstock collection. There are three possible passage routes through the adult ladder: trapped, shunted, and free passage. During the adult trapping season, all fish must swim through 0.305-m shunt pipes outfitted with PIT tag arrays that allow the selection of fish for trapping. Selected fish use the “trapped” route and are kept in a holding area for up to 20 h before being sampled and returned to the ladder. Unselected fish use the “shunted” route and immediately resume upstream migration after swimming through the pipes. When the trap is not in operation, the shunted route is inaccessible, and all fish use the “free passage” route to ascend the ladder without additional impediment. In 2016, a temporary change in ladder operations permitted free passage for a portion of the trapping season. Our study used this rare opportunity to evaluate how different passage routes affect in-ladder transit time and upstream homing success for five salmonid stocks: Sockeye Salmon Oncorhynchus nerka; steelhead O. mykiss; and spring-, summer-, and fall-run Chinook Salmon O. tshawytscha. In 2016, only Sockeye Salmon and spring- and summer-run Chinook Salmon were given access to free passage, and we found evidence that free passage increased subsequent detection at natal sites upstream. An expanded analysis of shunted versus trapped fish during the years 2012–2016 found no difference in rates of detection to home tributaries by route of passage for any of the five fish stocks examined.

Morrisett, Christina N.↗
DOE OSTI2025

A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints

Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.

PDE-constrained optimization↗
DOE OSTI2014

Cyclic steady states of treaded rolling bodies

The analysis of spinning axisymmetric bodies undergoing finite deformation is useful for understanding the behavior of a wide variety of engineering systems - for example, rolling tires. Furthermore, progress beyond the axisymmetric case has been lacking. In this work, we present a methodology for the treatment of treaded rolling bodies by devising a Newton-Krylov shooting scheme for the computation of cyclic steady states of motion. The scheme advocated permits one to determine the entire transient (cyclically steady) motion of a spinning treaded body with consideration of viscoelasticity as well as contact. We demonstrate the viability of the method through three examples, (1) a viscoelastic cylinder with eight sinusoidal tread blocks, (2) a viscoelastic cylinder with 16 square tread blocks, and (3) a viscoelastic oval in which repeated contact and separation from a rigid plane excites distinct transient vibrational modes in the body during each cycle. For practical values of the viscoelastic relaxation time, the new approach is found to converge faster than the naïve approach of evolving an initial guess over many cycles until a cyclic steady state is reached.

42 ENGINEERING↗
DOE OSTI2016

An n -material thresholding method for improving integerness of solutions in topology optimization

It is common in solving topology optimization problems to replace an integer-valued characteristic function design field with the material volume fraction field, a real-valued approximation of the design field that permits "fictitious" mixtures of materials during intermediate iterations in the optimization process. This is reasonable so long as one can interpolate properties for such materials and so long as the final design is integer valued. For this purpose, we present a method for smoothly thresholding the volume fractions of an arbitrary number of material phases which specify the design. This method is trivial for two-material design problems, for example, the canonical topology design problem of specifying the presence or absence of a single material within a domain, but it becomes more complex when three or more materials are used, as often occurs in material design problems. We take advantage of the similarity in properties between the volume fractions and the barycentric coordinates on a simplex to derive a thresholding, method which is applicable to an arbitrary number of materials. As we show in a sensitivity analysis, this method has smooth derivatives, allowing it to be used in gradient-based optimization algorithms. Finally, we present results, which show synergistic effects when used with Solid Isotropic Material with Penalty and Rational Approximation of Material Properties material interpolation functions, popular methods of ensuring integerness of solutions.

42 ENGINEERING↗
DOE OSTI2016

Field-gradient partitioning for fracture and frictional contact in the material point method: Field-gradient partitioning for fracture and frictional contact in the material point method [Fracture and frictional contact in material point method using damage-field gradients for velocity-field partitioning]

Contact and fracture in the material point method require grid-scale enrichment or partitioning of material into distinct velocity fields to allow for displacement or velocity discontinuities at a material interface. We present a new method which a kernel-based damage field is constructed from the particle data. The gradient of this field is used to dynamically repartition the material into contact pairs at each node. Our approach avoids the need to construct and evolve explicit cracks or contact surfaces and is therefore well suited to problems involving complex 3-D fracture with crack branching and coalescence. A straightforward extension of this approach permits frictional ‘self-contact’ between surfaces that are initially part of a single velocity field, enabling more accurate simulation of granular flow, porous compaction, fragmentation, and comminution of brittle materials. Finally, numerical simulations of self contact and dynamic crack propagation are presented to demonstrate the accuracy of the approach.

97 MATHEMATICS AND COMPUTING↗