Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Search for the Chiral Magnetic Effect via Charge-Dependent Azimuthal Correlations Relative to Spectator and Participant Planes in Au + Au Collisions at s N N = 200 GeV

The chiral magnetic effect (CME) refers to charge separation along a strong magnetic field due to imbalanced chirality of quarks in local parity and charge-parity violating domains in quantum chromodynamics. The experimental measurement of the charge separation is made difficult by the presence of a major background from elliptic azimuthal anisotropy. This background and the CME signal have different sensitivities to the spectator and participant planes, and could thus be determined by measurements with respect to these planes. We report such measurements in Au+Au collisions at a nucleon-nucleon center-of-mass energy of 200 GeV at the Relativistic Heavy-Ion Collider. It is found that the charge separation, with the flow background removed, is consistent with zero in peripheral (large impact parameter) collisions. Some indication of finite CME signals is seen with a significance of 1–3 standard deviations in mid-central (intermediate impact parameter) collisions. Furthermore, significant residual background effects may, however, still be present.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

SU-E-T-272: Direct Verification of a Treatment Planning System Megavoltage Linac Beam Photon Spectra Models, and Analysis of the Effects On Patient Plans

Purpose: For the photon dose calculation Philips Pinnacle Treatment Planning System (TPS) uses collapsed cone convolution algorithm, which relies on energy spectrum of the beam in computing the scatter component. The spectrum is modeled based on Linac’s standard commissioning data and typically is not independently verified. We explored a methodology of using transmission measurements in combination with regularization data processing to unfold Linac spectra. The measured spectra were compared to those modeled by the TPS, and the effect on patient plans was evaluated. Methods: Transmission measurements were conducted in narrow-beam geometry using a standard Farmer ionization chamber. Two attenuating materials and two build -up caps, having different atomic numbers, served to enhance discrimination between absorption of low and high-energy portions of the spectra, thus improving the accuracy of the results. The data was analyzed using a regularization technique implemented through spreadsheet-based calculations. Results: The unfolded spectra were found to deviate from the TPS beam models. The effect of such deviations on treatment planning was evaluated for patient plans through dose distribution calculations with either TPS modeled or measured energy spectra. The differences were reviewed through comparison of isodose distributions, and quantified based on maximum dose values for critical structures. While in most cases no drastic differences in the calculated doses were observed, plans with deviations of 4 to 8% in the maximum dose values for critical structures were discovered. The anatomical sites with large scatter contributions are the most vulnerable to inaccuracies in the modeled spectrum. Conclusion: An independent check of the TPS model spectrum is highly desirable and should be included as part of commissioning of a new Linac. The effect is particularly important for dose calculations in high heterogeneity regions. The developed approach makes acquisition of megavoltage Linac beam spectra achievable in a typical radiation oncology clinic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-F-BRD-13: A Phenomenological Relative Biological Effectiveness (RBE) Model for Proton Therapy Based On All Published in Vitro Cell Survival Data

Purpose: Proton therapy treatments are currently planned and delivered using the assumption that the proton relative biological effectiveness (RBE) relative to photons is 1.1. This assumption ignores strong experimental evidence that suggests the RBE varies. A recent review study (Paganetti 2014), collected over 72 experimental reports on proton RBE, providing a comprehensive dataset for predicting proton RBE. Using this dataset we develop a model to predict RBE based on dose, linear energy transfer (LET) and the tissue specific parameter (α/β)x. Methods: The relationship of the RBE on dose, dose average LET (LETd) and (α/β)x was explored using 290 experimental data points. The RBE was calculated for each experimental point for doses ranging from 1 to 10 Gy. An RBE model, based on the linear-quadratic model, was derived using a nonlinear regression fit to the experimental data. Results: The proposed RBE model predicts that the RBE increases with LET and decreases with (α/β)x. The model additionally predicts a decrease in RBE with increasing dose for low (α/β)x values (< 3 Gy). However the model predicts an increase in RBE with increasing dose for large (α/β)x (> 4 Gy). Previous phenomenological models were based on a small subset of experimental data and deviate significantly from our findings. Conclusion: The proposed RBE model is derived using the most comprehensive collection of proton RBE experimental data to date. The model agrees with previous theoretical predictions on the relationship between RBE, LETd and (α/β)x and also makes predictions on the relationship between RBE and dose. The proposed model shows a relationship between α and LETd but no significant relationship between β and LETd. Previously published phenomenological models based on a limited data set might have to be revised.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-EF-BRB-06: Implementation of a Modulated-Arc Total Body Irradiation (TBI) Technique Using the RayStation Treatment Planning System

Purpose: To develop a clinical workflow for delivering a modulated-arc total body irradiation (TBI) with RayStation scripting. This technique uses arc fields with the patient lying at floor level on a padded table and is validated through measurements taken on a custom-made TBI phantom. Methods: Treatment planning was performed for a retrospective cohort of eight patients with a diverse range of heights and body types. Each was replanned using an open-field dual arc method, with the patient in supine and prone positions on the floor of the vault. All plans were optimized using Raystation Planning 4.7.2.5 (RaySearch Laboratories, Stockholm, Sweden), with 200 cGy prescribed to the 95% of the body contour − 5mm. This results in an open-field beam that sweeps craniocaudally across the length of the patient. The technique is validated with measurements at 10 cm intervals in a custom-milled, 5 cm thick acrylic phantom. A centrally located CC13 ion chamber and a Mobile MOSFET (Best Medical Canada, Ottawa, ON) detector array were used to measure dose. Supine and prone arcs for each patient were consecutively delivered, and the aggregate dose at each point was compared to the planned dose calculated in the phantom. Results: The ion chamber measurements differed from the planned dose by an average of .5%, with a standard deviation of 2.1%. All measured data for the MOSFETS were within 10% of the corresponding planned dose except for two outlying points. The standard deviation of dose differences across the entire cohort was 4.0%. Most significant discrepancies occurred either in inhomogeneous regions with large gradients, or at inferior points where beam angle was steepest. Conclusion: We have confirmed that the planned dose is well matched to our measurements within 10% for this method of planning and delivery. We are currently incorporating this technique into our clinical workflow. This work is supported by RaySearch.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Tests of neutrino interaction models with the MicroBooNE detector

I measure a large set of observables in inclusive charged current muon neutrino scattering on argon with the MicroBooNE liquid argon time projection chamber operating at Fermilab. I evaluate three neutrino interaction models based on the widely used GENIE event generator using these observables. The measurement uses a data set consisting of neutrino interactions with a final state muon candidate fully contained within the MicroBooNE detector. These data were collected in 2016 with the Fermilab Booster Neutrino Beam, which has an average neutrino energy of $800$ MeV, using an exposure corresponding to $5.0\times10^{19}$ protons-on-target. The analysis employs fully automatic event selection and charged particle track reconstruction and uses a data-driven technique to separate neutrino interactions from cosmic ray background events. I find that GENIE models consistently describe the shapes of a large number of kinematic distributions for fixed observed multiplicity, but I show an indication that the observed multiplicity fractions deviate from GENIE expectations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

SPECKLE INTERFEROMETRY OF SECONDARY COMPONENTS IN NEARBY VISUAL BINARIES

Statistical characterization of secondary subsystems in binaries helps to distinguish between various scenarios of multiple-star formation. The Differential Speckle Survey Instrument was used at the Gemini-N telescope for several hours in 2015 July to probe the binarity of 25 secondary components in nearby solar-type binaries. Six new subsystems were resolved, with meaningful detection limits for the remaining targets. The large incidence of secondary subsystems agrees with other similar studies. The newly resolved subsystem HIP 115417 Ba,Bb causes deviations in the observed motion of the outer binary from which an astrometric orbit of Ba,Bb with a period of 117 years is deduced.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

Effect of wind turbulence and shear on landing performance of jet transports

Computer simulations of a Boeing 727 class aircraft landing in turbulence were developed by programming the longitudinal aircraft equations of motion into a digital computer with various input values of vertical and horizontal wind speeds. Turbulent wind data was fed to the computer in one-second intervals. The computer computed in one-second intervals the aircraft speed, altitude, horizontal distance traveled, rate-of-descent, pitch attitude, glide path angle (from edge of runway) and elevator angle. All computer runs were made in the 'stick-fixed' mode. The RMS values of altitude and velocity perturbations (from equilibrium) were found to be large when horizontal wind gusts had sinusoidal components at or near the phugoid (long period) frequency. Maximum RMS altitude deviations occurred when the vertical wind had sinusoidal components which were 1/10 to 1/5 of the phugoid frequency. When real wind data (obtained from NCAR Queen Air) were used as input winds good correlations were found to exist between RMS velocity perturbations and both horizontal and vertical wind shears.

Blick, E. F.↗
NASA NTRS1991

Jitter model and signal processing techniques for pulse width modulation optical recording

A jitter model and signal processing techniques are discussed for data recovery in Pulse Width Modulation (PWM) optical recording. In PWM, information is stored through modulating sizes of sequential marks alternating in magnetic polarization or in material structure. Jitter, defined as the deviation from the original mark size in the time domain, will result in error detection if it is excessively large. A new approach is taken in data recovery by first using a high speed counter clock to convert time marks to amplitude marks, and signal processing techniques are used to minimize jitter according to the jitter model. The signal processing techniques include motor speed and intersymbol interference equalization, differential and additive detection, and differential and additive modulation.

Liu, Max M.-K.↗
NASA NTRS1993

A zonal CFD method for three-dimensional wing simulations

The primary objective of this work is to demonstrate the feasibility of a 3D potential/viscous flow coupling procedure for reducing computational effort while maintaining solution accuracy. The closed-loop, overlapped, velocity-coupling concept has been developed in a new code, ZAP3D, that couples a potential flow panel code with a Navier-Stokes method. The current ZAP3D calculation for an aspect ratio 5 wing with an outer domain radius of about 1.2 chords represents a speed-up in CPU time over the ARC3D large domain calculation by about a factor of 2.5. This improvement is achieved for less than a 0.5 percent deviation in C(L), 10 counts change in C(D), and 0.0015 variation in C(My). Additional reductions in the required computational domain for ZAP3D are expected as the method is further developed and refined.

Summa, J. M.↗
NASA NTRS2022

KMT-2021-BLG-0912Lb: A Microlensing Super Earth Around A K-Type Star

Aims . The light curve of the microlensing event KMT-2021-BLG-0912 exhibits a very short anomaly relative to a single-lens single-source form. We investigate the light curve for the purpose of identifying the origin of the anomaly. Methods . We model the light curve under various interpretations. From this, we find four solutions, in which three solutions are found under the assumption that the lens is composed of two masses (2L1S models), and the other solution is found under the assumption that the source is comprised of binary stars (1L2S model). The 1L2S model is ruled out based on the contradiction that the faint source companion is bigger than its primary, and one of the 2L1S solutions is excluded from the combination of the poorer fit, blending constraint, and lower overall probability, leaving two surviving solutions with the planet/host mass ratios of q ∼ 2.8 × 10 −5 and ∼1.1 × 10 −5 . A subtle central deviation supports the possibility of a tertiary lens component, either a binary companion to the host with a very large or small separation, or a second planet lying near the Einstein ring, but it is difficult to claim a secure detection due to the marginal improvement of the fit, lack of consistency among different data sets, and difficulty in uniquely specifying the nature of the tertiary component. Results . With the observables of the event, it is estimated that the masses of the planet and host are ∼(6.9 M ⊕ , 0.75 M ⊙ ) according to one solution and ∼(2.8 M ⊕ , 0.80 M ⊙ ) according to the other, indicating that the planet is a super Earth around a K-type star, regardless of the solution. The fact that 16 (including the one reported in this work) out of 19 microlensing planets with M ≲ 10 M ⊕ were detected during the last 6 yr nicely demonstrates the importance of high-cadence global surveys in detecting very low-mass planets.

gravitational lensing↗
DOE OSTI2020

Characterizing Subsiding Shells in Shallow Cumulus Using Doppler Lidar and Large‐Eddy Simulation

Abstract The existence of subsiding shells on the periphery of shallow cumulus clouds has major implications concerning the parameterization of shallow convection, with the mass exchange between the shell and cloudy air representing a significant deviation from the commonly used bulk‐plume parameterization. We examine the structure and frequency of subsiding shells in shallow cumulus convection using Doppler lidars at the Atmospheric Radiation Measurement Southern Great Plains facility in the central United States and at the Jülich ObservatorY for Cloud Evolution in western Germany. Doppler lidar indicates that the vertical subsiding shell extent is asymmetric, while shell width is typically ~100 m. Large‐eddy simulation can reasonably simulate the observed shell structure using a grid spacing of 10 m and suggests that much of the observed asymmetry is not a result of transient cloud evolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A large interlaboratory electron diffraction study of monolayer graphene

Standardisation of data collection and analysis is essential to enable commercialisation of 2D materials in a wide range of technologies. Selected area electron diffraction (SAED) in the transmission electron microscope (TEM) is one of the key methods for distinguishing monolayer from bilayer and few-layer graphene by comparing the 1st and 2nd order diffraction spot intensities. Yet there are many factors that can affect the reliability of data collection and interpretation, causing the measurement of monolayer samples to deviate from the literature boundary condition of $I_{\{\bar{2}110\}}$$/$$I_{\{1\bar{1}00\}}$ < 1 for monolayer graphene (1LG). Here we present the results of a large interlaboratory SAED comparison study, where 15 international laboratories measured and analysed nominally identical samples of chemical vapour deposited graphene. Large variations were observed in the measured ratios of diffraction spot intensities, with the largest variance associated with poor quality SAED data resulting from inadequate specimen handling and storage. To inform the reliable determination of monolayer thickness from SAED patterns we provide a description of best practice for specimen handling, TEM operation, data collection and analysis. This work was undertaken within VAMAS Technical Working Area 41: Graphene and related 2D materials—Project 9, the results of which have been directly incorporated into ISO/TS 21356–2 for the characterisation of graphene sheets. We find that when this methodology is followed, 1LG can be distinguished from bilayer or thicker material with high confidence where analysis of a single SAED pattern gives $I_{\{\bar{2}110\}}$$/$$I_{\{1\bar{1}00\}}$ < 1.2, even in the absence of precise specimen tilting.

2D materials↗
DOE OSTI2016

Lensing corrections to the E {sub g} ( z ) statistics from large scale structure

We study the impact of the often neglected lensing contribution to galaxy number counts on the E {sub g} statistics which is used to constrain deviations from GR. This contribution affects both the galaxy-galaxy and the convergence-galaxy spectra, while it is larger for the latter. At higher redshifts probed by upcoming surveys, for instance at z = 1.5, neglecting this term induces an error of (25–40)% in the spectra and therefore on the E {sub g} statistics which is constructed from the combination of the two. Moreover, including it, renders the E {sub g} statistics scale and bias-dependent and hence puts into question its very objective.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Search for contact interactions and large extra dimensions in the dilepton channel using proton–proton collisions at √s = 8 TeV with the ATLAS detector

Research is conducted for non-resonant new phenomena in dielectron and dimuon final states, originating from either contact interactions or large extra spatial dimensions. The LHC 2012 proton–proton collision dataset recorded by the ATLAS detector is used, corresponding to 20 fb –1 at √s = 8 TeV. The dilepton invariant mass spectrum is a discriminating variable in both searches, with the contact interaction search additionally utilizing the dilepton forward-backward asymmetry. No significant deviations from the Standard Model expectation are observed. Lower limits are set on the ℓℓqq contact interaction scale Λ between 15.4 TeV and 26.3 TeV, at the 95% credibility level. For large extra spatial dimensions, lower limits are set on the string scale MS between 3.2 TeV to 5.0 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Variance reduction for Fokker–Planck based particle Monte Carlo schemes

Recently, Fokker–Planck based particle Monte Carlo schemes have been proposed and evaluated for simulations of rarefied gas flows [1–3]. In this paper, the variance reduction for particle Monte Carlo simulations based on the Fokker–Planck model is considered. First, deviational based schemes were derived and reviewed, and it is shown that these deviational methods are not appropriate for practical Fokker–Planck based rarefied gas flow simulations. This is due to the fact that the deviational schemes considered in this study lead either to instabilities in the case of two-weight methods or to large statistical errors if the direct sampling method is applied. Motivated by this conclusion, we developed a novel scheme based on correlated stochastic processes. The main idea here is to synthesize an additional stochastic process with a known solution, which is simultaneously solved together with the main one. By correlating the two processes, the statistical errors can dramatically be reduced; especially for low Mach numbers. To assess the methods, homogeneous relaxation, planar Couette and lid-driven cavity flows were considered. For these test cases, it could be demonstrated that variance reduction based on parallel processes is very robust and effective.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

PolyPole-1: An accurate numerical algorithm for intra-granular fission gas release

This paper describes the development of a new numerical algorithm (called PolyPole-1) to efficiently solve the equation for intra-granular fission gas release in nuclear fuel. The work was carried out in collaboration with Politecnico di Milano and Institute for Transuranium Elements. The PolyPole-1 algorithms is being implemented in INL's fuels code BISON code as part of BISON's fission gas release model. The transport of fission gas from within the fuel grains to the grain boundaries (intra-granular fission gas release) is a fundamental controlling mechanism of fission gas release and gaseous swelling in nuclear fuel. Hence, accurate numerical solution of the corresponding mathematical problem needs to be included in fission gas behaviour models used in fuel performance codes. Under the assumption of equilibrium between trapping and resolution, the process can be described mathematically by a single diffusion equation for the gas atom concentration in a grain. In this work, we propose a new numerical algorithm (PolyPole-1) to efficiently solve the fission gas diffusion equation in time-varying conditions. The PolyPole-1 algorithm is based on the analytic modal solution of the diffusion equation for constant conditions, with the addition of polynomial corrective terms that embody the information on the deviation from constant conditions. The new algorithm is verified by comparing the results to a finite difference solution over a large number of randomly generated operation histories. Furthermore, comparison to state-of-the-art algorithms used in fuel performance codes demonstrates that the accuracy of the PolyPole-1 solution is superior to other algorithms, with similar computational effort. Finally, the concept of PolyPole-1 may be extended to the solution of the general problem of intra-granular fission gas diffusion during non-equilibrium trapping and resolution, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Structure of β-AgGaO{sub 2}; ternary I–III–VI{sub 2} oxide semiconductor with a wurtzite-derived structure

The structure of the wurtzite-derived β-AgGaO{sub 2} was refined by Rietveld analysis of high-resolution powder diffraction data obtained using synchrotron X-ray radiation. The space group of the crystal is Pna2{sub 1} with lattice parameters of a{sub 0}=5.56175 Å, b{sub 0}=7.14749 Å, and c{sub 0}=5.46875 Å. The deviation of O–Ag–O and M–O–M bond angles from the regular tetrahedral angle of 109.5° was very large at ∼8° and ∼11°, respectively. The electronic structure of β-AgGaO{sub 2} is discussed based on its structure, and the indirect band gap of β-AgGaO{sub 2} was related to significant tetrahedral distortion. Although β-AgGaO{sub 2} decomposes into metallic silver and Ga{sub 2}O{sub 3} at a high temperature in any atmosphere, β-AgGaO{sub 2} is stable up to 690 °C under an O{sub 2} atmosphere. No direct transformation from the wurtzite-derived phase to a delafossite phase occurs in β-AgGaO{sub 2}. - Graphical abstract: Crystal structure of β-AgGaO{sub 2} was refined by Rietveld analysis. AgO{sub 4} and O(Ag,Ga){sub 4} tetrahedra are significantly distorted from ideal tetrahedron. - Highlights: • Orthorhombic β-AgGaO{sub 2} with a wurtzite-derived β-NaFeO{sub 2} structure was synthesized. • Its structure was refined by Rietveld analysis of high-resolution XRD data. • Silver and oxygen tetrahedra are significantly distorted from an ideal tetrahedron. • The extent of this tetrahedral distortion is related to the band gap nature. • β-AgGaO{sub 2} is a metastable phase but is stable up to 690 °C in an O{sub 2} atmosphere.

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

Measurement of top-quark pair production in association with charm quarks in proton–proton collisions at √s = 13 TeV with the ATLAS detector

Inclusive cross-sections or top-quark pair production in association with charm quarks are measured with proton-proton collision data at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 140 fb -1 , collected with the ATLAS experiment at LHC between 2015 and 2018. The measurements are performed by requiring one or two charged leptons (electrons and muons), two b-tagged jets, and at least one additional jet in the final state. A custom flavor-tagging algorithm is employed for the simultaneous identification of b-jets and c-jets. In a fiducial phase space that replicates the acceptance of the ATLAS detector, the cross-sections for $t\bar{t}$ + ≥ 2c and $t\bar{t}$ + 1c production are measured to be $1.28^{+0.27}_{-0.24}$ pb and $6.4^{+1.0}_{-0.9}$ pb, respectively. The measurements are primarily limited by uncertainties in the modeling of inclusive $t\bar{t}$ and $t\bar{t}$ + $b\bar{b}$ production, in the calibration of the flavor-tagging algorithm, and by data statistics. Cross-section predictions from various $t\bar{t}$ simulations are largely consistent with the measured cross-section values, though all underpredict the observed values by 0.5 to 2.0 standard deviations. In a phase-space volume without requirements on the $t\bar{t}$ decay products and the jet multiplicity, the cross-section ratios of $t\bar{t}$ + ≥ 2c and $t\bar{t}$ + 1c to total $t\bar{t}$ + jets production are determined to be (1.23 ± 0.25)% and (8.8 ± 1.3)%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗