SEARCH · Engineering Papers
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.
A method is proposed for determining the load in the gear engagements of a multiflux planetary gear with flexible carrier jaws. The influence of the design and manufacturing errors of the mechanism on the pliability of its components and the nonuniformity of the load distribution over the satellites is determined.
Rigorous testing methods are essential for ensuring the security and reliability of industrial controller software. Fuzzing, a technique that automatically discovers software bugs, has also proven effective in finding software vulnerabilities. Unsurprisingly, fuzzing has been applied to a wide range of platforms, including programmable logic controllers (PLCs). However, current approaches, such as coverage-guided evolutionary fuzzing implemented in the popular fuzzer American Fuzzy Lop Plus Plus (AFL++), are often inadequate for finding logical errors and bugs in PLC control logic applications. They primarily target generic programming languages like C/C++, Java, and Python, and do not consider the unique characteristics and behaviors of PLCs, which are often programmed using specialized programming languages like Structured Text (ST). Furthermore, these fuzzers are ill suited to deal with complex input structures encapsulated in ST, as they are not specifically designed to generate appropriate input sequences. This renders the application of traditional fuzzing techniques less efficient on these platforms. To address this issue, this paper presents a fuzzing framework designed explicitly for PLC software to discover logic bugs in applications written in ST specified by the IEC 61131-3 standard. The proposed framework incorporates a custom-tailored PLC runtime and a fuzzer designed for the purpose. We demonstrate its effectiveness by fuzzing a collection of ST programs that were crafted for evaluation purposes. We compare the performance against a popular fuzzer, namely, AFL++. The proposed fuzzing framework demonstrated its capabilities in our experiments, successfully detecting logic bugs in the tested PLC control logic applications written in ST. On average, it was at least 83 times faster than AFL++, and in certain cases, for example, it was more than 23,000 times faster.
This research project has had the goal of gaining a better understanding of the physics of electrons bound to the surface of superfluid helium from both experimental and theoretical perspectives. It has particularly been aimed at two areas which had not been well studied: the relaxation and decoherence of the spin of the electrons on the helium surface and how the properties of underlying metallic layers affect the behavior of the electrons when the helium covering the metal is thin. This work is motivated in part by interest in using the spin of these electrons as a quantum bit, or qubit. Low levels of decoherence are advantageous for qubits, and moving the electrons, as one might do in a quantum processor, will be easiest if thin helium films can be employed. It had been suggested that spin decoherence should be very weak for electrons bound to superfluid He, but before this work there have been no quantitative studies of spin relaxation and decoherence. It is especially important to know how moving the electrons across the helium surface would affect their spin coherence. Calculations performed as part of this project show that the Rashba effective magnetic field, the mechanism which limits the spin coherence of mobile electrons in silicon-based devices (an actively pursued qubit technology), is exceptionally weak for electrons bound to helium. This project has identified other decoherence mechanisms which are stronger, but still weak compared to analogous silicon-based structures. Calculated spin coherence times for mobile electrons approach one day, as compared to microseconds in silicon. With coherence times of this magnitude, the spin qubit errors on helium will be completely dominated by errors in the quantum gates. In related work, the possibility of using an artificial spin-orbit interaction (a gradient magnetic field) for quantum operations on the electrons spins was considered. The calculations show that a moderate gradient field, small enough to be generated by a narrow superconducting wire, will enable high-fidelity quantum operations on electrons held in lithographically-defined quantum dots by driving them with a microwave electric field. The spin and motional coherence of the electrons is sufficient to allow high-fidelity 2-qubit quantum operations between electrons in neighboring quantum dots. As an outgrowth of experiments aiming to measure electron spin coherence it was discovered that very high densities of electrons can be stably supported on thin helium films coating ultra-smooth amorphous metallic layers. The measured densities are high enough that the electron system has almost certainly transitioned from an ordered array of electrons, known as a Wigner crystal (ordered by the electrons’ mutual repulsion), to a quantum fluid known as a Fermi liquid. This transition has been a subject of intense interest for over 40 years, since the electron Wigner crystal was first observed with electrons bound to superfluid helium, but it has never been unambiguously observed. Experiments are still underway in these new structures to definitively determine whether true quantum melting of the Wigner crystal has been demonstrated. This work has also catalyzed the development of a new approach for measuring the transport of electrons across very thin helium films, as will be needed for some of the quantum computing applications. The high electron density experiments as well as experiments with electrons bound in quantum dots have led to new techniques which may enable spin coherence measurements.
In this work we provide a unified, comprehensive treatment of all operators that contribute to the antiferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We use the Hybrid Monte Carlo algorithm to perform a systematic, carefully controlled analysis in the temporal Trotter error and of the thermodynamic limit. We expect our findings to improve the consistency of Monte Carlo determinations of critical exponents. We perform a data collapse analysis and determine the critical exponent β = 0.898 (37) for the semimetal-Mott insulator transition in the hexagonal Hubbard Model. Our methods are applicable to a wide range of lattice theories of strongly correlated electrons.
Introduction: Measurement errors are ubiquitous in epidemiological studies, especially for doses arising from environmental and occupational exposures, which are difficult to assess. Correction methods for independent measurement errors, have been comprehensively described elsewhere, including regression calibration, simulation extrapolation (SIMEX), etc. However, when exposure estimates are constructed using complicated physical and biological models, uncertainties in the dosimetry system can become very complex. In such systems, uncertainties in some parameters may affect a large group of study participants simultaneously, and hence are “shared” as opposed to independent dosimetric uncertainties. Representing the structure of the uncertainty in the dosimetry by providing multiple realizations of possible dose has been described for several studies. The goal of this paper is to implement, and explore the properties of, a generalized estimating equation-based approach to shared uncertainties based on principles described in an earlier paper. Methods: The effects of shared and unshared dosimetry error on parameter estimates and confidence intervals in a generalized excess relative risk model are examined. A post-hoc correction method is implemented that constructs corrected confidence intervals based on an underlying variance structure described in an earlier paper. In this correction method it is assumed that the estimated slope parameter in a linear excess relative risk model is distributed as a mixture of normal and lognormal components. A simulation study of the behavior of this correction was developed based on multiple realizations provided by the Monte Carlo dosimetry system used for the plutonium exposures in the Mayak Worker Cohort and a risk model with close similarity to that used in analyses of lung cancer data in that cohort. Results: In the simulation study we found that shared dose errors severely degraded the coverage probabilities of confidence limits for the slope parameter in the linear excess relative risk model relating uncertain exposure to lung cancer risk. Using our correction method to adjust the confidence intervals provided much improved coverage probabilities for this parameter. A smaller amount of degradation in confidence interval coverage for other parameters (e.g. those in the baseline and dose modification terms) was noted, and the proposed correction method also improved coverage for these parameters as well. Discussion: Compared to other recently described methods our proposed correction method relies on relatively familiar software to fit dose response models. The basic modeling philosophy described here first uses existing Poisson regression methods to estimate parameters in the excess relative risk model; in this first stage mean dose (computed from the multiple realizations) is used as the independent variable in the regression. Once the model is fit we then correct (widen) the standard errors for all parameters using the proposed correction method. The simulation studies indicate that the methods described here are an important contribution to the available tools for dealing with shared measurement errors in complex risk analysis.
The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. Here, we apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value.
Using a mixed-ligand synthetic scheme, we create a family of quasi-two-dimensional antiferromagnets, namely, [Cu(HF 2 )(pyz) 2 ]ClO 4 [pyz = pyrazine], [CuL 2 (pyz) 2 ](ClO 4 ) 2 [L = pyO = pyridine-N-oxide and 4-phpy-O = 4-phenylpyridine-N-oxide. These materials are shown to possess equivalent two-dimensional [Cu(pyz) 2 ] 2+ nearly square layers, but exhibit interlayer spacings that vary from 6.5713 to 16.777 Å, as dictated by the axial ligands. We present the structural and magnetic properties of this family as determined via x-ray diffraction, electron-spin resonance, pulsed- and quasistatic-field magnetometry and muon-spin rotation, and compare them to those of the prototypical two-dimensional magnetic polymer Cu(pyz) 2 (ClO 4 ) 2 . We find that, within the limits of the experimental error, the two-dimensional, intralayer exchange coupling in our family of materials remains largely unaffected by the axial ligand substitution, while the observed magnetic ordering temperature (1.91 K for the material with the HF 2 axial ligand, 1.70 K for the pyO and 1.63 K for the 4-phpy-O) decreases slowly with increasing layer separation. Despite the structural motifs common to this family and Cu(pyz) 2 (ClO 4 ) 2 , the latter has significantly stronger two-dimensional exchange interactions and hence a higher ordering temperature. Here, we discuss these results, as well as the mechanisms that might drive the long-range order in these materials, in terms of departures from the ideal S = 1/2 two-dimensional square-lattice Heisenberg antiferromagnet. In particular, we find that both spin-exchange anisotropy in the intralayer interaction and interlayer couplings (exchange, dipolar, or both) are needed to account for the observed ordering temperatures, with the intralayer anisotropy becoming more important as the layers are pulled further apart.
The representative elementary volume (REV) is a fundamental property of a material, but no direct measurements exist for cementitious materials. In this paper, the REV of cement pastes with supplementary cementitious materials (GGBS, PFA, SF) was determined by analysing the three-dimensional pore structure (> 0.2 μm) using laser scanning confocal microscopy (LSCM). The effect of axial distortion inherent to LSCM on 3D pore structure was also investigated. A range of 3D pore parameters was measured using skeletonisation, maximal ball and random walker algorithms. Results show that axial distortion has insignificant effects on most parameters except Euler connectivity, average pore and throat volumes and directional diffusion tortuosities. Most pore parameters become independent of sampling volume at ≈ 60{sup 3} μm{sup 3} except diffusion tortuosities and formation factor. The REV for porosity calculated based on a statistical approach at eight realisations and 5% relative error was found to be ≈ 100{sup 3} μm{sup 3}.
Abstract High-fidelity analyses and forecasts of integrated vapor transport (VT) are central to the study of Earth’s hydrological cycle as well as high-impact phenomena such as monsoons and atmospheric rivers. The impact of the in-line analysis correction-based additive inflation (ACAI) on IVT biases and forecast errors is examined within the Navy Earth System Prediction Capability (Navy ESPC) global coupled system. The ACAI technique uses atmospheric analysis corrections from the data assimilation system to approximate model bias and as a representation of stochastic model error to simultaneously reduce systematic and random errors and improve ensemble performance. ACAI reduces the global average magnitude of the 7- and 14-day IVT bias by 16%–17% during Northern Hemisphere summer, reaching 70% reductions in some tropical regions. The global average IVT bias reduction is similar to the bias reduction for low-level wind speed bias and considerably smaller than the bias reduction in total precipitable water. The localized regions where ACAI increases IVT bias occur where the control IVT biases change sign and structure with increasing forecast lead time, such as the South Asian monsoon region. Substituting analyzed wind or moisture fields for the forecast fields when calculating the forecast IVT confirms that, on average, wind errors dominate the IVT error calculation in the tropics, although wind and moisture error contributions are comparable in the extratropics. The existence of regions where using either analyzed winds or analyzed moisture increases IVT bias or mean absolute error reveals areas with compensating errors.
The Global Ozone Monitoring Experiment (GOME) is a hyper-spectral satellite instrument measuring the ultraviolet backscatter at relatively high spectral resolution. GOME radiances have been slit averaged to emulate measurements of the Total Ozone Mapping Spectrometer (TOMS) made at discrete wavelengths and processed using the new TOMS Version 8 Ozone Algorithm. Compared to Differential Optical Absorption Spectroscopy (DOAS) techniques based on local structure in the Huggins Bands, the TOMS uses differential absorption between a pair of wavelengths including the local stiucture as well as the background continuum. This makes the TOMS Algorithm more sensitive to ozone, but it also makes the algorithm more sensitive to instrument calibration errors. While calibration adjustments are not needed for the fitting techniques like the DOAS employed in GOME algorithms, some adjustment is necessary when applying the TOMS Algorithm to GOME. Using spectral discrimination at near ultraviolet wavelength channels unabsorbed by ozone, the GOME wavelength dependent calibration drift is estimated and then checked using pair justification. In addition, the day one calibration offset is estimated based on the residuals of the Version 8 TOMS Algorithm. The estimated drift in the 2b detector of GOME is small through the first four years and then increases rapidly to +5% in normalized radiance at 331 nm relative to 385 nm by mid 2000. The lb detector appears to be quite well behaved throughout this time period.
Two fundamental approximations in macroscale solid-mechanics modeling are (1) the assumption of scale separation in homogenization theory and (2) the use of a macroscopic plasticity material model that represents, in a mean sense, the multitude of inelastic processes occurring at the microscale. With the goal of quantifying the errors induced by these approximations on engineering quantities of interest, we perform a set of direct numerical simulations (DNS) in which polycrystalline microstructures are embedded throughout a macroscale structure. The largest simulations model over 50,000 grains. The microstructure is idealized using a randomly close-packed Voronoi tessellation in which each polyhedral Voronoi cell represents a grain. An face centered cubic crystal-plasticity model is used to model the mechanical response of each grain. The overall grain structure is equiaxed, and each grain is randomly oriented with no overall texture. The detailed results from the DNS simulations are compared to results obtained from conventional macroscale simulations that use homogeneous isotropic plasticity models. The macroscale plasticity models are calibrated using a representative volume element of the idealized microstructure. Furthermore, we envision that DNS modeling will be used to gain new insights into the mechanics of material deformation and failure.
KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS G12C —sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS G12C . In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS G12C , we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.
Many radioactive materials within the nuclear fuel cycle present a significant hazard. Such materials include spent nuclear fuel, high-level waste, legacy waste and other nuclear materials. These materials are often held in long-term storage as an interim stage within their lifecycle, including reuse or disposal. Safety and security of spent nuclear fuel (SNF) interim storage installations and follow-on SNF transportations to its final repository sites are very important, due to a great concentration of fission products, actinides and activation products. Fracture mechanics approach in applying to SNF system reliability investigation, especially for the high burn up (HBU) SNF, during SNF long term dry storage or SNF transportation is necessary due to the inherited flaws and inhomogeneity structures existed in a SNF system from the reactor operation. For instance, such as hydride and oxide formation, surface flaw induced by SNF assembly contact interactions, or internal flaws induced by pellet-clad mechanical interaction (PCMI). However, none of the existing fracture toughness data deal with fuel cladding specific geometry or SNF material conditions, such as cladding structure with the segment fuel pellets configuration and the PCMI mechanism. Consequently, the application of existing fracture toughness data to SNF system reliability investigation can be error prone; thus, the development of an industry-wide consensus fracture mechanics approach is needed. The objective of this research is to develop an in-situ fracture mechanics testing protocol, including the associated analytical procedure, that is suitable for evaluating the SNF fracture toughness. Fracture toughness data was obtained under quasi-static fracture loading using Oak Ridge National Laboratory (ORNL) developed Spiral Notch Torsion Test (SNTT) technology carried out on a biaxial tension/torsion tester. These data will be used to support SNF reliability investigation during long-tern SNF dry storage stage or in the follow-on SNF system transport to the final repository sites. SNTT has been a recent breakthrough in measuring fracture toughness for different materials, including metals, ceramics, concrete and polymers composites. Due to its high geometry constraint and unique loading condition, SNTT can be used to measure the fracture toughness with smaller specimens without concern of size effects. The application of SNTT to brittle materials has been proved to be successful. The micro-cracks induced by original notches in brittle materials could ensure crack growth in SNTT samples. Therefore, no fatigue pre cracks are needed for brittle materials specimens. The application of SNTT to the ductile material to generate valid toughness data will require a test sample with sufficient pre-crack length to increase the sample localized constraint at crack front. Fatigue pre-crack growth techniques with the associated compliance function estimated was developed for estimating the crack penetration depth to monitor the fatigue crack growth evolution. In order to extend SNTT approach to a thin shell cladding structure material, a new testing protocol and the associated analytical procedure for estimating the SNF fracture toughness was developed. The detailed SNTT approach and its estimated fracture toughness for the Zr-4 cladding structure with segment alumina inserts are presented in this report. The SNTT test results indicate that SNTT method is a reliable test approach with good repeatability in applying to Zr-4 cladding material. For the medium and the short crack length, the estimate J Q upon fracture for the baseline Zr-4 cladding is at 285.7 lb/in with 2-sigma uncertainty of 18.59 lb/in., and the associated K Q is at 61.4 Ksi√in. For a long crack length, the crack initialization orientation is apparently deviated from the principle tensile stress contour – 45° spiral crack front (i.e., Mode I + Mode III failure mechanism); and the estimated J Q is at 108 lb/in., the associated K Q is at 37.7 Ksi√in. The estimate Mode I J Q (along the 45° pitch principle tensile stress contour) is at 200 lb/in.; which indicate a significant reduction in the fracture toughness in a tubing structure under mixed mode loading condition, Mode I + Mode III, compared to that of Mode I (tensile stress dominated) alone.
Preliminary design work on the New Generation Space Telescope (NGST) is currently under way. This telescope is envisioned as a lightweight, deployable Cassegrain reflector with an aperture of 8 meters, and an effective focal length of 80 meters. It is to be folded into a small-diameter package for launch by an Atlas booster, and unfolded in orbit. The primary is to consist of an octagon with a hole at the center, and with eight segments arranged in a flower petal configuration about the octagon. The comers of the petal-shaped segments are to be trimmed so that the package will fit atop the Atlas booster. This mirror, along with its secondary will focus the light from a point source into an image which is spread from a point by diffraction effects, figure errors, and scattering of light from the surface. The distribution of light in the image of a point source is called a point spread function (PSF). The obstruction of the incident light by the secondary mirror and its support structure, the trimmed corners of the petals, and the grooves between the segments all cause the diffraction pattern characterizing an ideal point spread function to be changed, with the trimmed comers causing the rings of the Airy pattern to become broken up, and the linear grooves causing diffraction spikes running radially away from the central spot, or Airy disk. Any figure errors the mirror segments may have, or any errors in aligning the petals with the central octagon will also spread the light out from the ideal point spread function. A point spread function for a mirror the size of the NGST and having an incident wavelength of 900 nm is considered. Most of the light is confined in a circle with a diameter of 0.05 arc seconds. The ring pattern ranges in intensity from 10(exp -2) near the center to 10(exp -6) near the edge of the plotted field, and can be clearly discerned in a log plot of the intensity. The total fraction of the light scattered from this point spread function is called the total integrated scattering (TIS), and the fraction remaining is called the Strehl ratio. The angular distribution of the scattered light is called the angle resolved scattering (ARS), and it shows a strong spike centered on a scattering angle of zero, and a broad , less intense distribution at larger angles. It is this scattered light, and its effect on the point spread function which is the focus of this study.
Radiation damage behavior of SiC depends on its lattice structure. Here, we use a pulsed-ion-beam method to study defect interaction dynamics in 6H-SiC and damage buildup in 4H- and 6H-SiC irradiated at 100 °C with 500 keV Ar ions. These results are compared with previously reported data for Ar-ion-irradiated 3C- and 4H-SiC. We find that, for these irradiation conditions, damage buildup in 3C- and 6H-SiC is statistically indistinguishable and is significantly more efficient than in 4H-SiC. Within a stimulated amorphization model, 4H-SiC is described by a reduced amorphization cross-section constant, while the point defect cluster production cross-section is the same (within experimental errors) for the three polytypes studied. Moreover, 4H-SiC exhibits slower defect relaxation dynamics than 3C- and 6H-SiC. These results clearly demonstrate the importance of the lattice structure in damage buildup and defect interaction dynamics in SiC.
Quantum annealing is a type of analog computation that aims to use quantum mechanical fluctuations in search of optimal solutions of QUBO (quadratic unconstrained binary optimization) or, equivalently, Ising problems. Since NP-hard problems can in general be mapped to Ising and QUBO formulations, the quantum annealing paradigm has the potential to help solve various NP-hard problems. Current quantum annealers, such as those manufactured by D-Wave Systems, Inc. have various practical limitations including the size (number of qubits) of the problem that can be solved, the qubit connectivity, and error due to the environment or system calibration, which can reduce the quality of the solutions. Typically, for an arbitrary problem instance, the corresponding QUBO (or Ising) structure will not natively embed onto the available qubit architecture on the quantum chip. Thus, in these cases, a minor embedding of the problem structure onto the device is necessary. However, minor embeddings on these devices do not always make use of the full sparse chip hardware graph, and a large portion of the available qubits stay unused during quantum annealing. In this work, we embed a disjoint random QUBO on the unused parts of the chip alongside the QUBO to be solved, which acts as an indicator of the solution quality of the device over time. Using experiments on three different D-Wave quantum annealers, we demonstrate that (i) long term trends in solution quality exist on the D-Wave device, and (ii) the unused qubits can be used to measure the current level of noise of the quantum system.
With the growing scale and complexity of high-performance computing (HPC) systems, resilience solutions that ensure continuity of service despite frequent errors and component failures must be methodically designed to balance the reliability requirements with the overheads to performance and power. Design patterns enable a structured approach to the development of resilience solutions, providing hardware and software designers with the building block elements for the rapid development of novel solutions and for adapting existing technologies for emerging, extreme-scale HPC environments. In this paper, we develop analytical models that enable designers to evaluate the reliability and performance characteristics of the design patterns. These models are particularly useful in building a unified framework that analyzes and compares various resilience solutions built using a combination of patterns.
Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.