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 91 records · Page 5

Improving Schrödinger Equation Implementations with Gray Code for Adiabatic Quantum Computers

We reformulate the continuous-space Schrödinger equation in terms of spin Hamiltonians. For the kinetic energy operator, the critical concept facilitating the reduction in model complexity is the idea of position encoding. A binary encoding of position produces a spin-1/2 Heisenberg-like model and yields exponential improvement in space complexity when compared to classical computing. Encoding with a binary reflected Gray code (BRGC), and a Hamming-distance-2 Gray code (H2GC) reduces the model complexity down to the 𝑋⁢𝑍 and transverse Ising model, respectively. For 𝐴 qubits BRGC yields 2 𝐴 positions and is reduced to its 2-local form with O⁡(𝐴) ancillary qubits. H2GC yields 2 𝐴/2+1 positions with O⁡(𝐴 2 ) three-local penalty terms. We also identify the bijective mapping between diagonal unitaries and the Walsh series, producing the mapping of any real potential to a series of 𝑘 -local Ising models through the fast Walsh transform. Finally, in a finite volume, we provide some numerical evidence to support the claim that the total time needed for adiabatic evolution is protected by the infrared cutoff of the system. As a result, initial state preparation from a free-field wave function to an interacting system is expected to exhibit polynomial time complexity with volume and constant scaling with respect to lattice discretization for all encodings. For H2GC, if the evolution starts with the transverse Hamiltonian due to hardware restrictions, then penalties are dynamically introduced such that the low-lying spectrum reproduces the energy levels of the Laplacian. The adiabatic evolution of the penalty Hamiltonian is therefore sensitive to the ultraviolet scale. It is expected to exhibit polynomial time complexity with lattice discretization, or exponential time complexity with respect to the number of qubits given a fixed volume.

97 MATHEMATICS AND COMPUTING↗

Hyperspectral Aquatic Radiative Transfer Modeling Using a High-Performance Cluster Computing-Based Approach

Abstract For aquatic studies, radiative transfer (RT) modeling can be used to compute hyperspectral above-surface remote sensing reflectance that can be utilized for inverse model development. Inverse models can provide bathymetry and inherent-and bottom-optical property estimation. Because measured oceanic field/organic datasets are often spatio-temporally sparse, synthetic data generation is useful in yielding sufficiently large datasets for inversion model development; however, these forward-modeled data are computationally expensive and time-consuming to generate. This study establishes the magnitude of wall-clock-time savings achieved for performing large, aquatic RT batch-runs using parallel computing versus a sequential approach. Given 2,600 simulations and identical compute-node characteristics, sequential architecture required ~100 hours until termination, whereas a parallel approach required only ~2.5 hours (42 compute nodes) a 40x speed-up. Tools developed for this parallel execution are discussed.

58 GEOSCIENCES↗

Uncertainty of SW Cloud Radiative Effect in Atmospheric Models Due to the Parameterization of Liquid Cloud Optical Properties

Clouds are largely responsible for the spread of climate models predictions. Here we focus on the uncertainties in cloud shortwave radiative effect due to the parameterization of liquid cloud single scattering properties (SSPs) from liquid water content (LWC) and droplet number concentration (N), named parameterization of cloud optical properties. Uncertainties arise from not accounting for the droplet size distribution (DSD)—which affects the estimation of the effective radius (r eff ) and modulates the r eff -dependency of the SSPs—and from averaging SSPs over wide spectral bands. To assess these uncertainties a series of r eff -dependent SSPs parameterizations corresponding to various DSDs and spectral averaging methods are derived and implemented in a radiative code. Combined with the DSD-dependent estimation of r eff they are used to compute the bulk radiative properties (reflectance, transmittance, absorptance) of various clouds (defined in terms of LWC and N), including a homogeneous cloud, more realistic case studies, and outputs of a climate model. The results show that the cloud radiative forcing can vary up to 20% depending on the assumed DSD. Likewise, differences up to 20% are obtained for heating rates. The estimation of r eff is the main source of uncertainty, while the SSPs parameterization contributes to around 20% of the total uncertainty. Spectral averaging is less an issue, except for atmospheric absorption. Overall, global shortwave cloud radiative effect can vary by 6 W m –2 depending on the assumed DSD shape, which is about 13% of the best observational estimate.

54 ENVIRONMENTAL SCIENCES↗

Careers in STEM: A Latina Perspective

Three Latina computing professionals at a large national laboratory reflect on the circumstances affecting the low representation of this segment of the population in STEM fields, and computing in particular. Here, the authors share highlights of their path to STEM careers and some of the efforts they are involved in for broadening participation in computing. They consider the roles of minority serving institutions, representation and mentoring, and advocacy.

99 GENERAL AND MISCELLANEOUS↗

A real-time applicator position monitoring system for gynecologic intracavitary brachytherapy

Purpose: To develop a real-time applicator position monitoring system (RAPS) for intracavitary brachytherapy using an infrared camera and reflective markers. Methods: 3D image-guided brachytherapy requires high accuracy of applicator localization; however, applicator displacement can happen during patient transfer for imaging and treatment delivery. No continuous applicator position monitoring system is currently available. The RAPS system was developed for real-time applicator position monitoring without additional radiation dose to patients. It includes an infrared camera, reflective markers, an infrared illuminator, and image processing software. After reflective markers are firmly attached to the applicator and the patient body, applicator displacement can be measured by computing the relative change in distance between the markers. The reflective markers are magnetic resonance imaging (MRI) compatible, which is suitable for MRI-guided HDR brachytherapy paradigm. In our prototype, a Microsoft Kinect sensor with a resolution of 640 by 480 pixels is used as an infrared camera. A phantom study was carried out to compare RAPS' measurements with known displacements ranging from −15 to +15 mm. A reproducibility test was also conducted. Results: The RAPS can achieve 4 frames/s using a laptop with Intel{sup ®} Core™2 Duo processor. When the pixel size is 0.95 mm, the difference between RAPS' measurements and known shift values varied from 0 to 0.8 mm with the mean value of 0.1 mm and a standard deviation of 0.44 mm. The system reproducibility was within 0.6 mm after ten reposition trials. Conclusions: This work demonstrates the feasibility of a real-time infrared camera based gynecologic intracavitary brachytherapy applicator monitoring system. Less than 1 mm accuracy is achieved when using an off-the-shelf infrared camera.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Reflections on the Shifting Experiences of Scientific Infrastructure

Infrastructure of all types is fundamental to modern work and life. Computing for scientific work, especially, extends from distributed local research sites, often at the edges of other major systems, outward into globally connected high-performance facilities and infrastructures. This commentary reviews longstanding research on the social characteristics of infrastructure. We reflect on social concerns that affect the ongoing development, use, and maintenance of a wide range of scientific computing and data resources. Reflecting on the social nature of infrastructure is timely for Computing in Science & Engineering readers, given continued emphasis on developing even more expansive platforms for data and artificial intelligence work in science (e.g., the United States’ Genesis Mission). We assert that, regardless of technological advances, the complex nature of scientific research and data will require continued understanding of longstanding and nascent social practices across varied communities. This is fundamentally necessary to build and sustain usable infrastructure or platforms that can productively advance scientific research.

Paine, Drew [Lawrence Berkeley National Laboratory↗

Final Reports of the 2020 Los Alamos National Laboratory Computational Physics Student Summer Workshop

For the past ten years, the workshop has been bringing a highly talented and diverse group of student every summer. Students work in teams of two, alongside typically two mentors, on research projects reflecting a broad range of topics within computational physics. In addition, students attend a series of lectures on topics within computational physics, facility tours, and networking events. The program lasts ten weeks, with this year’s workshop running from June 8 to August 14. At the end of the summer, students give a final presentation, along with a written report. Those reports are what make up the remaining sections of this document. Admission to the workshop is by a competitive process, with the mentors forming the selection committee. One of the important accomplishments of the workshop has been to create a student pipeline from diverse schools that sometimes are not normally tapped by LANL recruiting. Many workshop students maintain a continuing relationship with LANL, returning as students interns, post-doctoral researchers, and staff members. Additionally, workshop alumni act as ambassadors for LANL. The result is a wider awareness both of LANL as a potential employer, and of the technical work that happens at LANL. This year, the workshop format was changed in several ways, in order to accommodate the off-site, virtual format. Students worked on LANL virtual desktop systems remotely, also accessing LANL HPC resources. In order to facilitate communication, student were given accounts on both Webex, a video teleconferencing platform, and Mattermost, an online team collaboration and chat platform, similar to Slack. Daily communication between students and mentors was primarily on Mattermost, with Webex conferencing as needed. The lectures were all done on Webex. Given the difficulty of the virtual format, and a concern that students might have video teleconferencing burn-out after an academic semester largely moved to that format, all lectures were optional this year. In spite of this, the attendance was generally high. Lecturers were asked to try to move to a more high-level, ”What is it?,” format. Once again, the students did a tremendous job. Over the course of ten weeks, they did important research across a staggering array of disciplines. The following pages contain the final report for each team’s research efforts. We hope you will find reading them as exciting as it was for us to produce them.

36 MATERIALS SCIENCE↗

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗

Robust quantum gates via a photon triggering electromagnetically induced transparency

We propose a scheme to achieve a parallelized controlled-NOT (c-not) gate based on electromagnetically induced transparency (EIT) and cavity-assisted photon scattering. By using the flying photon as the control bit, we can implement a conditional transfer between two logical states of an atom trapped in the cavity with high fidelity and on a microsecond time scale. Thanks to the fact that the photon is well suited for scalable quantum computation, a quantum c-not gate between the single-photon pulses and the nonlocal gate on remote atoms are obtained by reflecting the photon pulse from an optical cavity with a single-trapped atom. Our protocol is robust for spontaneous emission and works quite well in the bad-cavity limit, which makes it more applicable in the laboratory with current experimental techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A review of computational tools for design and reconstruction of metabolic pathways

Metabolic pathways reflect an organism's chemical repertoire and hence their elucidation and design have been a primary goal in metabolic engineering. Various computational methods have been developed to design novel metabolic pathways while taking into account several prerequisites such as pathway stoichiometry, thermodynamics, host compatibility, and enzyme availability. The choice of the method is often determined by the nature of the metabolites of interest and preferred host organism, along with computational complexity and availability of software tools. In this paper, we review different computational approaches used to design metabolic pathways based on the reaction network representation of the database (i.e., graph or stoichiometric matrix) and the search algorithm (i.e., graph search, flux balance analysis, or retrosynthetic search). We also put forth a systematic workflow that can be implemented in projects requiring pathway design and highlight current limitations and obstacles in computational pathway design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ANALYTIC MODELS FOR ALBEDOS, PHASE CURVES, AND POLARIZATION OF REFLECTED LIGHT FROM EXOPLANETS

New observational facilities are becoming increasingly capable of observing reflected light from transiting and directly imaged extrasolar planets. In this study, we provide an analytic framework to interpret observed phase curves, geometric albedos, and polarization of giant planet atmospheres. We compute the observables for non-conservative Rayleigh scattering in homogeneous semi-infinite atmospheres using both scalar and vector formalisms. In addition, we compute phase curves and albedos for Lambertian, isotropic, and anisotropic scattering phase functions. We provide analytic expressions for geometric albedos and spherical albedos as a function of the scattering albedo for Rayleigh scattering in semi-infinite atmospheres. Given an observed geometric albedo our prescriptions can be used to estimate the underlying scattering albedo of the atmosphere, which in turn is indicative of the scattering and absorptive properties of the atmosphere. We also study the dependence of polarization in Rayleigh scattering atmospheres on the orbital parameters of the planet-star system, particularly on the orbital inclination. We show how the orbital inclination of non-transiting exoplanets can be constrained from their observed polarization parameters. We consolidate the formalism, solution techniques, and results from analytic models available in the literature, often scattered in various sources, and present a systematic procedure to compute albedos, phase curves, and polarization of reflected light.

79 ASTRONOMY AND ASTROPHYSICS↗

Final Reports of the 2021 Los Alamos National Laboratory Computational Physics Student Summer Workshop

Since 2011, the Los Alamos National Laboratory Computational Physics Student Summer Workshop has been bringing together a highly talented and diverse group of students every summer. Students work in teams of two, alongside typically two mentors, on research projects reflecting a broad range of topics within computational physics. In addition, students attend a series of lectures on topics within computational physics, facility tours, and networking events. The program lasts ten weeks, with this year’s workshop running from June 7 to August 13. At the end of the summer, students give a final presentation, along with a written report. Those reports are what make up the remaining sections of this document. Admission to the workshop is by a competitive process, with the mentors forming the selection committee. One of the important accomplishments of the workshop has been to create a student pipeline from diverse schools that sometimes are not normally tapped by LANL recruiting. Many workshop students maintain a continuing relationship with LANL, returning as student interns, post-doctoral researchers, and staff members. Additionally, workshop alumni act as ambassadors for LANL. The result is a wider awareness both of LANL as a potential employer, and of the technical work that happens at LANL. This year, the workshop was once again in an off-site, virtual format. Students worked on LANL virtual desktop systems remotely, also accessing LANL HPC resources. In order to facilitate communication, student were given accounts on both Webex, a video teleconferencing platform, and Mattermost, an online team collaboration and chat platform, similar to Slack. Daily communication between students and mentors was primarily on Mattermost, with Webex conferencing as needed. The lectures were all done on Webex. Given the difficulty of the virtual format, and a concern that students might have video teleconferencing burn-out, all lectures were optional this year. In spite of this, the attendance was generally high. Lecturers were asked to try to move to a more high-level, ”What is it?,” format. Once again, the students did a tremendous job. Over the course of ten weeks, they did important research across a staggering array of disciplines. The following pages contain the final report for each team’s research efforts. Enjoy!

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Rapid Reflective Facet Characterization Using Fringe Reflection Techniques

Reflective mirror facets for concentrating solar power (CSP) systems have stringent requirements on the surface slope accuracy in order to provide adequate system performance. This paper presents a tool that can fully characterize facets quickly enough for 100% inspection on a production line. A facet for a CSP system, specifically a dish concentrator has a parabolic design shape. This shape will concentrate near-parallel rays from the sun to a point (or a line for trough systems). Deviations of surface slope from the design shape impact the performance of the system, either losing power that misses the target or increasing peak fluxes to undesirable levels. During development or production, accurate knowledge of facet defects can lead to improvements to lower cost or improve performance. The reported characterization system, SOFAST (Sandia Optical Fringe Analysis Slope Tool), has a computer-connected camera that images the reflective surface, which is positioned so that it reflects an active target, such as an LCD screen, to the camera. A series of fringe patterns are displayed on the screen while images are captured. Using the captured information, the reflected target location of each pixel of mirror viewed can be determined, and thus through a mathematical transformation, a surface normal map can be developed. This is then fitted to the selected model equation, and the errors from design are characterized. While similar approaches have been explored, several key developments are presented here. The combination of the display, capture, and data reduction in one system allows rapid characterization. An “electronic boresight” approach is utilized to accommodate physical equipment positioning deviations, making the system insensitive to setup errors. Up to 1.5×10 6 points are characterized on each facet. Finally, while prior automotive industry commercial systems resolve the data to shape determination, SOFAST concentrates on slope characterization and reporting, which is tailored to solar applications. SOFAST can be used for facet analysis during development. However, the real payoff is in production, where complete analysis is performed in about 10s. With optimized coding, this could be further reduced.

14 SOLAR ENERGY↗

Classification of Dissolution Events Using Fusion of Effluents Measurements and Classifiers

Classifiers for dissolution events at a radiochemical processing facility are studied using gamma spectra measurements of effluents collected by a high purity germanium detector located at its off-gas stack. Data sets collected at the Oak Ridge National Laboratory’s Radiochemical Engineering Development Center under a Pu dissolution campaign spanning a three months period are utilized. Features corresponding to the activity levels of 15 radionuclides, including isotopes of iodine, krypton, and xenon, that are indicated by the target decay chains, are computed from the spectra at 1 hour intervals. A conceptualization diagram is developed to reflect the steps from the source to measurement to feature computation that depend on fission products indicated by decay chains, chemical processing, and effluents transport to the off-gas stack. A diverse set of eight classifiers based on different design principles are trained using the ground truth data for this campaign, and the outputs of top three classifiers, namely, classification trees, Ensemble of Trees (EOT), and k-nearest neighbor, are combined using EOT classifier-fuser. Our results show that for 5-fold cross validation, features associated with isotopes of xenon provide the lowest classification error among the different elements across the classifiers; the classification error is furthered improved when all 15 isotope features are used by each classifier, and it is again improved by the fusion of three classifiers. Further reduction in classification error is achieved by using a measurement window of 1-2 days which is identified based on half-life time estimates of the isotopes; it is long enough for the stabilization of feature estimates while being short enough not to be affected by the follow on dissolution events. As a net result of feature and classifier fusion, combined with the incorporation of decay chain and isotope half-life information, this approach achieves 98% detection rate while maintaining a false alarm rate under 2% for this data set.

Rao, Nageswara↗

Estimating the implied cost of carbon in future scenarios using a CGE model: The Case of Colorado

We develop a state-level computable general equilibrium (CGE) model that reflects the roles of coal, natural gas, wind, solar, and hydroelectricity in supplying electricity, using Colorado as a case study. Also, we focus on the economic impact of implementing Colorado's existing Renewable Portfolio Standard, updated in 2013. This requires that 25% of state generation come from qualifying renewable sources by 2020. We evaluate the policy under a variety of assumptions regarding wind integration costs and assumptions on the persistence of federal subsidies for wind. Specifically, we estimate the implied price of carbon as the carbon price at which a state-level policy would pass a state-level cost-benefit analysis, taking account of estimated greenhouse gas emission reductions and ancillary benefits from corresponding reductions in criteria pollutants. Our findings suggest that without the Production Tax Credit (federal aid), the state policy of mandating renewable power generation (RPS) is costly to state actors, with an implied cost of carbon of about $17 per ton of CO 2 with a 3% discount rate. Federal aid makes the decision between natural gas and wind nearly cost neutral for Colorado.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Insight into the Fast-Rechargeability of a Novel Mo 1.5 W 1.5 Nb 14 O 44 Anode Material for High-Performance Lithium-Ion Batteries

Wadsley–Roth phased niobates are promising anode materials for lithium-ion batteries, while their inherently low electrical conductivity still limits their rate-capability. Herein, a novel doped Mo 1.5 W 1.5 Nb 14 O 44 (MWNO) material is facilely prepared via an ionothermal-synthesis-assisted doping strategy. The detailed crystal structure of MWNO is characterized by neutron powder diffraction and aberration corrected scanning transmission electron microscope, unveiling the full occupation of Mo 6+ -dopant at the t1 tetrahedral site. In half-cells, MWNO exhibits enhanced fast-rechargeability. In this work, the origin of the improved performance is investigated by ultraviolet–visible diffuse reflectance spectroscopy, density functional theory (DFT) computation, and electrochemical impedance spectroscopy, revealing that bandgap narrowing improves the electrical conductivity of MWNO. Furthermore, operando X-ray diffraction elucidates that MWNO exhibits a typical solid-solution phase conversion-based lithium-ion insertion/extraction mechanism with reversible structural evolution during the electrochemical reaction. The boosted lithium-ion diffusivity of MWNO, due to the Mo 6+ /W 6+ doping effect, is confirmed by a galvanostatic intermittent titration technique and DFT. With the simultaneously enhanced electrical conductivity and lithium-ion diffusivity, MWNO successfully demonstrates its fast-rechargeability and practicality in the LiNi 0.5 Mn 1.5 O 4 -coupled full-cells. Therefore, this work illustrates the potential of ionothermal synthesis in energy storage materials and provides a mechanistic understanding of the doping effect on improving material's electrochemical performance.

25 ENERGY STORAGE↗

Anomalous lattice expansion in yttria stabilized zirconia under simultaneous applied electric and thermal fields: A time-resolved in situ energy dispersive x-ray diffractometry study with an ultrahigh energy synchrotron probe

Nonisothermal densification in 8% yttria doped zirconia (8YSZ) particulate matter of 250 nm median particle size was studied under 215 V/cm dc electric field and 9 Degree-Sign C/min heating rate, using time-resolved in-situ high temperature energy dispersive x-ray diffractometry with a polychromatic 200 keV synchrotron probe. Densification occurred in the 876-905 Degree-Sign C range, which resulted in 97% of the theoretical density. No local melting at particle-particle contacts was observed in scanning electron micrographs, implying densification was due to solid state mass transport processes. The maximum current draw at 905 Degree-Sign C was 3 A, corresponding to instantaneous absorbed power density of 570 W/cm{sup 3}. Densification of 8YSZ was accompanied by anomalous elastic volume expansions of the unit cell by 0.45% and 2.80% at 847 Degree-Sign C and 905 Degree-Sign C, respectively. The anomalous expansion at 905 Degree-Sign C at which maximum densification was observed is characterized by three stages: (I) linear stage, (II) anomalous stage, and (III) anelastic recovery stage. The densification in stage I (184 s) and II (15 s) was completed in 199 s, while anelastic relaxation in stage III lasted 130 s. The residual strains ({epsilon}) at room temperature, as computed from tetragonal (112) and (211) reflections, are {epsilon}{sub (112)} = 0.05% and {epsilon}{sub (211)} = 0.13%, respectively. Time dependence of (211) and (112) peak widths ({beta}) show a decrease with both exhibiting a singularity at 905 Degree-Sign C. An anisotropy in (112) and (211) peak widths of {l_brace} {beta}{sub (112)}/{beta}{sub (211)}{r_brace} = (3:1) magnitude was observed. No phase transformation occurred at 905 Degree-Sign C as verified from diffraction spectra on both sides of the singularity, i.e., the unit cell symmetry remains tetragonal. We attribute the reduction in densification temperature and time to ultrafast ambipolar diffusion of species arising from the superposition of mass fluxes due to Fickian diffusion, thermodiffusion (Soret effect), and electromigration, which in turn are a consequence of a superposition of chemical, temperature, and electrical potential gradients. On the other hand, we propose defect pile-up at particle-particle contacts and subsequent tunneling as a mechanism creating the 'burst-mode' discontinuous densification at the singularities observed at 847 and 905 Degree-Sign C.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗