Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational reflection”

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 865 records · Page 48

Planetary Rover Simulation for Lunar Exploration Missions

When planning planetary rover missions it is useful to develop intuition and skills driving in, quite literally, alien environments before incurring the cost of reaching said locales. Simulators make it possible to operate in environments that have the physical characteristics of target locations without the expense and overhead of extensive physical tests. To that end, NASA Ames and Open Robotics collaborated on a Lunar rover driving simulator based on the open source Gazebo simulation platform and leveraging ROS (Robotic Operating System) components. The simulator was integrated with research and mission software for rover driving, system monitoring, and science instrument simulation to constitute an end-to-end Lunar mission simulation capability. Although we expect our simulator to be applicable to arbitrary Lunar regions, we designed to a reference mission of prospecting in polar regions. The harsh lighting and low illumination angles at the Lunar poles combine with the unique reflectance properties of Lunar regolith to present a challenging visual environment for both human and computer perception. Our simulator placed an emphasis on high fidelity visual simulation in order to produce synthetic imagery suitable for evaluating human rover drivers with navigation tasks, as well as providing test data for computer vision software development.In this paper, we describe the software used to construct the simulated Lunar environment and the components of the driving simulation. Our synthetic terrain generation software artificially increases the resolution of Lunar digital elevation maps by fractal synthesis and inserts craters and rocks based on Lunar size-frequency distribution models. We describe the necessary enhancements to import large scale, high resolution terrains into Gazebo, as well as our approach to modeling the visual environment of the Lunar surface. An overview of the mission software system is provided, along with how ROS was used to emulate flight software components that had not been developed yet. Finally, we discuss the effect of using the high-fidelity synthetic Lunar images for visual odometry. We also characterize the wheel slip model, and find some inconsistencies in the produced wheel slip behaviour.

Allan, Mark↗

Systems and methods for plane-wave and fan-beam ultrasound-waveform tomography

A method for plane-wave or fan-beam ultrasound tomography includes: controlling, by a processor, a multi-channel ultrasound array controller to control a linear ultrasound transducer array to emit plane-wave or fan-beam ultrasound into a biological tissue, the linear ultrasound transducer array including a plurality of transducers; receiving, by the processor, plane-wave or fan-beam ultrasound waveform data from the transducers, the plane-wave or fan-beam ultrasound waveform data including reflections of the plane-wave or fan-beam ultrasound from the biological tissue; computing, by the processor, a synthetic ultrasound-waveform inversion from the received plane-wave or fan-beam ultrasound waveform data; and computing, by the processor, an image of the biological tissue from the synthetic ultrasound-waveform inversion.

Huang, Lianjie↗

Extracting Depth From Motion Parallax in Real-World and Synthetic Displays

In psychophysical studies on human sensitivity to visual motion parallax (MP), the use of computer displays is pervasive. However, a number of potential problems are associated with such displays: cue conflicts arise when observers accommodate to the screen surface, and observer head and body movements are often not reflected in the displays. We investigated observers' sensitivity to depth information in MP (slant, depth order, relative depth) using various real-world displays and their computer-generated analogs. Angle judgments of real-world stimuli were consistently superior to judgments that were based on computer-generated stimuli. Similar results were found for perceived depth order and relative depth. Perceptual competence of observers tends to be underestimated in research that is based on computer generated displays. Such findings cannot be generalized to more realistic viewing situations.

Hecht, Heiko↗

Theory of a Traveling Wave Feed for a Planar Slot Array Antenna

Planar arrays of waveguide-fed slots have been employed in many radar and remote sensing applications. Such arrays are designed in the standing wave configuration because of high efficiency. Traveling wave arrays can produce greater bandwidth at the expense of efficiency due to power loss in the load or loads. Traveling wave planar slot arrays may be designed with a long feed waveguide consisting of centered-inclined coupling slots. The feed waveguide is terminated in a matched load, and the element spacing in the feed waveguide is chosen to produce a beam squinted from the broadside. The traveling wave planar slot array consists of a long feed waveguide containing resonant-centered inclined coupling slots in the broad wall, coupling power into an array of stacked radiating waveguides orthogonal to it. The radiating waveguides consist of longitudinal offset radiating slots in a standing wave configuration. For the traveling wave feed of a planar slot array, one has to design the tilt angle and length of each coupling slot such that the amplitude and phase of excitation of each radiating waveguide are close to the desired values. The coupling slot spacing is chosen for an appropriate beam squint. Scattering matrix parameters of resonant coupling slots are used in the design process to produce appropriate excitations of radiating waveguides with constraints placed only on amplitudes. Since the radiating slots in each radiating waveguide are designed to produce a certain total admittance, the scattering (S) matrix of each coupling slot is reduced to a 2x2 matrix. Elements of each 2x2 S-matrix and the amount of coupling into the corresponding radiating waveguide are expressed in terms of the element S11. S matrices are converted into transmission (T) matrices, and the T matrices are multiplied to cascade the coupling slots and waveguide sections, starting from the load end and proceeding towards the source. While the use of non-resonant coupling slots may provide an additional degree of freedom in the design, resonant coupling slots simplify the design process. The amplitude of the wave going to the load is set at unity. The S11 parameter, r of the coupling slot closest to the load, is assigned an arbitrary value. A larger value of r will reduce the power dissipated in the load while increasing the reflection coefficient at the input port. It is now possible to obtain the excitation of the radiating waveguide closest to the load and the coefficients of the wave incident and reflected at the input port of this coupling slot. The next coupling slot parameter, r , is chosen to realize the excitation of that radiating waveguide. One continues this process moving towards the source, until all the coupling slot parameters r and hence the S11 parameter of the 4-port coupler, r, are known for each coupling slot. The goal is to produce the desired array aperture distribution in the feed direction. From an interpolation of the computed moment method data for the slot parameters, all the coupling slot tilt angles and lengths are obtained. From the excitations of the radiating waveguides computed from the coupling values, radiating slot parameters may be obtained so as to attain the desired total normalized slot admittances. This process yields the radiating slot parameters, offsets, and lengths. The design is repeated by choosing different values of r for the last coupling slot until the percentage of power dissipated in the load and the input reflection coefficient values are satisfactory. Numerical results computed for the radiation pattern, the tilt angles and lengths of coupling slots, and excitation phases of the radiating waveguides, are presented for an array with uniform amplitude excitation. The design process has been validated using computer simulations. This design procedure is valid for non-uniform amplitude excitations as well.

Rengarajan, Sembiam↗

Diverging polygon-based modeling (DPBM) of concentrated solar flux distributions

This paper presents an efficient and robust methodology for modeling concentrated solar flux distributions. Compared to ray tracing methods, which provide high accuracy but can be computationally intensive, this approach makes a number of simplifying assumptions in order to reduce complexity by modeling incident and reflected flux as a series of simple geometric diverging polygons, then applying shading and blocking effects. A reduction in processing time (as compared to ray tracing) allows for evaluating and visualizing numerous combinations of engineering and operational variables (easily exceeding 106 unique iterations) to ascertain instantaneous, transient, and annual system performance. The method is demonstrated on a linear Fresnel reflector array and a number of variable iteration examples presented. While some precision is sacrificed for computational speed, flux distributions were compared to ray tracing (SolTrace) and average concentration ratio generally found to agree within ~3%. Furthermore, this method presents a quick and very flexible coarse adjust method for concentrated solar power (CSP) field design, and can be used to both rapidly gain an understanding of system performance as well as to narrow variable constraint windows for follow-on high accuracy system optimization.

14 SOLAR ENERGY↗

Activation of Water by Pentavalent Actinide Dioxide Cations: Characteristic Curium Revealed by a Reactivity Turn after Americium

Swapping of an oxygen atom of water with that of a pentavalent actinide dioxide cation, AnO 2 + also called an “actinyl”, requires activation of an An–O bond. It was previously found that such oxo exchange in the gas phase occurs for the first two actinyls, PaO 2 + and UO 2 + , but not the next two, NpO 2 + and PuO 2 + . The An–O bond dissociation energies (BDEs) decrease from PaO 2 + to PuO 2 + , such that the observation of a parallel decrease in the An–O bond reactivity is intriguing. To elucidate oxo exchange, we here extend experimental studies to AmO 2 + , americyl(V), and CmO 2 + , curyl(V), which were produced in remarkable abundance by electrospray ionization of Am 3+ and Cm 3+ solutions. Like other AnO 2 + , americyl(V) and curyl(V) adsorb up to four H 2 O molecules to form tetrahydrates AnO 2 (H 2 O) 4 + with the actinide hexacoordinated by oxygen atoms. It was found that AmO 2 + does not oxo-exchange, whereas CmO 2 + does, establishing a “turn” to increasing the reactivity from americyl to curyl, which validates computational predictions. Because oxo exchange occurs via conversion of an actinyl(V) hydrate, AnO 2 (H 2 O) + , to an actinide(V) hydroxide, AnO(OH) 2 + , it reflects the propensity for actinyl(V) hydrolysis: PaO 2 + hydrolyzes and oxo-exchanges most easily, despite the fact that it has the highest BDE of all AnO 2 + . We note a reexamination of the computational results for actinyl(V) oxo exchange reveals distinctive properties and chemistry of curyl(V) species, particularly CmO(OH) 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vectorial approach of determining the wave propagation at metasurfaces

Vector approach often benefits optical engineers and physicists, and a vector formulation of the laws of reflection and refraction has been studied (Tkaczyk, 2012). However, the conventional reflection and refraction laws may be violated in the presence of a metasurface, and reflection and refraction at the metasurface obey generalized laws of reflection and refraction (Yu et al., 2011). In this letter, the vectorial laws of reflection and refraction at the metasurface were derived, and the matrix formulation of these vectorial laws are also obtained. These results enable highly efficient and unambiguous computations in ray-tracing problems that involve a metasurface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational study of critical mass curves for various ratios of {sup 238}Pu to {sup 239}Pu

A MCNP study for a water-reflected sphere of plutonium metal water pseudo-solution was performed for a various ratios of {sup 238}Pu to {sup 239}Pu to develop minimum critical mass curves. The {sup 238}Pu atom percent (at. %) was varied from 35 to 100 percent. It was demonstrated that at {sup 238}Pu ratios of ≥ 65 at. % the high absorption cross section of {sup 238}Pu in the thermal energy region dominates the system reactivity for water moderated systems resulting in a minimum critical mass in the metal rather than solution regime. However, due to the lack of available benchmark experiments to validated use of the {sup 238}Pu cross-section data for quantities above gram-level, the models fell outside the area of applicability (AOA) of the Los Alamos National Laboratory's (LANL) Nuclear Criticality Safety Division's (NCSD) MCNP validation report. An additional of margin of subcriticality (MOS) is needed for {sup 238}Pu systems. This additional MOS will be dependent on the amount and significance of the {sup 238}Pu, which varies in the MCNP study. Whisper is a computational software package designed to assist with validation studies for MCNP. Whisper uses sensitivity/uncertainty (S/U) methods to select relevant benchmarks to a particular application model. Therefore, Whisper 1.0.0 was used on the curve minima in the solution regime for the 65 at. % {sup 238}Pu curve to estimate an appropriate additional margin of subcriticality and create a subcritical curve based on the suggested upper subcritical limits (USL)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bidirectional reflectance of oceanic waters: A comparison of modeled and measured upward radiance fields

The bidirectional reflectance of oceanic waters is conveniently described in a normalized way by forming the ratio of the upwelling irradiance E(sub u) to any upwelling radiance L(sub u)(theta prime, phi). This ratio, Q (theta prime, theta(sub 0), (phi(sub 0) - phi), where theta prime, phi are the nadir and azimuth angles for the upward radiance and theta(sub 0), phi(sub 0) are the zenith and azimuth angles of the Sun, has been determined from measurements at sea and computed via Monte Carlo simulations using the inherent optical properties measured in the field and appropriate boundary conditions (clear sky, no wind, varying Sun angle). Experimental ad computed Q values are in excellent agreement. This successful comparison confirms the importance of the bidirectional character of ocean reflectance, already pointed out from a purely numerical approach without field validation, and corroborates the extended range of the Q variations. The later point is of importance when interpreting the marine signals detected by an ocean color satellite-borne sensor. The validation is extended by considering the historical data for the radiance distributions in Lake Pend Oreille determined at various depths. The closure issue in ocean optics is examined by solving the direct problem of radiative transfer and through a model-data comparison in terms of radiance field.

Morel, Andre↗

An Adaptive Tracking Algorithm for Convection in Simulated and Remote Sensing Data

A robust and computationally efficient object tracking algorithm is developed by incorporating various tracking techniques. Physical properties of the objects, such as brightness temperature or reflectivity, are not considered. Therefore, the algorithm is adaptable for tracking convection-like features in simulated data and remotely sensed two-dimensional images. In this algorithm, a first guess of the motion, estimated using the Fourier phase shift, is used to predict the candidates for matching. A disparity score is computed for each target–candidate pair. The disparity also incorporates overlapping criteria in the case of large objects. Then the Hungarian method is applied to identify the best pairs by minimizing the global disparity. The high-disparity pairs are unmatched, and their target and candidate are declared expired and newly initiated objects, respectively. They are tested for merger and split on the basis of their size and overlap with the other objects. The sensitivity of track duration is shown for different disparity and size thresholds. The report highlights the algorithm’s ability to study convective life cycles using radar and simulated data over Darwin, Australia. The algorithm skillfully tracks individual convective cells (a few pixels in size) and large convective systems. The duration of tracks and cell size are found to be lognormally distributed over Darwin. The evolution of size and precipitation types of isolated convective cells is presented in the Lagrangian perspective. This algorithm is part of a vision for a modular platform [viz., TINT is not TITAN (TINT) and Tracking and Object-Based Analysis of Clouds (tobac)] that will evolve into a sustainable choice to analyze atmospheric features.

54 ENVIRONMENTAL SCIENCES↗

Thin film concentrator panel development

The development and testing of a rigid panel concept that utilizes a thin film reflective surface for application to a low-cost point-focusing solar concentrator is discussed. It is shown that a thin film reflective surface is acceptable for use on solar concentrators, including 1500 F applications. Additionally, it is shown that a formed steel sheet substrate is a good choice for concentrator panels. The panel has good optical properties, acceptable forming tolerances, environmentally resistant substrate and stiffeners, and adaptability to low to mass production rates. Computer simulations of the concentrator optics were run using the selected reflector panel design. Experimentally determined values for reflector surface specularity and reflectivity along with dimensional data were used in the analysis. The simulations provided intercept factor and net energy into the aperture as a function of aperture size for different surface errors and pointing errors. Point source and Sun source optical tests were also performed.

Zimmerman, D. K.↗

Frequency-domain full-waveform inversion with non-linear descent directions

Full-waveform inversion (FWI) is a highly non-linear inverse problem, normally solved iteratively, with each iteration involving an update constructed through linear operations on the residuals. Incorporating a flexible degree of non-linearity within each update may have important consequences for convergence rates, determination of low model wavenumbers and discrimination of parameters. In this study, we examine one approach for doing so, wherein higher order scattering terms are included within the sensitivity kernel during the construction of the descent direction, adjusting it away from that of the standard Gauss–Newton approach. These scattering terms are naturally admitted when we construct the sensitivity kernel by varying not the current but the to-be-updated model at each iteration. Linear and/or non-linear inverse scattering methodologies allow these additional sensitivity contributions to be computed from the current data residuals within any given update. We show that in the presence of pre-critical reflection data, the error in a second-order non-linear update to a background of s 0 is, in our scheme, proportional to at most (Δs/s 0 ) 3 in the actual parameter jump Δs causing the reflection. In contrast, the error in a standard Gauss–Newton FWI update is proportional to (Δs/s 0 ) 2 . For numerical implementation of more complex cases, we introduce a non-linear frequency-domain scheme, with an inner and an outer loop. A perturbation is determined from the data residuals within the inner loop, and a descent direction based on the resulting non-linear sensitivity kernel is computed in the outer loop. We examine the response of this non-linear FWI using acoustic single-parameter synthetics derived from the Marmousi model. Lastly, the inverted results vary depending on data frequency ranges and initial models, but we conclude that the non-linear FWI has the capability to generate high-resolution model estimates in both shallow and deep regions, and to converge rapidly, relative to a benchmark FWI approach involving the standard gradient.

58 GEOSCIENCES↗

Review of Wind–Wave Coupling Models for Large-Eddy Simulation of the Marine Atmospheric Boundary Layer

Abstract We present a review of existing wind–wave coupling models and parameterizations used for large-eddy simulation of the marine atmospheric boundary layer. The models are classified into two main categories: (i) the wave-phase-averaged, sea surface–roughness models and (ii) the wave-phase-resolved models. Both categories are discussed from their implementation, validity, and computational efficiency viewpoints, with emphasis given on their applicability in offshore wind energy problems. In addition to the various models discussed, a review of laboratory-scale and field-measurement databases is presented thereafter. The majority of the presented data have been gathered over many decades of studying air–sea interaction phenomena, with the most recent ones compiled to reflect an offshore wind energy perspective. Both provide valuable data for model validation. We also discuss the modeling knowledge gaps and computational challenges ahead.

54 ENVIRONMENTAL SCIENCES↗

Pressure-induced phase transition and electrical properties of thermoelectric Al-doped Mg{sub 2}Si

A recent study has shown the thermoelectric performance of Al-doped Mg{sub 2}Si materials can be significantly enhanced at moderate pressure. To understand the cause of this phenomenon, we have performed in situ angle dispersive X-ray diffraction and infrared reflectivity measurements up to 17 GPa at room temperature. Contrary to previous experiment, using helium as a pressure transmission medium, no structural transformation was observed in pure Mg{sub 2}Si. In contrast, a phase transition from cubic anti-fluorite (Fm-3m) to orthorhombic anti-cotunnite (Pnma) was observed in the Al-doped sample at 10 GPa. Infrared reflectivity measurements show the electrical conductivity increases with pressure and is further enhanced after the phase transition. The electron density of states at the Fermi level computed form density functional calculations predict a maximum thermoelectric power factor at 1.9 GPa, which is in good agreement with the experimental observation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Design Through Manufacturing: The Solid Model - Finite Element Analysis Interface

State-of-the-art computer aided design (CAD) presently affords engineers the opportunity to create solid models of machine parts which reflect every detail of the finished product. Ideally, these models should fulfill two very important functions: (1) they must provide numerical control information for automated manufacturing of precision parts, and (2) they must enable analysts to easily evaluate the stress levels (using finite element analysis - FEA) for all structurally significant parts used in space missions. Today's state-of-the-art CAD programs perform function (1) very well, providing an excellent model for precision manufacturing. But they do not provide a straightforward and simple means of automating the translation from CAD to FEA models, especially for aircraft-type structures. The research performed during the fellowship period investigated the transition process from the solid CAD model to the FEA stress analysis model with the final goal of creating an automatic interface between the two. During the period of the fellowship a detailed multi-year program for the development of such an interface was created. The ultimate goal of this program will be the development of a fully parameterized automatic ProE/FEA translator for parts and assemblies, with the incorporation of data base management into the solution, and ultimately including computational fluid dynamics and thermal modeling in the interface.

Rubin, Carol↗

Atmospheric Correction Inter-comparison eXercise, ACIX-II Land: An Assessment of Amospheric Correction Processors for Landsat 8 and Sentinel-2 Over Land

The correction of the atmospheric effects on optical satellite images is essential for quantitative and multi-temporal remote sensing applications. In order to study the performance of the state-of-the-art methods in an integrated way, a voluntary and open-access benchmark Atmospheric Correction Inter-comparison eXercise (ACIX) was initiated in 2016 in the frame of Committee on Earth Observation Satellites (CEOS) Working Group on Calibration & Validation (WGCV). The first exercise was extended in a second edition wherein twelve atmospheric correction (AC) processors, a substantially larger testing dataset and additional validation metrics were involved. The sites for the inter-comparison analysis were defined by investigating the full catalogue of the Aerosol Robotic Network (AERONET) sites for coincident measurements with satellites' overpass. Although there were more than one hundred sites for Copernicus Sentinel-2 and Landsat 8 acquisitions, the analysis presented in this paper concerns only the common matchups amongst all processors, reducing the number to 79 and 62 sites respectively. Aerosol Optical Depth (AOD) and Water Vapour (WV) retrievals were consequently validated based on the available AERONET observations. The processors mostly succeeded in retrieving AOD for relatively light to medium aerosol loading (AOD < 0.2) with uncertainties <0.08, while the overall uncertainty values were typically 0.23 ± 0.15. Better performances were observed for WV retrievals with >90% of the results falling within the suggested empirical specifications and with the Root Mean Square Error (RMSE) being mostly <0.25 g/cm2. Regarding Surface Reflectance (SR) validation two main approaches were followed. For the first one, a simulated SR reference dataset was computed over all of the test sites by using the 6SV (Second Simulation of the Satellite Signal in the Solar Spectrum vector code) full radiative transfer modelling (RTM) and AERONET measurements for the required aerosol variables and water vapour content. The performance assessment demonstrated that the retrievals were not biased for most of the bands. The uncertainties ranged from approximately 0.003 to 0.01 (excluding B01) for the best performing processors in both sensors' analyses. For the second one, measurements from the radiometric calibration network RadCalNet over La Crau (France) and Gobabeb (Namibia) were involved in the validation. The performance of the processors was in general consistent across all bands for both sensors and with low standard deviations (<0.04) between on-site and estimated surface reflectance. Overall, our study provides a good insight of AC algorithms' performance to developers and users, pointing out similarities and differences for AOD, WV and SR retrievals. Such validation though still lacks of ground-based measurements of known uncertainty to better assess and characterize the uncertainties in SR retrievals.

Atmospheric correction↗

Orbital simulations of satellite escape/capture and the origin of satellites such as Triton

We investigate satellite escape/capture in the context of the restricted, circular three body problem as applied to the Sun, Neptune, and Triton. We have computed a large number of coplanar prograde and retrograde orbital simulations over a range of initial distances and velocities. The satellite starts at superior conjunction within approximately 2 Hill radii of Neptune and has a velocity orthogonal to the Sun-planet line. Orbits with these initial conditions can be reflected with respect to time, so an escape is simply the reverse of a capture. We numerically integrate the equations of motion to compute the satellite's position until it escapes, collides with Neptune, or after 100 planetary years fails to escape, when computations cease. The initial distance x and velocity v in the restricted problem uniquely define the Jacobi constant C, a conserved energy-like quantity. Plots of the simulation outcomes in the prograde and retrograde C, x phase spaces reveal distinct zones in which temporary satellites approach the planet closely enough that permanent capture can be effected by gas drag with a protoplanetary nebula or by collision with a pre-existing satellite. Single and double close-flybys constitute the most common possible capture orbits. Long term multiple flyby orbits occur near the stability limits between bound and unbound orbits, and are more common among retrograde captures.

Benner, Lance A. M.↗