Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conjunctions”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Dissipation-preserving spectral element method for damped seismic wave equations

This article describes the extension of the conformal symplectic method to solve the damped acoustic wave equation and the elastic wave equations in the framework of the spectral element method. The conformal symplectic method is a variation of conventional symplectic methods to treat non-conservative time evolution problems, which has superior behaviors in long-time stability and dissipation preservation. To reveal the intrinsic dissipative properties of the model equations, we first reformulate the original systems in their equivalent conformal multi-symplectic structures and derive the corresponding conformal symplectic conservation laws. We thereafter separate each system into a conservative Hamiltonian system and a purely dissipative ordinary differential equation system. Based on the splitting methodology, we solve the two subsystems respectively. The dissipative one is cheaply solved by its analytic solution. While for the conservative system, we combine a fourth-order symplectic Nyström method in time and the spectral element method in space to cover the circumstances in realistic geological structures involving complex free-surface topography. The Strang composition method is adopted thereby to concatenate the corresponding two parts of solutions and generate the completed conformal symplectic method. A relative larger Courant number than that of the traditional Newmark scheme is found in the numerical experiments in conjunction with a spatial sampling of approximately 5 points per wavelength. A benchmark test for the damped acoustic wave equation validates the effectiveness of our proposed method in precisely capturing dissipation rate. The classical Lamb problem is used to demonstrate the ability of modeling Rayleigh wave in elastic wave propagation. More comprehensive numerical experiments are presented to investigate the long-time simulation, low dispersion and energy conservation properties of the conformal symplectic methods in both the attenuating homogeneous and heterogeneous media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach

Here, a hybrid, design-order sliding mesh algorithm, which uses a control volume finite element method (CVFEM), in conjunction with a discontinuous Galerkin (DG) approach at non-conformal interfaces, is outlined in the context of a low-Mach fluid dynamics equation set. This novel hybrid DG approach is also demonstrated to be compatible with a classic edge-based vertex centered (EBVC) scheme. For the CVFEM, element polynomial, P , promotion is used to extend the low-order P = 1 CVFEM method to higher-order, i.e., P = 2. An equal-order low-Mach pressure-stabilized methodology, with emphasis on the non-conformal interface boundary condition, is presented. A fully implicit matrix solver approach that accounts for the full stencil connectivity across the non-conformal interface is employed. A complete suite of formal verification studies using the method of manufactured solutions (MMS) is performed to verify the order of accuracy of the underlying methodology. The chosen suite of analytical verification cases range from a simple steady diffusion system to a traveling viscous vortex across mixed-order non-conformal interfaces. Results from all verification studies demonstrate either second- or third-order spatial accuracy and, for transient solutions, second-order temporal accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Application of linear multifrequency-grey acceleration to preconditioned Krylov iterations for thermal radiation transport

The thermal radiative transfer (TRT) equations comprise a radiation equation coupled to the material internal energy equation. Linearization of these equations produces effective, thermally-redistributed scattering through absorption-reemission. In this paper, we investigate the effectiveness and efficiency of Linear-Multi-Frequency-Grey (LMFG) acceleration that has been reformulated for use as a preconditioner to Krylov iterative solution methods. We introduce two general frameworks, the scalar flux formulation (SFF) and the absorption rate formulation (ARF), and investigate their iterative properties in the absence and presence of true scattering. SFF has a group-dependent state size but may be formulated without inner iterations in the presence of scattering, while ARF has a group-independent state size but requires inner iterations when scattering is present. We compare and evaluate the computational cost and efficiency of LMFG applied to these two formulations using a direct solver for the preconditioners. Finally, this work is novel because the use of LMFG for the radiation transport equation, in conjunction with Krylov methods, involves special considerations not required for radiation diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Flow-driven spectral chaos (FSC) method for simulating long-time dynamics of arbitrary-order non-linear stochastic dynamical systems

Uncertainty quantification techniques such as the time-dependent generalized polynomial chaos (TD-gPC) use an adaptive orthogonal basis to better represent the stochastic part of the solution space (aka random function space) in time. However, because the random function space is constructed using tensor products, TD-gPC-based methods are known to suffer from the curse of dimensionality. Here, we introduce a new numerical method called the flow-driven spectral chaos (FSC) which overcomes this curse of dimensionality at the random-function-space level. The proposed method is not only computationally more efficient than existing TD-gPC-based methods but is also far more accurate. The FSC method uses the concept of enriched stochastic flow maps to track the evolution of a finite-dimensional random function space efficiently in time. To transfer the probability information from one random function space to another, two approaches are developed and studied herein. In the first approach, the probability information is transferred in the mean-square sense, whereas in the second approach the transfer is done exactly using a new theorem that was developed for this purpose. The FSC method can quantify uncertainties with high fidelity, especially for the long-time response of stochastic dynamical systems governed by ODEs of arbitrary order. Six representative numerical examples, including a nonlinear problem (the Van-der-Pol oscillator), are presented to demonstrate the performance of the FSC method and corroborate the claims of its superior numerical properties. Finally, a parametric, high-dimensional stochastic problem is used to demonstrate that when the FSC method is used in conjunction with Monte Carlo integration, the curse of dimensionality can be overcome altogether.

(nonlinear) stochastic dynamical systems↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability-based optimization of ACRT for the growth of CZT by the traveling heater method

Numerical modeling is employed to study the effects of several accelerated crucible rotation technique (ACRT) rotation schedules on the growth of cadmium zinc telluride (CZT) by the traveling heater method (THM). In conjunction with these analyses, a measure is developed to estimate the amplitude of disturbances associated with cellular interface growth morphologies that develop over time and across the melt-solid interface. Furthermore, this disturbance amplitude measure is used to assess the ability of ACRT to stabilize the growth interface and thereby minimize inclusion formation processes. Notably, our analyses reveal that interfacial instability is not well correlated with melt undercooling, explaining why classical ACRT approaches to reduce undercooling by mixing may not stabilize growth. Computations show that interface stability is best achieved via rotation schedules that accentuate outward Ekman flows during spin-up and, during spin-down, eliminate Taylor-Görtler flows and minimize inward Ekman flows.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Overview and technology opportunities for thermochemically-produced bio-blendstocks

Global demand for transportation fuels is projected to increase 40% by 2040, and biomass-derived fuels (biofuels) play a crucial role in substituting fossil fuels and mitigating greenhouse gas emissions. Currently, biofuels are mainly consumed as blendstocks combined with petroleum-based fuels, and effective conversion technologies can address the quality challenges for offering standalone biofuels. Thermochemical conversion process is one of the most promising pathways among existing technologies for biofuel production. However, the major barriers are unwanted characteristics (e.g., thermal instability) of intermediate products, such as bio-oil, and required upgrading treatments for producing compatible fuels. Here, this study highlights the merits and critical challenges of thermochemical conversion and physicochemical upgrading technologies for bio-blendstock production from lignocellulosic biomass. The novelty of this study lies in potential directions for future research through both critical and systematic literature reviews, and the proposed intensified process for lignocellulosic-based fuel blendstocks production. It is concluded that recovery and fractionation strategies (e.g., quenching and stripping) can maximize process yields and add values in the efficient conversion pathways. Effective quenching can stop secondary free radical reactions and improve liquid yields over gas and solid yields. Stripping process can improve process yield, catalyst lifespan, and thermal stability. It is further concluded that physicochemical treatments are not as effective as thermochemical treatments, but have advantages of mild operating conditions and potential for integrated solutions in conjunction with other treatments.

09 BIOMASS FUELS↗
DOE OSTI2022

Testing and signal identification for two-sample high-dimensional covariances via multi-level thresholding

The paper considers testing and signal identification for covariance matrices from two populations of marginally sub-Gaussian distributed. A multi-level thresholding procedure is proposed for testing the equality of two high-dimensional covariance matrices, which is designed to detect sparse and faint differences between the covariances. A novel -statistic composition is developed to establish the asymptotic distribution of the thresholding statistics in conjunction with the matrix blocking and the coupling techniques. It is shown that the proposed test is more powerful than the existing tests in detecting sparse and weak signals in covariances. Multiple testing procedures are constructed to discover different covariances and the sub-groups of variables with different covariance structures between the two populations. The proposed procedures are based on the multi-level thresholding test, which are able to control the false discovery proportion () with high power. Simulation experiments and a case study on the returns of the S&P 500 stocks before and after the COVID-19 pandemic are conducted to demonstrate and compare the utilities of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Improved performance comparisons of radioxenon systems for low level releases in nuclear explosion monitoring

The Comprehensive Nuclear-Test-Ban Treaty bans all nuclear tests and mandates development of verification measures to detect treaty violations. One verification measure is detection of radioactive xenon isotopes produced in the fission of actinides. The International Monitoring System (IMS) currently deploys automated radioxenon systems that can detect four radioxenon isotopes. Radioxenon systems with lower detection limits are currently in development. Historically, the sensitivity of radioxenon systems was measured by the minimum detectable concentration for each isotope. In this paper we analyze the response of radioxenon systems using rigorous metrics in conjunction with hypothetical representative releases indicative of an underground nuclear explosion instead of using only minimum detectable concentrations. Our analyses incorporate the impact of potential spectral interferences on detection limits and the importance of measuring isotopic ratios of the relevant radioxenon isotopes in order to improve discrimination from background sources particularly for low-level releases. To provide a sufficient data set for analysis, hypothetical representative releases are simulated every day from the same location for an entire year. The performance of three types of samplers are evaluated assuming they are located at 15 IMS radionuclide stations in the region of the release point. The performance of two IMS-deployed samplers and a next-generation system is compared with proposed metrics for detection and discrimination using representative releases from the nuclear test site used by the Democratic People’s Republic of Korea.

Haas, Derek A.↗
DOE OSTI2022

The impact of tritium phytoremediation on plant health as measured by fluorescence

Phytoremediation, using plants for soil, sediment, or water contaminant clean-up, is an established technology dependent on plant health. Tritium (3H), a radioactive isotope of hydrogen that is generally found in the environment as tritiated water (HTO), is a low-level beta emitter with a half-life of 12.32 years. Chlorophyll fluorescence (CF) for monitoring risk assessment of tritium to plant health was conducted at the Tritium Irrigation Facility (TIF) located on the US Department of Energy's Savannah River Site (SRS) near Aiken, SC. Two fluorometers were evaluated in conjunction with phytoremediation at the 25 –acre TIF where tritiated groundwater is being spray–irrigated on a mixed coniferous/deciduous forested watershed as a means of reducing tritium release to a nearby stream that serves as a tributary to the Savannah River. Tritium activity in irrigated water averaged 104 + 42 pCi mL -1 during the 2003 project. Fluorescence parameters measured by the two fluorometers were well correlated with each other (p < 0.0001). Tritium in water respired from oak leaves ranged up to 1845.13 pCi ml -1 and 2138.22 pCi ml -1 in pine needles. Trees in both the test and control sites were approximately 15 years old. In this report we demonstrated that fluorescence parameters provide an effective way to estimate the impact of HTO on plant health in a noninvasive, extremely rapid, and cost-effective manner. In the current study applying fluorometry, plants within the TIF phytoremediation site exposed to the site tritiated water were not significantly impacted by the tritium phytoremediation based on CF parameters as compared to the control, a nascent non-irrigated site.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Quantification of grain boundary effects on the geometrically necessary dislocation density evolution and strain hardening of polycrystalline Mg 4Al using in situ tensile testing in scanning electron microscope and HR-EBSD

In situ tensile testing in a scanning electron microscope (SEM) in conjunction with high-resolution electron backscatter diffraction (HR-EBSD) under load was used to characterize the evolution of geometrically necessary dislocation (GND) densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg—4Al alloy. The increase in GND density was investigated at plastic strains of 0 %, 0.6 %, 2.2 %, 3.3 % from the area including 76 grains and correlated with (i) geometric compatibility between slip systems across grain boundaries, and (ii) plastic incompatibility. We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility, from which uniaxial tensile stress-strain response of polycrystalline Mg—4Al are computed and compared with experimental measurement. The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.

Geometrically necessary dislocations↗
DOE OSTI2021

Origins of high ductility exhibited by an extruded magnesium alloy Mg-1.8Zn-0.2Ca: Experiments and crystal plasticity modeling

Low ductility and strength are major bottlenecks against Mg alloys' wide applications. In this work, we systematically design the composition and fabrication process for a low-alloyed Mg-Zn-Ca alloy, showing that it can be extruded at low temperatures (~ 250 degrees C) and high speeds (~2 mm/s). After the extrusion, this alloy exhibits a substantially weakened basal texture, relatively small grain size, very high tensile elongation (~ 30%), and good strength. The origin of the considerably improved ductility was studied using a combination of three-dimensional atom probe tomography (3D-APT), transmission electron microscopy (TEM), electron backscattered diffraction (EBSD) in conjunction with surface slip trace analysis, in-situ synchrotron X-ray diffraction, and elasto-plastic self-consistent (EPSC) modeling. Co-segregation of Zn and Ca atoms at a grain boundary is observed and associated with texture weakening and grain boundary mediated plasticity, both improving the ductility. While basal slip and prismatic slip are identified as the dominant deformation systems in the alloy, the ratio between their slip resistances is substantially reduced relative to pure Mg and most other Mg alloys, significantly contributing to the improved ductility of the alloy. This Mg-Zn-Ca alloy exhibiting excellent mechanical properties and low fabrication cost is a promising candidate for industrial productions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An effective numerical algorithm for intra-granular fission gas release during non-equilibrium trapping and resolution

Both fission gas release and gaseous swelling in nuclear fuel are driven by the transport of fission gas from within the fuel grains to grain boundaries (intra-granular fission gas release). The process involves gas atom diffusion in conjunction with trapping in and resolution from intra-granular bubbles, and is described mathematically by a system of two partial differential equations (PDE). Following Speight, under the assumption of equilibrium between trapping and resolution (quasi-stationary approximation) the system can be reduced to a single diffusion equation with an effective diffusion coefficient. Numerical solutions used in engineering fuel performance calculations invariably rely on this simplification. First, we investigate the validity of the quasi-stationary approximation compared to the solution of the general system of PDEs. Results demonstrate that the approximation is valid under most conditions of practical interest, but is inadequate to model intra-granular fission gas release during rapid transients to relatively high temperatures such as postulated reactivity initiated accidents (RIA). Then, we develop a novel numerical algorithm for the solution of the general problem in time-varying conditions. We verify the PolyPole-2 algorithm against a reference finite difference solution for a large number of randomly generated operation histories including prototypical RIA transients. Results demonstrate that PolyPole-2 captures the solution of the general PDE system with a high accuracy and a low computational cost. As a result, the PolyPole-2 algorithm overcomes the quasi-stationary approximation and the concept of an effective diffusion coefficient for the solution of the intra-granular fission gas release problem in nuclear fuel analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Measurement and modeling of the gas permeability of high burnup pressurized water reactor fuel rods

Here, the transmissibility of gas through the pellet stack of full-length pressurized water reactor (PWR) high burnup (HBU) spent nuclear fuel (SNF) rods is measured at room temperature and correlated to a permeability. Gas transmission and depressurization testing of eight high burn-up fuel rods revealed that gas communication from one end of the pellet stack to the other is unobstructed, but slow, at room temperature. Two traditional flow models, Darcy's Law and Muskat's application of Poiseuille's Law to compressible gas flow through porous media are explored. The derived permeability values for the two methods are compared and are considered in conjunction with rod cladding type and various indicators of rod operation such as burnup, High Duty Core Index (HDCI), and average fuel temperature. The derived permeability appears to be correlated to fuel temperature and maximum HDCI, but not to rod average burnup. Finally, three of the HBU rods were heat-treated to simulate dry storage vacuum drying conditions prior to the measurements and the heat-treatments may have resulted in increased permeability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Evolution of B2 and laves phases in a ferritic steel under Fe 2+ ion irradiation at 475 °C

Microstructural evolution in a novel ferritic steel (Fe–12Cr–3W–3Ni–3Al–1Nb, in wt.%) computationally designed to contain B2 and Laves phases after 4 MeV Fe 2+ ion irradiation up to 220 dpa at 475 °C was characterized using transmission electron microscopy in conjunction with x-ray energy dispersive spectroscopy. The ferritic matrix phase exhibited dislocation loops and tangled dislocations, but our focus was on stability of two types of intermetallic precipitates. The B2–NiAl precipitates ~13 nm in size remained crystalline and appeared to have slightly lower Al concentration after irradiation. The Laves phase, (Fe,Cr) 2 (Nb,W), were present in two size ranges: coarse micron-scale precipitates which were amorphized with a slight composition change at irradiation doses above~30 dpa, while the finer precipitate particles~100 nm in size were partially disintegrated with a noticeable composition change at doses above~70 dpa. Meanwhile, many Nb/Cr-enriched particles ~8 nm in size formed within a few hundreds of nanometers from the disintegrated particles. The understanding of the phase stability would help design advanced steels and engineer microstructures that are stable against high irradiation doses, while retaining good high temperature strength.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗