Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic time warping”

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.

45 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Global-Scale Comparison of Passive (SMOS) and Active (ASCAT) Satellite Based Microwave Soil Moisture Retrievals with Soil Moisture Simulations (MERRA-Land)

Global surface soil moisture (SSM) datasets are being produced based on active and passive microwave satellite observations and simulations from land surface models (LSM). This study investigates the consistency of two global satellite-based SSM datasets based on microwave remote sensing observations from the passive Soil Moisture and Ocean Salinity (SMOS;SMOSL3 version 2.5) and the active Advanced Scatterometer (ASCAT; version TUWien- WARP 5.5) with respect to LSM SSM from the MERRA-Land data product. The relationship between the global-scale SSM products was studied during the 2010-2012 period using (1) a time series statistics (considering both original SSM data and anomalies), (2) a space-time analysis using Hovmoller diagrams, and (3) a triple collocation error model. The SMOSL3 and ASCAT retrievals are consistent with the temporal dynamics of modeled SSM (correlation R (is) greater than 0.70 for original SSM) in the transition zones between wet and dry climates, including the Sahel, the Indian subcontinent, the Great Plains of North America, eastern Australia, and southeastern Brazil. Over relatively dense vegetation covers, a better consistency with MERRA-Land was obtained with ASCAT than with SMOSL3. However, it was found that ASCAT retrievals exhibit negative correlation versus MERRA-Land in some arid regions (e.g., the Sahara and the Arabian Peninsula). In terms of anomalies, SMOSL3 better captures the short term SSM variability of the reference dataset (MERRA-Land) than ASCAT over regions with limited radio frequency interference (RFI) effects (e.g., North America, South America, and Australia). The seasonal and latitudinal variations of SSM are relatively similar for the three products, although the MERRALand SSM values are generally higher and their seasonal amplitude is much lower than for SMOSL3 and ASCAT. Both SMOSL3 and ASCAT have relatively comparable triple collocation errors with similar spatial error patterns: (i) lowest errors in arid regions (e.g., Sahara and Arabian Peninsula), due to the very low natural variability of soil moisture in these areas, and Central America, and (ii) highest errors over most of the vegetated regions (e.g., northern Australia, India, central Asia, and South America). However, the ASCAT SSM product is prone to larger random errors in some regions (e.g., north-western Africa, Iran, and southern South Africa). Vegetation density was found to be a key factor to interpret the consistency with MERRA-Land between the two remotely sensed products (SMOSL3 and ASCAT) which provides complementary information on SSM. This study shows that both SMOS and ASCAT have thus a potential for data fusion into long-term data records.

ASCAT↗
NASA NTRS1994

Self-gravity and dissipation in polar rings

Studies of inclined rings inside galaxy potentials have mostly considered the influence of self-gravity and viscous dissipation separately. In this study, we construct models of highly inclined ('polar') rings in an external potential including both self-gravity and dissipation due to a drag force. We do not include pressure forces and thus ignore shock heating that dominates the evolution of gaseous rings inside strongly nonspherical potentials. We adopt an oblate spheroidal scale-free logarithmic potential with axis ratio q = 0.85 and an initial inclination of 80 deg for the self-gravitating rings. We find that stellar (dissipationless) rings suffer from mass loss during their evolution. Mass loss also drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than approximately 0.02 of the enclosed galaxy mass, rings remain bound and do not fall apart from differential precession. The rings precess at a constant rate for more than a precession period tau(sub p) finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self-gravity if the characteristic viscous inflow time scale tau(sub vi) is shorter than approximately 25(tau(sub p)). Rings with (tau(sub vi))/(tau(sub p)) less than or approximately equal to 25 collapse toward the nucleus of the potential within one precession period independent of the amount of self-gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.

Dubinski, John↗
DOE OSTI2020

Wobbly discs – corrugations seen in the dust lanes of edge-on galaxies

ABSTRACT We report the detection of small-scale bending waves, also known as corrugations, in the dust lanes of five nearby edge-on disc galaxies. This phenomenon, where the disc mid-plane bends to become wavy, just as in warps but on a smaller scale, is seen here for the first time, in the dust lanes running across the discs. Because they are seen in absorption, this feature must be present in the dust disc in the outskirts of these galaxies. We enhance the visibility of these features using unsharp masking, trace the dust mid-plane across the disc, measure the corrugation amplitude by eye and the corrugation wavelength using Fourier analysis. The corrugation amplitude is found to be in the range of 70–300pc and the wavelengths lie between 1 and 5 kpc. In this limited sample, we find that the amplitude of the corrugations tends to be larger for lower mass galaxies, whereas the wavelength of corrugation does not seem to depend on the mass of host galaxies. Linear stability analysis is performed to find out the dynamical state of these dust discs. Based on WKB analysis, we find that the dust corrugations in about half of our sample are stable. Further analysis, on a larger sample would be useful to strengthen the above results.

Narayan, Chaitra A.↗
DOE OSTI2020

X-ray spectral and flux variability of the microquasar GRS 1758-258 on timescales from weeks to years

We present the spectral and timing evolution of the persistent black hole X-ray binary GRS 1758-258 based on almost 12 years of observations using the Rossi X-ray Timing Explorer Proportional Counter Array. While the source was predominantly found in the hard state during this time, it entered the thermally dominated soft state seven times. In the soft state GRS 1758-258 shows a strong decline in flux above 3 keV rather than the pivoting flux around 10 keV more commonly shown by black hole transients. In its 3–20 keV hardness intensity diagram, GRS 1758-258 shows a hysteresis of hard and soft state fluxes typical for transient sources in outburst. The RXTE-PCA and RXTE-ASM long-term light curves do not show any orbital modulations in the range of 2–30 d. However, in the dynamic power spectra significant peaks drift between 18.47 and 18.04 d for the PCA data, while less significant signatures between 19 d and 20 d are seen for the ASM data as well as for the Swift/BAT data. We discuss different models for the hysteresis behavior during state transitions as well as possibilities for the origin of the long term variation in the context of a warped accretion disk.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

X-Ray Spectral and Flux Variability of the Microquasar GRS 1758-258 on Timescales from Weeks to Years

We present the spectral and timing evolution of the persistent black hole X-ray binary GRS 1758-258 based on almost 12 years of observations using the Rossi X-ray Timing Explorer Proportional Counter Array. While the source was predominantly found in the hard state during this time, it entered the thermally dominated soft state seven times. In the soft state GRS 1758-258 shows a strong decline in flux above 3 keV rather than the pivoting flux around 10 keV more commonly shown by black hole transients. In its 3-20 keV hardness intensity diagram, GRS 1758-258 shows a hysteresis of hard and soft state fluxes typical for transient sources in outburst. The RXTE-PCA and RXTE-ASM long-term light curves do not show any orbital modulations in the range of 2 to 30 d. However, in the dynamic power spectra significant peaks drift between 18.47d and 18.04d for the PCA data, while less significant signatures between 19d and 20d are seen for the ASM data as well as for the Swift/BAT data. We discuss different models for the hysteresis behavior during state transitions as well as possibilities for the origin of the long term variation in the context of a warped accretion disk.

Maria Hirsch↗
NASA NTRS1995

Radial alignment simulation of solar wind streams observed by Pioneers 10 and 11 in 1974

A particularly favorable radial lineup between spacecraft in Earth orbit and Pioneers 10 and 11 (near the ecliptic plane at approximately 5.5 AU and approximately 4.5 AU, respectively) occurred in mid-1974, when the solar corona was in a stable and well-defined warped-dipole configuration. The radial alignment study reported here differs from previous applications of the technique in two respects: (1) It is the first time a two-dimensional (2-D) MHD model has been tested over such a lengthy propagation interval; the 2-D capability is crucial for treatment of the nonradial shearing motions occurring across the stream interface. (2) The three-dimensional (3-D) structure observed in the white light corona at that time is related to the systematic patterns of nonradial flow deflections appearing at the Pioneer corotating interaction region (CIR) fronts. Comparison of predicted and observed flows for pairs of streams in two successive rotations reveals that when the parent coronal hole projects far across the solar equator in a predominantly north-south orientation (i.e., nearest the 2-D idealization), the mapping is accurate down to details of the flow structures. But where the spacecraft tracks along a latitudinal boundary of a hole or the associated stream front is inclined at a shallow angle to the equator, the numerical projections deviate systematically from the observations. Among the sources of error are 3-D dynamical interactions neglected in the model, differential rotation effects, and slow temporal evolution of the coronal structures. A characteristic pattern of north-south and east-west deflections indicative of the 3-D geometry appears across the CIRs, but these patterns evidently reflect primarily the local, as opposed to global, orientation of the 3-D stream fronts. Such patterns appear common in CIRs observed by Pioneer during this period. These findings thus hold special relevance for the analysis of Ulysses observations, since the present coronal configuration is similar to that of 1974.

Pizzo, V. J.↗
DOE OSTI2025

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

The Long-term Light Curves of X-ray Binaries Contain Simultaneous Periodic and Random Components

LMC X-3 and Cyg X-2 show large amplitude X-ray fluctuations that have been attributed to a warped accretion disk. Cyg X-3 displays high amplitude, apparently non-periodic oscillations. We reanalyze these systems using RXTE ASM data and time-frequency decomposition techniques. We find that the long-term variations in Cyg X-2 can be completely characterized by excursions whose durations are integer multiples of the orbital period, including one essentially identical to the reported "period" of 78 days. Cyg X-3 can be characterized in terms of integer multiples of a 71-day fundamental period unrelated to the 4.8 day orbital period, but suggestively close to the approximately equal to greater than 60 day reported precession period of the relativistic jet inferred from recent radio observations. The long-term excursions of LMC X-3 are related to each other by rational fractions, suggesting the characteristic time scale is 10.594 days, shorter than any observed excursion to date. We explore the phase space evolution of the light curves using a natural embedding and find that all three systems possess two rotation centers that organize the phase space trajectories, one of low luminosity and the other of high luminosity. The implications of this repeatable behavior on generic models of accretion disk dynamics and mass transfer variability are explored.

White, Nicholas E.↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗