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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

ACCURATE CHARACTERIZATION OF HIGH-DEGREE MODES USING MDI OBSERVATIONS

We present the first accurate characterization of high-degree modes, derived using the best Michelson Doppler Imager (MDI) full-disk full-resolution data set available. A 90 day long time series of full-disk 2 arcsec pixel{sup -1} resolution Dopplergrams was acquired in 2001, thanks to the high rate telemetry provided by the Deep Space Network. These Dopplergrams were spatially decomposed using our best estimate of the image scale and the known components of MDI's image distortion. A multi-taper power spectrum estimator was used to generate power spectra for all degrees and all azimuthal orders, up to l = 1000. We used a large number of tapers to reduce the realization noise, since at high degrees the individual modes blend into ridges and thus there is no reason to preserve a high spectral resolution. These power spectra were fitted for all degrees and all azimuthal orders, between l = 100 and l = 1000, and for all the orders with substantial amplitude. This fitting generated in excess of 5.2 Multiplication-Sign 10{sup 6} individual estimates of ridge frequencies, line widths, amplitudes, and asymmetries (singlets), corresponding to some 5700 multiplets (l, n). Fitting at high degrees generates ridge characteristics, characteristics that do not correspond to the underlying mode characteristics. We used a sophisticated forward modeling to recover the best possible estimate of the underlying mode characteristics (mode frequencies, as well as line widths, amplitudes, and asymmetries). We describe in detail this modeling and its validation. The modeling has been extensively reviewed and refined, by including an iterative process to improve its input parameters to better match the observations. Also, the contribution of the leakage matrix on the accuracy of the procedure has been carefully assessed. We present the derived set of corrected mode characteristics, which includes not only frequencies, but line widths, asymmetries, and amplitudes. We present and discuss their uncertainties and the precision of the ridge-to-mode correction schemes, through a detailed assessment of the sensitivity of the model to its input set. The precision of the ridge-to-mode correction is indicative of any possible residual systematic biases in the inferred mode characteristics. In our conclusions, we address how to further improve these estimates, and the implications for other data sets, like GONG+ and HMI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN EMPIRICAL FORMULA FOR THE DISTRIBUTION FUNCTION OF A THIN EXPONENTIAL DISC

An empirical formula for a Shu distribution function that reproduces a thin disc with exponential surface density to good accuracy is presented. The formula has two free parameters that specify the functional form of the velocity dispersion. Conventionally, this requires the use of an iterative algorithm to produce the correct solution, which is computationally taxing for applications like Markov Chain Monte Carlo model fitting. The formula has been shown to work for flat, rising, and falling rotation curves. Application of this methodology to one of the Dehnen distribution functions is also shown. Finally, an extension of this formula to reproduce velocity dispersion profiles that are an exponential function of radius is also presented. Our empirical formula should greatly aid the efficient comparison of disc models with large stellar surveys or N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Unification of binary star ephemeris solutions

Time-related binary system characteristics such as orbital period, its rate of change, apsidal motion, and variable light-time delay due to a third body, are measured in two ways that can be mutually complementary. The older way is via eclipse timings, while ephemerides by simultaneous whole light and velocity curve analysis have appeared recently. Each has its advantages, for example, eclipse timings typically cover relatively long time spans while whole curves often have densely packed data within specific intervals and allow access to systemic properties that carry additional timing information. Synthesis of the two information sources can be realized in a one step process that combines several data types, with automated weighting based on their standard deviations. Simultaneous light-velocity-timing solutions treat parameters of apsidal motion and the light-time effect coherently with those of period and period change, allow the phenomena to interact iteratively, and produce parameter standard errors based on the quantity and precision of the curves and timings. The logic and mathematics of the unification algorithm are given, including computation of theoretical conjunction times as needed for generation of eclipse timing residuals. Automated determination of eclipse type, recovery from inaccurate starting ephemerides, and automated data weighting are also covered. Computational examples are given for three timing-related cases—steady period change (XY Bootis), apsidal motion (V526 Sagittarii), and the light-time effect due to a binary's reflex motion in a triple system (AR Aurigae). Solutions for all combinations of radial velocity, light curve, and eclipse timing input show consistent results, with a few minor exceptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

redMaPPer. I. Algorithm and SDSS DR8 catalog

We describe redMaPPer, a new red sequence cluster finder specifically designed to make optimal use of ongoing and near-future large photometric surveys. The algorithm has multiple attractive features: (1) it can iteratively self-train the red sequence model based on a minimal spectroscopic training sample, an important feature for high-redshift surveys. (2) It can handle complex masks with varying depth. (3) It produces cluster-appropriate random points to enable large-scale structure studies. (4) All clusters are assigned a full redshift probability distribution P(z). (5) Similarly, clusters can have multiple candidate central galaxies, each with corresponding centering probabilities. (6) The algorithm is parallel and numerically efficient: it can run a Dark Energy Survey-like catalog in ∼500 CPU hours. (7) The algorithm exhibits excellent photometric redshift performance, the richness estimates are tightly correlated with external mass proxies, and the completeness and purity of the corresponding catalogs are superb. We apply the redMaPPer algorithm to ∼10, 000 deg{sup 2} of SDSS DR8 data and present the resulting catalog of ∼25,000 clusters over the redshift range z in [0.08, 0.55]. The redMaPPer photometric redshifts are nearly Gaussian, with a scatter σ {sub z} ≈ 0.006 at z ≈ 0.1, increasing to σ {sub z} ≈ 0.02 at z ≈ 0.5 due to increased photometric noise near the survey limit. The median value for |Δz|/(1 + z) for the full sample is 0.006. The incidence of projection effects is low (≤5%). Detailed performance comparisons of the redMaPPer DR8 cluster catalog to X-ray and Sunyaev-Zel'dovich catalogs are presented in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Full waveform inversion of solar interior flows

The inference of flows of material in the interior of the Sun is a subject of major interest in helioseismology. Here, we apply techniques of full waveform inversion (FWI) to synthetic data to test flow inversions. In this idealized setup, we do not model seismic realization noise, training the focus entirely on the problem of whether a chosen supergranulation flow model can be seismically recovered. We define the misfit functional as a sum of L {sub 2} norm deviations in travel times between prediction and observation, as measured using short-distance filtered f and p {sub 1} and large-distance unfiltered p modes. FWI allows for the introduction of measurements of choice and iteratively improving the background model, while monitoring the evolution of the misfit in all desired categories. Although the misfit is seen to uniformly reduce in all categories, convergence to the true model is very slow, possibly because it is trapped in a local minimum. The primary source of error is inaccurate depth localization, which, due to density stratification, leads to wrong ratios of horizontal and vertical flow velocities ({sup c}ross talk{sup )}. In the present formulation, the lack of sufficient temporal frequency and spatial resolution makes it difficult to accurately localize flow profiles at depth. We therefore suggest that the most efficient way to discover the global minimum is to perform a probabilistic forward search, involving calculating the misfit associated with a broad range of models (generated, for instance, by a Monte Carlo algorithm) and locating the deepest minimum. Such techniques possess the added advantage of being able to quantify model uncertainty as well as realization noise (data uncertainty).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A STATISTICAL RECONSTRUCTION OF THE PLANET POPULATION AROUND KEPLER SOLAR-TYPE STARS

Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1-4 R {sub ⊕}) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type (0.8 < R {sub *}/R {sub ☉} < 1.2) stars with 430 planets on 20-200 day orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an ''iterative simulation'' technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is 0.46 ± 0.03. The occurrence is ∼0.66 if planets interior to 20 days are included. We estimate the occurrence of Earth-size planets in the ''habitable zone'' (defined as 1-2 R {sub ⊕}, 0.99-1.7 AU for solar-twin stars) as 6.4{sub −1.1}{sup +3.4}%. Our results largely agree with those of Petigura et al., although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. stellar parameters, and also the exclusion of planets that may have been affected by proximity to the host star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THREE-DIMENSIONAL RADIO AND X-RAY MODELING AND DATA ANALYSIS SOFTWARE: REVEALING FLARE COMPLEXITY

Many problems in solar physics require analysis of imaging data obtained in multiple wavelength domains with differing spatial resolution in a framework supplied by advanced three-dimensional (3D) physical models. To facilitate this goal, we have undertaken a major enhancement of our IDL-based simulation tools developed earlier for modeling microwave and X-ray emission. The enhanced software architecture allows the user to (1) import photospheric magnetic field maps and perform magnetic field extrapolations to generate 3D magnetic field models; (2) investigate the magnetic topology by interactively creating field lines and associated flux tubes; (3) populate the flux tubes with user-defined nonuniform thermal plasma and anisotropic, nonuniform, nonthermal electron distributions; (4) investigate the spatial and spectral properties of radio and X-ray emission calculated from the model; and (5) compare the model-derived images and spectra with observational data. The package integrates shared-object libraries containing fast gyrosynchrotron emission codes, IDL-based soft and hard X-ray codes, and potential and linear force-free field extrapolation routines. The package accepts user-defined radiation and magnetic field extrapolation plug-ins. We use this tool to analyze a relatively simple single-loop flare and use the model to constrain the magnetic 3D structure and spatial distribution of the fast electrons inside this loop. We iteratively compute multi-frequency microwave and multi-energy X-ray images from realistic magnetic flux tubes obtained from pre-flare extrapolations, and compare them with imaging data obtained by SDO, NoRH, and RHESSI. We use this event to illustrate the tool's use for the general interpretation of solar flares to address disparate problems in solar physics.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗