Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PRESERVATION”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Consequences of fine-tuning for fifth force searches

Light bosonic fields mediate long range forces between objects. If these fields have self-interactions, i.e., non-quadratic terms in the potential, the experimental constraints on such forces can be drastically altered due to a screening (chameleon) or enhancement effect. We explore how technically natural values for such self-interaction coupling constants modify the existing constraints. We point out that assuming the existence of these natural interactions leads to new constraints, contrary to the usual expectation that screening leads to gaps in coverage. We discuss how screening can turn fundamentally equivalence principle (EP)-preserving forces into EP-violating ones. This means that when natural screening is present, searches for EP violation can be used to constrain EP-preserving forces. We show how this effect enables the recently discovered stellar triple system PSR J0337+1715 to place a powerful constraint on EP-preserving fifth forces. Finally, we demonstrate that technically natural cubic self-interactions modify the vacuum structure of the scalar potential, leading to new constraints from spontaneous and induced vacuum decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Stressful crystal histories recorded around melt inclusions in volcanic quartz

Abstract Magma ascent and eruption are driven by a set of internally and externally generated stresses that act upon the magma. We present microstructural maps around melt inclusions in quartz crystals from six large rhyolitic eruptions using synchrotron Laue X-ray microdiffraction to quantify elastic residual strain and stress. We measure plastic strain using average diffraction peak width and lattice misorientation, highlighting dislocations and subgrain boundaries. Quartz crystals across studied magma systems preserve similar and relatively small magnitudes of elastic residual stress (mean 53–135 MPa, median 46–116 MPa) in comparison to the strength of quartz (~ 10 GPa). However, the distribution of strain in the lattice around inclusions varies between samples. We hypothesize that dislocation and twin systems may be established during compaction of crystal-rich magma, which affects the magnitude and distribution of preserved elastic strains. Given the lack of stress-free haloes around faceted inclusions, we conclude that most residual strain and stress was imparted after inclusion faceting. Fragmentation may be one of the final strain events that superimposes stresses of ~ 100 MPa across all studied crystals. Overall, volcanic quartz crystals preserve complex, overprinted deformation textures indicating that quartz crystals have prolonged deformation histories throughout storage, fragmentation, and eruption.

58 GEOSCIENCES↗
DOE OSTI2019

Effects of cardiac resynchronization therapy on right ventricular function during rest and exercise, as assessed by radionuclide angiography, and on NT-proBNP levels

Aim: We carried out this study to investigate mid-term effects of cardiac resynchronization therapy (CRT) on right ventricular (RV) function and neurohormonal response, expressed by N-terminal pro-brain natriuretic peptide (NT-proBNP), in heart failure patients stratified by baseline RV ejection fraction (RVEF). Methods and Results: Thirty-six patients with nonischemic dilated cardiomyopathy underwent technetium-99m radionuclide angiography with bicycle exercise immediately after CRT implantation (during spontaneous rhythm and after CRT activation) and 3 months later. Plasma NT proBNP was assessed before implantation and after 3 months. At baseline, RVEF was impaired (≤35%) in 14 patients, preserved (>35%) in 22. At 3 months, RVEF improved during rest and exercise (P = .02) in patients with impaired RV function, while remaining unchanged in patients with preserved RV function. Rest and exercise RV dyssynchrony decreased in both groups at follow-up (P < .05). A similar mid-term improvement in left ventricular (LV) function and NT-proBNP was observed in patients with impaired and preserved RVEF. In the former, the decrease in NT-proBNP correlated with the improvements both in LV and RV dyssynchrony and functions. Conclusion: CRT may improve RV performance, during rest and exercise, and neurohormonal response in heart failure patients with nonischemic dilated cardiomyopathy and baseline RV dysfunction. RV dysfunction should not be considered per se a primary criterion for excluding candidacy to CRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Harada–Tsutsui gauge recovery procedure: From Abelian gauge anomalies to the Stueckelberg mechanism

Revisiting a path-integral procedure developed by Harada and Tsutsui for recovering gauge invariance from anomalous effective actions, it is shown that there are two ways to achieve gauge symmetry: one already presented by the authors, which is shown to preserve the anomaly in the sense of standard current conservation law, and another one which is anomaly-free, preserving current conservation. It is also shown that the application of the Harada–Tsutsui technique to other models which are not anomalous but do not exhibit gauge invariance allows the identification of the gauge invariant formulation of the Proca model, also done by the referred authors, with the Stueckelberg model, leading to the interpretation of the gauge invariant map as a generalization of the Stueckelberg mechanism. -- Highlights: • A gauge restoration technique from Abelian anomalous models is discussed. • It is shown that there is another way that leads to gauge symmetry restoration from such technique. • It is shown that the first gauge restoration preserves the anomaly, while the proposed second one is free from anomalies. • It is shown that the proposed gauge symmetry restoration can be identified with the Stueckelberg mechanism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Photon pumping in a weakly-driven quantum cavity–spin system

Highlights: • We study photon frequency conversion in a driven spin coupled to a cavity mode. • Quantized frequency conversion is excpected in the strong-drive adiabatic limit. • A new photon pumping effect is established in the accessible weak drive regime. • Pumping is linked to the delocalization of the corresponding Floquet states. • Quantum coherence is preserved in both the strong and ultraweak-drive limits. We investigate the photon pumping effect in a topological model consisting of a periodically driven spin-1/2 coupled to a quantum cavity mode out of the adiabatic limit. In the strong-drive adiabatic limit, a quantized frequency conversion of photons is expected as the temporal analog of the Hall current. We numerically establish a novel photon pumping phenomenon in the experimentally accessible nonadiabatic driving regime for a broad region of the parameter space. The photon frequency conversion efficiency exhibits strong fluctuations and high efficiency that can reach up 80% of the quantized value for commensurate frequency combinations. We link the pumping properties to the delocalization of the corresponding Floquet states which display multifractal behavior as the result of hybridization between localized and delocalized sectors. Finally we demonstrate that the quantum coherence properties of the initial state are preserved during the frequency conversion process in both the strong and ultra-weak-drive limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Regionally Correlated Oxygen and Carbon Isotope Zonation in Diagenetic Carbonates of the Bakken Formation

Diagenetic minerals preserve records of burial processes that overprint records of seawater chemistry and impact reservoir porosity. The Mississippian-Devonian aged Bakken Formation in the Williston Basin is a reservoir rock of economic importance whose productivity is affected by diagenetic carbonates, particularly dolomite-ankerite-series carbonates. To investigate how diagenetic carbonate alteration manifests in the Bakken and how that might change across the basin, a combined δ18O and δ13C isotope dataset for diagenetic carbonates was collected by in situ SIMS analysis of 10-μm spots from nine drill holes covering a ~250 km transect of the middle Bakken member. Observed core-to-rim isotopic variability in these small Fe-zoned dolomites and calcites frequently exceeds 10‰ in both δ 18 O and δ 13 C, indicating significant changes in thermal and chemical conditions during cementation. Individual ankeritic growth bands can be correlated across the basin by systematic similarities in minor-element compositions and isotope ratios. In the central part of the basin, δ 18 O and δ 13 C trends at sub-mm-scale decline consistently towards the rims of ankerite-series carbonates, which is interpreted to reflect mineral growth coincident with rising temperatures and an increasing organic contribution to inorganic carbon during burial. The most abrupt shifts in δ 18 O and δ 13 C (changes as large as 12‰ for δ 18 O and 6.5‰ for δ13C within distances < 15 μm) are observed along the basin margin, and are believed to preserve signals of hydrocarbon expulsion and fluid infiltration along basin margins. Based on the data, we conclude that ankerite-series carbonates preserve records of prolonged thermal and chemical processes that operated basinwide. The isotopic trends presented here show diagenetic rock evolution from shallow dolomite formation to deep burial ankerites, as well as the isotopic responses of carbonates to changing fluids and availability of organic δ 13 C. These results provide important constraints on how carbonate isotope records may be altered during burial in organic-rich sedimentary rocks and emphasize the need for caution in using bulk powder samples for geochemical analysis in carbonate systems.

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

Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations

In this work, we introduce a new residual distribution (RD) framework for the design of bound-preserving high-resolution finite element schemes. The continuous and discontinuous Galerkin discretizations of the linear advection equation are modified to construct local extremum diminishing (LED) approximations. To that end, we perform mass lumping and redistribute the element residuals in a manner which guarantees the LED property. The hierarchical correction procedure for high-order Bernstein finite element discretizations involves localization to subcells and definition of bound-preserving weights for subcell contributions. Using strong stability preserving (SSP) Runge–Kutta methods for time integration, we prove the validity of discrete maximum principles under CFL-like time step restrictions. The low-order version of our method has roughly the same accuracy as the one derived from a piecewise (multi)-linear approximation on a submesh with the same nodal points. In high-order extensions, we use an element-based flux-corrected transport (FCT) algorithm which can be interpreted as a nonlinear RD scheme. The proposed LED corrections are tailor-made for matrix-free implementations which avoid the rapidly growing cost of matrix assembly for high-order Bernstein elements. The results for 1D, 2D, and 3D test problems compare favorably to those obtained with the best matrix-based approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Early diagenetic processes in an iron-dominated marine depositional system

The early diagenetic interplay between reactive iron, sulfur, and organic matter in the bathymetrically isolated Santa Monica Basin (SMB) sediments are investigated in this study. Here, we explore solid-phase and porewater profiles from the basin, supplemented with a transect from 71 to 907 m water depth that includes oxygenated (>60 μM O 2 ) bottom waters near the coast and oxygen-deficient waters (~4 µM O 2 ) in the basin. The geochemical data of the basin sediments are further scrutinized by means of reactive transport modeling. The results show that the basin sediments do not follow the traditional geochemical signatures of oxygen-deficient settings. A lack of dissolved sulfide accumulation and sulfurized iron persists despite the sediments being deposited under reducing conditions (without bioturbation/bioirrigation), strong organic carbon input (TOC up to 5.0 wt%), and active dissimilatory sulfate reduction. Not only did we find an exceptional enrichment in highly reactive Fe in the surface sediments (~45 % of total Fe), but the enrichment of reactive Fe, including ferrihydrite, persists downcore and coexists with high levels of dissolved Fe. The enhanced preservation of Fe oxides and lack of iron-sulfide precipitation is in part explained by detection via Mössbauer spectra of iron oxides bounded to organic matter (Fe[III]-OM coprecipitates). The modeled Fe budget shows that most of the Fe oxides in the surface sediments are internally recycled by upward diffusion and subsequent oxidation of Fe 2+ . Sulfide oxidation coupled to Fe reduction effectively precludes sulfide accumulation while enhancing build-up of dissolved Fe, fueling the Fe cycle within the first 5 cm depth. Continuous reoxidation of Fe 2+ enhances the formation of Fe(III)-OM coprecipitates, limiting the amount of reactive organic matter. In the unavailability of labile organic matter, other than within the uppermost layers, the organic-rich sediment profiles are dominated by Fe cycling that limits the production and preservation of sulfides and enhances the preservation of Fe oxides and organic carbon. Finally, this study highlights key local controls on Fe availability in marginal basins and describes an intricate biogeochemical C-Fe-S cycling in modern and possibly ancient marine systems with important implications for Fe availability in the marine realm.

58 GEOSCIENCES↗
DOE OSTI2015

Long-term Results of the UCSF-LBNL Randomized Trial: Charged Particle With Helium Ion Versus Iodine-125 Plaque Therapy for Choroidal and Ciliary Body Melanoma

Purpose: Relevant clinical data are needed given the increasing national interest in charged particle radiation therapy (CPT) programs. Here we report long-term outcomes from the only randomized, stratified trial comparing CPT with iodine-125 plaque therapy for choroidal and ciliary body melanoma. Methods and Materials: From 1985 to 1991, 184 patients met eligibility criteria and were randomized to receive particle (86 patients) or plaque therapy (98 patients). Patients were stratified by tumor diameter, thickness, distance to disc/fovea, anterior extension, and visual acuity. Tumors close to the optic disc were included. Local tumor control, as well as eye preservation, metastases due to melanoma, and survival were evaluated. Results: Median follow-up times for particle and plaque arm patients were 14.6 years and 12.3 years, respectively (P=.22), and for those alive at last follow-up, 18.5 and 16.5 years, respectively (P=.81). Local control (LC) for particle versus plaque treatment was 100% versus 84% at 5 years, and 98% versus 79% at 12 years, respectively (log rank: P=.0006). If patients with tumors close to the disc (<2 mm) were excluded, CPT still resulted in significantly improved LC: 100% versus 90% at 5 years and 98% versus 86% at 12 years, respectively (log rank: P=.048). Enucleation rate was lower after CPT: 11% versus 22% at 5 years and 17% versus 37% at 12 years, respectively (log rank: P=.01). Using Cox regression model, likelihood ratio test, treatment was the most important predictor of LC (P=.0002) and eye preservation (P=.01). CPT was a significant predictor of prolonged disease-free survival (log rank: P=.001). Conclusions: Particle therapy resulted in significantly improved local control, eye preservation, and disease-free survival as confirmed by long-term outcomes from the only randomized study available to date comparing radiation modalities in choroidal and ciliary body melanoma.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

High-Dose-Rate Brachytherapy in the Management of Operable Rectal Cancer: A Systematic Review

Purpose: To evaluate the role of high-dose-rate endorectal brachytherapy (HDREBT) in the preoperative and definitive management of operable rectal cancer in terms of clinical outcomes and toxicities using a systematic review. Methods and Materials: A review of published articles from January 1990 to December 2016 was conducted using the PubMed, Embase, and Scopus databases using the search terms “rectal” or “rectum” in combination with “brachytherapy,” “high dose rate,” “HDR,” and “endorectal.” Additional publications were identified by scanning references. Only studies published in English reporting clinical outcomes with ≥30 patients treated with HDREBT were included. Results: The search identified 1688 articles, of which 22 met our inclusion criteria. Twelve studies were included in this systematic review. Following preoperative HDREBT with chemoradiation therapy (CRT), the pathologic complete response (pCR) rate ranged between 18% and 31% (weighted mean rate, 22.2%); R0 resection rate, between 80% and 99% (weighted mean rate, 95.5%); and sphincter-preservation rate, between 29% and 54% (weighted mean rate, 46.4%). The weighted mean 2-year progression-free survival and overall survival (OS) rates were 68.1% and 81.5%, respectively. After preoperative HDREBT alone, the pCR rate ranged between 10.4% and 27% (weighted mean rate, 23.8%), the R0 rate was 96.5% (1 study), and the sphincter-preservation rate ranged between 53.8% and 75.8% (weighted mean rate, 59.4%). The weighted mean 5-year progression-free survival and OS rates were 66.6% and 70.8%, respectively. There was only 1 study of HDREBT for nonsurgical management of rectal cancer, which reported a 2-year OS rate of 100%. Conclusions: Preoperative HDREBT either alone or in combination with CRT may result in a better pCR but may not necessarily translate into better survival, which is similar to outcomes seen following preoperative CRT alone. There were significant variations across studies in terms of patient selection, treatment approaches, and evaluation of clinical outcomes, suggesting the need for an international consensus on the dosimetric parameters and techniques of HDREBT, timing and methods of response assessment, definitions and assessment of toxicities, and optimal timing of surgery before further prospective studies. Future studies should include evaluation of the role of HDREBT in the nonsurgical curative treatment of screen-detected early cancers and organ preservation in lower rectal cancers.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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↗
DOE OSTI2014

A numerical method for the quasi-incompressible Cahn–Hilliard–Navier–Stokes equations for variable density flows with a discrete energy law

In this paper, we investigate numerically a diffuse interface model for the Navier–Stokes equation with fluid–fluid interface when the fluids have different densities [48]. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energy law preserving method will perform better for multiphase problems. Thus for the reformulated system, we design a C{sup 0} finite element method and a special temporal scheme where the energy law is preserved at the discrete level. Such a discrete energy law (almost the same as the continuous energy law) for this variable density two-phase flow model has never been established before with C{sup 0} finite element. A Newton method is introduced to linearise the highly non-linear system of our discretization scheme. Some numerical experiments are carried out using the adaptive mesh to investigate the scenario of coalescing and rising drops with differing density ratio. The snapshots for the evolution of the interface together with the adaptive mesh at different times are presented to show that the evolution, including the break-up/pinch-off of the drop, can be handled smoothly by our numerical scheme. The discrete energy functional for the system is examined to show that the energy law at the discrete level is preserved by our scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The construction of arbitrary order ERKN methods based on group theory for solving oscillatory Hamiltonian systems with applications

In general, extended Runge–Kutta–Nyström (ERKN) methods are more effective than traditional Runge–Kutta–Nyström (RKN) methods in dealing with oscillatory Hamiltonian systems. However, the theoretical analysis for ERKN methods, such as the order conditions, the symplectic conditions and the symmetric conditions, becomes much more complicated than that for RKN methods. Therefore, it is a bottleneck to construct high-order ERKN methods efficiently. In this paper, we first establish the ERKN group Ω for ERKN methods and the RKN group G for RKN methods, respectively. We then rigorously show that ERKN methods are a natural extension of RKN methods, that is, there exists an epimorphism η of the ERKN group Ω onto the RKN group G. This epimorphism gives a global insight into the structure of the ERKN group by the analysis of its kernel and the corresponding RKN group G. Meanwhile, we establish a particular mapping φ of G into Ω so that each image element is an ideal representative element of the congruence class in Ω. Furthermore, an elementary theoretical analysis shows that this map φ can preserve many structure-preserving properties, such as the order, the symmetry and the symplecticity. From the epimorphism η together with its section φ, we may gain knowledge about the structure of the ERKN group Ω via the RKN group G. In light of the theoretical analysis of this paper, we obtain high-order structure-preserving ERKN methods in an effective way for solving oscillatory Hamiltonian systems. Numerical experiments are carried out and the results are very promising, which strongly support our theoretical analysis presented in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates

In this paper, we present new Lagrangian discontinuous Galerkin (DG) hydrodynamic methods for compressible flows on unstructured meshes in axisymmetric coordinates. The physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with linear Taylor series polynomials. Two different approaches are used to discretize the evolution equations – the first one is the true volume approach and the second one is the area-weighted approach. For the true volume approach, the DG equations are derived using the true 3D volume that is consistent with the geometry conservation law (GCL). The Riemann velocity at the nodes on the surface of the element, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem using surfaces areas for axisymmetric coordinates. This true volume approach conserves mass, momentum, and total energy and satisfies the GCL. However, it can not preserve spherical symmetry on an equal-angle polar grid with 1D radial flows. For the area-weighted approach, the DG equations are based on the 2D Cartesian geometry that is rotated about the axis of symmetry using a single, element average radius. With this approach, the Riemann velocity at the nodes on the surface of the element, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem in 2D Cartesian geometry. This area-weighted approach, in the limit of an infinitesimal mesh size, conserves physical momentum, and physical total energy. The area-weighted approach preserves spherical symmetry on an equal-angle polar grid for 1D radial flows, but it does not satisfy the GCL. Finally, a suite of test problems are calculated to demonstrate stable mesh motion, the expected second order accuracy of these methods, and that the new area-weighted DG method preserves spherical symmetry on 1D radial flow problems with equal-angle polar meshes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Conservative remapping of material-dependent fields between possibly misaligned material regions

In this work, we propose an interpolation or remapping algorithm of material-dependent fields on polyhedral meshes where any source or target cell contains only one material. It is conservative and it preserves sharp material boundaries on the target mesh, even if the source and target regions delineating the same material are slightly misaligned. If those material regions are aligned, then the algorithm is also linearity-preserving and bounds-preserving. For a given material, it consists of a conservative field reconstruction on a target mesh part from a source mesh part associated with that material, followed by a repair step in case of misaligned boundaries. No assumption is made regarding the topology of the input meshes

36 MATERIALS SCIENCE↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Phase decomposition and bubble evolution in Xe implanted U 3 Si 2 at 450°C

We report on the microstructure investigations of a U 3 Si 2 specimen implanted by 84 MeV Xe ions at 450°C. In the region that corresponds to the highest irradiation dose, U 3 Si 2 was found to decompose into a Si-enriched nanocrystalline USi matrix phase and a U-enriched amorphous inclusion phase. Density functional theory (DFT) calculations were used to help understand the decomposition mechanism. Xe bubbles of different morphologies were observed in USi nano-grains, on USi grain boundaries, and phase boundaries. Original U 3 Si 2 micro-grains were preserved on both sides of the irradiation damage peak region within the sample, implying a radiation dose threshold for decomposition of approximately 150 dpa at 450°C. In the preserved U 3 Si 2 region beneath the irradiation damage peak, where critical amount of Xe ions were deposited, a monomodal size distribution of intragranular Xe bubbles with an average size of 2.68 nm formed. In both decomposed and preserved U 3 Si 2 , the size of Xe bubbles was found to be lower than 50 nm. Based on this observation, the fission gas behavior of U 3 Si 2 is controllable and free of run-away swelling despite the occurrence of decomposition at this irradiation temperature. (C) 2019 Elsevier B.V. All rights reserved.

36 MATERIALS SCIENCE↗