Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

STELLAR PARAMETERS FOR HD 69830, A NEARBY STAR WITH THREE NEPTUNE MASS PLANETS AND AN ASTEROID BELT

We used the CHARA Array to directly measure the angular diameter of HD 69830, home to three Neptune mass planets and an asteroid belt. Our measurement of 0.674 ± 0.014 mas for the limb-darkened angular diameter of this star leads to a physical radius of R {sub *} = 0.9058 ± 0.0190 R {sub ☉} and luminosity of L {sub *} = 0.622 ± 0.014 L {sub ☉} when combined with a fit to the spectral energy distribution of the star. Placing these observed values on an Hertzsprung-Russel diagram along with stellar evolution isochrones produces an age of 10.6 ± 4 Gyr and mass of 0.863 ± 0.043 M {sub ☉}. We use archival optical echelle spectra of HD 69830 along with an iterative spectral fitting technique to measure the iron abundance ([Fe/H] = –0.04 ± 0.03), effective temperature (5385 ± 44 K), and surface gravity (log g = 4.49 ± 0.06). We use these new values for the temperature and luminosity to calculate a more precise age of 7.5 ± 3 Gyr. Applying the values of stellar luminosity and radius to recent models on the optimistic location of the habitable zone produces a range of 0.61-1.44 AU; partially outside the orbit of the furthest known planet (d) around HD 69830. Finally, we estimate the snow line at a distance of 1.95 ± 0.19 AU, which is outside the orbit of all three planets and its asteroid belt.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Accurate gravities of F, G, and K stars from high resolution spectra without external constraints

We demonstrate a new procedure to derive accurate and precise surface gravities from high resolution spectra without the use of external constraints. Our analysis utilizes Spectroscopy Made Easy with robust spectral line constraints and uses an iterative process to mitigate degeneracies in the fitting process. We adopt an updated radiative transfer code, a new treatment for neutral perturber broadening, a line list with multiple gravity constraints and separate fitting for global stellar properties and abundance determinations. To investigate the sources of temperature dependent trends in determining log g noted in previous studies, we obtained Keck HIRES spectra of 42 Kepler asteroseismic stars. In comparison to asteroseismically determined log g our spectroscopic analysis has a constant offset of 0.01 dex with a rms scatter of 0.05 dex. We also analyzed 30 spectra which had published surface gravities determined using the a/R{sub ∗} technique from planetary transits and found a constant offset of 0.06 dex and rms scatter of 0.07 dex. The two samples covered effective temperatures between 5000 and 6700 K with log g between 3.7 and 4.6.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SPECTRAL CONFUSION FOR COSMOLOGICAL SURVEYS OF REDSHIFTED C II EMISSION

Far-infrared cooling lines are ubiquitous features in the spectra of star-forming galaxies. Surveys of redshifted fine-structure lines provide a promising new tool to study structure formation and galactic evolution at redshifts including the epoch of reionization as well as the peak of star formation. Unlike neutral hydrogen surveys, where the 21 cm line is the only bright line, surveys of redshifted fine-structure lines suffer from confusion generated by line broadening, spectral overlap of different lines, and the crowding of sources with redshift. We use simulations to investigate the resulting spectral confusion and derive observing parameters to minimize these effects in pencil-beam surveys of redshifted far-IR line emission. We generate simulated spectra of the 17 brightest far-IR lines in galaxies, covering the 150–1300 μm wavelength region corresponding to redshifts 0 < z < 7, and develop a simple iterative algorithm that successfully identifies the 158 μm [C ii] line and other lines. Although the [C ii] line is a principal coolant for the interstellar medium, the assumption that the brightest observed lines in a given line of sight are always [C ii] lines is a poor approximation to the simulated spectra once other lines are included. Blind line identification requires detection of fainter companion lines from the same host galaxies, driving survey sensitivity requirements. The observations require moderate spectral resolution 700 < R < 4000 with angular resolution between 20″ and 10′, sufficiently narrow to minimize confusion yet sufficiently large to include a statistically meaningful number of sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SHOCK WAVE STRUCTURE IN ASTROPHYSICAL FLOWS WITH AN ACCOUNT OF PHOTON TRANSFER

For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiative shock is calculated for a number of shock tube problems, including strong shock waves, and relativistic- and radiation-dominated cases. Calculations are performed using an iterative technique that consistently solves the equations of relativistic hydrodynamics and relativistic comoving radiative transfer. A comparison of radiative transfer solutions with the Eddington approximation and the M1 closure is made. A qualitative analysis of moment equations for radiation is performed and the conditions for the existence of jump discontinuity for non-relativistic cases are investigated numerically.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The dual-electron cyclotron emission based measurement of 3D structures on DIII-D tokamak

The measurements of the locked mode and other 3D magnetohydrodynamic modes are essential for real-time control of plasma stability in future fusion reactors. Current magnetic measurements present non-negligible challenges for the application on fusion pilot plants (FPPs) due to neutron radiation and pose challenging design issues. In this paper, we discuss an alternative solution with microwave-based measurements. The dual-electron cyclotron emission (ECE) method has been developed and applied on DIII-D. The dual-ECE measurement method demonstrates clear identification of magnetic islands compared to single ECE, which is highly beneficial for locked mode feedback control on the FPPs. Toroidally separated ECE measurements provide accurate three-dimensional structure of magnetic islands and tearing modes. This technique also offers deeper physical insights into the penetration of resonant magnetic perturbations. Numerical calculations indicate that small locked modes (magnetic island width > 2 cm) can be detected using the dual-ECE measurement method in ITER plasma. This will enhance measurement capabilities for the safe operation of FPP plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effect of noncircularity of experimental beam on CMB parameter estimation

Measurement of Cosmic Microwave Background (CMB) anisotropies has been playing a lead role in precision cosmology by providing some of the tightest constrains on cosmological models and parameters. However, precision can only be meaningful when all major systematic effects are taken into account. Non-circular beams in CMB experiments can cause large systematic deviation in the angular power spectrum, not only by modifying the measurement at a given multipole, but also introducing coupling between different multipoles through a deterministic bias matrix. Here we add a mechanism for emulating the effect of a full bias matrix to the PLANCK likelihood code through the parameter estimation code SCoPE. We show that if the angular power spectrum was measured with a non-circular beam, the assumption of circular Gaussian beam or considering only the diagonal part of the bias matrix can lead to huge error in parameter estimation. We demonstrate that, at least for elliptical Gaussian beams, use of scalar beam window functions obtained via Monte Carlo simulations starting from a fiducial spectrum, as implemented in PLANCK analyses for example, leads to only few percent of sigma deviation of the best-fit parameters. However, we notice more significant differences in the posterior distributions for some of the parameters, which would in turn lead to incorrect errorbars. These differences can be reduced, so that the errorbars match within few percent, by adding an iterative reanalysis step, where the beam window function would be recomputed using the best-fit spectrum estimated in the first step.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A description of the Galactic Center excess in the Minimal Supersymmetric Standard Model

Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range 84–92 GeV or 87–97 GeV annihilating into W bosons. A third solutions is found for a neutralino of mass 174–187 GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density 0.06<Ωh{sup 2}<0.13. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the pMSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Robust forecasts on fundamental physics from the foreground-obscured, gravitationally-lensed CMB polarization

Recent results from the BICEP, Keck Array and Planck Collaborations demonstrate that Galactic foregrounds are an unavoidable obstacle in the search for evidence of inflationary gravitational waves in the cosmic microwave background (CMB) polarization. Beyond the foregrounds, the effect of lensing by intervening large-scale structure further obscures all but the strongest inflationary signals permitted by current data. With a plethora of ongoing and upcoming experiments aiming to measure these signatures, careful and self-consistent consideration of experiments' foreground- and lensing-removal capabilities is critical in obtaining credible forecasts of their performance. We investigate the capabilities of instruments such as Advanced ACTPol, BICEP3 and Keck Array, CLASS, EBEX10K, PIPER, Simons Array, SPT-3G and SPIDER, and projects as COrE+, LiteBIRD-ext, PIXIE and Stage IV, to clean contamination due to polarized synchrotron and dust from raw multi-frequency data, and remove lensing from the resulting co-added CMB maps (either using iterative CMB-only techniques or through cross-correlation with external data). Incorporating these effects, we present forecasts for the constraining power of these experiments in terms of inflationary physics, the neutrino sector, and dark energy parameters. Made publicly available through an online interface, this tool enables the next generation of CMB experiments to foreground-proof their designs, optimize their frequency coverage to maximize scientific output, and determine where cross-experimental collaboration would be most beneficial. We find that analyzing data from ground, balloon and space instruments in complementary combinations can significantly improve component separation performance, delensing, and cosmological constraints over individual datasets. In particular, we find that a combination of post-2020 ground- and space-based experiments could achieve constraints such as σ(r)∼1.3×10{sup −4}, σ(n{sub t})∼0.03, σ( n{sub s} )∼1.8×10{sup −3}, σ(α{sub s})∼1.7×10{sup −3}, σ( M{sub ν} )∼31 meV, σ( w )∼0.09, σ( w{sub 0} )∼ 0.25, 0σ( w{sub a} )∼ 0.5, σ( N{sub eff} )∼0.024 and σ( Ω{sub k} )∼1.5×10{sup −3}, after component separation and iterative delensing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reconstructing cosmic polarization rotation with ResUNet-CMB

Cosmic polarization rotation, which may result from parity-violating new physics or the presence of primordial magnetic fields, converts E-mode polarization of the cosmic microwave background (CMB) into B-mode polarization. Anisotropic cosmic polarization rotation leads to statistical anisotropy in CMB polarization and can be reconstructed with quadratic estimator techniques similar to those designed for gravitational lensing of the CMB. At the sensitivity of upcoming CMB surveys, lensing-induced B-mode polarization will act as a limiting factor in the search for anisotropic cosmic polarization rotation, meaning that an analysis which incorporates some form of delensing will be required to improve constraints on the effect with future surveys. In this study we extend the ResUNet-CMB convolutional neural network to reconstruct anisotropic cosmic polarization rotation in the presence of gravitational lensing and patchy reionization, and we show that the network simultaneously reconstructs all three effects with variance that is lower than that from the standard quadratic estimator nearly matching the performance of an iterative reconstruction method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Speckle evolution with multiple steps of least-squares phase removal

We study numerically the evolution of speckle fields due to the annihilation of optical vortices after the least-squares phase has been removed. A process with multiple steps of least-squares phase removal is carried out to minimize both vortex density and scintillation index. Statistical results show that almost all the optical vortices can be removed from a speckle field, which finally decays into a quasiplane wave after such an iterative process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Condensate fraction of cold gases in a nonuniform external potential

Exact calculation of the condensate fraction in multidimensional inhomogeneous interacting Bose systems in a confining potential of arbitrary shape is a difficult computational problem. We have developed an iterative procedure which allows us to calculate the condensate fraction as well as the corresponding eigenfunction of the one-body density matrix. We successfully validate this procedure in diffusion Monte Carlo simulations of a Bose gas in an optical lattice at zero temperature. We also discuss the relation between different criteria used for testing coherence in cold Bose systems, such as the fraction of particles that are superfluid, condensed, or in the zero-momentum state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Basic prerequisites for limit-cycle oscillations within a synchronously pumped passive optical nonlinear fiber-ring resonator

By taking into account Kerr nonlinearity and negative second-order dispersion, we show by numerical investigations that in a passive optical nonlinear fiber ring resonator synchronously pumped with femtosecond pulses a temporal delay via a resonator length detuning or via third-order dispersion leads to limit cycle behavior, which is represented by an oscillatory evolution of the temporal as well as the spectral pulse shape. By introducing a two-dimensional iterative map, we show that limit cycles arise due to the mutual coupling of different spectral pulse components. The resulting system behavior can be used to introduce additional adjustable frequency components into frequency combs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Adaptive schemes for incomplete quantum process tomography

We propose an iterative algorithm for incomplete quantum process tomography with the help of quantum state estimation. The algorithm, which is based on the combined principles of maximum likelihood and maximum entropy, yields a unique estimator for an unknown quantum process when one has less than a complete set of linearly independent measurement data to specify the quantum process uniquely. We apply this iterative algorithm adaptively in various situations and so optimize the amount of resources required to estimate a quantum process with incomplete data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Machine learning assisted unfolding for neutrino cross-section measurements with the OmniFold technique

The choice of unfolding method for a cross-section measurement is tightly coupled to the model dependence of the efficiency correction and the overall impact of cross-section modeling uncertainties in the analysis. A key issue is the dimensionality used in unfolding, as the kinematics of all outgoing particles in an event typically affect the reconstruction performance in a neutrino detector. OmniFold is an unfolding method that iteratively reweights a simulated dataset, using machine learning to utilize arbitrarily high-dimensional information, that has previously been applied to proton-proton and proton-electron datasets. This paper demonstrates OmniFold’s application to a neutrino cross-section measurement for the first time using a public T2K near detector simulated dataset, comparing its performance with traditional approaches using a mock data study.

Machine learning↗
DOE OSTI2020

Relevance Vector Machine Image Reconstruction Algorithm for Electrical Capacitance Tomography With Explicit Uncertainty Estimates

We present a Relevance Vector Machine (RVM) based algorithm for electrical capacitance tomography (ECT) applications that can concurrently provide image reconstruction results and uncertainty estimates about the reconstruction. To illustrate the RVM operation in ECT, we simulate typical ECT scenarios, making explicit the connection between the reconstructed pixel values and the corresponding uncertainty estimates in each case. We compare the RVM reconstruction performance with that of the Iterative Landweber Method (ILM) and the least absolute shrinkage and selection operator (LASSO) in all the considered scenarios. The results show that, in addition to the key advantage of providing uncertainty measures, RVM can achieve similar reconstruction results with either lower or similar computational complexity.

42 ENGINEERING↗
DOE OSTI2021

Spread Spectrum Symbol Detection With Blind Interference Suppression in FBMC-SS

Recent works have demonstrated Filter Bank Multicarrier Spread Spectrum (FBMC-SS) to be a robust communication scheme in the presence of high-power interferers. Existing FBMC-SS symbol detector designs based on analysis filter banks (AFB) suggest using an optimal combining scheme to suppress the interferers, necessitating some noise/interference power estimation method. In this paper, we introduce a symbol detector with blind interference suppression by extending a recently developed packet detection method. We then provide an analysis to show that the existing AFB-based symbol detector and the one proposed in this paper are equivalent in typical usage scenarios. A fully-fledged receiver design is proposed utilizing this symbol detector, with specifics presented for estimation of the channel impulse response and carrier frequency offset (CFO). We also outline a method of iterating upon the channel and CFO estimations to improve the quality of both parameters. Moreover, a modification to allow improved performance of the symbol detector at high SNR is provided. Finally, simulated performance results are presented to corroborate these findings and demonstrate the efficiency of this receiver design.

99 GENERAL AND MISCELLANEOUS↗