Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adaptive 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 181 records · Page 10

Adaptive local learning in sampling based motion planning for protein folding

Background: Simulating protein folding motions is an important problem in computational biology. Motion planning algorithms, such as Probabilistic Roadmap Methods, have been successful in modeling the folding landscape. Probabilistic Roadmap Methods and variants contain several phases (i.e., sampling, connection, and path extraction). Most of the time is spent in the connection phase and selecting which variant to employ is a difficult task. Global machine learning has been applied to the connection phase but is inefficient in situations with varying topology, such as those typical of folding landscapes. Results: We develop a local learning algorithm that exploits the past performance of methods within the neighborhood of the current connection attempts as a basis for learning. It is sensitive not only to different types of landscapes but also to differing regions in the landscape itself, removing the need to explicitly partition the landscape. We perform experiments on 23 proteins of varying secondary structure makeup with 52–114 residues. We compare the success rate when using our methods and other methods. We demonstrate a clear need for learning (i.e., only learning methods were able to validate against all available experimental data) and show that local learning is superior to global learning producing, in many cases, significantly higher quality results than the other methods. Conclusions: We present an algorithm that uses local learning to select appropriate connection methods in the context of roadmap construction for protein folding. Our method removes the burden of deciding which method to use, leverages the strengths of the individual input methods, and it is extendable to include other future connection methods.

Mathematical & Computational Biology↗

System and method for controlling depth of imaging in tissues using fluorescence microscopy under ultraviolet excitation following staining with fluorescing agents

A method is disclosed for analyzing a thin tissue sample and adapted to be supported on a slide. The tissue sample may be placed on a slide and exposed to one or more different exogenous fluorophores excitable in a range of about 300 nm-200 nm, and having a useful emission band from about 350 nm-900 nm, and including one or more fluorescent dyes or fluorescently labeled molecular probes that accumulate in tissue or cellular components. The fluorophores may be excited with a first wavelength of UV light between about 200 nm-290 nm. An optical system collects emissions from the fluorophores at a second wavelength, different from the first wavelength, which are generated in response to the first wavelength of UV light, to produce an image for analysis.

47 OTHER INSTRUMENTATION↗

Rare Events via Cross-Entropy Population Monte Carlo

Rare events are events that happen with very low frequency. Estimating rare event probabilities using Monte Carlo techniques is computationally expensive, often to the point of intractability, and special methods are required. Importance sampling (IS) is a well known technique that uses a proposal distribution in place of a target distribution to lower the variance of the estimator. Key to the success of IS methods is the choice of a proposal distribution, or the parameters governing the distribution. Adaptive importance sampling improves the parameters of a family or population of proposal distributions iteratively through trials. We present a novel cross-entropy population Monte Carlo algorithm, which adapts the parameters of proposals through the cross-entropy method. The proposed method stands apart from previous work in that we are not optimizing a mixture distribution. Instead, we leverage deterministic mixture weights and optimize the distributions individually through a reinterpretation of the typical derivation of the cross-entropy method. Demonstrations on rare event examples show that the algorithm can outperform existing resampling based population Monte Carlo methods, especially for higher-dimensional problems. Finally, we also demonstrate efficacy on a conjunction analysis problem.

97 MATHEMATICS AND COMPUTING↗

Nb3Sn films for SRF cavities: genesis and RF properties

Understanding of Nb3Sn nucleation and growth is essential to the progress with Nb3Sn vapor diffusion coatings of SRF cavities. Samples representing different stages of Nb3Sn formation have been produced and examined to elucidate the effects of nucleation, growth, process conditions, and impurities. Nb3Sn films from few hundreds of nm up to ~15 µm were grown and characterized using AFM, SEM/EDS, XPS, EBSD, SIMS, and SAM. Microscopic examinations of samples suggest the mechanisms behind Nb3Sn thin film nucleation and growth. RF measurements of coated cavities were combined with material characterization of witness samples to adapt the coating process in "Siemens" coating configuration. Understanding obtained from sample studies, applied to cavities, resulted in Nb3Sn cavity with quality factor 2 ×1010 at 15 MV/m accelerating gradient at 4 K, without "Wuppertal" Q-slope. We discuss the genesis of the Nb3Sn thin film in a typical tin vapor diffusion process, and its consequences to the coating of SRF cavities.

Pudasaini, U.↗

Sequential Kalman tuning of the t -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems

Ensemble Kalman Inversion (EKI) has been proposed as an efficient method for the approximate solution of Bayesian inverse problems with expensive forward models. However, when applied to the Bayesian inverse problem EKI is only exact in the regime of Gaussian target measures and linear forward models. Here, in this work we propose embedding EKI and Flow Annealed Kalman Inversion, its normalizing flow (NF) preconditioned variant, within a Bayesian annealing scheme as part of an adaptive implementation of the t-preconditioned Crank-Nicolson (tpCN) sampler. The tpCN sampler differs from standard pCN in that its proposal is reversible with respect to the multivariate t-distribution. The more flexible tail behaviour allows for better adaptation to sampling from non-Gaussian targets. Within our Sequential Kalman Tuning (SKT) adaptation scheme, EKI is used to initialize and precondition the tpCN sampler for each annealed target. The subsequent tpCN iterations ensure particles are correctly distributed according to each annealed target, avoiding the accumulation of errors that would otherwise impact EKI. We demonstrate the performance of SKT for tpCN on three challenging numerical benchmarks, showing significant improvements in the rate of convergence compared to adaptation within standard SMC with importance weighted resampling at each temperature level, and compared to similar adaptive implementations of standard pCN. The SKT scheme applied to tpCN offers an efficient, practical solution for solving the Bayesian inverse problem when gradients of the forward model are not available. Code implementing the SKT schemes for tpCN is available at https://github.com/RichardGrumitt/KalmanMC.

97 MATHEMATICS AND COMPUTING↗

Trigger Detection for Adaptive Scientific Workflows Using Percentile Sampling

The increasing complexity of both scientific simulations and high-performance computing system architectures are driving the need for adaptive workflows, in which the composition and execution of computational and data manipulation steps dynamically depend on the evolutionary state of the simulation itself. Consider, for example, the frequency of data storage. Critical phases of the simulation should be captured with high frequency and with high fidelity for postanalysis; however, we cannot afford to retain the same frequency for the full simulation due to the high cost of data movement. We can instead look for triggers, indicators that the simulation will be entering a critical phase, and adapt the workflow accordingly. In this paper, we present a methodology for detecting triggers and demonstrate its use in the context of direct numerical simulations of turbulent combustion using S3D. We show that chemical explosive mode analysis (CEMA) can be used to devise a noise-tolerant indicator for rapid increase in heat release. However, exhaustive computation of CEMA values dominates the total simulation, and thus is prohibitively expensive. To overcome this computational bottleneck, we propose a quantile sampling approach. Our sampling-based algorithm comes with provable error/confidence bounds, as a function of the number of samples. Most importantly, the number of samples is independent of the problem size, and thus our proposed sampling algorithm offers perfect scalability. Our experiments on homogeneous charge compression ignition and reactivity controlled compression ignition simulations show that the proposed method can detect rapid increases in heat release, and its computational overhead is negligible. Our results will be used to make dynamic workflow decisions regarding data storage and mesh resolution in future combustion simulations.

97 MATHEMATICS AND COMPUTING↗

Three-Dimensional Mass Spectrometric Imaging of Biological Structures Using a Vacuum-Compatible Microfluidic Device

Three-dimensional (3D) molecular imaging of biological structures is important for a wide range of research. In recent decades, secondary ion mass spectrometry (SIMS) has been recognized as a powerful technique for both two-dimensional (2D) and 3D molecular imaging. Sample fixations (e. g., chemical fixation and cryogenic fixation methods) are necessary to adapt biological samples to the vacuum condition in the SIMS chamber, which has been demonstrated to be non-trivial and less controllable, thus limiting the wider application of SIMS on 3D molecular analysis of biological samples. Our group recently developed in situ liquid SIMS that offers great opportunities for the molecular study of various liquids and liquid interfaces. In this work, we demonstrate that a further development of the vacuum-compatible microfluidic device used in in situ liquid SIMS provides a convenient freeze-fixation of biological samples and leads to more controllable and convenient 3D molecular imaging. The special design of this new vacuum-compatible liquid chamber allows an easy determination of sputter rates of ice, which is critical for calibrating the depth scale of frozen biological samples. Sputter yield of a 20 keV Ar 1800 + ion on ice has been determined as 1500 (± 8%) water molecules per Ar 1800 + ion, consistent with our results from molecular dynamics simulations. Moreover, using the information of ice sputter yield, we successfully conduct 3D molecular imaging of frozen homogenized milk and observe network structures of interesting organic and inorganic species. Finally, taken together, our results will significantly benefit various research fields relying on 3D molecular imaging of biological structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural system reliability calculation using a probabilistic fault tree analysis method

The development of a new probabilistic fault tree analysis (PFTA) method for calculating structural system reliability is summarized. The proposed PFTA procedure includes: developing a fault tree to represent the complex structural system, constructing an approximation function for each bottom event, determining a dominant sampling sequence for all bottom events, and calculating the system reliability using an adaptive importance sampling method. PFTA is suitable for complicated structural problems that require computer-intensive computer calculations. A computer program has been developed to implement the PFTA.

Torng, T. Y.↗

Reducing Contingency through Sampling at the Luckey FUSRAP Site - 13186

Typically, the greatest risk in developing accurate cost estimates for the remediation of hazardous, toxic, and radioactive waste sites is the uncertainty in the estimated volume of contaminated media requiring remediation. Efforts to address this risk in the remediation cost estimate can result in large cost contingencies that are often considered unacceptable when budgeting for site cleanups. Such was the case for the Luckey Formerly Utilized Sites Remedial Action Program (FUSRAP) site near Luckey, Ohio, which had significant uncertainty surrounding the estimated volume of site soils contaminated with radium, uranium, thorium, beryllium, and lead. Funding provided by the American Recovery and Reinvestment Act (ARRA) allowed the U.S. Army Corps of Engineers (USACE) to conduct additional environmental sampling and analysis at the Luckey Site between November 2009 and April 2010, with the objective to further delineate the horizontal and vertical extent of contaminated soils in order to reduce the uncertainty in the soil volume estimate. Investigative work included radiological, geophysical, and topographic field surveys, subsurface borings, and soil sampling. Results from the investigative sampling were used in conjunction with Argonne National Laboratory's Bayesian Approaches for Adaptive Spatial Sampling (BAASS) software to update the contaminated soil volume estimate for the site. This updated volume estimate was then used to update the project cost-to-complete estimate using the USACE Cost and Schedule Risk Analysis process, which develops cost contingencies based on project risks. An investment of $1.1 M of ARRA funds for additional investigative work resulted in a reduction of 135,000 in-situ cubic meters (177,000 in-situ cubic yards) in the estimated base volume estimate. This refinement of the estimated soil volume resulted in a $64.3 M reduction in the estimated project cost-to-complete, through a reduction in the uncertainty in the contaminated soil volume estimate and the associated contingency costs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

Micro on a macroscale: relating microbial-scale soil processes to global ecosystem function

ABSTRACT Soil microorganisms play a key role in driving major biogeochemical cycles and in global responses to climate change. However, understanding and predicting the behavior and function of these microorganisms remains a grand challenge for soil ecology due in part to the microscale complexity of soils. It is becoming increasingly clear that understanding the microbial perspective is vital to accurately predicting global processes. Here, we discuss the microbial perspective including the microbial habitat as it relates to measurement and modeling of ecosystem processes. We argue that clearly defining and quantifying the size, distribution and sphere of influence of microhabitats is crucial to managing microbial activity at the ecosystem scale. This can be achieved using controlled and hierarchical sampling designs. Model microbial systems can provide key data needed to integrate microhabitats into ecosystem models, while adapting soil sampling schemes and statistical methods can allow us to collect microbially-focused data. Quantifying soil processes, like biogeochemical cycles, from a microbial perspective will allow us to more accurately predict soil functions and address long-standing unknowns in soil ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Triton Field Trials - Changes in Habitats, a Literature Review of Monitoring Technologies

Marine energy devices are installed in highly dynamic environments and have the potential to affect benthic and pelagic habitats around them. Regulatory bodies often require baseline characterization and/or post-installation monitoring to determine whether changes in these habitats are being observed. However, a great diversity of technologies is available for surveying and sampling marine habitats. Selecting the most suitable instrument to identify and measure changes in habitats at marine energy sites can become a daunting task. We conducted a thorough review of journal articles, survey reports, and grey literature to extract information about the technologies used, the data collection and processing methods, and the performance and effectiveness of these instruments. We examined documents related to marine energy development, offshore wind farms, oil and gas offshore sites, and other marine industries around the world over the last 20 years, as well as national and international guidelines for surveying habitats around offshore activities. A total of 120 different technologies were identified across six main habitat categories: seafloor, sediment, infauna, epifauna, pelagic, and biofouling. The technologies were organized into 12 broad technology classes: acoustic, corer, dredge, grab, hook and line, net and trawl, plate, remote sensing, scrape samples, trap, visual, and others. Visual was the most common and the most diverse technology class, with applications across all six habitat categories. Sampling designs varied considerably among the reviewed studies but transect was the predominant design for surveying seafloor, epifauna, and pelagic habitats. The most common data analyses were univariate and multivariate statistical analyses aimed at calculating and comparing biodiversity indices, characterizing faunal assemblages or sediment classes, or modeling the distribution of animals related to abiotic parameters. Technologies and sampling methods adaptable and designed to work efficiently in energetic environments have greater success at marine energy sites. In addition, sampling designs and statistical analyses should be carefully thought through to identify differences in faunal assemblages and spatiotemporal changes in habitats.

16 TIDAL AND WAVE POWER↗

Design of an L-band Microwave Radiometer with Active Mitigation of Interference

Radio frequency interference (RFI) impairs L-band radiometry outside the protected 20 MHz frequency band around 1413 MHz. However, bandwidths of 100 MHz or more are desired for certain remote sensing applications as well as certain astronomy applications. Because much of the RFI in this band is from radars with pulse lengths on the order of microseconds, traditional radiometers (i.e., those which directly measure total power or power spectral density integrated over time scales of milliseconds or greater) are poorly-suited to this task. Simply reducing integration time and discarding contaminated outputs may not be a practical answer due to the wide variety of modulations and pulse lengths observed in L-band RFI signals, the dynamic and complex nature of the associated propagation channels, and the logistical effort associated with post-measurement data editing. This motivates the design and development of radiometers capable of coherent sampling and adaptive, real-time mitigation of interference. Such a radiometer will be described in this presentation. This design is capable of coherently-sampling up to 100 MHz bandwidth at L-band. RFI mitigation is implemented in FPGA components so that real-time suppression is achieved. The system currently uses a cascade of basic time- and frequency-domain detection and blanking techniques; more advanced algorithms are under consideration. The modular FPGA-based architecture provides other benefits, such as the ability to implement extremely stable digital filters and the ability to reconfigure the system "on the fly". An overview of the basic design along with on-the-air results from an initial implementation will be provided in the presentation. Related L-band RFI surveys will be described to illustrate the relevance of this approach in a variety of operating conditions.

Hampson, G. A.↗

Design of an L-band Microwave Radiometer with Active Mitigation of Interference

Radio frequency interference (RFI) impairs L-band radiometry outside the protected 20 MHz frequency band around 1413 MHz. However, bandwidths of 100 MHz or more are desired for certain remote sensing applications as well as certain astronomy applications. Because much of the RFI in this band is from radars with pulse lengths on the order of microseconds, traditional radiometers (i.e., those which directly measure total power or power spectral density integrated over time scales of milliseconds or greater) are poorly-suited to this task. Simply reducing integration time and discarding contaminated outputs may not be a practical answer due to the wide variety of modulations and pulse lengths observed in L-band RFI signals, the dynamic and complex nature of the associated propagation channels, and the logistical effort associated with post-measurement data editing. This motivates the design and development of radiometers capable of coherent sampling and adaptive, real-time mitigation of interference. Such a radiometer will be described in this presentation. This design is capable of coherently-sampling up to 100 MHz bandwidth at L-band. RFI mitigation is implemented in FPGA components so that real-time suppression is achieved. The system currently uses a cascade of basic time- and frequency- domain detection and blanking techniques; more advanced algorithms are un- der consideration. The modular FPGA-based architecture provides other benefits, such as the ability to implement extremely stable digital filters and the ability to reconfigure the system "on the fly". An overview of the basic design along with on-the-air results from an initial implementation will be provided in the presentation. Related L-band RFI surveys will be described to illustrate the relevance of this approach in a variety of operating conditions.

Ellingson, S. W.↗

Organ sample generator for expected treatment dose construction and adaptive inverse planning optimization

Purpose: To create an organ sample generator (OSG) for expected treatment dose construction and adaptive inverse planning optimization. The OSG generates random samples of organs of interest from a distribution obeying the patient specific organ variation probability density function (PDF) during the course of adaptive radiotherapy. Methods: Principle component analysis (PCA) and a time-varying least-squares regression (LSR) method were used on patient specific geometric variations of organs of interest manifested on multiple daily volumetric images obtained during the treatment course. The construction of the OSG includes the determination of eigenvectors of the organ variation using PCA, and the determination of the corresponding coefficients using time-varying LSR. The coefficients can be either random variables or random functions of the elapsed treatment days depending on the characteristics of organ variation as a stationary or a nonstationary random process. The LSR method with time-varying weighting parameters was applied to the precollected daily volumetric images to determine the function form of the coefficients. Eleven h and n cancer patients with 30 daily cone beam CT images each were included in the evaluation of the OSG. The evaluation was performed using a total of 18 organs of interest, including 15 organs at risk and 3 targets. Results: Geometric variations of organs of interest during h and n cancer radiotherapy can be represented using the first 3 {approx} 4 eigenvectors. These eigenvectors were variable during treatment, and need to be updated using new daily images obtained during the treatment course. The OSG generates random samples of organs of interest from the estimated organ variation PDF of the individual. The accuracy of the estimated PDF can be improved recursively using extra daily image feedback during the treatment course. The average deviations in the estimation of the mean and standard deviation of the organ variation PDF for h and n cancer radiotherapy were less than 2 and 1 mm, respectively, for most organs after the second week of treatment. After the first three weeks of treatment, the mean discrepancy of the dose estimation accuracy was within 1% for most of organs, the corresponding standard deviation was within 2.5% for parotids, the brain stem and the cochleae, and within 1% for other organs. Conclusions: A patient specific OSG is feasible and can be used to generate random samples of organs of interest for the expected treatment dose construction and adaptive inverse planning. The accuracy of the OSG can be improved continuously and recursively during the adaptive treatment course using daily volumetric image feedback.

60 APPLIED LIFE SCIENCES↗