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 235 records · Page 13

Trying to predict the future: Resource planning and allocation in CMS

In the large LHC experiments the majority of computing resources are provided by the participating countries. These resource pledges account for more than three quarters of the total available computing. The experiments are asked to give indications of their requests three years in advance and to evolve these as the details and constraints become clearer. In this paper we will discuss the resource planning techniques used in CMS to predict the computing resources several years in advance. We will discuss how we attempt to implement the activities of the computing model in spread-sheets and formulas to calculate the needs. We will talk about how those needs are reflected in the 2012 running and how the planned long shutdown of the LHC in 2013 and 2014 impacts the planning process and the outcome. In the end we will speculate on the computing needs in the second major run of LHC.

Bloom, Kenneth↗

Conversion acoustic resonances in orthorhombic crystals

A classification of acoustic-beam reflection resonances in orthorhombic crystals under conditions where a proximity to conversion is implemented in the vicinity of total internal reflection is proposed. In this case, the energy from the incident pump beam falls almost entirely into a narrow intense reflected beam propagating at a small angle with respect to the surface. The crystal boundary is parallel to one of the elastic symmetry planes, and the excited beam propagates near one of axes 2 in this plane. Depending on the relations between the elastic moduli and the chosen propagation geometry, 18 types of resonances may occur, but no more than three in each crystal. The developed theory combines an approximate analytical description and accurate computer analysis. The relations between the elastic moduli providing minimum energy loss over the parasite reflected wave are determined. Some crystals with resonant excitation very close to conversion are revealed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A MULTIPLE SCATTERING POLARIZED RADIATIVE TRANSFER MODEL: APPLICATION TO HD 189733b

We present a multiple scattering vector radiative transfer model that produces disk integrated, full phase polarized light curves for reflected light from an exoplanetary atmosphere. We validate our model against results from published analytical and computational models and discuss a small number of cases relevant to the existing and possible near-future observations of the exoplanet HD 189733b. HD 189733b is arguably the most well observed exoplanet to date and the only exoplanet to be observed in polarized light, yet it is debated if the planet’s atmosphere is cloudy or clear. We model reflected light from clear atmospheres with Rayleigh scattering, and cloudy or hazy atmospheres with Mie and fractal aggregate particles. We show that clear and cloudy atmospheres have large differences in polarized light as compared to simple flux measurements, though existing observations are insufficient to make this distinction. Futhermore, we show that atmospheres that are spatially inhomogeneous, such as being partially covered by clouds or hazes, exhibit larger contrasts in polarized light when compared to clear atmospheres. This effect can potentially be used to identify patchy clouds in exoplanets. Given a set of full phase polarimetric measurements, this model can constrain the geometric albedo, properties of scattering particles in the atmosphere, and the longitude of the ascending node of the orbit. The model is used to interpret new polarimetric observations of HD 189733b in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Whole-Cell Scan using Automatic Variable-Angle and Variable-Illumination-Depth Pseudo—Total Internal Reflection Fluorescence Microscopy

An automatic calibration and angle-scanning prism-type total internal reflection fluorescence microscope (TIRFM) was modified to function in both TIRFM and pseudo-TIRFM modes. When the incident angle of the excitation laser beam was controlled to be larger than the critical angle, the instrument served as a variable-angle TIRFM. A homemade computer program automatically calibrates the laser illumination spot in the sample to overlap with the center of the microscope's field of view. Then, by measuring the fluorescence intensities at different incident angles, the z-positions of fluorescent nanospheres close to the cell basolateral membrane can be extracted. When the incident angle is reduced to be in the subcritical range, the instrument works as a pseudo-TIRFM. The whole cell body from bottom to top can be imaged in a vertical scan process. Furthermore, the illumination field depth in the pseudo-TIRFM can be controlled by changing the incident angle or the horizontal position of the laser spot.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exceptionally Complex Electronic Structures of Lanthanide Oxides and Small Molecules

Lanthanide (Ln) oxide clusters and molecular systems provide a bottom-up look at the electronic structures of the bulk materials because of close parallels in the patterns of Ln 4fN subshell occupancy between the molecular and bulk Ln 2 O 3 size limits. At the same time, these clusters and molecules offer a challenge to the theory community to find appropriate and robust treatments for the 4f N patterns across the Ln series. Anion photoelectron (PE) spectroscopy provides a powerful experimental tool for studying these systems, mapping the energies of the ground and low-lying excited states of the neutral relative to the initial anion state, providing spectroscopic patterns that reflect the Ln 4f N occupancy. In this Account, we review our anion PE spectroscopic and computational studies on a range of small lanthanide molecules and cluster species. The PE spectra of LnO – (Ln = Ce, Pr, Sm, Eu) diatomic molecules show spectroscopic signatures associated with detachment of an electron from what can be described as a diffuse Ln 6s-like orbital. While the spectra of all four diatomics share this common transition, the fine structure in the transition becomes more complex with increasing 4f occupancy. This effect reflects increased coupling between the electrons occupying the corelike 4f and diffuse 6s orbitals with increasing N. Understanding the PE spectra of these diatomics sets the stage for interpreting the spectra of polyatomic molecular and cluster species. In general, the results confirm that the partial 4f N subshell occupancy is largely preserved between molecular and bulk oxides and borides. However, they also suggest that surfaces and edges of bulk materials may support a low-energy, diffuse Ln 6s band, in contrast to bulk interiors, in which the 6s band is destabilized relative to the 5d band. We also identify cases in which the molecular Ln centers have 4f N+1 occupancy rather than bulklike 4f N , which results in weaker Ln–O bonding. Specifically, Sm centers in mixed Ce–Sm oxides or in Sm x O y – (y ≤ x) clusters have this higher 4f N+1 occupancy. The PE spectra of these particular species exhibit a striking increase in the relative intensities of excited-state transitions with decreasing photon energy (resulting in lower photoelectron kinetic energy). Furthermore, this is opposite of what is expected on the basis of the threshold laws that govern photodetachment. We relate this phenomenon to strong electron–neutral interactions unique to these complex electronic structures. The time scale of the interaction, which shakes up the electronic configuration of the neutral, increases with decreasing electron momentum. From a computational standpoint, we point out that special care must be taken when considering Ln cluster and molecular systems toward the center of the Ln series (e.g., Sm, Eu), where treatment of electrons explicitly or using an effective core potential can yield conflicting results on competing subshell occupancies. However, despite the complex electronic structures associated with partially filled 4f N subshells, we demonstrate that inexpensive and tractable calculations yield useful qualitative insight into the general electronic structural features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The c = 1 string theory S-matrix revisited

We revisit the perturbative S-matrix of c = 1 string theory from the worldsheet perspective. We clarify the origin of the leg pole factors, the non-analyticity of the string amplitudes, and the validity as well as limitations of earlier computations based on resonance momenta. We compute the tree level 4-point amplitude and the genus one 2-point reflection amplitude by numerically integrating Virasoro conformal blocks with DOZZ structure constants on the sphere and on the torus, with sufficiently generic complex Liouville momenta, and find agreement with known answers from the c = 1 matrix model

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling pulse-cleaning plasma mirrors from dielectric response to saturation: A particle-in-cell approach

Pulse-cleaning plasma mirrors are widely employed to improve ultraintense laser contrast, but much of the literature concerning their effect on the reflected pulse is empirical. A simulation study of pulse-cleaning plasma mirrors using the particle-in-cell code large scale plasma is presented here. The importance of capturing initial ionization from neutral atoms, collisional effects, and simulation dimensionality is considered. Excellent agreement with experimental data is obtained when a multiphoton ionization model is employed. Furthermore, a series of 2D simulations is shown to accurately replicate both the reflected light intensity and mode obtained from full 3D simulations at significantly reduced computation cost.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The atmospheres of earthlike planets after giant impact events

It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well-known example. In the aftermath of such collisions, the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO{sub 2} and H{sub 2}O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere, and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H{sub 2}O and CO{sub 2}, while the formation of CH{sub 4} and NH{sub 3} is quenched because of short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO{sub 2}. These are apparent in the emerging spectra and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of two lower brightness temperature in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30 m class telescopes. The 1-4 μm will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10{sup 5-6} yr on the basis of the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating Moving Target Defense with PLADD

This project evaluates the effectiveness of moving target defense (MTD) techniques using a new game we have designed, called PLADD, inspired by the game FlipIt. PLADD extends FlipIt by incorporating what we believe are key MTD concepts. We have analyzed PLADD and proven the existence of a defender strategy that pushes a rational attacker out of the game, demonstrated how limited the strategies available to an attacker are in PLADD, and derived analytic expressions for the expected utility of the game’s players in multiple game variants. We have created an algorithm for finding a defender’s optimal PLADD strategy. We show that in the special case of achieving deterrence in PLADD, MTD is not always cost effective and that its optimal deployment may shift abruptly from not using MTD at all to using it as aggressively as possible. We believe our effort provides basic, fundamental insights into the use of MTD, but conclude that a truly practical analysis requires model selection and calibration based on real scenarios and empirical data. We propose several avenues for further inquiry, including (1) agents with adaptive capabilities more reflective of real world adversaries, (2) the presence of multiple, heterogeneous adversaries, (3) computational game theory-based approaches such as coevolution to allow scaling to the real world beyond the limitations of analytical analysis and classical game theory, (4) mapping the game to real-world scenarios, (5) taking player risk into account when designing a strategy (in addition to expected payoff), (6) improving our understanding of the dynamic nature of MTD-inspired games by using a martingale representation, defensive forecasting, and techniques from signal processing, and (7) using adversarial games to develop inherently resilient cyber systems.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

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↗

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↗

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↗

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↗