Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observational methods”

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 1,135 records · Page 63

Refining Jets for CMS Run 3 using Fast Simulation

As the LHC moves into its high-luminosity phase, the CMS experiment must handle more complex data collected at much higher rates. While the Geant4-based simulation application (FullSim) provides highly accurate simulation to complement real data, FullSim’s intensive consumption of computing resources becomes an increasing liability as the rates increase, while faster tools offer an advantage. The fast MC production application (FastSim) delivers a complete simulation with a factor of 10 speedup over FullSim, but introduces inaccuracies in some observables. A specialized refinement method, Fast Perfekt, employs machine learning to improve the accuracy of FastSim. An initial report of this work focused on the refinement of jet flavor tagging observables. This article presents an update on the refinement, focusing on PUPPI jets with Run 3 data-taking conditions. Refinement is extended to include jet transverse momentum as well as its propagation to missing transverse momentum. A gridbased framework and real-time monitoring system have been developed to facilitate optimization and scaling of the refinement to a large number of target variables.

Güngördü, Açelya Deniz [Istanbul Tech. U.]↗

Spatio-Temporal Data Analysis at Scale Using Models Based on Gaussian Processes

Gaussian processes are the most commonly used statistical model for spatial and spatio-temporal processes that vary continuously. They are broadly applicable in the physical sciences and engineering and are also frequently used to approximate the output of complex computer models, deterministic or stochastic. We undertook research related to theory, computation, and applications of Gaussian processes as well as some work on estimating extremes of distributions for which a Gaussian process assumption might be inappropriate. Our theoretical contributions include the development of new classes of spatial-temporal covariance functions with desirable properties and new results showing that certain covariance models lead to predictions with undesirable properties. To understand how Gaussian process models behave when applied to deterministic computer models, we derived what we believe to be the first significant results on the large sample properties of estimators of parameters of Gaussian processes when the actual process is a simple deterministic function. Finally, we investigated some theoretical issues related to maxima of observations with varying upper bounds and found that, depending on the circumstances, standard large sample results for maxima may or may not hold. Our computational innovations include methods for analyzing large spatial datasets when observations fall on a partially observed grid and methods for estimating parameters of a Gaussian process model from observations taken by a polar-orbiting satellite. In our application of Gaussian process models to deterministic computer experiments, we carried out some matrix computations that would have been infeasible using even extended precision arithmetic by focusing on special cases in which all elements of the matrices under study are rational and using exact arithmetic. The applications we studied include total column ozone as measured from a polar-orbiting satellite, sea surface temperatures over the Pacific Ocean, and annual temperature extremes at a site in New York City. In each of these applications, our theoretical and computational innovations were directly motivated by the challenges posed by analyzing these and similar types of data.

97 MATHEMATICS AND COMPUTING↗

Coarse-Grained Reactive Molecular Dynamics Simulations of Heterogeneities in Shocked Energetic Materials (LDRD Final Report)

This report summarizes the result of LDRD project 16-0161, titled "Coarse-Grained Re- active Molecular Dynamics Simulations of Heterogeneities in Shocked Energetic Materials." The purpose of the project was to develop a coarse-grained reactive molecular dynamics capability in LAMMPS enabling simulations of initiation in energetic materials comparable in accuracy to what is currently possible using large-scale reactive molecular dynamics, but with greatly reduced computational cost. The starting point for this work was the reactive dissipative particle dynamics (DPD) approach, which has been implemented as the new USER-DPD package in LAMMPS by researchers at Army Research Laboratory. Using modified versions of the examples provided with the new package, we examined the computational efficiency of the method, as well as its ability to model energy release in energetic materials. We observed that the Shardlow splitting method provides a great speed and accuracy advantage over conventional velocity Verlet time integration. We observed that the generic model of an energetic material provided with the USER-DPD package exhibited incomplete reaction under constant volume conditions. This was caused by quenching of the internal temperature of the molecules due to a rapid build-up of repulsive interactions between product gas components under constant volume conditions. Under constant pressure conditions, complete reaction was observed, as volume expansion prevented the buildup of strong repulsive interactions. Finally, a more realistic model calibrated to reproduce the equation of state of the RDX molecular crystal was examined. This model exhibited much less quenching of the internal temperature under constant volume conditions and reacted very rapidly under constant pressure conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Low-pass spectral analysis of time-resolved serial femtosecond crystallography data

Low-pass spectral analysis (LPSA) is a recently developed dynamics retrieval algorithm showing excellent retrieval properties when applied to model data affected by extreme incompleteness and stochastic weighting. In this work, we apply LPSA to an experimental time-resolved serial femtosecond crystallography (TR-SFX) dataset from the membrane protein bacteriorhodopsin (bR) and analyze its parametric sensitivity. While most dynamical modes are contaminated by nonphysical high-frequency features, we identify two dominant modes, which are little affected by spurious frequencies. The dynamics retrieved using these modes shows an isomerization signal compatible with previous findings. We employ synthetic data with increasing timing uncertainty, increasing incompleteness level, pixel-dependent incompleteness, and photon counting errors to investigate the root cause of the high-frequency contamination of our TR-SFX modes. By testing a range of methods, we show that timing errors comparable to the dynamical periods to be retrieved produce a smearing of dynamical features, hampering dynamics retrieval, but with no introduction of spurious components in the solution, when convergence criteria are met. Using model data, we are able to attribute the high-frequency contamination of low-order dynamical modes to the high levels of noise present in the data. Finally, we propose a method to handle missing observations that produces a substantial dynamics retrieval improvement from synthetic data with a significant static component. Reprocessing of the bR TR-SFX data using the improved method yields dynamical movies with strong isomerization signals compatible with previous findings.

59 BASIC BIOLOGICAL SCIENCES↗

Inferring the gravitational potential of the Milky Way with a few precisely measured stars

The dark matter halo of the Milky Way is expected to be triaxial and filled with substructure. It is hoped that streams or shells of stars produced by tidal disruption of stellar systems will provide precise measures of the gravitational potential to test these predictions. We develop a method for inferring the Galactic potential with tidal streams based on the idea that the stream stars were once close in phase space. Our method can flexibly adapt to any form for the Galactic potential: it works in phase-space rather than action-space and hence relies neither on our ability to derive actions nor on the integrability of the potential. Our model is probabilistic, with a likelihood function and priors on the parameters. The method can properly account for finite observational uncertainties and missing data dimensions. We test our method on synthetic data sets generated from N-body simulations of satellite disruption in a static, multi-component Milky Way, including a triaxial dark matter halo with observational uncertainties chosen to mimic current and near-future surveys of various stars. We find that with just eight well-measured stream stars, we can infer properties of a triaxial potential with precisions of the order of 5%-7%. Without proper motions, we obtain 10% constraints on most potential parameters and precisions around 5%-10% for recovering missing phase-space coordinates. These results are encouraging for the goal of using flexible, time-dependent potential models combined with larger data sets to unravel the detailed shape of the dark matter distribution around the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of the Inclusive Electron Scattering Observables in the Resonance Region from the Experimental Data

The method for evaluation of inclusive electron scattering observables in the resonance region from the experimental data is presented. The computer codes in the web page allow user to evaluate inclusive electron scattering, virtual photon-proton cross sections and the F{sub 1} and F{sub 2} structure functions at the invariant masses of the final hadrons 1.0 < W < 4.0 GeV and at photon virtualities 0.5 < Q {sup 2} < 7.0 GeV{sup 2}. The developed approach is of interest for analyses of the experimental data on exploration of the hadron structure in the experiments with electromagnetic probes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrating AI with physics-based hydrological models and observations for insightinto changing climate and anthropogenic impacts

Focal Areas: Advanced computational methods that integrate AI, physics, and observations to provide predictive landscape hydrological modeling over large areas (regional, continental, worldwide) while incorporating increasingly available high-resolution data from drones, lidar and satellite. Science Challenge: Landscape data is available at finer scales than can be used in physics-based hydrological (PBH) models for regional or continental terrestrial water modeling. Thus, we throw away observable detail to achieve computability. We argue that integration of AI with PBH models and observed data can be used to provide upscaling for predictive models that are computable, retain physical conservation properties, and represent the fine-scale features that affect complex flow physics through both natural and urban environments. Developing such next-generation capabilities requires outside-the-box thinking that melds the different approaches of AI modeling, PBH modeling, and observation across multiple scales from local drones to satellites.

54 ENVIRONMENTAL SCIENCES↗

PREDICTING CORONAL MASS EJECTIONS USING MACHINE LEARNING METHODS

Of all the activity observed on the Sun, two of the most energetic events are flares and coronal mass ejections (CMEs). Usually, solar active regions that produce large flares will also produce a CME, but this is not always true. Despite advances in numerical modeling, it is still unclear which circumstances will produce a CME. Therefore, it is worthwhile to empirically determine which features distinguish flares associated with CMEs from flares that are not. At this time, no extensive study has used physically meaningful features of active regions to distinguish between these two populations. As such, we attempt to do so by using features derived from (1) photospheric vector magnetic field data taken by the Solar Dynamics Observatory ’s Helioseismic and Magnetic Imager instrument and (2) X-ray flux data from the Geostationary Operational Environmental Satellite’s X-ray Flux instrument. We build a catalog of active regions that either produced both a flare and a CME (the positive class) or simply a flare (the negative class). We then use machine-learning algorithms to (1) determine which features distinguish these two populations, and (2) forecast whether an active region that produces an M- or X-class flare will also produce a CME. We compute the True Skill Statistic, a forecast verification metric, and find that it is a relatively high value of ∼0.8 ± 0.2. We conclude that a combination of six parameters, which are all intensive in nature, will capture most of the relevant information contained in the photospheric magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

Observing programs in solar physics during the 1973 ATM Skylab program.

Description of the Apollo Telescope Mount (ATM), a set of spaceborne X-ray, ultraviolet, and visible-light telescopes intended to implement the extensive series of solar observations representing one of the major tasks assigned to the three-man crews who will perform three (respectively 28-, 56-, and 56-day long) manned orbital flights in mid-1973, known collectively as the Skylab mission. The ATM capabilities, its principal scientific goals, and planned observing program are reviewed, along with the methods to be used for coordinating the observations with ground-based and other kinds of solar observations.

Reeves, E. M.↗

Toward RADSCAT measurements over the sea and their interpretation

Investigations into several areas which are essential to the execution and interpretation of suborbital observations by composite radiometer - scatterometer sensor (RADSCAT) are reported. Experiments and theory were developed to demonstrate the remote anemometric capability of the sensor over the sea through various weather conditions. It is shown that weather situations found in extra tropical cyclones are useful for demonstrating the all weather capability of the composite sensor. The large scale fluctuations of the wind over the sea dictate the observational coverage required to correlate measurements with the mean surface wind speed. Various theoretical investigations were performed to establish a premise for the joint interpretation of the experiment data. The effects of clouds and rains on downward radiometric observations over the sea were computed. A method of predicting atmospheric attenuation from joint observations is developed. In other theoretical efforts, the emission and scattering characteristics of the sea were derived. Composite surface theories with coherent and noncoherent assumptions were employed.

Claassen, J. P.↗

Analysis of the Capability and Limitations of Relativistic Gravity Measurements Using Radio Astronomy Methods

The uses of radar observations of planets and very-long-baseline radio interferometric observations of extragalactic objects to test theories of gravitation are described in detail with special emphasis on sources of error. The accuracy achievable in these tests with data already obtained, can be summarized in terms of: retardation of signal propagation (radar), deflection of radio waves (interferometry), advance of planetary perihelia (radar), gravitational quadrupole moment of sun (radar), and time variation of gravitational constant (radar). The analyses completed to date have yielded no significant disagreement with the predictions of general relativity.

Shapiro, I. I.↗

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

Platform for immobilization and observation of subcellular processes

A method of immobilizing matter for imaging that includes providing an array of nanofibers and directing matter to the array of the nanofibers. The matter is immobilized when contacting at least three nanofibers of the array of nanofibers simultaneously. Adjacent nanofibers in the array of nanofibers may be separated by a pitch as great as 100 microns. The immobilized matter on the array of nanofibers may then be imaged. In some examples, the matter may be cell matter, such as protoplasts.

47 OTHER INSTRUMENTATION↗

Evaluation of Geologic Carbon Storage Resource Estimates (SREs) of Cambrian Ordovician Units within the CarbonSAFE Prefeasibility Study Region Subtasks 4.2 and 4.3 – Topical Report

This report by the Indiana Geological & Water Survey (IGWS) assessed the storage resource estimate (SRE) of Cambrian-Ordovician stratigraphic units located within the CarbonSAFE prefeasibility area within the Illinois Basin. A comprehensive data set of wireline logs and petrophysical information was used. These data include core analysis for porosity and permeability, and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for four potential storage reservoirs (Upper Ordovician Trenton Limestone/Black River Group, Middle Ordovician St. Peter Sandstone, Lower Ordovician and Upper Cambrian Knox Supergroup, including a separate assessment for the lower Knox Potosi Dolomite. The resultant SREs volumes are displayed using georeferenced maps and boxplots. The latter allow for comparing data statistics (mean values and variability) between five methods employed for the SRE calculations. Differences observed in SRE results from the five methods are mainly attributable to differences in the data used for porosity in each method. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗

Focused wave interactions with floating structures: a blind comparative study

The paper presents results from the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI) Blind Test Series 2. Without prior access to the physical data, participants, with numerical methods ranging from low-fidelity linear models to fully non-linear Navier-Stokes (NS) solvers, simulate the interaction between focused wave events and two separate, taut-moored, floating structures: a hemispherical-bottomed cylinder and a cylinder with a moonpool. The ‘blind’ numerical predictions for heave, surge, pitch and mooring load, are compared against physical measurements. Dynamic time warping is used to quantify the predictive capability of participating methods. In general, NS solvers and hybrid methods give more accurate predictions; however, heave amplitude is predicted reasonably well by all methods; and a WEC-Sim implementation, with CFD-informed viscous terms, demonstrates comparable predictive capability to even the stronger NS solvers. Large variations in the solutions are observed (even among similar methods), highlighting a need for standardisation in the numerical modelling of WSI problems.

42 ENGINEERING↗

INTERSTELLAR NEUTRAL HELIUM IN THE HELIOSPHERE FROM IBEX OBSERVATIONS. I. UNCERTAINTIES AND BACKGROUNDS IN THE DATA AND PARAMETER DETERMINATION METHOD

This paper is one of three companion papers presenting the results of our in-depth analysis of the interstellar neutral helium (ISN He) observations carried out using IBEX-Lo during the first six Interstellar Boundary Explorer (IBEX) observation seasons. We derive corrections for losses due to the limited throughput of the interface buffer and determine the IBEX spin-axis pointing. We develop an uncertainty system for the data, taking into account the resulting correlations between the data points. This system includes uncertainties due to Poisson statistics, background, spin-axis determination, systematic deviation of the boresight from the prescribed position, correction for the interface buffer losses, and the expected Warm Breeze (WB) signal. Subsequently, we analyze the data from 2009 to examine the role of various components of the uncertainty system. We show that the ISN He flow parameters are in good agreement with the values obtained by the original analysis. We identify the WB as the principal contributor to the global χ{sup 2} values in previous analyses. Other uncertainties have a much milder role and their contributions are comparable to each other. The application of this uncertainty system reduced the minimum χ{sup 2} value 4-fold. The obtained χ{sup 2} value, still exceeding the expected value, suggests that either the uncertainty system may still be incomplete or the adopted physical model lacks a potentially important element, likely an imperfect determination of the WB parameters. The derived corrections and uncertainty system are used in the accompanying paper by Bzowski et al. in an analysis of the data from six seasons.

79 ASTRONOMY AND ASTROPHYSICS↗

Synthetic galaxy clusters and observations based on Dark Energy Survey Year 3 Data

ABSTRACT We develop a novel data-driven method for generating synthetic optical observations of galaxy clusters. In cluster weak lensing, the interplay between analysis choices and systematic effects related to source galaxy selection, shape measurement, and photometric redshift estimation can be best characterized in end-to-end tests going from mock observations to recovered cluster masses. To create such test scenarios, we measure and model the photometric properties of galaxy clusters and their sky environments from the Dark Energy Survey Year 3 (DES Y3) data in two bins of cluster richness $\lambda \in [30; 45)$, $\lambda \in [45; 60)$ and three bins in cluster redshift ($z\in [0.3; 0.35)$, $z\in [0.45; 0.5)$ and $z\in [0.6; 0.65)$. Using deep-field imaging data, we extrapolate galaxy populations beyond the limiting magnitude of DES Y3 and calculate the properties of cluster member galaxies via statistical background subtraction. We construct mock galaxy clusters as random draws from a distribution function, and render mock clusters and line-of-sight catalogues into synthetic images in the same format as actual survey observations. Synthetic galaxy clusters are generated from real observational data, and thus are independent from the assumptions inherent to cosmological simulations. The recipe can be straightforwardly modified to incorporate extra information, and correct for survey incompleteness. New realizations of synthetic clusters can be created at minimal cost, which will allow future analyses to generate the large number of images needed to characterize systematic uncertainties in cluster mass measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Satellite observed thermodynamics during FGGE

During the First Global Atmospheric Research Program (GARP) Global Experiment (FGGE), determinations of temperature and moisture were made from TIROS-N and NOAA-6 satellite infrared and microwave sounding radiance measurements. The data were processed by two methods differing principally in their horizontal resolution. At the National Earth Satellite Service (NESS) in Washington, D.C., the data were produced operationally with a horizontal resolution of 250 km for inclusion in the FGGE Level IIb data sets for application to large-scale numerical analysis and prediction models. High horizontal resolution (75 km) sounding data sets were produced using man-machine interactive methods for the special observing periods of FGGE at the NASA/Goddard Space Flight Center and archived as supplementary Level IIb. The procedures used for sounding retrieval and the characteristics and quality of these thermodynamic observations are given.

Smith, W. L.↗