Search NASA⌕ Search

SEARCH · Search NASA

Results for “level sets”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗
DOE OSTI2022

CAMERA: A method for cost-aware, adaptive, multifidelity, efficient reliability analysis

Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distributions, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. Here, we propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential acquisition function based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest at the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set and probability estimation on synthetic test functions and two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier-Stokes equations. We show that our method predicts the failure boundary and probability more accurately and at a fraction of the computational cost compared with using just a single expensive high-fidelity model. Finally, we show that our sequential approach is guaranteed to asymptotically converge to the true failure boundary with high probability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On learning what to learn: Heterogeneous observations of dynamics and establishing possibly causal relations among them

Abstract Before we attempt to (approximately) learn a function between two sets of observables of a physical process, we must first decide what the inputs and outputs of the desired function are going to be. Here we demonstrate two distinct, data-driven ways of first deciding “the right quantities” to relate through such a function, and then proceeding to learn it. This is accomplished by first processing simultaneous heterogeneous data streams (ensembles of time series) from observations of a physical system: records of multiple observation processes of the system. We determine (i) what subsets of observables are common between the observation processes (and therefore observable from each other, relatable through a function); and (ii) what information is unrelated to these common observables, therefore particular to each observation process, and not contributing to the desired function. Any data-driven technique can subsequently be used to learn the input–output relation—from k-nearest neighbors and Geometric Harmonics to Gaussian Processes and Neural Networks. Two particular “twists” of the approach are discussed. The first has to do with the identifiability of particular quantities of interest from the measurements. We now construct mappings from a single set of observations from one process to entire level sets of measurements of the second process, consistent with this single set. The second attempts to relate our framework to a form of causality: if one of the observation processes measures “now,” while the second observation process measures “in the future,” the function to be learned among what is common across observation processes constitutes a dynamical model for the system evolution.

Sroczynski, David W.↗
NASA NTRS1991

Unstructured grid generation using the distance function

A new class of methods for obtaining level sets to generate unstructured grids is presented. The consecutive grid levels are computed using the distance functions, which corresponds to solving the Hamilton-Jacobi equations representing the equations of motion of fronts propagating with curvature-dependent speed. The relationship between the distance function and the governing equations will be discussed as well as its application to generating grids. Multi-ply connected domains and complex geometries are handled naturally, with a straightforward generalization to several space dimensions. The grid points for the unstructured grid are obtained simultaneously with the grid levels. The search involved in checking for overlapping triangles is minimized by triangulating the entire domain one level at a time.

Bihari, Barna L.↗
DOE OSTI2016

Topological classification of the Goryachev integrable case in rigid body dynamics

A topological analysis of the Goryachev integrable case in rigid body dynamics is made on the basis of the Fomenko-Zieschang theory. The invariants (marked molecules) which are obtained give a complete description, from the standpoint of Liouville classification, of the systems of Goryachev type on various level sets of the energy. It turns out that on appropriate energy levels the Goryachev case is Liouville equivalent to many classical integrable systems and, in particular, the Joukowski, Clebsch, Sokolov and Kovalevskaya-Yehia cases in rigid body dynamics, as well as to some integrable billiards in plane domains bounded by confocal quadrics -- in other words, the foliations given by the closures of generic solutions of these systems have the same structure. Bibliography: 15 titles.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Exospheric temperature and composition from satellite beacon measurements

Routine measurements of the slab thickness of the ionosphere, from 1965 to 1971, are used to infer the changes in neutral temperature and ion composition at a mean latitude of 40 S. Values of neutral temperature at solar maximum are 5 to 10% above Northern Hemisphere backscatter results. The diurnal and seasonal changes agree closely with satellite drag and backscatter measurements, except that the maximum temperature occurs after sunset in winter. Winter night-time values of the O(+)/H(+) transition height were 500 km in 1965-1966, 800 km in 1968-1969, and 700 km in 1971. Changes in the transition height lag about six months behind the changes in solar flux. Diurnal variations have a minimum just before sunrise and a maximum 1 to 3 hr after noon. On winter nights the transition height descends to the level set by chemical equilibrium. On summer nights the transition height is always above this level, giving a continual production of H(+) which serves as an additional source for maintaining the night-time ionosphere in the winter hemisphere.

Titheridge, J. E.↗
NASA NTRS1994

Oscillator strengths of selected resonance transitions in neutral sulfur

Mean lives and branching ratios for the 4s (3)S(sup 0 sub 1) and 4s(double prime) (3)P(sup 0 sub 1, 2) levels of neutral sulfur were determined at the Toledo Heavy Ion Accelerator using beam-foil spectroscopic techniques. The mean lives obtained for the 4s (3)S(sup 0 sub 1) 4s(double prime) (3)P(sup 0 sub 1), and (3)P(sup 0 sub 2) levels, 1.875 +/- 0.094 ns, 2.034 +/- 0.102 ns, and 2.146 +/- 0.129 ns, respectively, represent the most accurate results available to date. Oscillator strengths for the transitions 3p(sup 4)(3)P(sub 2, 1, 0)-4s (3)S(sup 0 sub 1) and 3p(sup 4)(3)P(sub 2, 1, 0)-4s(double prime) (3)P(sup 0 sub 1, 2) were derived from these mean lives together with the measured branching ratios. For comparison with published results which present only multiplet f-values, such f-values were computed from the measured individual line oscillator strengths. The value obtained for the 3p(sup 4)(3)P-4s(sup 3)S(sup 0) multiplet at 1814 A is 0.088 +/- 0.005, in good agreement with the mean of previous experimental measurements and theoretical calculations. With the mean life for the 4s(double prime) 3P(sup 0 sub 0) level set equal to the mean of our determination for the other two 4s(double prime) (3)P(sup 0) levels (a valid assumption for LS coupling conditions, and consistent with the observed equality of the mean lives within the accuracy of our measurement), an oscillator strength for the 1299 A multiplet of 0.121 +/- 0.004 was found. Both results are in good agreement with recent large-scale theoretical calculations that incorporate the effects of configuration interaction.

Beideck, D. J.↗
DOE OSTI2014

Assessment of dedicated low-dose cardiac micro-CT reconstruction algorithms using the left ventricular volume of small rodents as a performance measure

Purpose: Phase-correlated microcomputed tomography (micro-CT) imaging plays an important role in the assessment of mouse models of cardiovascular diseases and the determination of functional parameters as the left ventricular volume. As the current gold standard, the phase-correlated Feldkamp reconstruction (PCF), shows poor performance in case of low dose scans, more sophisticated reconstruction algorithms have been proposed to enable low-dose imaging. In this study, the authors focus on the McKinnon-Bates (MKB) algorithm, the low dose phase-correlated (LDPC) reconstruction, and the high-dimensional total variation minimization reconstruction (HDTV) and investigate their potential to accurately determine the left ventricular volume at different dose levels from 50 to 500 mGy. The results were verified in phantom studies of a five-dimensional (5D) mathematical mouse phantom. Methods: Micro-CT data of eight mice, each administered with an x-ray dose of 500 mGy, were acquired, retrospectively gated for cardiac and respiratory motion and reconstructed using PCF, MKB, LDPC, and HDTV. Dose levels down to 50 mGy were simulated by using only a fraction of the projections. Contrast-to-noise ratio (CNR) was evaluated as a measure of image quality. Left ventricular volume was determined using different segmentation algorithms (Otsu, level sets, region growing). Forward projections of the 5D mouse phantom were performed to simulate a micro-CT scan. The simulated data were processed the same way as the real mouse data sets. Results: Compared to the conventional PCF reconstruction, the MKB, LDPC, and HDTV algorithm yield images of increased quality in terms of CNR. While the MKB reconstruction only provides small improvements, a significant increase of the CNR is observed in LDPC and HDTV reconstructions. The phantom studies demonstrate that left ventricular volumes can be determined accurately at 500 mGy. For lower dose levels which were simulated for real mouse data sets, the HDTV algorithm shows the best performance. At 50 mGy, the deviation from the reference obtained at 500 mGy were less than 4%. Also the LDPC algorithm provides reasonable results with deviation less than 10% at 50 mGy while PCF and MKB reconstruction show larger deviations even at higher dose levels. Conclusions: LDPC and HDTV increase CNR and allow for quantitative evaluations even at dose levels as low as 50 mGy. The left ventricular volumes exemplarily illustrate that cardiac parameters can be accurately estimated at lowest dose levels if sophisticated algorithms are used. This allows to reduce dose by a factor of 10 compared to today's gold standard and opens new options for longitudinal studies of the heart.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Accumulated SET protein up-regulates and interacts with hnRNPK, increasing its binding to nucleic acids, the Bcl-xS repression, and cellular proliferation

Highlights: • hnRNPK is a new target of SET. • SET regulates hnRNPK. • SET and hnRNPK accumulation promotes tumorigenesis. • SET accumulation is a potential model to study genes regulated by SET-hnRNPK. - Abstract: SET and hnRNPK are proteins involved in gene expression and regulation of cellular signaling. We previously demonstrated that SET accumulates in head and neck squamous cell carcinoma (HNSCC); hnRNPK is a prognostic marker in cancer. Here, we postulate that SET and hnRNPK proteins interact to promote tumorigenesis. We performed studies in HEK293 and HNSCC (HN6, HN12, and HN13) cell lines with SET/hnRNPK overexpression and knockdown, respectively. We found that SET and/or hnRNPK protein accumulation increased cellular proliferation. SET accumulation up-regulated hnRNPK mRNA and total/phosphorylated protein, promoted hnRNPK nuclear location, and reduced Bcl-x mRNA levels. SET protein directly interacted with hnRNPK, increasing both its binding to nucleic acids and Bcl-xS repression. We propose that hnRNPK should be a new target of SET and that SET–hnRNPK interaction, in turn, has potential implications in cell survival and malignant transformation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Two-Level Convergence Theory for Multigrid Reduction in Time (MGRIT)

Here we develop a two-grid convergence theory for the parallel-in-time scheme known as multigrid reduction in time (MGRIT), as it is implemented in the open-source package [XBraid: Parallel Multigrid in Time, http://llnl.gov/casc/xbraid]. MGRIT is a scalable and multilevel approach to parallel-in-time simulations that nonintrusively uses existing time-stepping schemes, and in a specific two-level setting it is equivalent to the widely known parareal algorithm. The goal of this paper is twofold. First, we present a two-level MGRIT convergence analysis for linear problems where the spatial discretization matrix can be diagonalized, and then apply this analysis to our two basic model problems, the heat equation and the advection equation. One important assumption is that the coarse and fine time-grid propagators can be diagaonalized by the same set of eigenvectors, which is often the case when the same spatial discretization operator is used on the coarse and fine time grids. In many cases, the MGRIT algorithm is guaranteed to converge, and we demonstrate numerically that the theoretically predicted convergence rates are sharp in practice for our model problems. Second, we explore how the convergence of MGRIT compares to the stability of the chosen time-stepping scheme. In particular, we demonstrate that a stable time-stepping scheme does not necessarily imply convergence of MGRIT, although MGRIT with FCF-relaxation always converges for the diffusion dominated problems considered here.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1.

Chemistry↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2006

Space and Ground Trades for Human Exploration and Wearable Computing

Human exploration of the Moon and Mars will present unique trade study challenges as ground system elements shift to planetary bodies and perhaps eventually to the bodies of human explorers in the form of wearable computing technologies. This presentation will highlight some of the key space and ground trade issues that will face the Exploration Initiative as NASA begins designing systems for the sustained human exploration of the Moon and Mars, with an emphasis on wearable computing. We will present some preliminary test results and scenarios that demonstrate how wearable computing might affect the trade space noted below. We will first present some background on wearable computing and its utility to NASA's Exploration Initiative. Next, we will discuss three broad architectural themes, some key ground and space trade issues within those themes and how they relate to wearable computing. Lastly, we will present some preliminary test results and suggest guidance for proceeding in the assessment and creation of a value-added role for wearable computing in the Exploration Initiative. The three broad ground-space architectural trade themes we will discuss are: 1. Functional Shift and Distribution: To what extent, if any, should traditional ground system functionality be shifted to, and distributed among, the Earth, Moon/Mars, and the human. explorer? 2. Situational Awareness and Autonomy: How much situational awareness (e.g. environmental conditions, biometrics, etc.) and autonomy is required and desired, and where should these capabilities reside? 3. Functional Redundancy: What functions (e.g. command, control, analysis) should exist simultaneously on Earth, the Moon/Mars, and the human explorer? These three themes can serve as the axes of a three-dimensional trade space, within which architectural solutions reside. We will show how wearable computers can fit into this trade space and what the possible implications could be for the rest of the ground and space architecture(s). We intend this to be an example of explorer-centric thinking in a fully integrated explorer paradigm, where integrated explorer refers to a human explorer having instant access to all relevant data, knowledge of the environment, science models, health and safety-related events, and other tools and information via wearable computing technologies. The trade study approach will include involvement from the relevant stakeholders (Constellation Systems, CCCI, EVA Project Office, Astronaut office, Mission Operations, Space Life Sciences, etc.) to develop operations concepts (and/or operations scenarios) from which a basic high-level set of requirements could be extracted. This set of requirements could serve as a foundation (along with stakeholder buy-in) that would help define the trade space and assist in identifying candidate technologies for further study and evolution to higher-level technology readiness levels.

Lupisella, Mark↗
DOE OSTI2015

WE-G-BRD-04: BEST IN PHYSICS (JOINT IMAGING-THERAPY): An Integrated Model-Based Intrafractional Organ Motion Tracking Approach with Dynamic MRI in Head and Neck Radiotherapy

Purpose: In-treatment dynamic cine images, provided by the first commercially available MRI-guided radiotherapy system, allow physicians to observe intrafractional motion of head and neck (H&N) internal structures. Nevertheless, high anatomical complexity and relatively poor cine image contrast/resolution have complicated automatic intrafractional motion evaluation. We proposed an integrated model-based approach to automatically delineate and analyze moving structures from on-board cine images. Methods: The H&N upper airway, a complex and highly deformable region wherein severe internal motion often occurs, was selected as the target-to-be-tracked. To reliably capture its motion, a hierarchical structure model containing three statistical shapes (face, face-jaw, and face-jaw-palate) was first built from a set of manually delineated shapes using principal component analysis. An integrated model-fitting algorithm was then employed to align the statistical shapes to the first to-be-detected cine frame, and multi-feature level-set contour propagation was performed to identify the airway shape change in the remaining frames. Ninety sagittal cine MR image sets, acquired from three H&N cancer patients, were utilized to demonstrate this approach. Results: The tracking accuracy was validated by comparing the results to the average of two manual delineations in 20 randomly selected images from each patient. The resulting dice similarity coefficient (93.28+/−1.46 %) and margin error (0.49+/−0.12 mm) showed good agreement with the manual results. Intrafractional displacements of anterior, posterior, inferior, and superior airway boundaries were observed, with values of 2.62+/−2.92, 1.78+/−1.43, 3.51+/−3.99, and 0.68+/−0.89 mm, respectively. The H&N airway motion was found to vary across directions, fractions, and patients, and highly correlated with patients’ respiratory frequency. Conclusion: We proposed the integrated computational approach, which for the first time allows to automatically identify the H&N upper airway and quantify in-treatment H&N internal motion in real-time. This approach can be applied to track other structures’ motion, and provide guidance on patient-specific prediction of intra-/inter-fractional structure displacements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Crickets{sup TM} and Curies: Lessons Learned from the Steel Industry - 18221

The global nuclear power industry, and to a lesser extent, the decommissioning and demolition of old nuclear weapons plants, obviously present challenges to those in the radioactive waste management field. Some if not many of these challenges have already been dealt with, and resolved, in the steel industry, including its recycling component. This paper presents proven state-of-the-art radiation detection systems used in steel and related recycling facilities, emphasizing current, reliable and flexible technologies for characterization, volume reduction, and remote handling of recyclable solid wastes contaminated and, just as important, not contaminated, with a wide range of radioactive isotopes. It will also cover features of remotely operated waste handling equipment effective in reducing or eliminating exposure of workers engaged in decontamination and demolition work on active and/or inactive nuclear installations. The industry group identified above has a long-standing policy regarding any radioactivity in potentially recyclable as well as waste materials it receives: if it is detected, it is rejected. The US Steel Manufacturers Association (SMA) put it this way in one of their Public Policy Statements: 'SMA members have not, and will not, accept scrap that is known to be radioactively contaminated' (emphasis included in the original) (1). This policy, adopted around the globe, obviously has major implications for decommissioning and demolition projects. To put it bluntly, if rad-contaminated metal, concrete, or other materials are sent to a facility for possible recycling, reuse or even disposal, they may well be rejected. If so, the materials will be returned to the generator--at the generator's cost-and a new solution will have to be found. This also results in additional costs for the generator. Facilities in the steel making and metal recycling industries have access to, and routinely count on, a wide range of proven detection technologies and techniques for monitoring incoming shipments and sorting the materials in them. The approaches range from using sensitive handheld instruments for pinpointing and isolating small devices such as isotope-containing gauges and radium dials all the way up to technologies effective for large volumes and weights of contaminated construction and demolition material. In this latter realm, the most suitable, robust, and field-proven system for decommissioning and demolition work is a radiation detection system fitted to a grappling and/or magnetic device, known as the Crickets{sup TM} technology, which automatically scans bulk material with no added effort or time in material handling operations. Highly versatile, a Crickets{sup TM} system can be put to work on any material-handling system, including grapples handling volumes ranging from 1 up to 18 cubic yards, while detecting radioactivity in real-time over a very wide dose range, from nSv to Sv levels. Other field-proven real-time measuring radiation detection devices can be fixed to conveyor belts with traverse speeds up 3 meters/second, also measuring in real time and over the same dose range, allowing both recyclable and waste materials with a pre-set level of radioactivity to be diverted off the belt. Finally, there are highly sensitive radiation portal installations for detecting radiation from materials being transported by vehicles, including those with roll-offs, dump trailers, and other containers, again in real-time and for the same dose ranges. For vehicles, the speed must be carefully controlled-the vehicles need to travel from 5 to 10 km/h through the scanning zone to ensure that maximum detection capability is maintained. When radiation detection systems at recycling and steel making facilities detect radioactive material in an incoming shipment, these facilities have a zero tolerance for accepting the cargo. Therefore, it is imperative that material from a demolition operation that is slated to be sent for recycling or other reuse be checked at the point of demolition at similar detection levels to avoid possible rejection later of a shipment because of its radioactivity. In addition to measuring overall radiation levels (for example, total gamma), when necessary, combinations of different radiation detection media can be used to detect the presence and intensities of specific isotopes such as U-235, Cs-137, Co-60, Am-241, etc., again in real-time and over a similar dose range. These detection technologies can be fixed or transportable. More important for the rad-waste management profession, this overall systems approach can be tailored and adapted during use on-site in any dismantling, demolition or remediation project, to allow materials to be thoroughly scanned before being released to the recycling stream for reuse, or to appropriate long-term storage and disposal facilities. Potentially recyclable materials with contamination such as steel beams and rebar will set off alarms at recycling facilities, resulting in unnecessary added expenses to the generator, primarily due to increased transportation costs, time delays and additional rad-waste handling and disposal charges. Additionally, there are situations where relatively lightly contaminated materials, once identified and appropriately characterized, can be reused in new construction applications, such as concrete from demolition being used in the base layer for roads or recycled steel used for repair or construction of bridges. In other cases, specific smaller volumes can be sent to specialized rad-waste disposal sites. These radiation detection systems can also be custom tailored to ensure proper handling and disposition of contaminated materials destined for reuse in specific applications where, for example, acceptable levels of radioactivity in NORM or TENORM can be reused or safely disposed of. (authors)

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

A stabilized Runge–Kutta–Legendre method for explicit super-time-stepping of parabolic and mixed equations

Parabolic partial differential equations appear in several physical problems, including problems that have a dominant hyperbolic part coupled to a sub-dominant parabolic component. Explicit methods for their solution are easy to implement but have very restrictive time step constraints. Implicit solution methods can be unconditionally stable but have the disadvantage of being computationally costly or difficult to implement. Super-time-stepping methods for treating parabolic terms in mixed type partial differential equations occupy an intermediate position. In such methods each superstep takes “s” explicit Runge–Kutta-like time-steps to advance the parabolic terms by a time-step that is s{sup 2} times larger than a single explicit time-step. The expanded stability is usually obtained by mapping the short recursion relation of the explicit Runge–Kutta scheme to the recursion relation of some well-known, stable polynomial. Prior work has built temporally first- and second-order accurate super-time-stepping methods around the recursion relation associated with Chebyshev polynomials. Since their stability is based on the boundedness of the Chebyshev polynomials, these methods have been called RKC1 and RKC2. In this work we build temporally first- and second-order accurate super-time-stepping methods around the recursion relation associated with Legendre polynomials. We call these methods RKL1 and RKL2. The RKL1 method is first-order accurate in time; the RKL2 method is second-order accurate in time. We verify that the newly-designed RKL1 and RKL2 schemes have a very desirable monotonicity preserving property for one-dimensional problems – a solution that is monotone at the beginning of a time step retains that property at the end of that time step. It is shown that RKL1 and RKL2 methods are stable for all values of the diffusion coefficient up to the maximum value. We call this a convex monotonicity preserving property and show by examples that it is very useful in parabolic problems with variable diffusion coefficients. This includes variable coefficient parabolic equations that might give rise to skew symmetric terms. The RKC1 and RKC2 schemes do not share this convex monotonicity preserving property. One-dimensional and two-dimensional von Neumann stability analyses of RKC1, RKC2, RKL1 and RKL2 are also presented, showing that the latter two have some advantages. The paper includes several details to facilitate implementation. A detailed accuracy analysis is presented to show that the methods reach their design accuracies. A stringent set of test problems is also presented. To demonstrate the robustness and versatility of our methods, we show their successful operation on problems involving linear and non-linear heat conduction and viscosity, resistive magnetohydrodynamics, ambipolar diffusion dominated magnetohydrodynamics, level set methods and flux limited radiation diffusion. In a prior paper (Meyer, Balsara and Aslam 2012 [36]) we have also presented an extensive test-suite showing that the RKL2 method works robustly in the presence of shocks in an anisotropically conducting, magnetized plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗