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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Dynamic phasing of multichannel cw laser radiation by means of a stochastic gradient algorithm

The phasing of a multichannel laser beam by means of an iterative stochastic parallel gradient (SPG) algorithm has been numerically and experimentally investigated. The operation of the SPG algorithm is simulated, the acceptable range of amplitudes of probe phase shifts is found, and the algorithm parameters at which the desired Strehl number can be obtained with a minimum number of iterations are determined. An experimental bench with phase modulators based on lithium niobate, which are controlled by a multichannel electronic unit with a real-time microcontroller, has been designed. Phasing of 16 cw laser beams at a system response bandwidth of 3.7 kHz and phase thermal distortions in a frequency band of about 10 Hz is experimentally demonstrated. The experimental data are in complete agreement with the calculation results. (control of laser radiation parameters)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Scheme of a hydrogen-molecule quantum simulator based on two ultracold rubidium atoms

A scheme is proposed for implementing a hydrogen-molecule quantum simulator based on two ultracold rubidium atoms trapped into spatially separated optical dipole traps. The scheme includes the adiabatic preparation of the initial quantum state of two atoms and the iterative quantum phase estimation. The accuracy of measuring the ground state energy of a molecule is numerically calculated as a function of the number of iterations. The simulation is performed using two-qubit gates based on the dipole blockade effect under short-term excitation of atoms into the Rydberg states with allowance for the finite lifetime of Rydberg states and the finite energies of the van der Waals interaction. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

An iterative approach to non-overdetermined inverse scattering at fixed energy

We present an iterative approximate reconstruction algorithm for non-overdetermined inverse scattering at fixed energy E with incomplete data, where the dimension d≥2. In particular, we obtain rapidly converging approximate reconstructions for this inverse scattering for E→+∞. Bibliography: 38 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Weighted inequalities for quasilinear integral operators on the semi-axis and applications to Lorentz spaces

A precise characterization of inequalities in weighted Lebesgue spaces with positive quasilinear integral operators of iterative type on the half-axis is given. All cases of positive integration parameters are treated, including the case of supremum. Applications to the solution of the well-known problem of the boundedness of the Hardy-Littlewood maximal operator in weighted Lorentz Γ-spaces are given. Bibliography: 41 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advanced Reactor Control and Operations (ARCO): A University Research Facility for Developing Optimized Digital Control Rooms

The Advanced Reactor Control and Operations (ARCO) facility was constructed in January 2018 to serve as a test bed for advanced reactor control rooms and operator support systems. Since then, it has supported human-machine interface user experience research, fault detection and mitigation technology development, control room concept of operations development, and remote operations research. ARCO serves as the control room for the Compact Integral Effects Test (CIET) facility, which replicates the primary-side flow paths and thermal-hydraulic behavior of a fluoride-salt-cooled high-temperature reactor (FHR) using simulant fluids and scaling principles. New reactor designs feature different operating conditions and scenarios than those in existing reactors. ARCO supports the research and development of digital tools for operator communications, intuitive real-time data analysis, online health monitoring and prognostics, and control room cybersecurity. By integrating these different technologies, ARCO acts as a prototypical control system to iteratively develop methods and tools of operation in advanced small modular nuclear reactors. This paper describes the features of and challenges to operating advanced small modular reactors underlying the design basis for ARCO and its operator support systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

THE EFFECTS OF CLOSE COMPANIONS (AND ROTATION) ON THE MAGNETIC ACTIVITY OF M DWARFS

We present a study of close white dwarf and M dwarf (WD+dM) binary systems and examine the effect that a close companion has on the magnetic field generation in M dwarfs. We use a base sample of 1602 white dwarf main-sequence binaries from Rebassa-Mansergas et al. to develop a set of color cuts in GALEX, SDSS, UKIDSS, and 2MASS color space. Then using the SDSS Data Release 8 spectroscopic database, we construct a sample of 1756 WD+dM high-quality pairs from our color cuts and previous catalogs. We separate the individual WD and dM from each spectrum using an iterative technique that compares the WD and dM components to best-fit templates. Using the absolute height above the Galactic plane as a proxy for age, and the H{alpha} emission line as an indicator for magnetic activity, we investigate the age-activity relation for our sample for spectral types {<=} M7. Our results show that early-type M dwarfs ({<=}M4) in close binary systems are more likely to be active and have longer activity lifetimes compared to their field counterparts. However, at a spectral type of M5 (just past the onset of full convection in M dwarfs), the activity fraction and lifetimes of WD+dM binary systems become more comparable to that of the field M dwarfs. One of the implications of having a close binary companion is presumed to be increased stellar rotation through disk disruption, tidal effects, or angular momentum exchange. Thus, we interpret the similarity in activity behavior between late-type dMs in WD+dM pairs and late-type field dMs to be due to a decrease in sensitivity in close binary companions (or stellar rotation), which has implications for the nature of magnetic activity in fully convective stars. Using the WD components of the pairs, we find WD cooling ages to use as an additional constraint on the age-activity relation for our sample. We find that, on average, active early-type dMs tend to be younger and that active late-type dMs span a much broader age regime making them indistinguishable from the inactive late-type population. We also show that magnetic strength, as measured by H{alpha}, is comparable between paired and field M dwarfs until a spectral type of M6/M7 where M dwarf activity for stars with close companions becomes much stronger. In addition, we present 37 very close candidate pairs with fast-moving orbits that display radial velocity changes over hour timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE McGill PLANAR HYDROGEN ATMOSPHERE CODE (McPHAC)

The McGill Planar Hydrogen Atmosphere Code (McPHAC) v1.1 calculates the hydrostatic equilibrium structure and emergent spectrum of an unmagnetized hydrogen atmosphere in the plane-parallel approximation, at surface gravities appropriate for neutron stars. McPHAC incorporates several improvements over previous codes for which tabulated model spectra are available: (1) Thomson scattering is treated anisotropically, which is shown to result in a 0.2%-3% correction in the emergent spectral flux across the 0.1-5 keV passband; (2) the McPHAC source code is made available to the community, allowing it to be scrutinized and modified by other researchers wishing to study or extend its capabilities; and (3) the numerical uncertainty resulting from the discrete and iterative solution is studied as a function of photon energy, indicating that McPHAC is capable of producing spectra with numerical uncertainties <0.01%. The accuracy of the spectra may at present be limited to {approx}1%, but McPHAC enables researchers to study the impact of uncertain inputs and additional physical effects, thereby supporting future efforts to reduce those inaccuracies. Comparison of McPHAC results with spectra from one of the previous model atmosphere codes (NSA) shows agreement to {approx}<1% near the peaks of the emergent spectra. However, in the Wien tail a significant deficit of flux in the spectra of the previous model is revealed, determined to be due to the previous work not considering large enough optical depths at the highest photon frequencies. The deficit is most significant for spectra with T{sub eff} < 10{sup 5.6} K, though even there it may not be of much practical importance for most observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

DERIVING PLASMA DENSITIES AND ELEMENTAL ABUNDANCES FROM SERTS DIFFERENTIAL EMISSION MEASURE ANALYSIS

We use high-resolution spectral emission line data obtained by the SERTS instrument during three rocket flights to demonstrate a new approach for constraining electron densities of solar active region plasma. We apply differential emission measure (DEM) forward-fitting techniques to characterize the multithermal solar plasma producing the observed EUV spectra, with constraints on the high-temperature plasma from the Yohkoh Soft X-ray Telescope. In this iterative process, we compare line intensities predicted by an input source distribution to observed line intensities for multiple iron ion species, and search a broad range of densities to optimize {chi}{sup 2} simultaneously for the many available density-sensitive lines. This produces a density weighted by the DEM, which appears to be useful for characterizing the bulk of the emitting plasma over a significant range of temperature. This 'DEM-weighted density' technique is complementary to the use of density-sensitive line ratios and less affected by uncertainties in atomic data and ionization fraction for any specific line. Once the DEM shape and the DEM-weighted density have been established from the iron lines, the relative elemental abundances can be determined for other lines in the spectrum. We have also identified spectral lines in the SERTS wavelength range that may be problematic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ON THE GRAVITATIONAL FIELDS OF MACLAURIN SPHEROID MODELS OF ROTATING FLUID PLANETS

Hubbard recently derived an important iterative equation for calculating the gravitational coefficients of a Maclaurin spheroid that does not require an expansion in a small distortion parameter. We show that this iterative equation, which is based on an incomplete solution of the Poisson equation, diverges when the distortion parameter is not sufficiently small. We derive a new iterative equation that is based on a complete solution of the Poisson equation and, hence, always converges when calculating the gravitational coefficients of a Maclaurin spheroid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

POINT-SPREAD FUNCTIONS FOR THE EXTREME-ULTRAVIOLET CHANNELS OF SDO/AIA TELESCOPES

We present the stray-light point-spread functions (PSFs) and their inverses we characterized for the Atmospheric Imaging Assembly (AIA) EUV telescopes on board the Solar Dynamics Observatory (SDO) spacecraft. The inverse kernels are approximate inverses under convolution. Convolving the original Level 1 images with them produces images with improved stray-light characteristics. We demonstrate the usefulness of these PSFs by applying them to two specific cases: photometry and differential emission measure (DEM) analysis. The PSFs consist of a narrow Gaussian core, a diffraction component, and a diffuse component represented by the sum of a Gaussian-truncated Lorentzian and a shoulder Gaussian. We determined the diffraction term using the measured geometry of the diffraction pattern identified in flare images and the theoretically computed intensities of the principal maxima of the first few diffraction orders. To determine the diffuse component, we fitted its parameterized model using iterative forward-modeling of the lunar interior in the SDO/AIA images from the 2011 March 4 lunar transit. We find that deconvolution significantly improves the contrast in dark features such as miniature coronal holes, though the effect was marginal in bright features. On a percentage-scattering basis, the PSFs for SDO/AIA are better by a factor of two than that of the EUV telescope on board the Transition Region And Coronal Explorer mission. A preliminary analysis suggests that deconvolution alone does not affect DEM analysis of small coronal loop segments with suitable background subtraction. We include the derived PSFs and their inverses as supplementary digital materials.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

FOCAL PLANE WAVEFRONT SENSING USING RESIDUAL ADAPTIVE OPTICS SPECKLES

Optical imperfections, misalignments, aberrations, and even dust can significantly limit sensitivity in high-contrast imaging systems such as coronagraphs. An upstream deformable mirror (DM) in the pupil can be used to correct or compensate for these flaws, either to enhance the Strehl ratio or suppress the residual coronagraphic halo. Measurement of the phase and amplitude of the starlight halo at the science camera is essential for determining the DM shape that compensates for any non-common-path (NCP) wavefront errors. Using DM displacement ripples to create a series of probe and anti-halo speckles in the focal plane has been proposed for space-based coronagraphs and successfully demonstrated in the lab. We present the theory and first on-sky demonstration of a technique to measure the complex halo using the rapidly changing residual atmospheric speckles at the 6.5 m MMT telescope using the Clio mid-IR camera. The AO system's wavefront sensor measurements are used to estimate the residual wavefront, allowing us to approximately compute the rapidly evolving phase and amplitude of speckle halo. When combined with relatively short, synchronized science camera images, the complex speckle estimates can be used to interferometrically analyze the images, leading to an estimate of the static diffraction halo with NCP effects included. In an operational system, this information could be collected continuously and used to iteratively correct quasi-static NCP errors or suppress imperfect coronagraphic halos.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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↗