Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpolation models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗
NASA NTRS2002

Determining Greenland Ice Sheet Accumulation Rates from Radar Remote Sensing

An important component of NASA's Program for Arctic Regional Climate Assessment (PARCA) is a mass balance investigation of the Greenland Ice Sheet. The mass balance is calculated by taking the difference between the areally Integrated snow accumulation and the net ice discharge of the ice sheet. Uncertainties in this calculation Include the snow accumulation rate, which has traditionally been determined by interpolating data from ice core samples taken from isolated spots across the ice sheet. The sparse data associated with ice cores juxtaposed against the high spatial and temporal resolution provided by remote sensing , has motivated scientists to investigate relationships between accumulation rate and microwave observations as an option for obtaining spatially contiguous estimates. The objective of this PARCA continuation proposal was to complete an estimate of surface accumulation rate on the Greenland Ice Sheet derived from C-band radar backscatter data compiled in the ERS-1 SAR mosaic of data acquired during, September-November, 1992. An empirical equation, based on elevation and latitude, is used to determine the mean annual temperature. We examine the influence of accumulation rate, and mean annual temperature on C-band radar backscatter using a forward model, which incorporates snow metamorphosis and radar backscatter components. Our model is run over a range of accumulation and temperature conditions. Based on the model results, we generate a look-up table, which uniquely maps the measured radar backscatter, and mean annual temperature to accumulation rate. Our results compare favorably with in situ accumulation rate measurements falling within our study area.

Jezek, Kenneth C.↗
DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal photon sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal electron sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2018

Assessing the Altimetric Measurement from CYGNSS Data

The Cyclone Global Navigation Satellite System (CYGNSS) mission was designed to study hurricane intensification by measuring wind speeds in tropical cyclones. However, the delay-Doppler maps (DDM) that are produced can be used to estimate the sea surface height (SSH) at the specular reflection point on the ocean surface. Proofof-concept studies that DDMs are suitable to solve for SSH have been recently reported (Clarizia et al., 2016; Mashburn et al., 2018), based on data acquired by the demonstration satellite experiment Tech Demo Sat - 1 (TDS-1) carrying a GNSS-R receiver similar to the ones onboard CYGNSS. Although the precision of each 1sec averaged SSH is considerably lower than that of the existing satellite altimeters, by virtue of the dense coverage and frequent revisit time exhibited by the constellation of 8 microsats, the error may be smoothed down considerably by optimal interpolation (Li et al., 2016). Hence the CYGNSS dataset presents a potential opportunity to sample the tropical oceans, and investigate the sensitivity of the SSH measurements to mesoscale eddies. Our objective is to analyze the CYGNSS data and test a suite of retrieval algorithms, including that of Mashburn et al., 2018 to obtain SSH, to understand the error sources and possible corrections when available. Because of the limited antenna gain of the CYGNSS instrument, the power SNR for general wind conditions is expected to be too low to accurately track the reflection point contribution, resulting in altimetry errors of the order of several meters. Other significant error sources are the CYGNSS satellites orbit knowledge, the ionosphere, the troposphere, the mean sea surface and the tides; some of these errors can be reduced by using models. The corrected SSH data will be assimilated into the ROMS high-resolution ocean model to investigate their ability to reproduce mesoscale eddies in the tropical oceans.

Voronovich, Alexander↗
NASA NTRSDate not supplied

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗
NASA NTRS1999

Estimated Benefits of Variable-Geometry Wing Camber Control for Transport Aircraft

Analytical benefits of variable-camber capability on subsonic transport aircraft are explored. Using aerodynamic performance models, including drag as a function of deflection angle for control surfaces of interest, optimal performance benefits of variable camber are calculated. Results demonstrate that if all wing trailing-edge surfaces are available for optimization, drag can be significantly reduced at most points within the flight envelope. The optimization approach developed and illustrated for flight uses variable camber for optimization of aerodynamic efficiency (maximizing the lift-to-drag ratio). Most transport aircraft have significant latent capability in this area. Wing camber control that can affect performance optimization for transport aircraft includes symmetric use of ailerons and flaps. In this paper, drag characteristics for aileron and flap deflections are computed based on analytical and wind-tunnel data. All calculations based on predictions for the subject aircraft and the optimal surface deflection are obtained by simple interpolation for given conditions. An algorithm is also presented for computation of optimal surface deflection for given conditions. Benefits of variable camber for a transport configuration using a simple trailing-edge control surface system can approach more than 10 percent, especially for nonstandard flight conditions. In the cruise regime, the benefit is 1-3 percent.

Bolonkin, Alexander↗
NASA NTRS2012

Development of a Robust and Efficient Parallel Solver for Unsteady Turbomachinery Flows

The traditional design and analysis practice for advanced propulsion systems relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. In the design of advanced propulsion systems, CFD plays a major role in defining the required performance over the entire flight regime, as well as in testing the sensitivity of the design to the different modes of operation. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment. The turbomachinery simulation capability presented here is being developed in a computational tool called Loci-STREAM [1]. It integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci [2] which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective is to be able to routinely simulate problems involving complex geometries requiring large unstructured grids and complex multidisciplinary physics. An immediate application of interest is simulation of unsteady flows in rocket turbopumps, particularly in cryogenic liquid rocket engines. The key components of the overall methodology presented in this paper are the following: (a) high fidelity unsteady simulation capability based on Detached Eddy Simulation (DES) in conjunction with second-order temporal discretization, (b) compliance with Geometric Conservation Law (GCL) in order to maintain conservative property on moving meshes for second-order time-stepping scheme, (c) a novel cloud-of-points interpolation method (based on a fast parallel kd-tree search algorithm) for interfaces between turbomachinery components in relative motion which is demonstrated to be highly scalable, and (d) demonstrated accuracy and parallel scalability on large grids (approx 250 million cells) in full turbomachinery geometries.

West, Jeff↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
DOE OSTI2020

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.

Non-NASA Center↗
NASA NTRS2015

Seasonal Variability of Saturn's Tropospheric Temperatures, Winds and Para-H2 from Cassini Far-IR Spectroscopy

Far-IR 16-1000 micrometer spectra of Saturn's hydrogen-helium continuum measured by Cassini's Composite Infrared Spectrometer (CIRS) are inverted to construct a near-continuous record of upper tropospheric (70-700 mbar) temperatures and para-H2 fraction as a function of latitude, pressure and time for a third of a saturnian year (2004-2014, from northern winter to northern spring). The thermal field reveals evidence of reversing summertime asymmetries superimposed onto the belt/zone structure. The temperature structure is almost symmetric about the equator by 2014, with seasonal lag times that increase with depth and are qualitatively consistent with radiative climate models. Localised heating of the tropospheric hazes (100-250 mbar) create a distinct perturbation to the temperature profile that shifts in magnitude and location, declining in the autumn hemisphere and growing in the spring. Changes in the para-H2 (f(sub p)) distribution are subtle, with a 0.02-0.03 rise over the spring hemisphere (200-500 mbar) perturbed by (i) low-f(sub p) air advected by both the springtime storm of 2010 and equatorial upwelling; and (ii) subsidence of high-f(sub p) air at northern high latitudes, responsible for a developing north-south asymmetry in f(sub p). Conversely, the shifting asymmetry in the para-H2 disequilibrium primarily reflects the changing temperature structure (and hence the equilibrium distribution of f(sub p)), rather than actual changes in f(sub p) induced by chemical conversion or transport. CIRS results interpolated to the same point in the seasonal cycle as re-analysed Voyager-1 observations (early northern spring) show qualitative consistency from year to year (i.e., the same tropospheric asymmetries in temperature and f(sub p)), with the exception of the tropical tropopause near the equatorial zones and belts, where downward propagation of a cool temperature anomaly associated with Saturn's stratospheric oscillation could potentially perturb tropopause temperatures, para-H2 and winds. Quantitative differences between the Cassini and Voyager epochs suggest that the oscillation is not in phase with the seasonal cycle at these tropospheric depths (i.e., it should be described as quasi-periodic rather than 'semi annual'). Variability in the zonal wind field derived from latitudinal thermal gradients is small (less than 10 m/s per scale height near the tropopause) and mostly affects the broad retrograde jets, with the notable exception of large variability on the northern flank of the equatorial jet. The meridional potential vorticity (PV) gradient, and hence the 'staircase of PV' associated with spatial variations in the vigour of vertical mixing, has varied over the course of the mission but maintained its overall shape. PV gradients in latitude and altitude are used to estimate the atmospheric refractive index for the propagation of stationary planetary (Rossby) waves, predicting that such wave activity would be confined to regions of real refractivity (tropical regions plus bands at 35-45 in both hemispheres). The penetration depth of these regions into the upper troposphere is temporally variable (potentially associated with stratification changes), whereas the latitudinal structure is largely unchanged over time (associated with the zonal jet system).

Fletcher, Leigh N.↗
DOE OSTI2021

DSGAN

This study develops a neural network-based approach for emulating high-resolution modeled precipitation data with comparable statistical properties but at greatly reduced computational cost. The key idea is to use combination of low-and high- resolution simulations to train a neural network to map from the former to the latter. Specifically, we define two types of CNNs, one that stacks variables directly and one that encodes each variable before stacking, and we train each CNN type both with a conventional loss function, such as mean square error (MSE), and with a conditional generative adversarial network (CGAN), for a total of four CNN variants. We compare the four new CNN-derived high-resolution precipitation results with precipitation generated from original high resolution simulations, a bilinear interpolater and the state-of-the-art CNN-based super-resolution (SR) technique. Results show that the SR technique produces results similar to those of the bilinear interpolator with smoother spatial and temporal distributions and smaller data variabilities and extremes than the original high resolution simulations. While the new CNNs trained by MSE generate better results over some regions than the interpolator and SR technique do, their predictions are still not as close as the original high resolution simulations. The CNNs trained by CGAN generate more realistic and physically reasonable results, better capturing not only data variability in time and space but also extremes such as intense and long-lasting storms. The new proposed CNN-based downscaling approach can downscale precipitation from 50~km to 12~km in 14~min for 30~years

LIU, ZHENGCHUN↗
NASA NTRS1987

Baseline estimation from simultaneous satellite laser tracking

Simultaneous Range Differences (SRDs) to Lageos are obtained by dividing the observing stations into pairs with quasi-simultaneous observations. For each of those pairs the station with the least number of observations is identified, and at its observing epochs interpolated ranges for the alternate station are generated. The SRD observables are obtained by subtracting the actually observed laser range of the station having the least number of observations from the interpolated ranges of the alternate station. On the basis of these observables semidynamic single baseline solutions were performed. The aim of these solutions is to further develop and implement the SRD method in the real data environment, to assess its accuracy, its advantages and disadvantages as related to the range dynamic mode methods, when the baselines are the only parameters of interest. Baselines, using simultaneous laser range observations to Lageos, were also estimated through the purely geometric method. These baselines formed the standards the standards of comparison in the accuracy assessment of the SRD method when compared to that of the range dynamic mode methods. On the basis of this comparison it was concluded that for baselines of regional extent the SRD method is very effective, efficient, and at least as accurate as the range dynamic mode methods, and that on the basis of a simple orbital modeling and a limited orbit adjustment. The SRD method is insensitive to the inconsistencies affecting the terrestrial reference frame and simultaneous adjustment of the Earth Rotation Parameters (ERPs) is not necessary.

Dedes, George C.↗
NASA NTRS2003

Characterization of the Crab Pulsar's Timing Noise

We present a power spectral analysis of the Crab pulsar's timing noise, mainly using radio measurements from Jodrell Bank taken over the period 1982-1989, an interval bounded by sparse data sampling and a large glitch. The power spectral analysis is complicated by nonuniform data sampling and the presence of a steep red power spectrum that can distort power spectra measurement by causing severe power 'leakage'. We develop a simple windowing method for computing red noise power spectra of uniformly sampled data sets and test it on Monte Carlo generated sample realizations of red power-law noise. We generalize time-domain methods of generating power-law red noise with even integer spectral indices to the case of noninteger spectral indices. The Jodrell Bank pulse phase residuals are dense and smooth enough that an interpolation onto a uniform time series is possible. A windowed power spectrum is computed revealing a periodic or nearly periodic component with a period of 568 +/- 10 days and a l/f(exp 3) power-law noise component in pulse phase with a noise strength S(sub infinity)=(1.24 +/- 0.067) x 10(exp 16) cycles(exp 2)/sec(exp 2) over the analysis frequency range f=0.003- 0.1 cycles/day. This result deviates from past analyses which characterized the pulse phase timing residuals as either l/f(sub 4) power-law noise or a quasiperiodic process. The analysis was checked using the Deeter polynomial method of power spectrum estimation that was developed for the case of nonuniform sampling, but has lower spectral resolution. The timing noise is consistent with a torque noise spectrum rising with analysis frequency as f implying blue torque noise, a result not predicted by current models of pulsar timing noise. If the periodic or nearly periodic component is due to a binary companion, we find a mass function f(M) = (6.8 +/- 2.4) x 10(exp -16) solar mass and a companion mass, M(sub c) is greater than or equal to 3.2 solar mass assuming a Crab pulsar mass of 1.4 solar mass.

Scott, D. M.↗
DOE OSTI2025

Thermal bootstrap of matrix quantum mechanics

We implement a bootstrap method that combines stationary state conditions, thermal inequalities, and semidefinite relaxations of matrix logarithm in the ungauged one-matrix quantum mechanics, at finite rank N as well as in the large N limit, and determine finite temperature observables that interpolate between available analytic results in the low and high temperature limits respectively. We also obtain bootstrap bounds on thermal phase transition as well as preliminary results in the ungauged two-matrix quantum mechanics.

1/N Expansion↗
DOE OSTI2023

Vehicle Powertrain Simulation Accuracy for Various Drive Cycle Frequencies and Upsampling Techniques

As connected and automated vehicle technologies emerge and proliferate, lower frequency vehicle trajectory data is becoming more widely available. In some cases, entire fleets are streaming position, speed, and telemetry at sample rates of less than 10 seconds. This presents opportunities to apply powertrain simulators such as the National Renewable Energy Laboratory's Future Automotive Systems Technology Simulator to model how advanced powertrain technologies would perform in the real world. However, connected vehicle data tends to be available at lower temporal frequencies than the 1-10 Hz trajectories that have typically been used for powertrain simulation. Higher frequency data, typically used for simulation, is costly to collect and store and therefore is often limited in density and geography. This paper explores the suitability of lower frequency, high availability, connected vehicle data for detailed powertrain simulation. A large data set of 1 Hz trajectories is used to quantify the accuracy loss when simulating energy consumption for conventional, hybrid, and battery electric powertrains using less than 1 Hz data. Techniques to upsample lower frequency drive cycle data in order to increase accuracy are also explored. Median energy consumption errors when simulating energy consumption for a 1/10 Hz trajectory are found to be 3-6% when compared to 1 Hz trajectories. Applying upsampling and interpolation techniques are shown to reduce the simulation errors by roughly 50%. The findings in this work can guide connected vehicle data collection specifications and processing techniques applied when using collected data for powertrain simulation.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Coupling laser physics to radiation-hydrodynamics

In order to accurately model implosion hydrodynamics in a radiation-hydrodynamics code, it is essential to include accurate accounting for energy deposition physics. In inertial confinement fusion (ICF), where capsules are driven by lasers or laser-driven x-rays, energy deposition profiles and energy transport have a strong impact on the development and evolution of capsule dynamics and hydrodynamic instabilities. Nevertheless, accurately modeling laser beam propagation in radiation-hydrodynamics codes presents unique challenges associated with disparate resolution requirements, the potential to seed spurious noise in highly unstable systems, and computational expense. We discuss a new method for coupling laser ray-tracing physics to a radiation hydrodynamics code, developed in the process of implementing the Mazinisin laser ray-trace into the xRAGE radiation hydrodynamics code. In contrast to previous approaches, in which laser ray-tracing is performed on the radiation-hydrodynamics mesh, our method involves a mesh generation and evolution strategy that addresses the unique requirements of the laser ray-trace in a separate mesh, enabling performance enhancements and strategies to reduce noise seeded by the discretization of beams into computational rays. In addition, we have employed several methods to ensure that spurious mesh imprinting is minimized. These involved optimizing the laser and radiation-hydrodynamics meshes as well as interpolation between them and requires the use of an exact initialization method for the radiation-hydrodynamics mesh. Furthermore, these techniques have enabled efficient computation of laser-driven implosions and other experiments with minimal introduction of spurious noise.

47 OTHER INSTRUMENTATION↗