Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ray trace”

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 361 records · Page 20

Predicting Weak Lensing Statistics from Halo Mass Reconstructions (Final Paper)

As dark matter does not absorb or emit light, its distribution in the universe must be inferred through indirect effects such as the gravitational lensing of distant galaxies. While most sources are only weakly lensed, the systematic alignment of background galaxies around a foreground lens can constrain the mass of the lens which is largely in the form of dark matter. In this paper, I have implemented a framework to reconstruct all of the mass along lines of sight using a best-case dark matter halo model in which the halo mass is known. This framework is then used to make predictions of the weak lensing of 3,240 generated source galaxies through a 324 arcmin 2 field of the Millennium Simulation. The lensed source ellipticities are characterized by the ellipticity-ellipticity and galaxy-mass correlation functions and compared to the same statistic for the intrinsic and ray-traced ellipticities. In the ellipticity-ellipticity correlation function, I and that the framework systematically under predicts the shear power by an average factor of 2.2 and fails to capture correlation from dark matter structure at scales larger than 1 arcminute. The model predicted galaxy-mass correlation function is in agreement with the ray-traced statistic from scales 0.2 to 0.7 arcminutes, but systematically underpredicts shear power at scales larger than 0.7 arcminutes by an average factor of 1.2. Optimization of the framework code has reduced the mean CPU time per lensing prediction by 70% to 24 ± 5 ms. Physical and computational shortcomings of the framework are discussed, as well as potential improvements for upcoming work.

79 ASTRONOMY AND ASTROPHYSICS↗

Polarization Considerations for the Laser Interferometer Space Antenna (LISA)

A polarization ray trace model of the Laser Interferometer Space Antenna's (LISA) optical path is being created. The model will be able to assess the effects of various polarizing elements and the optical coatings on the required picometer level interferometry. All of the computational steps are described in detail. This should eliminate any ambiguities associated with polarization ray trace modeling of interferometers and provide a basis for determining its limitations and serve as a clearly defined starting point for future improvements.

Waluschka, Eugene↗

SU-D-213-03: Towards An Optimized 3D Scintillation Dosimetry Tool for Quality Assurance of Dynamic Radiotherapy Techniques

Purpose: The purpose of this work is to simulate a multi-focus plenoptic camera used as the measuring device in a real-time three-dimensional scintillation dosimeter. Simulating and optimizing this realistic optical system will bridge the technological gap between concept validation and a clinically viable tool that can provide highly efficient, accurate and precise measurements for dynamic radiotherapy techniques. Methods: The experimental prototype, previously developed for proof of concept purposes, uses an off-the-shelf multi-focus plenoptic camera. With an array of interleaved microlenses of different focal lengths, this camera records spatial and angular information of light emitted by a plastic scintillator volume. The three distinct microlens focal lengths were determined experimentally for use as baseline parameters by measuring image-to-object magnification for different distances in object space. A simulated plenoptic system was implemented using the non-sequential ray tracing software Zemax: this tool allows complete simulation of multiple optical paths by modeling interactions at interfaces such as scatter, diffraction, reflection and refraction. The active sensor was modeled based on the camera manufacturer specifications by a 2048×2048, 5 µm-pixel pitch sensor. Planar light sources, simulating the plastic scintillator volume, were employed for ray tracing simulations. Results: The microlens focal lengths were determined to be 384, 327 and 290 µm. A realistic multi-focus plenoptic system, with independently defined and optimizable specifications, was fully simulated. A f/2.9 and 54 mm-focal length Double Gauss objective was modeled as the system’s main lens. A three-focal length hexagonal microlens array of 250-µm thickness was designed, acting as an image-relay system between the main lens and sensor. Conclusion: Simulation of a fully modeled multi-focus plenoptic camera enables the decoupled optimization of the main lens and microlens specifications. This work leads the way to improving the 3D dosimeter’s achievable resolution, efficiency and build for providing a quality assurance tool fully meeting clinical needs. M.R. is financially supported by a Master’s Canada Graduate Scholarship from the NSERC. This research is also supported by the NSERC Industrial Research Chair in Optical Design.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Polarization Considerations for the Laser Interferometer Space Antenna

A polarization ray trace model of the Laser Interferometer Space Antenna s (LISA) optical path is being created. The model will be able to assess the effects of various polarizing elements and the optical coatings on the required, very long path length, picometer level dynamic interferometry. The computational steps are described. This should eliminate any ambiguities associated with polarization ray tracing of interferometers and provide a basis for determining the computer model s limitations and serve as a clearly defined starting point for future work.

Waluschka, Eugene↗

Optimizing Lunar Map Partitioning for Multipath Fade Loss Analyses

A methodology is presented pertaining to the run-time optimization of quasi-optical ray-tracing simulations. These shooting-and-bouncing ray methods are used for assessing multipath fade loss in cislunar communication links. Our approach employs partitioning lunar terrain maps in shapes that approximate two-dimensional projections of Fresnel zone cross-sections. This technique drastically reduces computational run-time without causing significant loss of simulation accuracy.

Digital elevation models↗

New simulation capability for gamma ray mirror experiments

This report provides a description of the simulation toolkit developed at Lawrence Livermore National Laboratory to support the design of nuclear safeguards experiments using grazing incidence multilayer mirrors in the energy band of uranium (U) and plutonium (Pu) emission lines. This effort was motivated by the data analysis of a scoping experiment at the Irradiated Fuels Examination Facility (IFEL) at Oak Ridge National Laboratory in FY13 and of a benchmark experiment at the Idaho National Laboratory (INL) in FY14 that highlighted the need for predictive tools built around a ray-tracing capability. This report presents the simulation toolkit and relevant results such as the simulated spectra for TMI, MOX, and ATM106 fuel rods based on spent fuel models provided by Los Alamos National Laboratory and for a virgin high 240 Pu-content fuel plate, as well as models of the IFEL and INL experiments implemented in the ray tracing tool. The beam position and height were validated against the INL ~60 keV americium data. Examples of alternate configurations of the optics or experimental set-up illustrate the future use of the simulation suite to guide the next IFEL experimental campaign.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗

Solar Thermal Concept Evaluation

Concentrated solar thermal energy can be utilized in a variety of high temperature applications for both terrestrial and space environments. In each application, knowledge of the collector and absorber's heat exchange interaction is required. To understand this coupled mechanism, various concentrator types and geometries, as well as, their relationship to the physical absorber mechanics were investigated. To conduct experimental tests various parts of a 5,000 watt, thermal concentrator, facility were made and evaluated. This was in anticipation at a larger NASA facility proposed for construction. Although much of the work centered on solar thermal propulsion for an upper stage (less than one pound thrust range), the information generated and the facility's capabilities are applicable to material processing, power generation and similar uses. The numerical calculations used to design the laboratory mirror and the procedure for evaluating other solar collectors are presented here. The mirror design is based on a hexagonal faceted system, which uses a spherical approximation to the parabolic surface. The work began with a few two dimensional estimates and continued with a full, three dimensional, numerical algorithm written in FORTRAN code. This was compared to a full geometry, ray trace program, BEAM 4, which optimizes the curvatures, based on purely optical considerations. Founded on numerical results, the characteristics of a faceted concentrator were construed. The numerical methodologies themselves were evaluated and categorized. As a result, the three-dimensional FORTRAN code was the method chosen to construct the mirrors, due to its overall accuracy and superior results to the ray trace program. This information is being used to fabricate and subsequently, laser map the actual mirror surfaces. Evaluation of concentrator mirrors, thermal applications and scaling the results of the 10 foot diameter mirror to a much larger concentrator, were studied. Evaluations, recommendations and pit falls regarding the structure, materials and facility design are presented.

Hawk, Clark W.↗

Extension of the flux fit method for estimating power deposition profiles

The flux fit method is used to self-consistently estimate the power deposition profile and heat transport profiles from temperature measurements originating from perturbative experiments with a modulated source. Furthermore, this letter improves on this method by addressing the limitations and assumptions. The most crucial improvement is the additional freedom in the source deposition profile. Allowing for a variable central deposition location and height and including a skewness parameter produces deposition profiles more consistent with the measurement data, but still wider than equilibrium ray tracing in two different DIII-D discharges. Moreover, we show that the quality of the estimated deposition profile is key to the accuracy of diffusivity and connectivity estimates, but inversely, the estimated transport parameters hardly affect the quality of the power deposition estimate. Using this method, we show that the power deposition profile estimate is broadened with respect to ray-tracing by about 1.7-1.8 times in two DIII-D discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport

Fast-pairwise neutrino oscillations potentially affect many aspects of core-collapse supernovae (CCSNe): the explosion mechanism, neutrino signals, and nucleosynthesis in the ejecta. This particular mode of collective neutrino oscillations has a deep connection to the angular structure of neutrinos in momentum space; for instance, the appearance of electron-neutrino lepton number (ELN) angular crossings in momentum space is a good indicator of the occurrence of flavor conversions. However, many multidimensional CCSN simulations are carried out with approximate neutrino transport (such as two-moment methods), which limits access to the angular distributions of neutrinos, i.e., inhibits ELNcrossing searches. In this paper, we develop a new method of searching for ELN crossing in these CCSN simulations. The required data is the zeroth and first angular moments of neutrinos and the matter profile, all of which are available in CCSN models with the two-moment method. One of the novelties of our new method is the use of a ray-tracing neutrino transport to determine ELNs in the direction of the stellar center. It is designed to compensate for shortcomings of the crossing searches with only two angular moments. We assess the capability of the method by carrying out a detailed comparison with results of full Boltzmann neutrino transport in 1D and 2D CCSN models. In this work, we find that the ray-tracing neutrino transport improves the accuracy of crossing searches; indeed, the appearance/disappearance of the crossings is accurately detected even in the region of forward-peaked angular distributions. The new method is computationally cheap and has the benefit of efficient parallelization; hence, it will be useful for ELN-crossing searches in any CCSN models that employ two-moment neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗

Formula for the rms blur circle radius of Wolter telescope based on aberration theory

A formula for the rms blur circle for Wolter telescopes has been derived using the transverse ray aberration expressions of Saha (1985), Saha (1984), and Saha (1986). The resulting formula for the rms blur circle radius over an image plane and a formula for the surface of best focus based on third-, fifth-, and seventh-order aberration theory predict results in good agreement with exact ray tracing. It has also been shown that one of the two terms in the empirical formula of VanSpeybroeck and Chase (1972), for the rms blur circle radius of a Wolter I telescope can be justified by the aberration theory results. Numerical results are given comparing the rms blur radius and the surface of best focus vs the half-field angle computed by skew ray tracing and from analytical formulas for grazing incidence Wolter I-II telescopes and a normal incidence Cassegrain telescope.

Shealy, David L.↗

PICS: SIMULATIONS OF STRONG GRAVITATIONAL LENSING IN GALAXY CLUSTERS

Gravitational lensing has become one of the most powerful tools available for investigating the “dark side” of the universe. Cosmological strong gravitational lensing, in particular, probes the properties of the dense cores of dark matter halos over decades in mass and offers the opportunity to study the distant universe at flux levels and spatial resolutions otherwise unavailable. Studies of strongly lensed variable sources offer even further scientific opportunities. One of the challenges in realizing the potential of strong lensing is to understand the statistical context of both the individual systems that receive extensive follow-up study, as well as that of the larger samples of strong lenses that are now emerging from survey efforts. Motivated by these challenges, we have developed an image simulation pipeline, Pipeline for Images of Cosmological Strong lensing (PICS), to generate realistic strong gravitational lensing signals from group- and cluster-scale lenses. PICS uses a low-noise and unbiased density estimator based on (resampled) Delaunay Tessellations to calculate the density field; lensed images are produced by ray-tracing images of actual galaxies from deep Hubble Space Telescope observations. Other galaxies, similarly sampled, are added to fill in the light cone. The pipeline further adds cluster member galaxies and foreground stars into the lensed images. The entire image ensemble is then observed using a realistic point-spread function that includes appropriate detector artifacts for bright stars. Noise is further added, including such non-Gaussian elements as noise window-paning from mosaiced observations, residual bad pixels, and cosmic rays. The aim is to produce simulated images that appear identical—to the eye (expert or otherwise)—to real observations in various imaging surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

PROBING PRIMORDIAL NON-GAUSSIANITY WITH WEAK-LENSING MINKOWSKI FUNCTIONALS

We study the cosmological information contained in the Minkowski functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local-type primordial non-Gaussianity parameter f {sub NL}. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing to generate 100 independent convergence maps of a 25 deg{sup 2} field of view for f {sub NL} = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude {sigma}{sub 8}. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For upcoming wide-field observations such as those from the Subaru Hyper Suprime-Cam survey with a proposed survey area of 1500 deg{sup 2}, the primordial non-Gaussianity can be constrained with a level of f {sub NL} {approx} 80 and w {approx} 0.036 by weak-lensing MFs. If simply scaled by the effective survey area, a 20,000 deg{sup 2} lensing survey using the Large Synoptic Survey Telescope will yield constraints of f {sub NL} {approx} 25 and w {approx} 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.

79 ASTRONOMY AND ASTROPHYSICS↗

Design and analysis of optical systems for the Stanford/MSFC Multi-Spectral Solar Telescope Array

This paper reports on the design and the theoretical ray trace analysis of the optical systems which will comprise the primary imaging components for the Stanford/MSFC Multi-Spectral Solar Telescope Array (MSSTA). This instrument is being developed for ultra-high resolution investigations of the sun from a sounding rocket. Doubly reflecting systems of sphere-sphere, ellipsoid-sphere (Dall-Kirkham), paraboloid-hyperboloid (Cassegrain), and hyperboloid-hyperboloid (Ritchey-Chretien) configurations were analyzed. For these mirror systems, ray trace analysis was performed and through-focus spot diagrams, point spread function plots, and geometrical and diffraction MTFs were generated. The results of these studies are presented along with the parameters of the Ritchey-Chretien optical system selected for the MSSTA flight. The payload, which incorporates seven of these Ritchey-Chretien systems, is now being prepared for launch in late September 1989.

Hadaway, James B.↗

Modeling and Simulation of Radiative Compressible Flows in Aerodynamic Heating Arc-Jet Facility

Numerical simulations of an arc heated flow inside NASA's 20 [MW] Aerodynamics heating facility (AHF) are performed in order to investigate the three-dimensional swirling flow and the current distribution inside the wind tunnel. The plasma is considered in Local Thermodynamics Equilibrium(LTE) and is composed of Air-Argon gas mixture. The governing equations are the Navier-Stokes equations that include source terms corresponding to Joule heating and radiative cooling. The former is obtained by solving an electric potential equation, while the latter is calculated using an innovative massively parallel ray-tracing algorithm. The fully coupled system is closed by the thermodynamics relations and transport properties which are obtained from Chapman-Enskog method. A novel strategy was developed in order to enable the flow solver and the radiation calculation to be preformed independently and simultaneously using a different number of processors. Drastic reduction in the computational cost was achieved using this strategy. Details on the numerical methods used for space discretization, time integration and ray-tracing algorithm will be presented. The effect of the radiative cooling on the dynamics of the flow will be investigated. The complete set of equations were implemented within the COOLFluiD Framework. Fig. 1 shows the geometry of the Anode and part of the constrictor of the Aerodynamics heating facility (AHF). Fig. 2 shows the velocity field distribution along (x-y) plane and the streamline in (z-y) plane.

Pasma flows↗

Steering of Solar Sails Using Optical Lift Force

Optical wing structures were theoretically and numerically analyzed, and prototype arrays of wings called optical flying carpets were fabricated with solar sail material clear polyimide (CP1). This material was developed at NASA Langley to better withstand damaging ultraviolet radiation found in outer space. Various optical wing sizes and shapes were analyzed to develop design strategies for thrust and torque applications. The developed ray-tracing model has undergone continual advancement, and stands as an effective tool for modeling most types of solar sails. To our understanding, such a model does not exist else where. The distributed forces and torques have been reduced to a simple theoretical whereby the fundamental mechanics may be understood in terms of the numerically determined center of pressure off set from the center of mass. This description applies to any type of solar sail, affording our ray-tracing model a general utility. This research has established a foundation for understanding the force and torque afforded by optical wings. The study began by considering transparent wings and ended by considering wings having a reflecting face. The latter was found to afford the advantages of high thrust and both intrinsic and extrinsic torque. Our discovery of the intrinsic torque on optical wings (meaning that a moment arm is not required) has no analogy for a flat reflective solar sail, and therefore provides an extra degree of control that may be useful for sail craft attitude and navigation purposes.

Swartzlander, Gover↗

Aperture alignment in autocollimator-based deflectometric profilometers

During the last ten years, deflectometric profilometers have become indispensable tools for the precision form measurement of optical surfaces. They have proven to be especially suitable for characterizing beam-shaping optical surfaces for x-ray beamline applications at synchrotrons and free electron lasers. Deflectometric profilometers use surface slope (angle) to assess topography and utilize commercial autocollimators for the contactless slope measurement. To this purpose, the autocollimator beam is deflected by a movable optical square (or pentaprism) towards the surface where a co-moving aperture limits and defines the beam footprint. In this paper, we focus on the precise and reproducible alignment of the aperture relative to the autocollimator's optical axis. Its alignment needs to be maintained while it is scanned across the surface under test. The reproducibility of the autocollimator's measuring conditions during calibration and during its use in the profilometer is of crucial importance to providing precise and traceable angle metrology. In the first part of the paper, we present the aperture alignment procedure developed at the Advanced Light Source, Lawrence Berkeley National Laboratory, USA, for the use of their deflectometric profilometers. In the second part, we investigate the topic further by providing extensive ray tracing simulations and calibrations of a commercial autocollimator performed at the Physikalisch-Technische Bundesanstalt, Germany, for evaluating the effects of the positioning of the aperture on the autocollimator's angle response. The investigations which we performed are crucial for reaching fundamental metrological limits in deflectometric profilometry.

47 OTHER INSTRUMENTATION↗

GEOS satellite tracking corrections for refraction in the ionosphere

The analytic formulations at different elevation angles and at a frequency of 2-GHz for the ionospheric refraction corrections used on the GEOS satellite tracking data are compared. The formulas and ray trace results for elevations greater than 10 deg, where most satellite tracking is done, differ in elevation, range, and range rate by less than 0.4 millidegrees (1.4 arc-seconds), 12 meters, and 12 cm/sec, respectively. In comparison to most operational requirements, this is insignificant. However, for the GEOS Observation Systems Intercomparison Investigation, these differences are equivalent in size to observed differences in system biases for some of the best electronic geodetic tracking systems and are probably contributing to the observed biases. The ray trace results and most of the more detailed analytic correction formulas show that the ionospheric refraction correction for range rate on an overhead pass is a maximum for elevation angles between 15 deg and 30 deg and falls off rapidly for both higher and lower elevation angles, contrary to the effect of the troposphere and to some reports in the literature.

Berbert, J. H.↗