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 757 records · Page 42

Cross-Measurement Comparisons for a CFD Validation Dataset on Mach 2.5 Axisymmetric Turbulent Shock-Wave/Boundary-Layer Interactions

Experimental data for a shock-wave/boundary-layer interaction has been collected using multiple measurement techniques. Unfortunately, diversity of methods for acquisition begets an aggregation of data which does not directly quantify the same properties of the flowfield. The objective in this paper was to (a) present the collection of flowfield measurements in one venue in a format more usable for those validating CFD models against the data, (b) evaluate the degree to which the experimental data support each other, and (c) highlight the differences and relative advantages/shortcomings of each measurement techniques. To present the various measured quantities into common format, CFD results from a companion paper also submitted for presentation at this meeting are utilized. For the case where the boundary layer remains attached, there is agreement between the various measurements as well as with Reynolds-averaged Navier-Stokes simulation solutions. As the impinging shock strength is increased beyond the point of separating the boundary layer, the congruity of the data wanes. Generally, agreement among the measurements exceeds the degree to which the CFD solutions agree with experiments. This suggests that unmodeled physical phenomena, such as transient motion of the reflected shock and separation bubble, give rise to the discrepancies observed in the computational results.

Shock-wave Boundary-layer interaction↗

Neutronics Analysis for ITER Diagnostic Generic Upper Port Plug

This paper presents new results of neutronics analysis performed for the updated design of the ITER Diagnostic Generic Upper Port Plug (DGUPP). Such analysis is performed regularly during the development of the ITER Ports design. Neutronics analysis is an important step in support of the development of ITER components irradiated by 14 MeV neutrons generated in D-T fusion thermonuclear reactions. The neutrons are born in the ITER tokamak plasma chamber and then impinge on the surrounding blanket and port plugs, after that they diffuse through shielding material of vacuum vessel components or could stream through the diagnostic channels and assembly gaps far away from plasma up to the bio-shield and in a small amount beyond it. While ITER machine is in operation it is inaccessible, after the shut-down it is planned to provide personnel access to the Upper Port Inter-Space Structure (UP ISS) behind the Port Plug via the corridors. Neutrons activate the ITER materials and those radioactive materials emit decay gammas during the shut-down causing a Shut-Down Dose Rate (SDDR) field distributed all around the activated materials. Decay gammas could also stream through the void channels. That emphasizes the need for calculation of decay gamma transport in evaluations of SDDR, which is reflected in the use of transport - activation interface code described in Section of Computational Modelling. The assessment of radiation environment (fluxes and SDDR) at DGUPP was focused on the port back-side area of UP ISS, where personnel access is planned along the ISS maintenance corridors. Minimization of SDDR at the front part of the ISS maintenance corridors was the target value for optimization of the DGUPP shielding and analysis of related neutronics effects. Following the dose ALARA (As Low As Reasonably Achievable) principle, every possibility to reduce SDDR should be investigated, emphasizing the importance of this analysis. The challenges of such analysis are inextricably connected with the engineering complexity of the ITER machine, with integration of many diagnostics systems having channels and gaps inside the ITER ports. Those gaps result in radiation streaming and negatively affect the shielding performance of the port plugs. The computational complication is caused by the need to model precisely the ITER geometry and then to run deep penetrating radiation transport in the substantially heterogenic toroidal geometry of ITER large-size tokamak, with 30 m height and diameter of its bio-shield. The whole mass of Upper Port Plug (UPP) could be up to 25 tons and it must be very precisely positioned inside the port opening in blanket and Vacuum Vessel. Facing the engineering challenges for the assembly of DGUPP, there were proposed two Project Change Requests (PCRs) No. 439 and 709 concerning the design updates of UPP/blanket interface and its remote handling.

61 RADIATION PROTECTION AND DOSIMETRY↗

AVIRIS Reflectance Retrievals: UCSB Users Manual

The following write-up is designed to help students and researchers take Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) radiance data and retrieve surface reflectance. In the event that the software is not available, but a user has access to a reflectance product, this document is designed to provide a better understanding of how AVIRIS reflectance was retrieved. This guide assumes that the reader has both a basic understanding of the UNIX computing environment, and that of spectroscopy. Knowledge of the Interactive Data Language (IDL) and the Environment for Visualizing Images (ENVI) is helpful. This is a working document, and many of the fine details described in the following pages have been previously undocumented. After having read this document the reader should be able to process AVIRIS to reflectance, provided access to all of the code is possible. The AVIRIS radiance data itself is pre-processed at the Jet Propulsion Laboratory (JPL) in Pasadena, California. The first section of this paper describes how to read data from tape and byte-swap the data. Section 2 describes the procedure in preparing support files before running the 'h2o' suite of programs. Section 3 describes the four programs used in the process, h2olut9.f, h2ospl9.f, vlsfit9.f and rfl9.f.

Roberts, Dar A.↗

A study of the reflection effect in close spectroscopic binary stars

Four W Ursae Majoris-type binaries were investigated using the simple model of blackbody solid spheres. Numerical computations were carried out for the rotationally broadened profiles of spectral lines Sr II (4077.7) and Ca I (4226.7) and the mean radial velocities. These mean velocities at different phases outside eclipses were compared with sinusoidal values of the radial velocity curve for point masses. Within the present accuracy of observations, the reflection effect is small and may be neglected in the analysis of radial velocity curves.

Chen, K.-Y.↗

A computationally efficient algorithm for computing convex hull prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may be required to execute side payments to participants whose costs are not covered through energy sales as determined via traditional locational marginal pricing schemes. Convex hull pricing minimizes this and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. Computing convex hull prices involves solving either a large-scale linear program or the Lagrangian dual of the corresponding non-convex scheduling problem. Further, the former approach requires explicit descriptions of market participants' convex hulls. While linear programs for computing convex hull prices are large, their structure is naturally decomposable by generators. Here, in this work, we propose and empirically analyze a Benders decomposition approach to computing convex hull prices that leverages recent advances in convex hull formulations for thermal generating units. We demonstrate across a large set of test instances that our decomposition approach only requires modest computational effort, obtaining solutions at least an order of magnitude faster than the equivalent large-scale linear programming approach. Overall, we provide a computationally feasible method for computing convex hull prices for industrial scale market clearing problems, enabling the possibility of practical adoption of this advanced pricing mechanism.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Convection in a two-layer fluid system

Experimental results are presented, and preliminary computations are performed on a system of two immiscible liquid layers with a temperature gradient applied parallel to the interface. The experiments reflect the combined contribution of buoyancy and surface-tension-induced (Marangoni) convection. It is concluded that buoyancy effects appear to be dominant and mask any surface-tension-induced convection present. Numerical computations show significant modification of pure buoyant convection by surface-tension gradients. The results are of interest in connection with the liquid encapsulation of GaAs melts in a microgravity environment.

Prakash, A.↗

The Dark Target aerosol retrieval algorithm applied to Low Earth Orbit and GEOstationary imagers: progress towards an integrated LEO-GEO view of global aerosol

The relatively simple dark-target (DT) aerosol retrieval algorithm provides products of spectral aerosol optical depth (AOD) from measurements of multi-spectral reflectance in visible, near-infrared and shortwave infrared wavelength bands. Originally developed for Moderate-resolution Imaging Spectroradiometer (MODIS aboard Terra and Aqua) in Low-Earth Orbit (LEO), DT has been ported to Visible Infrared Imaging Suite (VIIRS aboard Suomi-NPP and NOAA-20, also in LEO), to enhanced-MODIS Airborne Simulator (eMAS, on an airborne platform), and now to sensors in GEOstationary orbit (Advanced Himawari Imager - AHI aboard Himawari-8 and Advanced Baseline Imagers – ABI aboard GOES-16 and 17). Together, these new datasets not only extend upon the 20+ year MODIS aerosol record, but also expand the temporal sampling and/or spatial resolution. Between July and October of 2019, NASA participated in two field experiments on opposite sides of the globe. These included FIREX-AQ which focused on fire and smoke in the Western U.S., and then CAMP2EX which targeted aerosol/cloud interactions around the Philippines. We have performed DT aerosol retrievals on all images from all sensors during these three months, validated against ground observations from stationary and mobile sunphotometer sites, and have begun to develop a synergy that represents semi-global observations every half hour. The resulting aerosol products are being used as context and for model assimilation, thus providing the framework for more complete characterization of global aerosol transport and lifecycle. Here, we report on progress, as well as remaining challenges such as data management, computer processing, and accounting for differences between GEO and LEO observation geometry and surface reflectance parameterization.

dark target↗

Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals

The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities vL, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants c16 or c26 - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry which is two orders of magnitude faster than the previous method.

36 MATERIALS SCIENCE↗

Radiation transfer in plant canopies - Scattering of solar radiation and canopy reflectance

The one-dimensional vertical model of radiation transfer in a plant canopy described by Verstraete (1987) is extended to account for the transfer of diffuse radiation. This improved model computes the absorption and scattering of both visible and near-infrared radiation in a multilayer canopy as a function of solar position and leaf orientation distribution. Multiple scattering is allowed, and the spectral reflectance of the vegetation stand is predicted. The results of the model are compared to those of other models and actual observations.

Verstraete, Michel M.↗

Extent of hydrogen coverage of Si(001) under chemical vapor deposition conditions from ab initio approaches

The extent of hydrogen coverage of the Si(001) c(4 × 2) surface in the presence of hydrogen gas has been studied with dispersion corrected density functional theory. Electronic energy contributions are well described using a hybrid functional. The temperature dependence of the coverage in thermodynamic equilibrium was studied computing the phonon spectrum in a supercell approach. As an approximation to these demanding computations, an interpolated phonon approach was found to give comparable accuracy. The simpler ab initio thermodynamic approach is not accurate enough for the system studied, even if corrections by the Einstein model for surface vibrations are considered. The on-set of H{sub 2} desorption from the fully hydrogenated surface is predicted to occur at temperatures around 750 K. Strong changes in hydrogen coverage are found between 1000 and 1200 K in good agreement with previous reflectance anisotropy spectroscopy experiments. These findings allow a rational choice for the surface state in the computational treatment of chemical reactions under typical metal organic vapor phase epitaxy conditions on Si(001).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

StressNet - Deep learning to predict stress with fracture propagation in brittle materials

Abstract Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a deep learning model, StressNet, is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 h, with an average MAPE of 2% relative to test data.

36 MATERIALS SCIENCE↗

Implementing a Laser Stabilization System for Trapping Ca+ Ions: an Internship Reflection

At Lawrence Livermore National Laboratory, I contributed to a project developing 3D printed micro ion traps for quantum computing. I designed, implemented, and assessed a laser stabilization system that locked lasers to the frequencies required for calibrating our High Finesse WS8-10 wavelength meter and for laser cooling and trapping of Ca+ ions. I also programmed a Python interface for hardware communication, data collection, and statistical analysis. Additionally, I optimized and aligned laser beam paths, and I implemented a closed digital feedback loop using Proportional, Integral, and Derivative (PID) control parameters. I analyzed both the long-term and short-term behavior of our locked lasers and adjusted PID parameters to enhance performance. Furthermore, I used COMSOL to simulate the capacitance of a linear Paul trap design and predict our trap’s performance. The procedures I developed for the interface, analysis, and simulations will continue to support the ion trapping experiment after my appointment. I strengthened my skills in data analysis, Python coding, and optical alignment for laser systems. My confidence as a researcher grew, particularly in communicating my research. This experience taught me the importance of careful planning and consideration in research and solidified my desire to continue exploring novel quantum technology as an undergraduate

42 ENGINEERING↗

Jupiter's atmospheric composition and cloud structure deduced from absorption bands in reflected sunlight

The spectrum of sunlight reflected by Jupiter is analyzed by comparing observations of Woodman (1979) with multiple-scattering computations. The analysis yields information on the vertical cloud structure at several latitudes and on the abundance of CH4 and NH3 in the atmosphere of Jupiter. The abundances of CH4 and NH3 suggest that all ices and rocks are overabundant on Jupiter by a factor of 2 or more, providing an important constraint on models for the formation of Jupiter from the primitive solar nebula. The pressure level of the clouds, the gaseous NH3 abundance, the mean temperature profile, and the Clausius-Clapeyron relation suggest that these clouds are predominantly ammonia crystals with the cloud bottom at 600-700 mb. A diffuse distribution of aerosols exists between 150 and 500 mb, and the spectral variation of albedo reflects a changing bulk absorption coefficient of the material composing the aerosols and is diagnostic of the aerosol composition.

Sato, M.↗

The Odd Power of Dispersion

In ancient China, water has been regarded as one of the five vital components of life. It has been observed that water has many fascinating properties: water is ‘soft’ yet it can penetrate a hard rock; water is ‘pure’ yet it can tolerate other beings. Because of its unique properties, water is often associated with good quality and has been given the highest praise by Laozi in his book Tao Te Ching saying: the highest/best quality that one can have is being like water. However, little did people understand why and how water possesses such fascinating properties. Modern scientific developments made people realize that the macroscopic liquid water is made of a large number of water molecules held together via a network of hydrogen bonds. And those wonderful properties of water are merely the macroscopic manifestations of the interactions between water molecules and other molecules. For example, the dissolving ability of water is due to the fact that the interaction between a water molecule and the other molecular species is stronger than the interactions among their own molecular species. In fact the interactions between any two molecules are governed by the same physics and are termed intermolecular interaction (or intermolecular forces in some literature, although technically ‘force’ is incorrect usage here). Although the very existence of the intermolecular interactions is easily proved, e.g. the mere presence of the solid phase of matter, and scientists today have recognized that the seemingly weak intermolecular interactions essentially hold the world together through a delicate and cooperative process, the theoretical understanding of various intermolecular interactions is still far from satisfactory. On the practical side, theoreticians need to balance computational cost and accuracy. Because of the relatively small magnitudes of the intermolecular interactions, errors that appear tiny compared to the usual chemical (covalent) bonding may change conclusions qualitatively. High-level ab initio methods including explicit description of electron correlation can achieve the desired accuracy at very high computational cost. (Chapter 5 and 6) However the cooperative network of hundreds of thousands of molecules that reflects the true power of intermolecular interactions cannot be modeled easily by ab initio methods. Deeper understanding of intermolecular interactions yields better theoretical models; better theoretical models facilitate and even deepen the understanding of intermolecular interactions. With the aforementioned motivation in mind, a significant portion of this dissertation is dedicated to developing a method to describe the intermolecular interactions accurately with affordable computational resources.

74 ATOMIC AND MOLECULAR PHYSICS↗

Interrelationships between information and energy using knowledge management tools

Edward Fredkin was an enthusiastic advocate of information-based theoretical physics, who, in the early 1980s, proposed a new theory of physics based on the idea that the universe is ultimately composed of software. According to Fredkin, reality should be considered as being composed not of particles, matter and forces or energy but of bits of data or information modified according to computational rules. Fredkin went on to demonstrate that, while energy is necessary for storing and retrieving information, it can be arbitrarily reduced in order to carry out any particular instance of information processing, and this operation does not have a lower bound. This implies that it is information rather than matter or energy that should be considered at the ultimate fundamental constituent of reality. This possibility had already been suggested by other scientists. Norbert Wiener heralded a fundamental shift from energy to information and suggested that the universe was founded essentially on the transformation of information, not energy. However, Konrad Zuse was the first, back in 1967, to defend the idea that a digital computer is computing the universe. Richard P. Feynman showed this possibility in a similar light in his reflections on how information related to matter and energy. Other pioneering research on the theory of digital physics was published by Kantor in 1977 and more recently by Stephen Wolfram in 2002, who thereby joined the host of voices upholding that it is patterns of information, not matter and energy, that constitute the cornerstones of reality. In this paper, we introduce the use of knowledge management tools for the purpose of analysing this topic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-E-J-254: Utility of Pinnacle Dynamic Planning Module Utilizing Deformable Image Registration in Adaptive Radiotherapy

Purpose For certain highly conformal treatment techniques, changes in patient anatomy due to weight loss and/or tumor shrinkage can result in significant changes in dose distribution. Recently, the Pinnacle treatment planning system added a Dynamic Planning module utilizing Deformable Image Registration (DIR). The objective of this study was to evaluate the effectiveness of this software in adapting to altered anatomy and adjusting treatment plans to account for it. Methods We simulated significant tumor response by changing patient thickness and altered chin positions using a commercially-available head and neck (H and N) phantom. In addition, we studied 23 CT image sets of fifteen (15) patients with H and N tumors and eight (8) patients with prostate cancer. In each case, we applied deformable image registration through Dynamic Planning module of our Pinnacle Treatment Planning System. The dose distribution of the original CT image set was compared to the newly computed dose without altering any treatment parameter. Result was a dose if we did not adjust the plan to reflect anatomical changes. Results For the H and N phantom, a tumor response of up to 3.5 cm was correctly deformed by the Pinnacle Dynamic module. Recomputed isodose contours on new anatomies were within 1 mm of the expected distribution. The Pinnacle system configuration allowed dose computations resulting from original plans on new anatomies without leaving the planning system. Original and new doses were available side-by-side with both CT image sets. Based on DIR, about 75% of H and N patients (11/15) required a re-plan using new anatomy. Among prostate patients, the DIR predicted near-correct bladder volume in 62% of the patients (5/8). Conclusions The Dynamic Planning module of the Pinnacle system proved to be an accurate and useful tool in our ability to adapt to changes in patient anatomy during a course of radiotherapy.

60 APPLIED LIFE SCIENCES↗

Exact wave solver for nonparaxial laser beam propagation

Simulations of inertial confinement fusion (ICF) experiments require high-fidelity models for laser beam propagation in a nonuniform plasma with varying index of refraction. We describe a new numerical wave solver that is applicable to centimeter-scale length plasmas encountered in indirect drive ICF applications. The one-way Helmholtz equation (OHE) generalizes the time-harmonic paraxial wave equation to large angles. Here, we present a methodology to numerically evaluate the exact solution to the OHE. This solution is computed by analytically advancing eigenfunctions of the one-way Helmholtz operator along a propagation direction and is applicable to any given index of a refraction profile. We compare our exact method with a commonly used approximate split-step technique for solving the OHE. As a test problem, we consider nonparaxial propagation of Gaussian and speckled beams in a plasma density channel with internal reflection. We find that the split-step approach incurs significant errors compared to the exact solution computed using the novel algorithm.

Belyaev, Mikhail A. (ORCID:0000000224908887)↗