Search NASA⌕ Search

SEARCH · Search NASA

Results for “Taylor series”

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 55 records · Page 3

Exponential Time Differencing Schemes for Fuel Depletion and Transport in Molten Salt Reactors: Theory and Implementation

A numerical framework for modeling depletion and mass transport in liquid-fueled molten salt reactions is presented based on exponential time differencing. The solution method involves using the finite volume method to transform the system of partial differential equations (PDEs) into a much larger system of ordinary differential equations. The key part of this method involves solving for the exponential of a matrix. We explore six different algorithms to compute the exponential in a series of progression problems that explore physical transport phenomena in molten salt reactors. This framework shows good results for solving linear parabolic PDEs with each of the six matrix exponential algorithms. For large problems, the series solvers such as Padé and Taylor have large run times, which can be mitigated by using the Krylov subspace.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Report Viewgraphs for IC project: Is high resolution or a more realistic bathymetry necessary to simulate Maud Rise polynyas in the Southern Ocean?

Maud Rise and Weddell Sea polynyas in the Southern Ocean form because of a series of feedbacks between upper ocean high salinity anomalies, Taylor columns found near steep bathymetric features such the Maud Rise seamounts and the Astrid Ridge complex, and subsurface accumulation of ocean heat that can sustain ice openings once deep convection is triggered.

58 GEOSCIENCES↗

L4-3 Developing DOE M441.1-1 Compliant Container Components: Filter Cup Experimental Fixture Test Data Verification

This document serves to provide technical documentation and justification for using the Filter Cup Experimental Fixture (FCEP) for filter efficiency and pressure drop testing in the Enhanced Filter Test System (EFTS). This document assumes that the reader has prior knowledge related to both of these systems, as well as the SAVY-4000® nuclear material container series filter technology.

42 ENGINEERING↗

The Consequences of Rayleigh-Taylor Instabilities on Implosion Dynamics in Gas-puff Z-pinch Experiments (Final Report)

The Laboratory of Plasma Studies (LPS), Cornell University, proposed to carry out a comprehensive series of gas-puff z-pinch experiments on a 1 MA pulsed power machine for the purpose of understanding the dynamics of imploding magnetized high energy density plasmas in the presence of Rayleigh-Taylor Instabilities (RTI). Cornell’s pulsed power facilities, including the associated diagnostic instruments, are largely supported under the auspices of the Center for Pulsed-power-driven High-energy-density Plasmas (the Center), an NNSA-sponsored Center of Excellence. The proposed experiments were carried out predominantly by LPS graduate student, Sophia Rocco (the “primary grad student”), in collaboration with other graduate students and the scientific and technical staff of the Center using Cornell’s 1 MA COBRA pulsed power facility and all of its associated diagnostic instruments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementing Spatio-Temporal 3D-Convolution Neural Networks and UAV Time Series Imagery to Better Predict Lodging Damage in Sorghum

Unmanned aerial vehicle (UAV)-based remote sensing is gaining momentum in a variety of agricultural and environmental applications. Very-high-resolution remote sensing image sets collected repeatedly throughout a crop growing season are becoming increasingly common. Analytical methods able to learn from both spatial and time dimensions of the data may allow for an improved estimation of crop traits, as well as the effects of genetics and the environment on these traits. Multispectral and geometric time series imagery was collected by UAV on 11 dates, along with ground-truth data, in a field trial of 866 genetically diverse biomass sorghum accessions. We compared the performance of Convolution Neural Network (CNN) architectures that used image data from single dates (two spatial dimensions, 2D) versus multiple dates (two spatial dimensions + temporal dimension, 3D) to estimate lodging detection and severity. Lodging was detected with 3D-CNN analysis of time series imagery with 0.88 accuracy, 0.92 Precision, and 0.83 Recall. This outperformed the best 2D-CNN on a single date with 0.85 accuracy, 0.84 Precision, and 0.76 Recall. The variation in lodging severity was estimated by the best 3D-CNN analysis with 9.4% mean absolute error (MAE), 11.9% root mean square error (RMSE), and goodness-of-fit (R2) of 0.76. This was a significant improvement over the best 2D-CNN analysis with 11.84% MAE, 14.91% RMSE, and 0.63 R2. The success of the improved 3D-CNN analysis approach depended on the inclusion of “before and after” data, i.e., images collected on dates before and after the lodging event. The integration of geometric and spectral features with 3D-CNN architecture was also key to the improved assessment of lodging severity, which is an important and difficult-to-assess phenomenon in bioenergy feedstocks such as biomass sorghum. This demonstrates that spatio-temporal CNN architectures based on UAV time series imagery have significant potential to enhance plant phenotyping capabilities in crop breeding and Precision agriculture applications.

3D-convolution neural networks↗

Searching for the QCD Critical Point Along the Pseudo-critical/Freeze-out Line Using Padé-resummed Taylor Expansions of Cumulants of Conserved Charge Fluctuations

Using high-statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we construct estimators for the radius of convergence from an eighth order series expansion of the pressure as well as the number density. We show that the estimator for pressure and number density will be identical in the asymptotic limit. In the vicinity of the pseudo-critical temperature, T pc ≃ 156.5 MeV, we find the estimator of the radius of convergence to be µ B /T ≳ 3 for strangeness-neutral matter. We also present results for the pole structure of the Pad´e approximants for the pressure at non-zero values of the baryon chemical potential and show that the pole structure of the [4,4] Pad´e is consistent with not having a critical point at temperatures larger than 135 MeV and a baryon chemical potential smaller than µ B /T ~ 2.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

cppTPSA/pyTPSA: a C++/Python package for truncated power series algebra

The truncated power series algebra (TPSA), also referred to as differential algebra (DA), is a well-established and widely used method in particle accelerator physics and astronomy. The most straightforward usage of TPSA/DA is to calculate the Taylor expansion of a given function at a specific point up to order ?. In recent years, as the application of TPSA/TA has been extended to other fields, a reusable implementation of TPSA/DA as a modern C++ library or other high level programming language like Python has become desirable. The cppTPSA package implements TPSA/DA in C++11 and provides developers a convenient library with which to build advanced TPSA/DA-based methods. A Python 3 library, pyTPSA, has also been developed based on the C++ lib.

97 MATHEMATICS AND COMPUTING↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗

Optimization of TAM16, a Benzofuran That Inhibits the Thioesterase Activity of Pks13; Evaluation toward a Preclinical Candidate for a Novel Antituberculosis Clinical Target

With increasing drug resistance in tuberculosis (TB) patient populations, there is an urgent need for new drugs. Ideally, new agents should work through novel targets so that they are unencumbered by preexisting clinical resistance to current treatments. Benzofuran 1 was identified as a potential lead for TB inhibiting a novel target, the thioesterase domain of Pks13. Although, having promising activity against Mycobacterium tuberculosis, its main liability was inhibition of the hERG cardiac ion channel. This article describes the optimization of the series toward a preclinical candidate. Despite improvements in the hERG liability in vitro, when new compounds were assessed in ex vivo cardiotoxicity models, they still induced cardiac irregularities. Further series development was stopped because of concerns around an insufficient safety window. However, the demonstration of in vivo activity for multiple series members further validates Pks13 as an attractive novel target for antitubercular drugs and supports development of alternative chemotypes.

60 APPLIED LIFE SCIENCES↗

Analysis as a Key Guiding Tool for Waste Carbon Utilization

Over a billion metric tons of waste and biomass are projected to be available in a future mature market in the United States. These resources represent an opportunity to decouple chemical and polymer production from conventional fossil fuel feedstocks, but such a broad solution space can also make for challenging decision-making. This talk will provide researchers with an introduction to key analysis techniques such as techno-economic analysis, life cycle assessment, and material flow analysis: how they are conducted and how they can be used to benchmark the costs, environmental impacts, and circularity of new innovations as well as to identify opportunities for prioritization and optimization. Using a series of examples related to plastic recycling and chemical manufacturing, we will explore how analysis can guide where and how to leverage waste carbon in supply chains towards a future circular economy.

chemical↗

Particle transport in Markov stochastic mixtures with spatial gradients

Markov geometries are a prototype class of stochastic media that are widely used to model complex disordered systems. In a series of recent works, we have shown that spatially homogeneous (but possibly non-isotropic) Poisson tessellations can be conveniently used in order to generate three-dimensional realizations of Markov media and thus characterize the features of particle transport e.g. in fragmented fuel elements following severe nuclear accidents or Rayleigh-Taylor turbulent mixing for inertial confinement fusion. In this paper we show that Poisson tessellations can be extended to take into account the presence of spatial gradients, which commonly occur in real-world applications due to material stratification. For this purpose, we will provide an explicit construction of spatially-non-homogeneous Poisson tessellations, fully characterize their statistical properties, and finally illustrate the behaviour of particle transport in such random media. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electron vacancy-level dependent hybrid photoionization of the F–@C60+molecule: A novel effect

Our previous studies [J. Phys. B 53, 125101 (2020); Euro. Phys. J. D 74, 191 (2020)] have predicted that the atom-fullerene hybrid photoionization properties for X = Cl, Br and I endohedrally confined in C60 are different before and after an electron transfers from C60 to the halogen. It was further found as a rule that the ionization dynamics is insensitive to the C60level the electron originates from to produce X–@C60+ . In the current study, we report an exception to this rule in F@C60. It is found that when the electron vacancy is situated in the C60level that participates in the hybridization in F–@C60+, the mixing becomes dramatically large leading to strong modifications in the photoionization of the hybrid levels. This novel effect is fundamentally based on a level-crossing phenomenon driven by the electron transfer in F@C60. But when the vacancy is at any other pure level of C60, the level-invariance is retained showing weak hybridization. Even though this case of F@C60 is an exception in the halogen@C60 series, the phenomenon can be more general and can occur with compounds of other atoms caged in a variety of fullerenes. Possible experimental studies are suggested to benchmark the present results.

74 ATOMIC AND MOLECULAR PHYSICS↗

Building Interfaces with Echo Using R

The following report details the functionality that was developed to directly query the MongoDB and read in echo records from within R, and the ability to create a cinema database directly from echo records. Both of these efforts used the data collected under the MINOS project for development. The nuclear facility under consideration is the High-Flux Isotope Reactor (HFIR) and co-located Radiochemical Engineering Development Center (REDC) at Oak Ridge National Laboratory in Oak Ridge, TN. HFIR is an 85 MW research reactor and is used primarily for production of a medical radioisotopes, material irradiation experiments, neutron activation, and neutron scattering. Targets for the reactor are fabricated and processed and dissolved at REDC, in addition to other glove-box and hot-cell type activities. An overhead view of the facility is provided. The rest of the report is organized as follows: Section 2 describes how to query the MongoDB and load in echo records from within R, with background information on the database type, MongoDB, and the file format of echo records, hdf5, as well as lessons learned during the development phase. Section 3 details the functionality and usage of creating a cinema database directly from echo records within R. A series of appendixes display all the R code (Section 5), some example scripts where the code is used (Section 6), and useful links for further exploration (Section 7).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗