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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Propagation of a Strong Shock Over a Random Bed of Spherical Particles

Propagation of a strong shock through a bed of particles results in complex wave dynamics such as a reflected shock, a transmitted shock, and highly unsteady flow inside the particle bed. In this paper we present three-dimensional numerical simulations of shock propagation in air over a random bed of particles. We assume the flow is inviscid and governed by the Euler equations of gas dynamics. Simulations are carried out by varying the volume fraction of the particle bed at a fixed shock Mach number. We compute the unsteady inviscid streamwise and transverse drag coefficients as a function of time for each particle in the random bed as a function of volume fraction. We show that (i) there are significant variations in the peak drag for the particles in the bed, (ii) the mean peak drag as a function of streamwise distance through the bed decreases with a slope that increases as the volume fraction increases, and (iii) the deviation from the mean peak drag does not correlate with local volume fraction. We also present the local Mach number and pressure contours for the different volume fractions to explain the various observed complex physical mechanisms occurring during the shock-particle interactions. Since the shock interaction with the random bed of particles leads to transmitted and reflected waves, we compute the average flow properties to characterize the strength of the transmitted and reflected shock waves and quantify the energy dissipation inside the particle bed. Finally, to better understand the complex wave dynamics in a random bed, we consider a simpler approximation of a planar shock propagating in a duct with a sudden area change. We obtain Riemann solutions to this problem, which are used to compare with fully resolved numerical simulations.

42 ENGINEERING↗
DOE OSTI2022

Cross Inference of Throughput Profiles Using Micro Kernel Network Method

Dedicated network connections are being increasingly deployed in cloud, centralized and edge computing and data infrastructures, whose throughput profiles are critical indicators of the underlying data transfer performance. Due to the cost and disruptions to physical infrastructures, network emulators, such as Mininet, are often used to generate measurements needed to estimate throughput profiles, typically expressed as a function of the connection round trip time. The profiles estimated using measurements from such emulated networks are usually inaccurate for high bandwidth and high latency connections, since they do not accurately reflect the critical network transport dynamics mainly due to computing and memory constraints of the host. We present a machine learning (ML) method to estimate the throughput profiles using emulation measurements to closely match the testbed and production network profiles. In particular, we propose a micro Kernel Network (mKN) that provides baseline throughput measurements on the host running Mininet emulations, which are used to learn a regression map that converts them to the corresponding testbed measurement estimates. Once initially learned, this map is applied to measurements from subsequent network emulations on the same host. We present experimental measurements to illustrate this approach, and derive generalization equations for the proposed mKN-ML method. Using a four-site scenario emulation, we show the effectiveness of this method in providing accurate concave throughput profiles from inaccurate convex or non-smooth ones indicated by Mininet emulation.

Rao, Nageswara↗
DOE OSTI2017

Design Rules for Tailoring Antireflection Properties of Hierarchical Optical Structures

Hierarchical structures consisting of small sub-wavelength features stacked atop larger structures have been demonstrated as an effective means of reducing the reflectance of surfaces. However, optical devices require different antireflective properties depending on the application, and general unifying guidelines on hierarchical structures' design to attain a desired antireflection spectral response are still lacking. The type of reflectivity (diffuse, specular, or total/hemispherical) and its angular- and spectral-dependence are all dictated by the structural parameters. Through computational and experimental studies, guidelines have been devised to modify these various aspects of reflectivity across the solar spectrum by proper selection of the features of hierarchical structures. In this wavelength regime, micrometer-scale substructures dictate the long-wavelength spectral response and effectively reduce specular reflectance, whereas nanometer-scale substructures dictate primarily the visible wavelength spectral response and reduce diffuse reflectance. Coupling structures having these two length scales into hierarchical arrays impressively reduces surfaces' hemispherical reflectance across a broad spectrum of wavelengths and angles. Furthermore, such hierarchical structures in silicon are demonstrated having an average total reflectance across the solar spectrum of 1.1% (average weighted reflectance of 1% in the 280–2500 nm range of the AM 1.5 G spectrum) and specular reflectance <1% even at angles of incidence as high as 67°.

14 SOLAR ENERGY↗
DOE OSTI2015

SHELXT – Integrated space-group and crystal-structure determination

SHELXT automates routine small-molecule structure determination starting from single-crystal reflection data, the Laue group and a reasonable guess as to which elements might be present. The new computer program SHELXT employs a novel dual-space algorithm to solve the phase problem for single-crystal reflection data expanded to the space group P1. Missing data are taken into account and the resolution extended if necessary. All space groups in the specified Laue group are tested to find which are consistent with the P1 phases. After applying the resulting origin shifts and space-group symmetry, the solutions are subject to further dual-space recycling followed by a peak search and summation of the electron density around each peak. Elements are assigned to give the best fit to the integrated peak densities and if necessary additional elements are considered. An isotropic refinement is followed for non-centrosymmetric space groups by the calculation of a Flack parameter and, if appropriate, inversion of the structure. The structure is assembled to maximize its connectivity and centred optimally in the unit cell. SHELXT has already solved many thousand structures with a high success rate, and is optimized for multiprocessor computers. It is, however, unsuitable for severely disordered and twinned structures because it is based on the assumption that the structure consists of atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Prediction of perovskite and other ternary oxide multilayers as mirrors for soft X-rays

Highlights: • Theoretical framework to identify optimal multilayers for high x-ray reflectivity. • Perovskite oxide multilayers are proposed as soft x-ray mirrors. • High soft x-ray reflectivity up to 33% are expected in perovskite multilayers. - Abstract: We propose multilayers of perovskite and related ternary oxides with a general formula of ABO{sub 3} as high reflectivity mirrors, especially for use in water window region (2.3–4.4 nm). The high reflectivity combinations of oxide multilayers were deduced using evolutionary search algorithms such as genetic algorithms, and the calculation speed was accelerated via parallel computing methods We also propose a figure of merit for X-ray reflectivity in periodic multilayers systems, which can simplify the future efforts on identifying material combinations, and the search through this multi-dimensional parameter space. The highest reflectivity value was found to be over 33% at 3.1 nm in the water window region. The effect of interface roughness was simulated and the decrease in reflectivity was found to be modest for practically achievable roughness values. This work establishes the foundation for future experimental and theoretical studies towards achieving high reflectivity x-ray mirrors of complex oxide multilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2019

An ab-initio Investigation: The Physical Properties of ScIr{sub 2} Superconductor

Using ab initio technique the physical properties of ScIr{sub 2} superconductor have been investigated with T{sub c} 1.03 K with a MgCu{sub 2}-type structure. We have carried out the plane-wave pseudopotential approach within the framework of the first-principles density functional theory (DFT) implemented within the CASTEP code. The calculated structural parameters confirm a good agreement with the experimental and other theoretical results. Using the Voigt-Reuss-Hill (VRH) averaging scheme the most important elastic properties including the bulk modulus B, shear modulus G, Young’s modulus E and Poisson’s ratio ν of ScIr{sub 2} are determined. At ambient condition, the values of Cauchy pressure and Pugh’s ratio exhibit ductile nature of ScIr{sub 2}. The electronic and optical properties of ScIr{sub 2} were investigated for the first time. The electronic band structure reveals metallic conductivity and the major contribution comes from Ir-5d states. In the ultraviolet region the reflectivity is high up to 50 eV as evident from the reflectivity spectrum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Uncertainty quantification for evaluating impacts of caprock and reservoir properties on pressure buildup and ground surface displacement during geological CO2 sequestration

A series of numerical test cases reflecting broad and realistic ranges of geological formation properties was developed to systematically evaluate and compare the impacts of those properties on geomechanical responses to CO2 injection. A coupled hydro-geomechanical subsurface transport simulator, STOMP (Subsurface Transport over Multiple Phases), was adopted to simulate the CO2 migration process and geomechanical behaviors of the surrounding geological formations. A quasi-Monte Carlo sampling method was applied to efficiently sample a high-dimensional parameter space consisting of injection rate and 14 subsurface formation properties, including porosity, permeability, entry pressure, irreducible gas and aqueous saturation, Young’s modulus, and Poisson’s ratio for both reservoir and caprock. Generalized cross-validation and analysis of variance methods were used to quantitatively measure the significance of the 15 input parameters. Reservoir porosity, permeability, and injection rate were found to be among the most significant factors affecting the geomechanical responses to the CO2 injection. We used a quadrature generalized linear model to build a reduced-order model that can estimate the geomechanical response instantly instead of running computationally expensive numerical simulations.

CO2 geological sequestration↗
DOE OSTI2025

An evaluation of multi-fidelity methods for quantifying uncertainty in projections of ice-sheet mass change

Abstract. This study investigated the computational benefits of using multi-fidelity statistical estimation (MFSE) algorithms to quantify uncertainty in the mass change of Humboldt Glacier, Greenland, between 2007 and 2100 using a single climate change scenario. The goal of this study was to determine whether MFSE can use multiple models of varying cost and accuracy to reduce the computational cost of estimating the mean and variance of the projected mass change of a glacier. The problem size and complexity were chosen to reflect the challenges posed by future continental-scale studies while still facilitating a computationally feasible investigation of MFSE methods. When quantifying uncertainty introduced by a high-dimensional parameterization of the basal friction field, MFSE was able to reduce the mean-squared error in the estimates of the statistics by well over an order of magnitude when compared to a single-fidelity approach that only used the highest-fidelity model. This significant reduction in computational cost was achieved despite the low-fidelity models used being incapable of capturing the local features of the ice-flow fields predicted by the high-fidelity model. The MFSE algorithms were able to effectively leverage the high correlation between each model's predictions of mass change, which all responded similarly to perturbations in the model inputs. Consequently, our results suggest that MFSE could be highly useful for reducing the cost of computing continental-scale probabilistic projections of sea-level rise due to ice-sheet mass change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Perfectly Matched Layers and Characteristic Boundaries in Lattice Boltzmann: Accuracy vs Cost

Artificial boundary conditions (BCs) play a ubiquitous role in numerical simulations of transport phenomena in several diverse fields, such as fluid dynamics, electromagnetism, acoustics, geophysics, and many more. They are essential for accurately capturing the behavior of physical systems whenever the simulation domain is truncated for computational efficiency purposes. Ideally, an artificial BC would allow relevant information to enter or leave the computational domain without introducing artifacts or unphysical effects. Boundary conditions designed to control spurious wave reflections are referred to as nonreflective boundary conditions (NRBCs). Another approach is given by the perfectly matched layers (PMLs), in which the computational domain is extended with multiple dampening layers, where outgoing waves are absorbed exponentially in time. Here, in this work, the definition of PML is revised in the context of the lattice Boltzmann method. The impact of adopting different types of BCs at the edge of the dampening zone is evaluated and compared, in terms of both accuracy and computational costs. It is shown that for sufficiently large buffer zones, PMLs allow stable and accurate simulations even when using a simple zeroth-order extrapolation BC. Moreover, employing PMLs in combination with NRBCs potentially offers significant gains in accuracy at a modest computational overhead, provided the parameters of the BC are properly tuned to match the properties of the underlying fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Making use of x-ray optical effects in photoelectron-, Auger electron-, and x-ray emission spectroscopies: Total reflection, standing-wave excitation, and resonant effects

We present a general theoretical methodology and related open-access computer program for carrying out the calculation of photoelectron, Auger electron, and x-ray emission intensities in the presence of several x-ray optical effects, including total reflection at grazing incidence, excitation with standing-waves produced by reflection from synthetic multilayers and at core-level resonance conditions, and the use of variable polarization to produce magnetic circular dichroism. Calculations illustrating all of these effects are presented, including in some cases comparisons to experimental results. Sample types include both semi-infinite flat surfaces and arbitrary multilayer configurations, with interdiffusion/roughness at their interfaces. These x-ray optical effects can significantly alter observed photoelectron, Auger, and x-ray intensities, and in fact lead to several generally useful techniques for enhancing surface and buried-layer sensitivity, including layer-resolved densities of states and depth profiles of element-specific magnetization. The computer program used in this study should thus be useful for a broad range of studies in which x-ray optical effects are involved or are to be exploited in next-generation surface and interface studies of nanoscale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2022

Capacitance of Carbon Nanotube/Graphene Composite Electrodes with [BMIM + ][BF$_4^–$]/Acetonitrile: Fixed Voltage Molecular Dynamics Simulations

The capacitance of nanoporous carbon electrode materials is dictated by the physical interactions with electrolyte ions and molecules at the accessible interior electrode surface. While significant progress has been made in designing and synthesizing carbon materials with well-defined and relatively homogeneous nanoporosity, the majority of materials remain heterogeneous, and the resulting properties reflect an average over this distribution. In this regard, computer simulations can be a valuable tool to predict or validate structure/property relationships by systematically investigating well-defined, model electrode morphologies. In this work, we utilize fixed-voltage molecular dynamics simulations to predict structure/capacitance relationships for five different model morphologies of carbon nanotube/graphene (CNT/G) composite electrodes with 1-butyl-3-methylimidazolium tetrafluoroborate/acetonitrile electrolyte. The CNT/G electrode models are inspired by experimental “layer-by-layer” syntheses and strike a balance between realism and computational tractability. Comparison between different model electrode architectures elucidates important structure/property relationships. We find that CNT/graphene contact points serve as “hot spots” with significantly enhanced charge separation relative to the rest of the electrode. Furthermore, we demonstrate a specific nanoconfinement motif that provides substantial 3–4$\times$ enhancement of local capacitance, resulting in a ~40% increase of the total electrode differential capacitance. Because the accessible surface area of the model CNT/G electrodes is precisely determined, a comparison of per-area capacitance across systems is unambiguous. Our results thus complement a prior computational demonstration of capacitance enhancement in nanoconfinement while elucidating additional interaction motifs at CNT/G electrochemical interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mid- to long-wave infrared computational spectroscopy with a graphene metasurface modulator

In recent years there has been much interest concerning the development of modulators in the mid- to long-wave infrared, based on emerging materials such as graphene. These have been frequently pursued for optical communications, though also for other specialized applications such as infrared scene projectors. Here we investigate a new application for graphene modulators in the mid- to long-wave infrared. We demonstrate, for the first time, computational spectroscopy in the mid- to long-wave infrared using a graphene-based metasurface modulator. Furthermore, our metasurface device operates at low gate voltage. To demonstrate computational spectroscopy, we provide our algorithm with the measured reflection spectra of the modulator at different gate voltages. We also provide it with the measured reflected light power as a function of the gate voltage. The algorithm then estimates the input spectrum. We show that the reconstructed spectrum is in good agreement with that measured directly by a Fourier transform infrared spectrometer, with a normalized mean-absolute-error (NMAE) of 0.021.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Excess entropy and crystallization in Stillinger-Weber and Lennard-Jones fluids

Molecular dynamics simulations are used to contrast the supercooling and crystallization behaviour of monatomic liquids that exemplify the transition from simple to anomalous, tetrahedral liquids. As examples of simple fluids, we use the Lennard-Jones (LJ) liquid and a pair-dominated Stillinger-Weber liquid (SW{sub 16}). As examples of tetrahedral, water-like fluids, we use the Stillinger-Weber model with variable tetrahedrality parameterized for germanium (SW{sub 20}), silicon (SW{sub 21}), and water (SW{sub 23.15} or mW model). The thermodynamic response functions show clear qualitative differences between simple and water-like liquids. For simple liquids, the compressibility and the heat capacity remain small on isobaric cooling. The tetrahedral liquids in contrast show a very sharp rise in these two response functions as the lower limit of liquid-phase stability is reached. While the thermal expansivity decreases with temperature but never crosses zero in simple liquids, in all three tetrahedral liquids at the studied pressure, there is a temperature of maximum density below which thermal expansivity is negative. In contrast to the thermodynamic response functions, the excess entropy on isobaric cooling does not show qualitatively different features for simple and water-like liquids; however, the slope and curvature of the entropy-temperature plots reflect the heat capacity trends. Two trajectory-based computational estimation methods for the entropy and the heat capacity are compared for possible structural insights into supercooling, with the entropy obtained from thermodynamic integration. The two-phase thermodynamic estimator for the excess entropy proves to be fairly accurate in comparison to the excess entropy values obtained by thermodynamic integration, for all five Lennard-Jones and Stillinger-Weber liquids. The entropy estimator based on the multiparticle correlation expansion that accounts for both pair and triplet correlations, denoted by S{sub trip}, is also studied. S{sub trip} is a good entropy estimator for liquids where pair and triplet correlations are important such as Ge and Si, but loses accuracy for purely pair-dominated liquids, like LJ fluid, or near the crystallization temperature (T{sub thr}). Since local tetrahedral order is compatible with both liquid and crystalline states, the reorganisation of tetrahedral liquids is accompanied by a clear rise in the pair, triplet, and thermodynamic contributions to the heat capacity, resulting in the heat capacity anomaly. In contrast, the pair-dominated liquids show increasing dominance of triplet correlations on approaching crystallization but no sharp rise in either the pair or thermodynamic heat capacities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Theoretical study of triaxial shapes of neutron-rich Mo and Ru nuclei

Here, whether atomic nuclei can possess triaxial shapes at their ground states is still a subject of ongoing debate. According to theory, good prospects for low-spin triaxiality are in the neutron-rich Mo-Ru region. Recently, transition quadrupole moments in rotational bands of even-mass neutron-rich isotopes of molybdenum and ruthenium nuclei have been measured. The new data have provided a challenge for theoretical descriptions invoking stable triaxial deformations. The purpose of this study is to understand experimental data on rotational bands in the neutron-rich Mo-Ru region, we carried out theoretical analysis of moments of inertia, shapes, and transition quadrupole moments of neutron-rich even-even nuclei around 110 Ru using self-consistent mean-field and shell model techniques. Methods: To describe yrast structures in Mo and Ru isotopes, we use nuclear density functional theory (DFT) with the optimized energy density functional UNEDF0. We also apply triaxial projected shell model (TPSM) to describe yrast and positive-parity, near-yrast band structures. As a result, our self-consistent DFT calculations predict triaxial ground-state deformations in 106,108 Mo and 108,110,112 Ru and reproduce the observed low-frequency behavior of moments of inertia. As the rotational frequency increases, a negative-gamma structure, associated with the aligned ν(h 11/2 ) 2 pair, becomes energetically favored. The computed transition quadrupole moments vary with angular momentum, which reflects deformation changes with rotation; those variations are consistent with experiment. The TPSM calculations explain the observed band structures assuming stable triaxial shapes. Lastly, the structure of neutron-rich even-even nuclei around Ru-110 is consistent with triaxial shape deformations. Our DFT and TPSM frameworks provide a consistent and complementary description of experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An Improved Algorithm for Estimating Surface Shortwave Radiation: Preliminary Evaluation With MODIS Products

Cloud parameters, as key inputs in radiative transfer algorithms, have a critical impact on surface shortwave radiation (SSR) computation. By introducing a parameterization of cloud transmittance and reflectance, based on radiative transfer simulations, this study improves the accuracy of an existing physically based model which severely underestimates SSR under thick cloud conditions. The cloud parameterization adopts the single-layer cloud model and simulates cloud transmittances and reflectances by varying cloud optical thickness, cloud particle size, and solar zenith angle. The revised model is applied to estimate instantaneous SSR using Moderate-resolution Imaging Spectroradiometer (MODIS) atmospheric and land products. The retrieved SSR is evaluated against observation data from 41 Baseline Surface Radiation Network (BSRN) stations and is also compared with the MODIS official SSR product. The root mean square error (RMSE) of the estimated instantaneous radiation is approximately 52 and 98 W m -2 under clear-sky and all-sky conditions, respectively. The accuracy of the improved parameterization is higher than that of the original model, and there is no obvious underestimation of SSR in the case of high cloud optical thickness. Therefore, the new algorithm improves the accuracy of SSR estimates in the presence of thick clouds. Retrievals with the improved model also achieve higher accuracy than the MODIS official SSR product (MCD18A1). To conclude, the reliable performance of the scheme at most BSRN stations illustrates that the improved model can be used to map SSR on a global scale.

42 ENGINEERING↗
DOE OSTI2019

Nearly quantized conductance plateau of vortex zero mode in an iron-based superconductor

Majorana zero modes (MZMs) are spatially-localized zero-energy fractional quasiparticles with non-Abelian braiding statistics that hold promise for topological quantum computing. Owing to the particle-antiparticle equivalence, MZMs reflect quantized conductance at low temperature. By utilizing variable-tunnel-coupled scanning tunneling spectroscopy, we study tunneling conductance of vortex bound states on FeTe 0.55 Se 0.45 superconductors. Here, we report observations of conductance plateaus as a function of tunnel coupling for zero-energy vortex bound states with values close to or even reaching the 2 e 2 / h quantum conductance (here e is the electron charge and h is Planck’s constant). In contrast, no plateaus were observed on either finite energy vortex bound states or in the continuum of electronic states outside the superconducting gap. This behavior of the zero-mode conductance confirms the existence of MZMs in FeTe 0.55 Se 0.45 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Method for In-Situ Q₀ Measurements of High-Quality SRF Resonators

Accelerator projects such as LCLS-II naturally require low-loss superconducting (SRF) cavities. Due to strong demand for improving intrinsic quality factor (Q₀), importance of accurate cavity characterization increases. We propose a method to measure Q₀ in situ for an SRF resonator installed in its cryogenic module and connected with a RF feed source via a fixed RF coupler. The method exploits measurements of a response for an SRF resonator fed by an amplitude-modulated signal. Such a signal can be synthesized as a beat-wave composed of two frequencies that are close to the resonant frequency. Analyzing the envelope of the reflected signal, one can find the difference in reflection for the chosen frequencies and use them to compute the intrinsic Q. We also develop the methodology to carry out measurements of Q₀ at the nominal cavity operating voltage. We verified our method in experiments with a room temperature copper resonator and with two SRF resonators including Fermilab’s 650 MHz cavity and JLab’s 1500 MHz cavity.

43 PARTICLE ACCELERATORS↗