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 289 records · Page 16

Final Report for Industrial Base Cybersecurity Initiative (IBCI) Project (FY2019)

Information security issues impact all manufacturing sectors and are particularly troublesome when designs for advanced weapons systems are lost via attacks that propagate through the relatively unprotected lower tiers of the nation's supply chains. The US Manufacturing Industrial Base (MIB) consists of suppliers and distributors that range in size from large prime contractors with extensive cybersecurity resources to small shops with little or no in-house information technology capabilities. The acquisition and maintenance of the technologies, systems, and financial resources needed to assure secure collaborations can be an overwhelming challenge for many manufacturers yet the cost of poor cybersecurity practices can lead to the compromise of national security information worth hundreds of millions of dollars. The Department of Defense Manufacturing Industrial Base Policy organization and the Department of Energy, Energy Efficiency and Renewable Energy Advanced Manufacturing Office have chartered the Los Alamos National Laboratory (LANL) and the Oak Ridge Y-12 National Security Complex (Y-12) to design, develop, and demonstrate a Manufacturing Threat Evaluation Wargaming System (MTEWS). The MTEWS concept is based on the use of a virtual model of the computer networks and cybersecurity protocols used across the MIB that will be continuously updated to reflect the dynamic threat environment. The MTEWS cyber-physical test range will demonstrate a prototype wargaming system that analyzes the impact of new and existing cybersecurity threats while also providing the information needed to sustain the integrity of the MIB supply chains. A long-term goal is the establishment of an enduring national manufacturing operations center that proactively responds to the continuously changing attacks that threaten the MIB and US national security. LANL and Y-12 have formed an Industrial Base Cybersecurity Initiative Team (IBCI), consisting of partners from the public, private and academic sectors, that is developing the MTEWS. To date, most of the project tasks have been centered on the technologies required to create a representative model of the MIB because this is a foundational element for the entire project. Modeling activities at the University of North Carolina at Charlotte (UNC Charlotte) have demonstrated a VMware-based prototype system that can accommodate a network model with up to 350 nodes and includes features such as the option to incorporate network scanning information when building the model and the ability to down select the elements of duplicative networks to a statistically significant subset that avoids redundantly duplicating network configurations. Broad Ordering Agreements have been established with UNC Charlotte and the University of West Florida (UWF) that support a variety of cybersecurity for manufacturing tasks and wide-spread endorsement has been received from government, industry and academia sectors concerning the IBCI goals and technical direction.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A Novel Sensor Platform Matching the Improved Version of IPMVP Option C for Measuring Energy Savings

It is easy to measure energy consumption with a power meter. However, energy savings cannot be directly computed by the powers measured using existing power meter technologies, since the power consumption only reflects parts of the real energy flows. The International Performance Measurement and Verification Protocol (IPMVP) was proposed by the Efficiency Valuation Organization (EVO) to quantify energy savings using four different methodologies of A, B, C and D. Although energy savings can be estimated following the IPMVP, there are limitations on its practical implementation. Moreover, the data processing methods of the four IPMVP alternatives use multiple sensors (thermometer, hygrometer, Occupant information) and power meter readings to simulate all facilities, in order to determine an energy usage benchmark and the energy savings. This study proposes a simple sensor platform to measure energy savings. Using usually the Electronic Product Code (EPC) global standard, an architecture framework for an information system is constructed that integrates sensors data, power meter readings and occupancy conditions. The proposed sensor platform is used to monitor a building with a newly built vertical garden system (VGS). A VGS shields solar radiation and saves on energy that would be expended on air-conditioning. With this platform, the amount of energy saved in the whole facility is measured and reported in real-time. The data are compared with those obtained from detailed measurement and verification (M&V) processes. The discrepancy is less than 1.565%. Using measurements from the proposed sensor platform, the energy savings for the entire facility are quantified, with a resolution of ±1.2%. The VGS gives an 8.483% daily electricity saving for the building. Thus, the results show that the simple sensor platform proposed by this study is more widely applicable than the four complicated IPMVP alternatives and the VGS is an effective tool in reducing the carbon footprint of a building.

vertical garden system (VGS)↗

CACTI CSAPR2 Taranis Retrievals

Taranis is an end-to-end processing chain for radar data written in Python with C extension for computation performance. Features include: masking for quality control, specific differential phase (Kdp), attenuation correction for reflectivity factor (Z) and differential reflectivity (Zdr) in rain, and additional geophysical retrievals. Retrievals are mostly drawn from literature or open-source software when appropriate, and have been tested, tuned, and modified to work with one another cohesively rather than using isolated off-the-shelf algorithms. Incorporated algorithms include hydrometeor (echo) identification, rain water content, raindrop mass-weighted mean diameter (gamma size distribution assumption), and rainfall rate (QPE). Taranis data sets exist for CSAPR2 PPI, HSRHI, and sector RHI scans. Cartesian-gridded data sets were also produced as well as a near-surface rain rate retrieval. More details can be found in the README.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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↗

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↗

A block-tridiagonal solver with two-level parallelization for finite element-spectral codes

Two-level parallelization is introduced to solve a massive block-tridiagonal matrix system. One-level is used for distributing blocks whose size is as large as the number of block rows due to the spectral basis, and the other level is used for parallelizing in the block row dimension. The reason for the added parallelization dimension is to retard the saturation of the scaling due to communication overhead and inefficiencies in the single-level parallelization only distributing blocks. As a technique for parallelizing the tridiagonal matrix, the combined method of “Partitioned Thomas method” and “Cyclic Odd–Even Reduction” is implemented in an MPI-Fortran90 based finite element-spectral code (TORIC) that calculates the propagation of electromagnetic waves in a tokamak. The two-level parallel solver using thousands of processors shows more than 5 times improved computation speed with the optimized processor grid compared to the single-level parallel solver under the same conditions. Three-dimensional RF field reconstructions in a tokamak are reflected as examples of the physics simulations that have been enabled by this algorithmic advance.

97 MATHEMATICS AND COMPUTING↗

Reflection and transmission of electromagnetic pulses at a planar dielectric interface -- theory and quantum lattice simulations

There is considerable interest in the application of quantum information science to advance computations in plasma physics. A particular point of curiosity is whether it is possible to take advantage of quantum computers to speed up numerical simulations relative to conventional computers. Many of the topics in fusion plasma physics are classical in nature. In order to implement them on quantum computers it will require couching a classical problem in the language of quantum mechanics. Electromagnetic waves are routinely used in fusion experiments to heat a plasma or to generate currents in the plasma. The propagation of electromagnetic waves is described by Maxwell equations with an appropriate description of the plasma as a dielectric medium. Before advancing to the tensor dielectric of a magnetized plasma, this paper considers electromagnetic wave propagation in a one-dimensional inhomogeneous scalar dielectric. The classic theory of scattering of plane electromagnetic waves at a planar interface, separating two different dielectric media, leads to Fresnel equations for reflection and transmission coefficients. In contrast to plane waves, this paper is on the reflection and transmission of a spatially confined electromagnetic pulse. Following an analytical formulation for the scattering of a Gaussian pulse, it is deduced that the maximum transmission coefficient for a pulse is $\sqrt{n_2/n_1}$ times that for a plane wave; the incident and transmitted pulses propagating in dielectric media with refractive indices $n_1$ and $n_2$, respectively. The analytical theory is complemented by numerical simulations using a quantum lattice algorithm for Maxwell equations. The algorithm, based on the Riemann-Silberstein-Weber representation of the electromagnetic fields and expressed in term of qubits, is an interleaved sequence of entangling operators at each lattice site and unitary streaming operators which transmit information from one site to an adjacent lattice site. Besides substantiating results from the theory for Gaussian pulses, numerical simulations show their validity for non-Gaussian pulses. Apart from their time-asymptotic forms, the simulations display an interplay between the incident, reflected, and transmitted pulses in the vicinity of the transition region between two dielectric media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of laser shock experiments on precompressed samples using a quartz reference and application to warm dense hydrogen and helium

Megabar (1 Mbar = 100 GPa) laser shocks on precompressed samples allow reaching unprecedented high densities and moderately high ~10 3 –10 4 K temperatures. We describe in this paper a complete analysis framework for the velocimetry (VISAR) and pyrometry (SOP) data produced in these experiments. Since the precompression increases the initial density of both the sample of interest and the quartz reference for pressure-density, reflectivity, and temperature measurements, we describe analytical corrections based on available experimental data on warm dense silica and density-functional-theory based molecular dynamics computer simulations. Finally, using our improved analysis framework, we report a re-analysis of previously published data on warm dense hydrogen and helium, compare the newly inferred pressure, density, and temperature data with most advanced equation of state models and provide updated reflectivity values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computation of glint, glare, and solar irradiance distribution

Described herein are technologies pertaining to computing the solar irradiance distribution on a surface of a receiver in a concentrating solar power system or glint/glare emitted from a reflective entity. At least one camera captures images of the Sun and the entity of interest, wherein the images have pluralities of pixels having respective pluralities of intensity values. Based upon the intensity values of the pixels in the respective images, the solar irradiance distribution on the surface of the entity or glint/glare corresponding to the entity is computed.

14 SOLAR ENERGY↗

Computation of glint, glare, and solar irradiance distribution

Described herein are technologies pertaining to computing the solar irradiance distribution on a surface of a receiver in a concentrating solar power system or glint/glare emitted from a reflective entity. At least one camera captures images of the Sun and the entity of interest, wherein the images have pluralities of pixels having respective pluralities of intensity values. Based upon the intensity values of the pixels in the respective images, the solar irradiance distribution on the surface of the entity or glint/glare corresponding to the entity is computed.

14 SOLAR ENERGY↗