Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase object approximation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

ZPRD Database: ZPPR-15 Monte Carlo Results

The ZPPR-15 experiments [1] were mockups of a 330 MWe Integral Fast Reactor (IFR). The ZPPR-15 assembly consisted of a clean, two zone, approximately circular core surrounded by a thin depleted uranium (DU) blanket with sodium (Na) cooling and a thick stainless steel reflector. The ZPPR-15 program was conducted in four phases: A, B, C, and D. Each phase was marked by a particular composition of the reference assembly, with the last three being representative of the three stages of the IFR fuel cycle. This report documents all MCNP [2] runs of the ZPPR-15 loadings that were included in the ZPRD database on GitLab. In the present work, eigenvalues and associated standard deviations were computed for each ZPPR-15 loading with the use of the three data libraries (ENDF/B-VII.0 [3], ENDF/B-VII.1 [4] and ENDF/B-VIII.0 [5]). Relevant results from the obtained outputs are also discussed in this report with two main objectives: a) provide and evaluate updated calculated values with respect to previous reports, notably Refs. [6] and [7] that were based on the use of the ENDF/B-VII.0 library. b) address any change in the observed reactivity effects relevant to the analysis of the experimental data when a different data library is used (previously reported results were mostly based on the use of ENDF/B-VII.0 data only).

Aliberti, Gerardo↗
DOE OSTI2019

Differential Many-Body Effects for Initial and Core Ionic States: Impact on XPS Spectra

In this paper, the contribution of many-body effects to the X-ray photoelectron spectroscopy, XPS, of an NO molecule is studied using wavefunction theory where the specific consequences of different many-body terms are examined and contrasted. It is shown that there is a differential importance of the many-body effects for the different configurations involved in the XPS. These are the ground, initial state configuration and final, N(1s) and O(1s) core hole ionic configurations. The consequences of the many-body effects are examined for the binding energies, BEs, to the two final state multiplets, triplet, and singlet, for each of the core ions and for the relative intensities of the XPS transitions to these multiplets. The many-body effects examined are those described as static effects that arise for individual terms that are important. The objective is to understand the chemical and physical origins that determine the importance of the correlation effects for the XPS, rather than to obtain very accurate predictions of the BEs. An important theoretical construct that is tested and justified is the equivalent core approximation where the core-ionized atom is replaced by the next higher element in the periodic table. This construct allows us to establish a correlation for the relative importance of the many-body effects in terms of effective charges of the different atoms. This is a correlation that has not been considered before and that we expect may have general relevance. The potential of the effects that we have identified for the XPS of NO to be relevant for the XPS of more complex, condensed phase systems is considered.

XPS, Static Correlation, Equivalent core model↗
DOE OSTI2017

Magnetar Broadband X-Ray Spectra Correlated with Magnetic Fields: Suzaku Archive of SGRs and AXPs Combined with NuSTAR, Swift, and RXTE

The 1–70 keV persistent spectra of 15 magnetars, observed with Suzaku from 2006 to 2013, were studied as a complete sample. Combined with early NuSTAR observations of four hard X-ray emitters, nine objects showed a hard power-law emission dominating at ≳10 keV with the 15–60 keV flux of ∼1–11×10{sup −11} erg s{sup −1} cm{sup −2}. The hard X-ray luminosity L{sub h}, relative to that of a soft-thermal surface radiation L{sub s}, tends to become higher toward younger and strongly magnetized objects. Their hardness ratio, updated from a previous study and defined as ξ=L{sub h}/L{sub s}, is correlated with the measured spin-down rate P-dot as ξ=0.62×( P-dot /10{sup −11}ss{sup −1}){sup 0.72}, corresponding to positive and negative correlations with the dipole field strength B{sub d} (ξ∝B{sub d}{sup 1.41}) and the characteristic age τ{sub c} (ξ∝τ{sub c}{sup −0.68}), respectively. Among our sample, five transients were observed during X-ray outbursts, and the results are compared with their long-term 1–10 keV flux decays monitored with Swift/XRT and RXTE/PCA. Fading curves of three bright outbursts are approximated by an empirical formula used in the seismology, showing a ∼10–40 day plateau phase. Transients show the maximum luminosities of L{sub s} ∼ 10{sup 35} erg s{sup −1}, which are comparable to those of persistently bright ones, and fade back to ≲10{sup 32} erg s{sup −1}. Spectral properties are discussed in the framework of the magnetar hypothesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Kinetic theory for strongly coupled Coulomb systems

The calculation of dynamical properties for matter under extreme conditions is a challenging task. The popular Kubo-Greenwood model exploits elements from equilibrium density-functional theory (DFT) that allow a detailed treatment of electron correlations, but its origin is largely phenomenological; traditional kinetic theories have a more secure foundation but are limited to weak ion-electron interactions. The objective here is to show how a combination of the two evolves naturally from the short-time limit for the generator of the effective single-electron dynamics governing time correlation functions without such limitations. This provides a theoretical context for the current DFT-related approach, the Kubo-Greenwood model, while showing the nature of its corrections. The method is to calculate the short-time dynamics in the single-electron subspace for a given configuration of the ions. This differs from the usual kinetic theory approach in which an average over the ions is performed as well. In this way the effective ion-electron interaction includes strong Coulomb coupling and is shown to be determined from DFT. The correlation functions have the form of the random-phase approximation for an inhomogeneous system but with renormalized ion-electron and electron-electron potentials. The dynamic structure function, density response function, and electrical conductivity are calculated as examples. The static local field corrections in the dielectric function are identified in this way. The current analysis is limited to semiclassical electrons (quantum statistical potentials), so important quantum conditions are excluded. However, a quantization of the kinetic theory is identified for broader application while awaiting its detailed derivation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Pilot Scale Operation and Testing of Syngas Chemical Looping For Hydrogen Production

The Syngas Chemical looping (SCL) process provides efficient and economic means to utilize the abundant fossil reserve of coal. The main problem associated with coal utilization is the CO 2 emissions resulting from its combustion. Even though CO 2 regulation or carbon tax is currently not in place, its enforcement is expected in the near future. Such a greenhouse gas emission control, if adopted in current power plant systems, will drive the efficiency down and increase the cost of electricity, due to the energy and capital-intensive nature of current CO 2 separation techniques. This highlights the need to develop technologies that present a solution to the rising cost of electricity in the future. The integrated gasification combined cycle (IGCC) presents an improvement to the above predicament, but it is capital intensive due to the extensive unit operations involved. The efficiency for a CO 2 capture incorporated IGCC system is around ~32%, resulting in approximately 45% increase in cost of electricity. This is definitely a better path to traverse than the conventional pulverized coal (PC) power plants as it has improved economics and efficiencies accompanied by product flexibility. The SCL process advances the benefits even further when integrated with the IGCC process. The SCL process removes the expensive WGS system and CO 2 separation columns, thereby saving on capital and operational expenses. This integration results in a 12 – 21% increase in efficiency accompanied by a reduction in cost of electricity by 15 – 28%, over conventional IGCC systems. Therefore, SCL process provides the best route to harness the energy from coal. The overall project objective was to construct and operate a syngas chemical looping pilot scale test unit at the NCCC. The project scope of work was divided into 3 phases. In Phase I, cold flow model studies using a 1:1 scale acrylic test unit constructed at Particulate Solids Research, Inc. (PSRI) was successfully completed confirming robust solid flow control is achievable in the non-mechanical system design. In Phase II, the high pressure, high temperature, chemical looping reactor was successfully designed with all necessary equipment and safety instrumentation/controls specified to allow for fabrication, site construction, and assembly to commence in Phase III. In Phase III, the SCL pilot plant was successfully assembled and all necessary functional checks and pre-startup safety reviews completed. During unit commission over the course of 200 hours of testing, operational issues were observed with premix burner used for system startup. With the grants awarded by the National Energy Laboratory (NETL) and the Ohio Development Services Agency (ODSA) under awards DE-FE0023915 and D-14-18, the SCL pilot unit underwent three auxiliary equipment modifications to resolve all of the startup operational issues encountered. The unit was successfully demonstrated with 300+hr continuous operation. Key results obtained include high syngas conversion of 97.95% with 16.03% oxygen carrier conversion in the moving bed reducer and >99% purity H 2 produced from the moving bed oxidizer.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2016

Assessment of a Smoothed Particle Hydrodynamics Method for Thin Film on Textured Surfaces

A large number of previous studies suggest that the behavior of the near-wall liquid micro-layer plays a key role in boiling heat transfer and boiling crisis; see examples. The thin liquid layer interacts with heater surface and 'external' two-phase flow. The net effects of these interactions define the micro-hydrodynamics that leads to dryout. As a result, a good understanding of the thin-layer hydrodynamics is significant for the accurate prediction of CHF. Recently, with the advent of nanotechnology including micro-scale manufacturing, it becomes possible to perform high heat-flux boiling experiments with controlled surface conditions. Although many reports show rich observations about the effects of surface characteristics on CHF, physical mechanisms remain not fully understood. In order to have better understanding of the above phenomena, several computational models are developed, ranging from the lubrication theory to CFD methods. The former models are limited to one-dimensional formulation and unable to deal with the dry-out phenomena. It turns out that CFD methods have a greater potential for the mechanistic study of boiling. Towards the study objective, a particle-based CFD method called Smoothed Particle Hydrodynamics (SPH) is adopted. SPH is a Lagrangian particle method for solving systems of partial differential equations. Originally SPH is proposed for fluid dynamics in the context of astrophysical problems. Because mass and density are introduced during the particle approximation, SPH is also applied in fluid dynamic problems, where mass and density are important parameters to determine, without mass loss (a limitation of interface capturing schemes in mesh-based methods). Besides, geometrically complex and/or dynamic boundaries and interfaces can be handled without undue difficulty in SPH. Currently, most SPH simulations focus on large-scale scenarios (meter scale), like flooding and solitary wave. Also some works have been done with SPH in pore-scale. As for micro- or nano-scale simulation, people usually apply Dissipative Particle Dynamics (DPD) and Molecular Dynamics (MD) for particle-based simulation. Almost no work has been done on micro-scale with SPH, so assessments are needed to investigate the capability of SPH in application of micro-hydrodynamics. The Phenomena Identification and Ranking Table (PIRT) is applied to align application's technical requirements with associated physics and materials, models, simulation code capabilities and verification, validation and uncertainty quantification activities. The figures of merit are determined as: (1) Timing of micro-layer rewetting dry surface; (2) Effects of surface features. The corresponding key phenomena are determined as: (1) Hydrodynamics; (2) Interface Dynamics; (3) Energy; (4) Surface Effect. This study deals with isothermal scenario and no energy term is considered. Based on the current needs, a new SPH package, Micro-SPH, is compiled based on SPHysics with surface tension and wettability built in. Test cases are performed with Micro-SPH to assess the capability of SPH to each phenomenon identified. NEUTRINO with incompressible SPH solver is employed to investigate micro-fluid dynamics on a surface with uniformly distributed cylindrical pillars for application purposes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

The evolution of precipitate crystal structures in an Al-Mg-Si(-Cu) alloy studied by a combined HAADF-STEM and SPED approach

Highlights: • Scanning electron diffraction and machine learning complement lattice imaging • Developed approach gives increased objectivity in precipitate phase analysis • Precipitate phase fractions and solute content obtained with superior statistics • 0.01 at.% Cu changes the Al-Mg-Si system precipitation during overageing • Q’-phase eventually dominates precipitate structures and shows varying Cu-occupancy - Abstract: This work presents a detailed investigation into the effect of a low Cu addition (0.01 at.%) on precipitation in an Al-0.80Mg-0.85Si alloy during ageing. The precipitate crystal structures were assessed by scanning transmission electron microscopy combined with a novel scanning precession electron diffraction approach, which includes machine learning. The combination of techniques enabled evaluation of the atomic arrangement within individual precipitates, as well as an improved estimate of precipitate phase fractions at each ageing condition, through analysis of a statistically significant number of precipitates. Based on the obtained results, the total amount of solute atoms locked inside precipitates could be approximated. It was shown that even with a Cu content close to impurity levels, the Al-Mg-Si system precipitation was significantly affected with overageing. The principal change was due to a gradually increasing phase fraction of the Cu-containing Q′-phase, which eventually was seen to dominate the precipitate structures. The structural overtake could be explained based on a continuous formation of the thermally stable Q′-phase, with Cu atomic columns incorporating less Cu than what could potentially be accommodated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2022

Quasinormal g-modes of neutron stars with quarks

Quasinormal oscillation modes of neutron stars provide a means to probe their interior composition using gravitational wave astronomy. We compute the frequencies and damping times of composition-dependent core 𝑔-modes of neutron stars containing quark matter employing linearized perturbative equations of general relativity. We find that ignoring background metric perturbations due to the oscillating fluid, as in the Cowling approximation, underestimates the 𝑔-mode frequency by up to 10% for higher mass stars, depending on the parameters of the nuclear equation of state and how the mixed phase is constructed. The 𝑔-mode frequencies are well described by a linear scaling with the central lepton (or combined lepton and quark) fraction for nucleonic (hybrid) stars. Furthermore, our findings suggest that neutron stars with and without quarks are manifestly different with regards to their quasinormal 𝑔-mode spectrum, and may thus be distinguished from one another in future observations of gravitational waves from merging neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The short-term comprehensive impact of the phase-out of global coal combustion on air pollution and climate change

With the continuous intensification of global warming, the reduction and ultimate phase-out of coal combustion is an inevitable trend in the future global energy transformation. This study comprehensively analyzed the impact of phasing out coal combustion on global emissions and concentrations of air pollutants, radiative fluxes, meteorology and climate using Community Earth System Model 2 (CESM2). The results indicate that after the global phase-out of coal combustion, there is a marked decrease in the concentrations of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ) and fine particulate matter (PM 2.5 ), with some regions experiencing a reduction of exceeding 50%. There is no significant change in global ozone (O 3 ) concentration. There are decreasing AOD and positive radiative fluxes globally in the short term, though the cloud contributes minor negative radiative fluxes. The global air temperature may increase by approximately (0.02 ± 0.15) °C on average with regional and seasonal variations, and the precipitation may potentially increase by approximately (2.7 ± 40.6) mm yr −1 globally and over 20% in equatorial regions in the short term. But combined with the decreasing trend of cloud water content in the Northern Hemisphere, it indicates a potential increase in the extremity of precipitation events. In conclusion, this study provides references for global control of air pollution, mitigation strategies of climate change, and transformation of energy structures under the objective of “carbon neutrality”, such as focusing on the negative climate impacts of exacerbating regional warming and increasing extreme precipitation resulting from the rapid reduction of aerosols in the short term.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Integrated Solar Thermochemical Reaction System (Final Report)

This report describes progress accomplished in a three-phase project that has advanced the concept for a Solar Thermochemical Advanced Reactor System (STARS) – a system that is designed to convert solar energy into storable/useful chemical energy – from Technology Readiness Level 3 to Technology Readiness Level 6. STARS accomplishes its objective by directing the heat from a dish concentrator onto a chemical reaction system, increasing the energy content of the feedstock, methane, a chemical fuel that is available from multiple sources including natural gas, landfills and anaerobic digesters. The system uses a parabolic dish solar concentrator, previously developed for electrical power generation, and a compact, process-intensive chemical reaction system based on micro- and meso-channel process technology (MMPT). The combined Dish-STARSTM system provides a solar augment to the incoming methane stream, increasing its chemical energy content by 20-30% while decreasing its carbon intensity. When the immediate chemical product, known as syngas, is directly used for power generation the result is electricity with approximately 20% less CO 2 emissions. Alternately, if the syngas is further reacted to produce valuable chemical products, such as hydrogen, the reduction in carbon emissions is retained and reduced carbon intensities can be attained for the chemical products. In a co-production mode, lowcarbon hydrogen or electricity plus various hydrocarbons (for example, methanol, olefins or plastics) can be produced. Over the course of this project, the efficiency of Dish-STARSTM was improved from a previous value of 63% to slightly over 70%, with results and analyses indicating a potential for values greater than 80%; manufacturing investigations identified fabrication methods for the MMPT components, proofed-out additive manufacturing for microchannel heat exchangers, developed equipment cost models and, together with on-sun testing, supported technoeconomic evaluations. If developed into commercial products, project results indicate the potential for the efficient use of concentrated solar energy in nearterm applications for electrical power generation and the production of chemicals, including hydrogen.

14 SOLAR ENERGY↗
DOE OSTI2026

Beam Performance of the Positron Transport Line for CEBAF Positron Upgrade

The Low Energy Recirculator Facility (LERF) at Jefferson Lab, formerly operated for the Free-Electron Laser program, has been proposed as the injector complex for the planned 12 GeV CEBAF positron upgrade (Ce+BAF), with an additional pathway to support a potential 22 GeV CEBAF electron upgrade. In this configuration, LERF would generate and pre-accelerate positrons to 123 MeV, matching the present injection energy into the North Linac. Due to the relatively large emittance expected from the positron source, a comprehensive acceptance study has been performed from LERF through the CEBAF recirculating linacs and beam transport lines to the experimental halls. The objective is to establish the positron phase-space acceptance and provide design feedback to the positron production and capture systems. Furthermore, given CEBAF’s capability to deliver highly polarized beams, spin-tracking simulations have been carried out including magnet imperfections, alignment errors, and synchrotron-radiation–induced energy spread. Particular attention is given to the evolution of the spin tune and the corresponding depolarization mechanisms along the beam delivery path, especially for providing longitudinal polarization at the experimental halls. These results inform injector design choices and assess the overall feasibility of delivering high-polarization positron beams in CEBAF.

Ogur, S. [Thomas Jefferson National Accelerator Fa↗
DOE OSTI2022

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗
DOE OSTI2022

Universal relations for neutron star 𝑓-mode and 𝑔-mode oscillations

Among the various oscillation modes of neutron stars, 𝑓- and 𝑔- modes are the most likely to be ultimately observed in binary neutron star mergers due to their relatively large coupling and shared frequencies with tidal excitations. The 𝑓-mode frequency and damping time are known to correlate in normal neutron stars with their compactness, and previous fits to hadronic stars are extended and shown to be valid for an extremely broad sampling of equations of state using a piecewise polytropic parametrization scheme for hadrons and a constant sound-speed parametrization for quark matter. Separate fits applicable to quark (self-bound) stars are improved. Much more significant correlations exist with tidal deformability, and therefore with moment of inertia and quadrupole moment. It is conclusively demonstrated that these correlations are the same for all types of stars, whether hadronic, hybrid, or pure quark, and its accuracy is quantified. A novel 1-node branch of the 𝑓-mode that occurs in low-mass hybrid stars in a narrow mass range just beyond the critical mass necessary for a phase transition to appear is identified. This 1-node branch shows the largest, but still small, deviations from the universal correlation for any configuration. It is characterized by a nonmonotonic relation between neutron star mass and 𝑓-mode frequency, in contrast to the behavior otherwise observed in normal, quark and hybrid stars. The 𝑔-mode only exists in matter with a nonbarotropic equation of state involving temperature, chemical potential or composition (such as being out of beta equilibrium), or a phase transition in barotropic matter. Here, the 𝑔-mode therefore could serve as a probe for studying phase transitions in hybrid stars. In contrast with the 𝑓-mode, 𝑔-mode frequencies do not correlate well with tidal deformability, but depend strongly on properties of the transition (the density and the magnitude of the discontinuity) at the transition. Imposing causality and maximum mass constraints, a fit involving neutron star and phase transition properties is found and the 𝑔-mode frequency is determined to have an upper bound of about 1.25 kHz. However, if the sound speed 𝑐𝑠 in the inner core at densities above the phase transition density is restricted to 𝑐$^{2}_{𝑠}$ ≤1/3, 𝑔-mode frequencies can only reach about 0.8 kHz, which are significantly lower than 𝑓-mode frequencies (1.3–2.8 kHz). 𝑔-mode gravitational wave damping times are found to be extremely long, >10 4 s (102 s) in the inner core with 𝑐$^{2}_{𝑠}$ ≤1/3 (1), in comparison with 𝑓-mode damping times (0.1–1 s).

Composition of astronomical objects↗
DOE OSTI2021

Scalable, In-situ Data Clustering Data Analysis for Extreme Scale Scientific Computing (Final Report)

The objective of this project is to address challenges in the design and development of scalable in-situ data clustering and analytics algorithms and software. Our goal is to develop parallel software consisting of a set of spatio-temporal data clustering and anomaly detection functions, both of which are very important for large-scale analysis and have wide applicability for in-situ runs as well as post-processing analysis. Our design principles for in-situ analysis consider the following: (1) identify parts of the computation can be done close to the data within the nodes, while it is still in memory; (2) extract analysis components can (and should) be performed in remote staging and analysis nodes; (3) develop error-bound approximation methods for applications tolerable for small errors; (4) identify the type of derived distributions and statistics, for spatio-temporal data, that can be kept locally in order to both accelerate computations and meet energy constraints in subsequent iterations and phases; (5) use a self-describing data format so that data can be consistent and understood among local storage (memory and SSDs) and at staging and analysis nodes, thereby providing portability and flexibility; (6) develop service-oriented functions that can schedule in-situ and post-hoc analysis tasks based on the dynamic requirements of applications. Our development focus is to produce the parallel data analysis software/library that will be scalable, reusable, extensible, and generic for applications in different disciplines. The software will be able to run in-situ with the simulations as well as post-hoc analysis. This approach will satisfy many synergistic requirements for data intensive applications executed on data coming from instruments and experiments. In particular, the proposed multilevel approach is directly applicable to perform design tradeoffs for running part of the algorithms near the instruments and the rest on remote (analysis) systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems

The objective of this paper is to present a local bounds preserving stabilized finite element scheme for hyperbolic systems on unstructured meshes based on continuous Galerkin (CG) discretization in space. A CG semi-discrete scheme with low order artificial dissipation that satisfies the local extremum diminishing (LED) condition for systems is used to discretize a system of conservation equations in space. The low order artificial diffusion is based on approximate Riemann solvers for hyperbolic conservation laws. In this case we consider both Rusanov and Roe artificial diffusion operators. In the Rusanov case, two designs are considered, a nodal based diffusion operator and a local projection stabilization operator. The result is a discretization that is LED and has first order convergence behavior. To achieve high resolution, limited antidiffusion is added back to the semi-discrete form where the limiter is constructed from a linearity preserving local projection stabilization operator. The procedure follows the algebraic flux correction procedure usually used in flux corrected transport algorithms. To further deal with phase errors (or terracing) common in FCT type methods, high order background dissipation is added to the antidiffusive correction. The resulting stabilized semi-discrete scheme can be discretized in time using a wide variety of time integrators. Numerical examples involving nonlinear scalar Burgers equation, and several shock hydrodynamics simulations for the Euler system are considered to demonstrate the performance of the method. In conclusion, for time discretization, Crank–Nicolson scheme and backward Euler scheme are utilized.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Successful In Situ Remediation of Cr(VI)-Impacted Deeper Vadose Zone Soils for Groundwater Protection - 18603

Shallow soil at a former industrial wood treating facility was impacted with chromium, arsenic, and copper above cleanup goals and shallow soil has been excavated for protection of direct contact exposure pathways. Deeper vadose zone soil was predominantly impacted with chromium and was undergoing remediation for long-term groundwater protection. The 920-meter (m) long Cr(VI) groundwater plume has undergone full-scale remediation with groundwater recirculation and in situ bioremediation. Following approximately 4 years of operation, the Cr(VI) groundwater plume, with former source area groundwater concentrations over 10,000 micrograms per liter (μg/L), has been reduced to below the laboratory method reporting limit (0.50 μg/L). Approximately 23,000 cubic meters (m{sup 3}) of deeper (4.5 m - 11 m below ground surface, bgs) vadose zone soil required Cr(VI) remediation. The remediation strategy relied on reduction of Cr(VI) and geochemical fixation as trivalent chromium (Cr(III)) and iron hydroxide precipitates. Chromium reduction was facilitated by a carbon source (corn syrup and methanol) to establish and maintain biogeochemical reducing conditions and ferrous iron to support direct Cr(VI) chemical reduction. Combined, the injected soil moisture, carbon source and ferrous iron supported fast and sustained ferrous/ferric iron cycling and Cr(VI) reduction. Remedy performance monitoring included time-series water quality data from a network of vacuum lysimeters and confirmation soil borings. Reagent was applied either through liquid injection wells and direct-push injection points, or by deep soil mixing methodology in areas with difficult geology (i.e., fine-grained soils and discontinuous iron-silica hardpan layer). Following both phases of deep soil remediation, soils met long-term protection of groundwater remediation action objectives. In situ treatment of soils provides a cost-effective remedial solution compared to excavation, specifically when impacts are deep and long-term protection of groundwater is a primary remediation goal. The delivery methods employed in this case study can be applied to sites impacted with a broad range of contaminants amendable to geochemical fixation (e.g., chromium, uranium, strontium, nitrate, etc.). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗