Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse sampling”

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 217 records · Page 12

Stochastic AC optimal power flow: A data-driven approach

There is an emerging need for efficient solutions to stochastic AC Optimal Power Flow (AC-OPF) to ensure optimal and reliable grid operations in the presence of increasing demand and generation uncertainty. Herein this paper presents a highly scalable data-driven algorithm for stochastic AC-OPF that has extremely low sample requirement. The novelty behind the algorithm’s performance involves an iterative scenario design approach that merges information regarding constraint violations in the system with data-driven sparse regression. Compared to conventional methods with random scenario sampling, our approach is able to provide feasible operating points for realistic systems with much lower sample requirements. Furthermore, multiple sub-tasks in our approach can be easily paralleled and based on historical data to enhance its performance and application. We demonstrate the computational improvements of our approach through simulations on different test cases in the IEEE PES PGLib-OPF benchmark library.

42 ENGINEERING↗

Scalable Techniques for Stochastic Power Flow Problems (Final Report)

The proposed research focuses on developing scalable algorithms for two-stage security-constrained OPF problems with AC power flow constraints, a class of problems complicated by (i) scale arising from a scenario representation; and (ii) the presence of nonlinearity, nonconvexity, and possibly second-stage discreteness or complementarity. Unfortunately, most existing solvers cannot contend with both challenges simultaneously; accordingly, the proposed research focuses on developing solution techniques that can both scale with the number of scenarios and contend with nonconvexity and second-stage complementarity. We consider three avenues for addressing such problems: (i) Variable sample-size SQP (VS-SQP) methods that combine sparse Quasi-Newton updates with a scalable variance-reduced stochastic gradient scheme for stochastic QP subproblems, allowing for contending with second-stage complementarity via regularization; (ii) Variable sample-size stochastic Interior-point (VS-sIP) schemes that propose a sampling-based regularized (to allow for contending with complementarity) interior-point schemes in which a Schur-complement technique is employed for decomposing the Newton direction computation step; (iii) Variable sample-size tractable ADMM (VS-tADMM) schemes combine variable sample-sizes with carefully designed techniques for resolving each of the nonconvex updates (by leveraging the QCQP structures). We intend to compare the three schemes using performance profiles in terms of solution quality, scalability, etc. and then select one scheme which will then be developed and further refined in Python for purposes of the GO competition.

42 ENGINEERING↗

LUMINOSITY FUNCTIONS OF SPITZER-IDENTIFIED PROTOSTARS IN NINE NEARBY MOLECULAR CLOUDS

We identify protostars in Spitzer surveys of nine star-forming (SF) molecular clouds within 1 kpc: Serpens, Perseus, Ophiuchus, Chamaeleon, Lupus, Taurus, Orion, Cep OB3, and Mon R2, which combined host over 700 protostar candidates. These clouds encompass a variety of SF environments, including both low-mass and high-mass SF regions, as well as dense clusters and regions of sparsely distributed star formation. Our diverse cloud sample allows us to compare protostar luminosity functions in these varied environments. We combine near- and mid-infrared photometry from the Two Micron All Sky Survey and Spitzer to create 1-24 {mu}m spectral energy distributions (SEDs). Using protostars from the c2d survey with well-determined bolometric luminosities, we derive a relationship between bolometric luminosity, mid-IR luminosity (integrated from 1-24 {mu}m), and SED slope. Estimations of the bolometric luminosities for protostar candidates are combined to create luminosity functions for each cloud. Contamination due to edge-on disks, reddened Class II sources, and galaxies is estimated and removed from the luminosity functions. We find that luminosity functions for high-mass SF clouds (Orion, Mon R2, and Cep OB3) peak near 1 L{sub Sun} and show a tail extending toward luminosities above 100 L{sub Sun }. The luminosity functions of the low-mass SF clouds (Serpens, Perseus, Ophiuchus, Taurus, Lupus, and Chamaeleon) do not exhibit a common peak, however the combined luminosity function of these regions peaks below 1 L{sub Sun }. Finally, we examine the luminosity functions as a function of the local surface density of young stellar objects. In the Orion molecular clouds, we find a significant difference between the luminosity functions of protostars in regions of high and low stellar density, the former of which is biased toward more luminous sources. This may be the result of primordial mass segregation, although this interpretation is not unique. We compare our luminosity functions to those predicted by models and find that our observed luminosity functions are best matched by models that invoke competitive accretion, although we do not find strong agreement between the high-mass SF clouds and any of the models.

79 ASTRONOMY AND ASTROPHYSICS↗

Stacked reverberation mapping of high-redshift quasars in DESI. I. Feasibility analysis

The broad-line region of quasars has long been probed by reverberation mapping techniques that measure time lags between continuum and broad emission-line variations. Stacked reverberation mapping has been proposed as a less observationally expensive alternative to traditional methods. This ensemble approach also reduces biases from small-number statistics. The Dark Energy Spectroscopic Instrument (DESI) is conducting the most extensive spectroscopic survey of quasars to date. We create mock light curves emulating expected DESI quasar observations at redshifts $1.48\lt z\lt 5.2$ and luminosities $44.68 \le \log \lambda L_{1350 \mathring{\rm A}{}} / \mathrm{erg\, s^{-1}} \le 45.99$ to test stacked reverberation mapping feasibility using sparse spectroscopic data paired with well-sampled photometric data. The pipeline, using the lag estimation code JAVELIN (Just Another Vehicle for Estimating Lags In Nuclei), successfully recovers the simulated C IV lags within 1σ of the true values using spectroscopic light curves composed of only a few spectral epochs (2–10) with irregular cadences. We investigate how observational factors, including C IV flux error magnitude, number of stacked quasars, and spectral epoch count, affect performance. This work motivates a pathway for future stacked reverberation mapping projects with large-scale spectroscopic surveys of quasars having $\ge 2$ spectroscopic observations. Our results suggest an economical alternative for constraining and extending the radius–luminosity relation to higher redshifts and luminosities. Subsequently, this relation can be employed more reliably in single-epoch black hole mass measurements and quasar cosmology in these distant regimes.

quasars: general, quasars: supermassive black hole↗

The plug-based nanovolume Microcapillary Protein Crystallization System (MPCS)

The Microcapillary Protein Crystallization System (MPCS) embodies a new semi-automated plug-based crystallization technology which enables nanolitre-volume screening of crystallization conditions in a plasticware format that allows crystals to be easily removed for traditional cryoprotection and X-ray diffraction data collection. Protein crystals grown in these plastic devices can be directly subjected to in situ X-ray diffraction studies. The MPCS integrates the formulation of crystallization cocktails with the preparation of the crystallization experiments. Within microfluidic Teflon tubing or the microfluidic circuitry of a plastic CrystalCard, {approx}10-20 nl volume droplets are generated, each representing a microbatch-style crystallization experiment with a different chemical composition. The entire protein sample is utilized in crystallization experiments. Sparse-matrix screening and chemical gradient screening can be combined in one comprehensive 'hybrid' crystallization trial. The technology lends itself well to optimization by high-granularity gradient screening using optimization reagents such as precipitation agents, ligands or cryoprotectants.

59 BASIC BIOLOGICAL SCIENCES↗

A Bayesian Approach for Extracting Fluorescence Lifetimes from Sparse Data Sets and Its Significance for Imaging Experiments

The measurement of fluorescence lifetimes, especially in small sample volumes, presents the dual challenge of probing a small number of fluorophores and fitting the concomitant sparse data set to the appropriate excited-state decay function. A common method of analysis, such as the maximum likelihood (ML) technique, assumes a uniform probability distribution of the parameters describing the fluorescence decay function. An improvement is thus suggested by implementing a suitable nonuniform distribution, as is provided by a Bayesian framework, where the distribution of parameters is obtained from both their prior knowledge and the evidence-based likelihood of an event for a given set of parameters. We have also considered the Dirichlet prior distribution, whose mathematical form enables analytical solutions of the fitting parameters to be rapidly obtained. If Gaussian and exponential prior distributions are judiciously chosen, they reproduce the experimental target lifetime to within 20% with as few as 20 total photon counts for the data set, as does the Dirichlet prior distribution. But because of the analytical solutions afforded by the Dirichlet prior distribution, it is proposed to employ a Dirichlet prior to search parameter space rapidly to provide, if necessary, appropriate parameters for subsequent employment of a Gaussian or exponential prior distribution

97 MATHEMATICS AND COMPUTING↗

Assessing the limitations of commercial sensors and models for supporting marine carbon dioxide removal monitoring: a case study

Several unknowns remain surrounding marine Carbon Dioxide Removal (mCDR) monitoring, reporting, and verification (MRV) practices and capabilities. Current in-situ sensor technology is limited (primarily pH and pCO 2 ), requiring calculations and assumptions to estimate changes in carbonate chemistry parameters, including total alkalinity (TA). Considering that cost, energy consumption, and accuracy of commercial sensors can vary by orders of magnitude, understanding how well existing sensors perform in an mCDR context is important for this emerging community. Likewise, documenting sensor limitations and how relatively simple models can optimize sensor deployments will improve MRV efforts and support protocol development. Here we (1) compare performance a variety of commercially available sensors in a blind mesocosm experiment simulating ocean alkalinity enhancement (OAE), and how sensor performance impacted carbonate chemistry estimates; (2) evaluate if sensors can distinguish the OAE signal from natural variability during a small scale OAE field test in Sequim Bay, WA, USA, and (3) use an idealized ocean biogeochemistry model to explore optimal sensor network design based on (1) and (2). Our mesocosm results indicate that correctly constraining pH uncertainty will be critical for accurate TA estimates with current sensor technology compared to the less impactful variation caused by uncertainty in pCO 2 (pH data that are presented throughout are reported on the total scale (pH T ) unless otherwise noted). Our pilot field test demonstrated that sensors were capable of distinguishing mCDR signatures from natural variability under optimal real-world conditions. Idealized modeling simulations of the field test showed that a range of sparse and dense (3 to 100) sensors sampling areas of detectable increases will underestimate the net change in surface pH by at least 35–55%, at both realistic and highly elevated alkalinity input levels. We also highlight the limitations of current sensing technology for MRV, and the importance of ocean biogeochemistry models as critical tools for predicting when and where mCDR signals will be detectable using available sensors. Overall, our findings suggest that commercially available pCO 2 sensors and some pH sensors will form an important backbone for mCDR MRV tasks, though complete MRV characterization will require these data to be used in combination with other tools.

OAE↗

Sensing analytical instrument parameters, specimen characteristics, or both from sparse datasets

Disclosed are methods for sensing conditions of an electron microscope system and/or a specimen analyzed thereby. Also disclosed are sensor systems and electron microscope systems able to sense system conditions, and/or conditions of the specimen being analyzed by such systems. In one embodiment, a sparse dataset can be acquired from a random sub-sampling of the specimen by an electron beam probe of the electron microscope system. Instrument parameters, specimen characteristics, or both can be estimated from the sparse dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

YOUNG STELLAR POPULATIONS IN MYStIX STAR-FORMING REGIONS: CANDIDATE PROTOSTARS

The Massive Young Star-Forming Complex in Infrared and X-ray (MYStIX) project provides a new census on stellar members of massive star-forming regions within 4 kpc. Here the MYStIX Infrared Excess catalog and Chandra -based X-ray photometric catalogs are mined to obtain high-quality samples of Class I protostars using criteria designed to reduce extragalactic and Galactic field star contamination. A total of 1109 MYStIX Candidate Protostars (MCPs) are found in 14 star-forming regions. Most are selected from protoplanetary disk infrared excess emission, but 20% are found from their ultrahard X-ray spectra from heavily absorbed magnetospheric flare emission. Two-thirds of the MCP sample is newly reported here. The resulting samples are strongly spatially associated with molecular cores and filaments on Herschel far-infrared maps. This spatial agreement and other evidence indicate that the MCP sample has high reliability with relatively few “false positives” from contaminating populations. But the limited sensitivity and sparse overlap among the infrared and X-ray subsamples indicate that the sample is very incomplete with many “false negatives.” Maps, tables, and source descriptions are provided to guide further study of star formation in these regions. In particular, the nature of ultrahard X-ray protostellar candidates without known infrared counterparts needs to be elucidated.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity study of a full-scale industrial spray-injected fluidized bed reactor

The industrial fluidized bed reactor (FBR) described here is designed to convert an aqueous solid laden stream into a consistent solid product material. Additionally, the FBR is heated to 650 °C and the aqueous reactant is spray-atomized into the fluidized bed, converting it into a granular solid through a series of heterogeneous chemical reactions. The objective of this study is to use the MFiX two-fluid model as a screening tool to investigate the effects of off-nominal conditions for the purpose of establishing performance guidelines. This information will be helpful to assess whether the identified boundaries of the system are steep or gradual. Initially, 13 independent operating parameters were sampled in a one-at-a-time manner using high, low, and base values and ranked against a set of performance criteria. Five operating parameters were found to have the largest effect (bed particle size, bed particle density, coal particle size, spray feed flow rate, and fluidizing gas flow rate) on three quantities of interest—bed differential temperature, solids velocity, and bed voidage. Latin hypercube sampling was used to generate a minimal number of random values for various combinations of input parameters. The spray feed flow rate was the most significant parameter for the temperature differences, and the coal particle and bed particle size were significant contributors. The bed particle size had the largest effect on the low velocity of the bed particles, with the coal particle size as a contributing second effect. The high solid packing shows the coal particle size with the largest effect, and significant secondary contributions from the feed flow rate and fluidizing gas flow rate. Finally, the fits to the five-dimensional Gaussian Process models were 0.7797, 0.8664, and 0.9440 for the temperature, velocity, and solids packing respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Disk Density Tuning of a Maximal Random Packing

We introduce an algorithmic framework for tuning the spatial density of disks in a maximal random packing, without changing the sizing function or radii of disks. Starting from any maximal random packing such as a Maximal Poisson-disk Sampling (MPS), we iteratively relocate, inject (add), or eject (remove) disks, using a set of three successively more-aggressive local operations. We may achieve a user-defined density, either more dense or more sparse, almost up to the theoretical structured limits. The tuned samples are conflict-free, retain coverage maximality, and, except in the extremes, retain the blue noise randomness properties of the input. We change the density of the packing one disk at a time, maintaining the minimum disk separation distance and the maximum domain coverage distance required of any maximal packing. These properties are local, and we can handle spatially-varying sizing functions. Using fewer points to satisfy a sizing function improves the efficiency of some applications. We apply the framework to improve the quality of meshes, removing non-obtuse angles; and to more accurately model fiber reinforced polymers for elastic and failure simulations.

97 MATHEMATICS AND COMPUTING↗

Background-free dual-mode optical and 13 C magnetic resonance imaging in diamond particles

Multimodal imaging—the ability to acquire images of an object through more than one imaging mode simultaneously—has opened additional perspectives in areas ranging from astronomy to medicine. In this paper, we report progress toward combining optical and magnetic resonance (MR) imaging in such a “dual” imaging mode. They are attractive in combination because they offer complementary advantages of resolution and speed, especially in the context of imaging in scattering environments. Our approach relies on a specific material platform, microdiamond particles hosting nitrogen vacancy (NV) defect centers that fluoresce brightly under optical excitation and simultaneously “hyperpolarize” lattice C 13 nuclei, making them bright under MR imaging. We highlight advantages of dual-mode optical and MR imaging in allowing background-free particle imaging and describe regimes in which either mode can enhance the other. Leveraging the fact that the two imaging modes proceed in Fourier-reciprocal domains (real and k-space), we propose a sampling protocol that accelerates image reconstruction in sparse-imaging scenarios. Our work suggests interesting possibilities for the simultaneous optical and low-field MR imaging of targeted diamond nanoparticles.

36 MATERIALS SCIENCE↗

The Value of Information From Horizontal Distributed Acoustic Sensing Compared to Multicomponent Geophones Via Machine Learning

Abstract Faults play an important role in recharging many geothermal reservoirs, and seismic information can image the locations of these faults. The value of information (VOI) metric is used to objectively quantify and compare the value of two types of seismic receiver data via a machine learning approach. The demonstrated VOI methodology is novel by including spatial models from seismic data and obtaining the information statistics from machine learning. Our two-dimensional numerical experiments compare images created from sparsely spaced (80 m), two-component geophone sampling to high spatial resolution (1 m), single-component DAS. We used a three-fold cross validation of a U-Net convolutional neural networks to achieve average classification statistics. The results suggest that when horizontal sources are utilized, geophones and DAS identify reflectors and non-reflectors at roughly the same rate. The average F1 score for horizontal DAS is 0.939 and 0.931 for geophones. For images created from a vertical source, DAS performed marginally better (F1 = 0.919) than geophones (F1 = 0.877). Our transferrable methodology can provide guidance on which acquisition scenarios can improve images of important structures in the subsurface and present an efficient method for obtaining reliability statistics from high-dimensional, spatial data.

Energy & Fuels↗

Resolution-adapted recombination of structural features significantly improves sampling in restraint-guided structure calculation

Recent work has shown that NMR structures can be determined by integrating sparse NMR data with structure prediction methods such as Rosetta. The experimental data serve to guide the search for the lowest energy state towards the deep minimum at the native state which is frequently missed in Rosetta de novo structure calculations. However, as the protein size increases, sampling again becomes limiting; for example, the standard Rosetta protocol involving Monte Carlo fragment insertion starting from an extended chain fails to converge for proteins over 150 amino acids even with guidance from chemical shifts (CS-Rosetta) and other NMR data. The primary limitation of this protocol—that every folding trajectory is completely independent of every other—was recently overcome with the development of a new approach involving resolution-adapted structural recombination (RASREC). Here we describe the RASREC approach in detail and compare it to standard CS-Rosetta. We show that the improved sampling of RASREC is essential in obtaining accurate structures over a benchmark set of 11 proteins in the 15-25 kDa size range using chemical shifts, backbone RDCs and H N -H N NOE data; in a number of cases the improved sampling methodology makes a larger contribution than incorporation of additional experimental data. Experimental data are invaluable for guiding sampling to the vicinity of the global energy minimum, but for larger proteins, the standard Rosetta fold-from-extended-chain protocol does not converge on the native minimum even with experimental data and the more powerful RASREC approach is necessary to converge to accurate solutions.

59 BASIC BIOLOGICAL SCIENCES↗

Predicting SARS-CoV-2 Variant Using Non-Invasive Hand Odor Analysis: A Pilot Study

The adaptable nature of the SARS-CoV-2 virus has led to the emergence of multiple viral variants of concern. This research builds upon a previous demonstration of sampling human hand odor to distinguish SARS-CoV-2 infection status in order to incorporate considerations of the disease variants. This study demonstrates the ability of human odor expression to be implemented as a non-invasive medium for the differentiation of SARS-CoV-2 variants. Volatile organic compounds (VOCs) were extracted from SARS-CoV-2-positive samples using solid phase microextraction (SPME) coupled with gas chromatography–mass spectrometry (GC–MS). Sparse partial least squares discriminant analysis (sPLS-DA) modeling revealed that supervised machine learning could be used to predict the variant identity of a sample using VOC expression alone. The class discrimination of Delta and Omicron BA.5 variant samples was performed with 95.2% (±0.4) accuracy. Omicron BA.2 and Omicron BA.5 variants were correctly classified with 78.5% (±0.8) accuracy. Lastly, Delta and Omicron BA.2 samples were assigned with 71.2% (±1.0) accuracy. This work builds upon the framework of non-invasive techniques producing diagnostics through the analysis of human odor expression, all in support of public health monitoring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Streamline-Based Time-Lapse-Seismic-Data Integration Incorporating Pressure and Saturation Effects

In this paper, we present an efficient history-matching technique that simultaneously integrates 4D repeat seismic surveys with well-production data. This approach is particularly well-suited for the calibration of the reservoir properties of high-resolution geologic models because the seismic data are a really dense but sparse in time, whereas the production data are finely sampled in time but spatially averaged. The joint history matching is performed by use of streamline-based sensitivities derived from either finite-difference or streamline-based flow simulation. For the most part, earlier approaches have focused on the role of saturation changes, but the effects of pressure have largely been ignored. Here, we present a streamline-based semianalytic approach for computing model-parameter sensitivities, accounting for both pressure and saturation effects. The novelty of the method lies in the semianalytic sensitivity computations, making it computationally efficient for high-resolution geologic models. The approach is implemented by use of a finite-difference simulator incorporating the detailed physics. Its efficacy is demonstrated by use of both synthetic and field applications. Lastly, for both the synthetic and the field cases, the advantages of incorporating the time-lapse variations are clear, seen through the improved estimation of the permeability distribution, the pressure profile, the evolution of the fluid saturation, and the swept volumes.

58 GEOSCIENCES↗

Sparse Multi-Static Arrays for Near-Field Millimeter-Wave Imaging

This paper describes a novel design technique for sparse multi-static linear arrays. The methods described allow the development of densely sampled linear arrays suitable for high-resolution near-field imaging that require dramatically fewer antenna and switch elements than the previous state of the art. The techniques used are related to sparse array techniques used in radio astronomy applications, but differ significantly in design due to the transmit-receive nature of the arrays, and the application to linear arrays that achieve dense uniform sampling suitable for high-resolution near-field imaging. As many as 3 to 5 or more samples per antenna can be obtained, compared to 1 sample per antenna for the current state of the art. This could dramatically reduce cost and improve performance over current active millimeter-wave imaging systems.

Millimeter waves, imaging, personnel screening, co↗

Linking Associations of Rare Low-Abundance Species to Their Environments by Association Networks

Studies of microbial communities by targeted sequencing of rRNA genes lead to recovering numerous rare low-abundance taxa with unknown biological roles. We propose to study associations of such rare organisms with their environments by a computational framework based on transformation of the data into qualitative variables. Namely, we analyze the sparse table of putative species or OTUs (operational taxonomic units) and samples generated in such studies, also known as an OTU table, by collecting statistics on co-occurrences of the species and on shared species richness across samples. Based on the statistics we built two association networks, of the rare putative species and of the samples respectively, using a known computational technique, Association networks (Anets) developed for analysis of qualitative data. Clusters of samples and clusters of OTUs are then integrated and combined with metadata of the study to produce a map of associated putative species in their environments. We tested and validated the framework on two types of microbiomes, of human body sites and that of the Populus tree root systems. We show that in both studies the associations of OTUs can separate samples according to environmental or physiological characteristics of the studied systems.

54 ENVIRONMENTAL SCIENCES↗