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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sum-of-Fractions Methodology for Actinides in Water- and Polyethylene-Moderated and -Reflected Systems

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232 U, 233 U, 234 U, 235 U, 237 Np, 236 Pu, 238 Pu, 239 Pu, 240 Pu, 241 Pu, 242 Pu, 241 Am, 242 m Am, 243 Am, 242 Cm, 243 Cm, 244 Cm, 245 Cm, 246 Cm, 247 Cm, 249 Cf, and 251 Cf. The methodology uses available benchmark data for mixtures of 233 U, 235 U, and 239 Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Sum-of-Fractions Method

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232U, 233U, 234U, 235U, 237Np, 236Pu, 238Pu, 239Pu, 240Pu, 241Pu, 242Pu, 241Am, 242mAm, 243Am, 242Cm, 243Cm, 244Cm, 245Cm, 246Cm, 247Cm, 249Cf, and 251Cf. The methodology uses available benchmark data for mixtures of 233U, 235U, and 239Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

criticality safety, Actinide↗
DOE OSTI2017

Memory-optimized shift operator alternating direction implicit finite difference time domain method for plasma

Through introducing the alternating direction implicit (ADI) technique and the memory-optimized algorithm to the shift operator (SO) finite difference time domain (FDTD) method, the memory-optimized SO-ADI FDTD for nonmagnetized collisional plasma is proposed and the corresponding formulae of the proposed method for programming are deduced. In order to further the computational efficiency, the iteration method rather than Gauss elimination method is employed to solve the equation set in the derivation of the formulae. Complicated transformations and convolutions are avoided in the proposed method compared with the Z transforms (ZT) ADI FDTD method and the piecewise linear JE recursive convolution (PLJERC) ADI FDTD method. The numerical dispersion of the SO-ADI FDTD method with different plasma frequencies and electron collision frequencies is analyzed and the appropriate ratio of grid size to the minimum wavelength is given. The accuracy of the proposed method is validated by the reflection coefficient test on a nonmagnetized collisional plasma sheet. The testing results show that the proposed method is advantageous for improving computational efficiency and saving computer memory. The reflection coefficient of a perfect electric conductor (PEC) sheet covered by multilayer plasma and the RCS of the objects coated by plasma are calculated by the proposed method and the simulation results are analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

In Pursuit of 2D Materials for Maximum Optical Response

Despite being only a few atoms thick, single-layer two-dimensional (2D) materials display strong electron–photon interactions that could be utilized in efficient light modulators on extreme subwavelength scales. In various applications involving light modulation and manipulation, materials with strong optical response at different wavelengths are required. Using qualitative analytical modeling and first-principles calculations, we determine the theoretical limit of the maximum optical response such as absorbance ($A$) and reflectance ($R$) in 2D materials and also conduct a computational survey to seek out those with best $A$ and $R$ in various frequency ranges, from mid-infrared to deep-ultraviolet. We find that 2D boron has broadband reflectance R > 99% for >100 layers, surpassing conventional thin films of bulk metals such as silver. Moreover, we identify 2D monolayer semiconductors with maximum response, for which we obtain quantitative estimates by calculating quasiparticle energies and accounting for excitonic effects by solving the Bethe–Salpeter equation. We found several monolayer semiconductors with absorbances ≳30% in different optical ranges, which are more than half of the maximum possible value, Alim = 1/2, for a freestanding 2D material. Our study predicts 2D materials which can potentially be used in ultrathin reflectors and absorbers for optoelectronic application in various frequency ranges.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solving the Grid Optimization Competition Challenge 3 Problem

The Grid Optimization Competition Challenge 3 Problem posed a multiperiod security-constrained unit commitment problem with base-case AC power flow. The problem formulation includes binary unit commitment decisions, nonlinear AC power flow and balance, dispatchable loads, and linearized contingency real power flow, among other features. This talk will present a modified consensus ADMM algorithm, which splits the problem into mixed-integer linear and nonlinear components, as a heuristic solution method for this large-scale mixed integer nonlinear program. We will present some computational results from the competition for our implementation and reflect on the challenges of participating the grid optimization competition.

AC power flow↗
DOE OSTI2019

A Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT): Part II. The cloudy-sky model

The Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT) reported in Part I of this report is enhanced to include the requirements for cloudy-sky conditions. Surface radiances in 2002 narrow-wavelength bands from 0.28 to 4.0 um are analytically computed by solving the radiative transfer equation for five independent photon paths accounting for clear-sky absorption, Rayleigh scattering, and cloud absorption and scattering. The Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS) is used to show the optical thickness of the clear-sky atmosphere. Unlike Part I, which approximates the computation of aerosol scattering using the single-scattering phase function, the cloud transmittance and reflectance are efficiently retrieved from a comprehensive look-up table pre-computed by a 32-stream DIScrete Ordinates Radiative Transfer (DISORT) model for possible cloud conditions as well as solar and viewing geometries. A resolution assessment is performed to observe the optimal balance between the computational efficiency and accuracy in the development of the look-up table. Model simulations by DISORT and TMYSPEC are used to evaluate the performance of FARMS-NIT under cloudy-sky conditions. Compared to DISORT, FARMS-NIT yields 2-3% uncertainties on average, but it substantially reduces the computational time because of the independent computation of cloud properties and the implementation of the look-up table. In contrast to TMYSPEC, which uses successive steps to empirically compute plane-of-array (POA) irradiances and spectral irradiances, FARMS-NIT directly solves spectral radiances from the radiative transfer equation, which profoundly increases the accuracy in surface irradiances, especially over inclined photovoltaics (PV) panels.

14 SOLAR ENERGY↗
DOE OSTI2013

X-RAY REFLECTED SPECTRA FROM ACCRETION DISK MODELS. III. A COMPLETE GRID OF IONIZED REFLECTION CALCULATIONS

We present a new and complete library of synthetic spectra for modeling the component of emission that is reflected from an illuminated accretion disk. The spectra were computed using an updated version of our code XILLVER that incorporates new routines and a richer atomic database. We offer in the form of a table model an extensive grid of reflection models that cover a wide range of parameters. Each individual model is characterized by the photon index {Gamma} of the illuminating radiation, the ionization parameter {xi} at the surface of the disk (i.e., the ratio of the X-ray flux to the gas density), and the iron abundance A{sub Fe} relative to the solar value. The ranges of the parameters covered are 1.2 {<=} {Gamma} {<=} 3.4, 1 {<=} {xi} {<=} 10{sup 4}, and 0.5 {<=} A{sub Fe} {<=} 10. These ranges capture the physical conditions typically inferred from observations of active galactic nuclei, and also stellar-mass black holes in the hard state. This library is intended for use when the thermal disk flux is faint compared to the incident power-law flux. The models are expected to provide an accurate description of the Fe K emission line, which is the crucial spectral feature used to measure black hole spin. A total of 720 reflection spectra are provided in a single FITS file (http://hea-www.cfa.harvard.edu/{approx}javier/xillver/) suitable for the analysis of X-ray observations via the atable model in XSPEC. Detailed comparisons with previous reflection models illustrate the improvements incorporated in this version of XILLVER.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Toward Molecular Catalysts by Computer

Rational design of molecular catalysts requires a systematic approach to designing ligands with specific functionality and precisely tailored electronic and steric properties. It then becomes possible to devise computer protocols to predict accurately the required properties and ultimately to design catalysts by computer. In this account we first review how thermodynamic properties such as oxidation-reduction potentials (E0), acidities (pKa), and hydride donor abilities (ΔGH-) form the basis for a systematic design of molecular catalysts for reactions that are critical for a secure energy future (hydrogen evolution and oxidation, oxygen and nitrogen reduction, and carbon dioxide reduction). We highlight how density functional theory allows us to determine and predict these properties within “chemical” accuracy (~ 0.06 eV for redox potentials, ~ 1 pKa unit for pKa values, and ~ 1.5 kcal/mol for hydricities). These quantities determine free energy maps and profiles associated with catalytic cycles, i.e. the relative energies of intermediates, and help us distinguish between desirable and high-energy pathways and mechanisms. Good catalysts have flat profiles that avoid high activation barriers due to low and high energy intermediates. We illustrate how the criterion of a flat energy profile lends itself to the prediction of design points by computer for optimum catalysts. This research was carried out in the Center for Molecular Electro-catalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory (PNNL) is operated for the DOE by Battelle.

acidity↗
DOE OSTI2019

Adsorption mechanism of alkyl hydroxamic acid onto bastnäsite: Fundamental steps toward rational collector design for rare earth elements

Rare earth element (REE) production is limited in part by inefficient strategies for beneficiation. Hydroxamic acid ligands are promising reagents for the selective flotation of bastnäsite [(Ce,La)FCO 3 ], a major REE ore mineral, but the mechanism and energetics of adsorption are not understood, interfering with the design of new, more efficient reagents. In this work, the adsorption of octyl hydroxamic acid onto bastnäsite was calculated using a combination of experimental and computational methods. In-situ attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy revealed changes in the hydroxamate functional group vibrational frequencies, corresponding to chelation with cerium cations at the bastnäsite surface. The results indicate a monodentate chemisorption mechanism at low surface loading that changes to bidentate chemisorption at higher concentrations. This interpretation is supported by molecular vibrational frequency shifts calculated using density functional theory (DFT), and orientation of the hydrocarbon chain measured by sum frequency generation (SFG) vibrational spectroscopy. The binding enthalpies of octyl hydroxamic acid interacting with La and Ce-bastnäsite surfaces were measured using isothermal titration calorimetry (ITC) revealing a stronger coordinating ability with bastnäsite than with a common gangue mineral, calcite (CaCO 3 ). Because octyl hydroxamate favors monodentate adsorption at low surface coverages, the relative chelating strength of metal ions could be a poor predictor for selectivity under monolayer adsorption conditions. At higher surface loadings, where the bidentate mode of adsorption is active, selectivity is likely to be limited by increased flotation of gangue ores.

58 GEOSCIENCES↗
DOE OSTI2023

Open-Shell Tensor Hypercontraction

The extension of least-squares tensor hypercontracted second- and third-order Møller–Plesset perturbation theory (LS-THC-MP2 and LS-THC-MP3) to open-shell systems is an important development due to the scaling reduction afforded by THC and the ubiquity of molecular ions, radicals, and other open-shell reactive species. The complexity of wavefunction-based quantum-chemical methods such as Møller–Plesset and coupled cluster theory is reflected in the steep scaling of the computational costs with the molecular size. The least-squares tensor hypercontraction (LS-THC) method is an efficient, single-step factorization for the two-electron integral tensor but can also be used to factorize the double excitation amplitudes, leading to significant scaling reduction. Herein we extend this promising method to open-shell variants of LS-THC-MP2 and -MP3 by using diagrammatic techniques and explicit spin summation. The accuracy of the resulting methods for open-shell species is benchmarked on standard test systems such as regular alkanes as well as realistic systems involving bond breaking, radical stabilization, and other effects. We find that open-shell LS-THC-MPn methods exhibit errors highly comparable to those produced by closed-shell LS-THC-MPn and are highly insensitive to particular chemical interactions, geometries, or even moderate spin contamination.

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

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗
DOE OSTI2014

Using Hessian update formulae to construct modified Shepard interpolated potential energy surfaces: Application to vibrating surface atoms

Modified Shepard interpolation based on second order Taylor series expansions has proven to be a flexible tool for constructing potential energy surfaces in a range of situations. Extending this to gas–surface dynamics where surface atoms are allowed to move represents a substantial increase in the dimensionality of the problem, reflected in a dramatic increase in the computational cost of the required Hessian (matrix of second derivatives) evaluations. This work demonstrates that using approximate Hessians derived from well known Hessian update formulae and a single accurate Hessian can provide an effective way to avoid this expensive accurate Hessian determination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Laboratory-based x-ray reflectometer for multilayer characterization in the 15–150 keV energy band

A laboratory-based X-ray reflectometer has been developed to measure the performance of hard X-ray multilayer coatings at their operational X-ray energies and incidence angles. The instrument uses a sealed-tube X-ray source with a tungsten anode that can operate up to 160 kV to provide usable radiation in the 15–150 keV energy band. Two sets of adjustable tungsten carbide slit assemblies, spaced 4.1 m apart, are used to produce a low-divergence white beam, typically set to 40 μm × 800 μm in size at the sample. Multilayer coatings under test are held flat using a vacuum chuck and are mounted at the center of a high-resolution goniometer used for precise angular positioning of the sample and detector; additionally, motorized linear stages provide both vertical and horizontal adjustments of the sample position relative to the incident beam. A CdTe energy-sensitive detector, located behind a third adjustable slit, is used in conjunction with pulse-shaping electronics and a multi-channel analyzer to capture both the incident and reflected spectra; the absolute reflectance of the coating under test is computed as the ratio of the two spectra. The instrument’s design, construction, and operation are described in detail, and example results are presented obtained with both periodic, narrow-band and depth-graded, wide-band hard X-ray multilayer coatings.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Sci—Thur PM: Imaging — 06: Canada's National Computed Tomography (CT) Survey

The value of computed tomography (CT) in medical imaging is reflected in its' increased use and availability since the early 1990's; however, given CT's relatively larger exposures (vs. planar x-ray) greater care must be taken to ensure that CT procedures are optimised in terms of providing the smallest dose possible while maintaining sufficient diagnostic image quality. The development of CT Diagnostic Reference Levels (DRLs) supports this process. DRLs have been suggested/supported by international/national bodies since the early 1990's and widely adopted elsewhere, but not on a national basis in Canada. Essentially, CT DRLs provide guidance on what is considered good practice for common CT exams, but require a representative sample of CT examination data to make any recommendations. Canada's National CT Survey project, in collaboration with provincial/territorial authorities, has collected a large national sample of CT practice data for 7 common examinations (with associated clinical indications) of both adult and pediatric patients. Following completion of data entry into a common database, a survey summary report and recommendations will be made on CT DRLs from this data. It is hoped that these can then be used by local regions to promote CT practice optimisation and support any dose reduction initiatives.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Cyclops I: Modeling Error Study for Inertially Confined Detonation of Insensitive High Explosive (IHE) PBX-9502

A validation study is performed using the hydrodynamics code PAGOSA to assess the accuracy of reflected shock, burn front propagation, and Photon Doppler Velocimetry (PDV) predictions using several burn models and strength models. Experimental data collected from proton radiography (pRad) is compared to burn front propagation times and reflected shock patterns. Burn front isochrones are computed and the lag time for the burn front is quantified for a given burn model. The effect of strength models and damage are assessed by comparison of simulations to experimental PDV data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Conservative Estimation of Tail Probabilities from Limited Sample Data

Several sparse-sample uncertainty quantification (UQ) methods are compared for conservative but not overly conservative estimation of small tail probabilities involving responses that lay beyond specified thresholds in the tails of probability distributions. Sixteen very differently shaped distributions (or probability density functions, PDFs) and tail probability magnitudes ranging from 10 -5 to 10 -1 are considered in order for the study to be relevant to a wide range of risk analysis and quantification of margins and uncertainty (QMU) problems. The emphasis of the study is on limited data regimes ranging from N = 2 to 20 samples, reflective of most experimental and some expensive computational situations. Relatively simple sparse-sample UQ methods tested for this regime involve statistical tolerance interval "Equivalent Normal and related "Ensemble of Normals" and "Superdistribution (SD) approaches. (The independently derived SD is effectively equivalent to the Bayesian posterior predictive distribution given the assumptions of the derivation.) The performance of the methods was generally improved for N ≥ 5 samples with a generalized Jackknife resampling technique, which determines a tail probability estimate by averaging estimates from smaller sub-samples. Several quantitative metrics for method conservatism and accuracy of tail probability estimation are used to assess and rank the methods' performance over many random trials for each test PDF and probability magnitude. A variant of Bootstrap resampling was also tried, but did not significantly improve tail probability estimates in most cases. Detailed results are presented from over 100-million tests over the above factors that provide useful granular information on which methods or combination of methods perform best in various areas of the factor space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

The O{sub 2} A-Band in the Fluxes and Polarization of Starlight Reflected by Earth-Like Exoplanets

Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmospheres and surfaces can be derived by analyzing the light of the parent star reflected by the planet. We investigate the influence of the surface albedo A {sub s}, the optical thickness b {sub cloud}, the altitude of water clouds, and the mixing ratio of biosignature O{sub 2} on the strength of the O{sub 2} A-band (around 760 nm) in the flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ( η < 0.4) will yield moderately deep bands in flux and moderate-to-small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on properties such as optical thickness, top altitude, particle phase, coverage fraction, and horizontal distribution. Depending on the surface albedo and cloud properties, different O{sub 2} mixing ratios η can give similar absorption-band depths in flux and band strengths in polarization, especially if the clouds have moderate-to-high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, especially not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O{sub 2} or any other absorbing gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Incredibility of a Postulated Criticality Event for Enriched Uranyl Nitrate Storage

At the H-Canyon outside facilities, highly enriched uranyl nitrate solutions collected from leaks or during deliberate maintenance activities is collected into 5-gallon polyethylene bottles. The current criticality safety evaluation for the facility postulates that a criticality event is possible if the solutions are allowed to concentrate through evaporation, freezing, storing multiple bottles in a geometrically unfavorable array, or any combination of these conditions. Using chemical data on uranyl nitrate freezing from Savannah River National Laboratory, protected limits on concentration and enrichment, and normal operational restrictions, criticality was shown to be incredible under all normal and credible abnormal conditions. A bounding normal solution is analyzed under normal conditions and shown to be substantially subcritical. A bounding credible abnormal solution is defined as 9.2 g U{sup 235}/L and 73 wt. % enriched and 0 M excess nitric acid. The concentration and enrichment are limitations already protected by the facility for other processes. Because acidity is not specifically restricted, the most reactive lower bound is selected. The credible abnormal worst placement is in a concrete basin reflected on three sides. There the bottle(s) would be subject in the worst cases to flooding (full water reflection) and postulated concentration through evaporation or freezing. Computational modeling was performed using KENOVI in the SCALE 6.1 code package. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗