Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set 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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Active Power Control of Wind Turbines for Ancillary Services: A Comparison of Pitch and Torque Control Methodologies

As wind energy generation becomes more prevalent in some regions, there is increased demand for wind power plants to provide ancillary services, which are essential for grid reliability. This paper compares two different wind turbine control methodologies to provide active power control (APC) ancillary services, which include derating or curtailing power generation, providing automatic generation control (AGC), and providing primary frequency control (PFC). The torque APC controller provides all power control through the power electronics whereas the pitch APC controller uses the blade pitch actuators as the primary means of power control. These controllers are simulated under various wind conditions with different derating set points and AGC participation levels. The metrics used to compare their performance are the damage equivalent loads (DELs) induced on the structural components and AGC performance metrics, which are used to determine the payments for AGC services by system operators in the United States. The simulation results show that derating the turbine reduces structural loads for both control methods, with the APC pitch control providing larger reductions in DELs, lower AGC performance scores, and higher root-mean-square pitch rates. Providing AGC increases the structural loads when compared to only derating the turbine, but even the AGC DELs are generally lower than those of the baseline control system. The torque APC control methodology also allows for more sustained PFC responses under certain derating conditions.

17 WIND ENERGY↗
DOE OSTI2022

Behavioral and Population Data-Driven Distribution System Load Modeling

Distribution system residential load modeling and analysis for different geographic areas within a utility or an independent system operator territory are critical for enabling small-scale, aggregated distributed energy resources to participate in grid services under Federal Energy Regulatory Commission Order No. 2222 [1]. In this study, we develop a methodology of modeling residential load profiles in different geographic areas with a focus on human behavior impact. First, we construct a behavior-based load profile model leveraging state-of-the-art appliance models. We simulate human activity and occupancy using Markov chain Monte Carlo methods calibrated with the American Time Use Survey data set. Second, we link our model with cleaned Current Population Survey data from the U.S. Census Bureau. Finally, we populate two sets of 500 households using California and Texas census data, respectively, to perform an initial analysis of the load in different geographic areas with various group features (e.g., different income levels). To distinguish the effect of population behavior differences on aggregated load, we simulate load profiles for both sets assuming fixed physical household parameters and weather data. Analysis shows that average daily load profiles vary significantly by income and income dependency varies by locality.

American Time Use Survey↗
DOE OSTI2016

WE-FG-BRA-08: Potential Role of the Glycolytic Oscillator in Acute Hypoxia in Tumors

Purpose: Oscillatory dynamics in acute hypoxia have been observed, but poorly understood. They have mostly been attributed to vascular perturbations, but no link has yet been made to metabolic causes. We set out to determine the fundamental frequencies and test for coherence in tumor oxygen dynamics and spatial properties. Methods: Severe combined immunodeficient (SCID) mice were inoculated onto bilateral flanks with human derived head and neck carcinoma (UW-SCC22) cell line xenografts. Oxygen dynamics were monitored in the tumor every minute for an hour using three modalities: blood oxygen level dependent - magnetic resonance imaging (BOLD-MRI), hemoglobin oxygen saturation photoacoustic, and locally manufactured optical probes for spectral fitting. A statistical test was used to separate fluctuating from non-fluctuating voxels and pixels in BOLD-MRI and photoacoustic data respectively. The power spectrum density (PSD) and the autocorrelation functions were calculated for the time series of each voxel, pixel and region, of the BOLD-MRI, photoacoustic or fiber optic data respectively. Results: Using all three techniques, intermittent oxygen dynamics with both coherent and incoherent signatures was observed in the tumors. Upon averaging the PSDs of fluctuating voxels and pixels, it was found that these oscillations occurred with periods of minutes to tens of minutes from all three approaches. Observations from the BOLD-MRI and photoacoustic data showed that clusters of voxels oscillated in a synchronized manner. Conclusion: We were able to use three different modalities to show that fluctuation in tumor oxygen is both coherent and incoherent, with periods of minutes to tens of minutes. These periods are very similar to those from the well-established metabolic, non-linear biomechanical phenomenon called the glycolytic oscillator. This may provide an additional explanation to the cause of cyclic hypoxia. Such dynamics could have profound implications in hypofractionated radiotherapy regiments and could help guide treatment and make it more patient specific. The authors would like to thank the University of Wisconsin Carbone Cancer Center (UWCCC) for the funds to complete this project. This work is also supported in part by NIH/NCI P30 CA014520- UW Comprehensive Cancer Center Support”.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Diagnosing United States hurricane landfall risk: An alternative to count-based methodologies: Diagnosing U.S. Hurricane Landfall Risk

Assessing hurricane landfall risk is of immense public utility, yet extant methods of diagnosing annual tropical cyclone (TC) activity demonstrate no skill in diagnosing U.S. hurricane landfalls. Atlantic TC count itself has limited skill, explaining less than 20% of interannual variance in landfall incidence. Using extended landfall activity and reanalysis data sets, we employed empirical Poisson modeling to produce a landfall diagnostic index (LDI), incorporating spatially and temporally averaged upper level divergence, relative sea surface temperature, meridional wind, and zonal shear vorticity. LDI captures 31% of interannual variability of U.S. hurricane landfalls and offers physical insight into why indices that successfully capture TC activity fail to diagnose landfalls: there is inherent tension between conditions likely to steer hurricanes toward the U.S. and conditions favorable for TC development. Given this tension, attempting to diagnose, predict, or understand TC count is inadequate for quantifying societal impacts due to landfalling hurricanes.

Geology↗
DOE OSTI2014

Model-based PSF and MTF estimation and validation from skeletal clinical CT images

Purpose: A method was developed to correct for systematic errors in estimating the thickness of thin bones due to image blurring in CT images using bone interfaces to estimate the point-spread-function (PSF). This study validates the accuracy of the PSFs estimated using said method from various clinical CT images featuring cortical bones. Methods: Gaussian PSFs, characterized by a different extent in the z (scan) direction than in the x and y directions were obtained using our method from 11 clinical CT scans of a cadaveric craniofacial skeleton. These PSFs were estimated for multiple combinations of scanning parameters and reconstruction methods. The actual PSF for each scan setting was measured using the slanted-slit technique within the image slice plane and the longitudinal axis. The Gaussian PSF and the corresponding modulation transfer function (MTF) are compared against the actual PSF and MTF for validation. Results: The differences (errors) between the actual and estimated full-width half-max (FWHM) of the PSFs were 0.09 ± 0.05 and 0.14 ± 0.11 mm for the xy and z axes, respectively. The overall errors in the predicted frequencies measured at 75%, 50%, 25%, 10%, and 5% MTF levels were 0.06 ± 0.07 and 0.06 ± 0.04 cycles/mm for the xy and z axes, respectively. The accuracy of the estimates was dependent on whether they were reconstructed with a standard kernel (Toshiba's FC68, mean error of 0.06 ± 0.05 mm, MTF mean error 0.02 ± 0.02 cycles/mm) or a high resolution bone kernel (Toshiba's FC81, PSF FWHM error 0.12 ± 0.03 mm, MTF mean error 0.09 ± 0.08 cycles/mm). Conclusions: The method is accurate in 3D for an image reconstructed using a standard reconstruction kernel, which conforms to the Gaussian PSF assumption but less accurate when using a high resolution bone kernel. The method is a practical and self-contained means of estimating the PSF in clinical CT images featuring cortical bones, without the need phantoms or any prior knowledge about the scanner-specific parameters.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Bias-free lateral terahertz emitters—A simulation study

The design and performance of bias-free InN-based THz emitters that exploit lateral photocurrents is studied by means of numerical simulations. We use a drift-diffusion model with adjusted carrier temperatures and mobilities. The applicability of this approach is demonstrated by a comparison with results from Monte-Carlo simulations. We consider a simple but robust lateral emitter concept using metal stripes with two different thicknesses with one of them being thin enough to be transparent for THz radiation. This arrangement can be easily multiplexed and the efficiency of this concept has already been demonstrated by experiment for GaAs substrates. In the present study, we consider InN, which is known to be an efficient photo-Dember emitter because of its superior transport properties. Our main focus is on the impact of the emitter design on the emission efficiency assuming different operation principles. Both the lateral photo-Dember (LPD) effect and built-in lateral field effects are considered. The appropriate choice of the metal stripe and window geometry as well as the impact of surface Fermi level pinning are investigated in detail, and design guidelines for efficient large area emitters using multiplexed structures are provided. We find that InN LPD emitters do not suffer from Fermi level pinning at the InN surface. The optimum emission efficiency is found for LPD emitter structures having 200 nm wide illumination windows and mask stripes. Emitter structures in which lateral electric fields are induced by the metal mask contacts can have a considerably higher efficiency than pure LPD emitters. In the best case, the THz emission of such structures is increased by one order of magnitude. Their optimum window size is 1 μm without the necessity of a partially transparent set of mask stripes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Assessment of interpatient heterogeneity in tumor radiosensitivity for nonsmall cell lung cancer using tumor-volume variation data

Purpose: In our previous work, the authors showed that a distribution of cell surviving fractionsS{sub 2} in a heterogeneous group of patients could be derived from tumor-volume variation curves during radiotherapy for head and neck cancer. In this research study, the authors show that this algorithm can be applied to other tumors, specifically in nonsmall cell lung cancer. This new application includes larger patient volumes and includes comparison of data sets obtained at independent institutions. Methods: Our analysis was based on two data sets of tumor-volume variation curves for heterogeneous groups of 17 patients treated for nonsmall cell lung cancer with conventional dose fractionation. The data sets were obtained previously at two independent institutions by using megavoltage computed tomography. Statistical distributions of cell surviving fractionsS{sub 2} and clearance half-lives of lethally damaged cells T{sub 1/2} have been reconstructed in each patient group by using a version of the two-level cell population model of tumor response and a simulated annealing algorithm. The reconstructed statistical distributions of the cell surviving fractions have been compared to the distributions measured using predictive assays in vitro. Results: Nonsmall cell lung cancer presents certain difficulties for modeling surviving fractions using tumor-volume variation curves because of relatively large fractional hypoxic volume, low gradient of tumor-volume response, and possible uncertainties due to breathing motion. Despite these difficulties, cell surviving fractionsS{sub 2} for nonsmall cell lung cancer derived from tumor-volume variation measured at different institutions have similar probability density functions (PDFs) with mean values of 0.30 and 0.43 and standard deviations of 0.13 and 0.18, respectively. The PDFs for cell surviving fractions S{sub 2} reconstructed from tumor volume variation agree with the PDF measured in vitro. Conclusions: The data obtained in this work, when taken together with the data obtained previously for head and neck cancer, suggests that the cell surviving fractionsS{sub 2} can be reconstructed from the tumor volume variation curves measured during radiotherapy with conventional fractionation. The proposed method can be used for treatment evaluation and adaptation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

The benzene radical anion: A computationally demanding prototype for aromatic anions

The benzene radical anion is studied with ab initio coupled-cluster theory in large basis sets. Unlike the usual assumption, we find that, at the level of theory investigated, the minimum energy geometry is non-planar with tetrahedral distortion at two opposite carbon atoms. The anion is well known for its instability to auto-ionization which poses computational challenges to determine its properties. Despite the importance of the benzene radical anion, the considerable attention it has received in the literature so far has failed to address the details of its structure and shape-resonance character at a high level of theory. Here, we examine the dynamic Jahn-Teller effect and its impact on the anion potential energy surface. We find that a minimum energy geometry of C{sub 2} symmetry is located below one D{sub 2h} stationary point on a C{sub 2h} pseudo-rotation surface. The applicability of standard wave function methods to an unbound anion is assessed with the stabilization method. The isotropic hyperfine splitting constants (A{sub iso}) are computed and compared to data obtained from experimental electron spin resonance experiments. Satisfactory agreement with experiment is obtained with coupled-cluster theory and large basis sets such as cc-pCVQZ.

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

The ADAQ framework: An integrated toolkit for data acquisition and analysis with real and simulated radiation detectors

Here, the ADAQ framework is a collection of software tools that is designed to streamline the acquisition and analysis of radiation detector data produced in modern digital data acquisition (DAQ) systems and in Monte Carlo detector simulations. The purpose of the framework is to maximize user scientific pro-ductivity by minimizing the effort and expertise required to fully utilize radiation detectors in a variety of scientific and engineering disciplines. By using a single set of tools to span the real and simulation domains, the framework eliminates redundancy and provides an integrated workflow for high-fidelity comparison between experimental and simulated detector performance. Built on the ROOT data analysis framework, the core of the ADAQ framework is a set of Cþþ and Python libraries that enable high-level control of digital DAQ systems and detector simulations with data stored into standardized binary ROOT files for further analysis. Two graphical user interface programs utilize the libraries to create powerful tools: ADAQAcquisition handles control and readout of real-world DAQ systems and ADAQAnalysis provides data analysis and visualization methods for experimental and simulated data. At present, the ADAQ framework supports digital DAQ hardware from CAEN S.p.A. and detector simulations performed in Geant4; however, the modular design will facilitate future extension to other manufacturers and simulation platforms.

detector simulation↗
DOE OSTI2013

Polyad quantum numbers and multiple resonances in anharmonic vibrational studies of polyatomic molecules

In the theory of anharmonic vibrations of a polyatomic molecule, mixing the zero-order vibrational states due to cubic, quartic and higher-order terms in the potential energy expansion leads to the appearance of more-or-less isolated blocks of states (also called polyads), connected through multiple resonances. Such polyads of states can be characterized by a common secondary integer quantum number. This polyad quantum number is defined as a linear combination of the zero-order vibrational quantum numbers, attributed to normal modes, multiplied by non-negative integer polyad coefficients, which are subject to definition for any particular molecule. According to Kellman's method [J. Chem. Phys. 93, 6630 (1990)], the corresponding formalism can be conveniently described using vector algebra. In the present work, a systematic consideration of polyad quantum numbers is given in the framework of the canonical Van Vleck perturbation theory (CVPT) and its numerical-analytic operator implementation for reducing the Hamiltonian to the quasi-diagonal form, earlier developed by the authors. It is shown that CVPT provides a convenient method for the systematic identification of essential resonances and the definition of a polyad quantum number. The method presented is generally suitable for molecules of significant size and complexity, as illustrated by several examples of molecules up to six atoms. The polyad quantum number technique is very useful for assembling comprehensive basis sets for the matrix representation of the Hamiltonian after removal of all non-resonance terms by CVPT. In addition, the classification of anharmonic energy levels according to their polyad quantum numbers provides an additional means for the interpretation of observed vibrational spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗
DOE OSTI2013

PHOTOMETRIC REDSHIFTS AND SYSTEMATIC VARIATIONS IN THE SPECTRAL ENERGY DISTRIBUTIONS OF LUMINOUS RED GALAXIES FROM SDSS DR7

We describe the construction of a template set of spectral energy distributions (SEDs) for the estimation of photometric redshifts of luminous red galaxies (LRGs) with a Bayesian template fitting method. By examining the color properties of several publicly available SED sets within a redshift range of 0 < z {approx}< 0.5 and comparing them to Sloan Digital Sky Survey (SDSS) Data Release 7 data, we show that only some of the investigated SEDs approximately match the colors of the LRG data throughout the redshift range, however not at the quantitative level required for precise photometric redshifts. This is because the SEDs of galaxies evolve with time (and redshift) and because at fixed redshift the LRG colors have an intrinsic spread such that they cannot be matched by one SED only. We generate new SEDs by superposing model SEDs of composite stellar populations with a burst model, allowing both components to be reddened by dust, in order to match the data in five different redshift bins. We select a set of SEDs which represents the LRG data in color space within five redshift bins, thus defining our new SED template set for photometric redshift estimates. The results we obtain with the new template set and our Bayesian template fitting photometric redshift code (PhotoZ) are nearly unbiased, with a scatter of {sigma}{sub {Delta}z} = 0.027 (including outliers), a fraction of catastrophic outliers (|z{sub phot} - z{sub spec}|/(1 + z{sub spec}) > 0.15) of {eta} = 0.12%, and a normalized median absolute rest frame deviation (NMAD) of {sigma}{sub NMAD} = 1.48 Multiplication-Sign MAD = 0.017 for non-outliers. We show that templates that optimally describe the brightest galaxies (-24.5 {<=} M{sub R} {<=} -22.7) indeed vary from z = 0.1 to z = 0.5, consistent with aging of the stellar population. Furthermore, we find that templates that optimally describe galaxies at z < 0.1 strongly differ as a function of the absolute magnitude of the galaxies, indicating an increase in star formation activity for less luminous galaxies. Our findings based on the photometry of the SDSS LRGs and our SED template fitting are supported by comparison to the average SDSS LRG spectra in different luminosity and redshift bins.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Transfer Learning Trained LSTM Models for Household Load Profile Forecasting

Grid edge renewable energy resources, such as rooftop solar photovoltaics, closely interact with consumer load profiles. Therefore, forecasting future electricity demand, ideally at the individual household level, is indispensable. In this paper, we present a transfer learning enhanced household load profile forecasting method. First, we tune a long short-term memory forecasting model to perform day-ahead prediction of household electricity load profiles. Then we improve these individualized models using transfer learning, and we use k-means clustering to create optimal source data sets. We find average improvements of 4.38% (largest improvement of 10.71%) when the entire data set was used to train the source model and 2.45% (largest improvement of 11.57%) in the mean absolute error when households were first clustered and used to train separate source models for each cluster. We find that transfer learning with clustered data can effectively boost the forecasting performance of the LSTM models. We use realistic household power measurements for 148 real residential households in Austin, Texas.

deep learning↗
DOE OSTI2011

Effects of exposure equalization on image signal-to-noise ratios in digital mammography: A simulation study with an anthropomorphic breast phantom

Purpose: The scan equalization digital mammography (SEDM) technique combines slot scanning and exposure equalization to improve low-contrast performance of digital mammography in dense tissue areas. In this study, full-field digital mammography (FFDM) images of an anthropomorphic breast phantom acquired with an anti-scatter grid at various exposure levels were superimposed to simulate SEDM images and investigate the improvement of low-contrast performance as quantified by primary signal-to-noise ratios (PSNRs). Methods: We imaged an anthropomorphic breast phantom (Gammex 169 ''Rachel,'' Gammex RMI, Middleton, WI) at various exposure levels using a FFDM system (Senographe 2000D, GE Medical Systems, Milwaukee, WI). The exposure equalization factors were computed based on a standard FFDM image acquired in the automatic exposure control (AEC) mode. The equalized image was simulated and constructed by superimposing a selected set of FFDM images acquired at 2, 1, 1/2, 1/4, 1/8, 1/16, and 1/32 times of exposure levels to the standard AEC timed technique (125 mAs) using the equalization factors computed for each region. Finally, the equalized image was renormalized regionally with the exposure equalization factors to result in an appearance similar to that with standard digital mammography. Two sets of FFDM images were acquired to allow for two identically, but independently, formed equalized images to be subtracted from each other to estimate the noise levels. Similarly, two identically but independently acquired standard FFDM images were subtracted to estimate the noise levels. Corrections were applied to remove the excess system noise accumulated during image superimposition in forming the equalized image. PSNRs over the compressed area of breast phantom were computed and used to quantitatively study the effects of exposure equalization on low-contrast performance in digital mammography. Results: We found that the highest achievable PSNR improvement factor was 1.89 for the anthropomorphic breast phantom used in this study. The overall PSNRs were measured to be 79.6 for the FFDM imaging and 107.6 for the simulated SEDM imaging on average in the compressed area of breast phantom, resulting in an average improvement of PSNR by {approx}35% with exposure equalization. We also found that the PSNRs appeared to be largely uniform with exposure equalization, and the standard deviations of PSNRs were estimated to be 10.3 and 7.9 for the FFDM imaging and the simulated SEDM imaging, respectively. The average glandular dose for SEDM was estimated to be 212.5 mrad, {approx}34% lower than that of standard AEC-timed FFDM (323.8 mrad) as a result of exposure equalization for the entire breast phantom. Conclusions: Exposure equalization was found to substantially improve image PSNRs in dense tissue regions and result in more uniform image PSNRs. This improvement may lead to better low-contrast performance in detecting and visualizing soft tissue masses and micro-calcifications in dense tissue areas for breast imaging tasks.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Characterization of a constrained paired-view technique in iterative reconstruction for breast tomosynthesis

Purpose: The order in which the projection views are employed in the reconstruction of tomosynthesis by iterative algorithms, such as simultaneous algebraic reconstruction technique and maximum likelihood, has a strong effect on the rate of convergence, accuracy, and the edge-blurring artifacts in the reconstructed image. The purpose of this investigation was to characterize and evaluate the effects of ordering schemes on image quality for breast tomosynthesis reconstruction and to explore a new constrained paired-view technique that could provide reduction of reconstruction artifacts. In this work, the authors compared several different ordering schemes and characterized the image quality and the formation of out-of-plane artifacts. Furthermore, a new normalization method is presented. It produces more accurate reconstructions with reduced artifacts comparing to the standard method of sequential ordering.Methods: In addition to visual assessment of image quality, several indices such as the signal-difference-to-noise ratio, the artifact-spread function, and the lesion detectability (d{sup ′}) were computed to quantitatively evaluate the effect of ordering scheme. The sets of breast tomosynthesis projection images were simulated for reconstruction; one set had uniform background (white noise only) and the other two contained both anatomic background and quantum noise. Clinical breast images were also studied for comparison.Results: The authors have quantified the image quality in reconstructed slices for a range of tumor sizes. The authors’ proposed method provides better performance for all of the metrics tested (contrast, d{sup ′}, and the level of artifacts) both for the uniform phantom case and in the presence of anatomical structure.Conclusions: The paired projection normalization provides better performance in the image quality of the reconstructed slices, and results in a lower level of artifacts in the Z direction. This implies that even a relatively simple method like the “side-to-side” sequence, which pairs two symmetrical projections with equal angular distance from the central projection, would achieve better reconstructed image quality than the conventional “step-by-step” method, which uses sequential projections one after another.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Comparison Study of Conventional and Adaptive Mesh Refinement in Organic Material Decomposition Models

This study compares conventional mesh refinement techniques, specifically Uniform Mesh Refinement (UMR), with a new Adaptive Mesh Refinement (AMR) method, applied to Organic Material Decomposition (OMD) models. The proposed benefit of AMR is that only areas that require refinement, based on minimizing a specific field gradient, are refined thus decreasing model wall time compared to conventional UMR methods. This work specifically focuses on comparing UMR and AMR methods on decomposing (both No-Flow and Porous-Flow material models) Polymeric Methylene Diisocyanate (PMDI) polyurethane foam. Throughout the work, the geometry increased in complexity to assess the refinement methods performance at varying levels geometric intricacy. While AMR has been shown to work well in a variety of applications, the UMR approach proved to be computationally faster, for many of the geometries and foam decomposition models, than AMR. However, it was observed that at higher levels of refinement, greater than 3 UMR, AMR begins to be computationally better. Additionally, the settings used to perform AMR greatly impact its performance, and lessons learned, in terms of OMD models, are shared. Due to physics involved in material decomposition, specifically the evolution of state variables, these problems don’t fully benefit from the advantages of AMR.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Method for Measuring Dry Density of Compacted Bentonite Mixture Soil Using 3D Laser Scanner data - 19119

The bentonite mixture soil is planned to be used for the low-permeability layer of low-level radioactive waste (LLW) disposal facilities. The construction of the bentonite mixture soil is similar to that of normal soil construction, and the rolling pressure and/or vibration compaction are used. The quality control of construction is the same as normal soil work, and the construction management has to confirm that the compaction degree (that is ratio of the dry density to the modified proctor curve) is above specified value. As a method of measuring the dry density after construction, the conventional methods are the core cutter method and the sand replacement method. The measurement principle of these methods is simple and the reliability of the measurement value is high, but during measurements digging holes disturb the construction place, and the measurement duration is a disadvantage in these methods. Therefore, we developed a non-destructive method that can make measurements in a short time. The developed method uses the correlation between the density change and the layer thickness change before and after compaction, and estimates the density after compaction from the height change obtained from the survey results. However, since the material is compressed while the cross-sectional area spreads, this cross-sectional area change is also taken into consideration. By using the point cloud data of the 3D laser scanner surveying, it is possible to obtain almost continuous data in a short period of time. Corrected dry density values calculated by the above procedure ρ'{sub 3D} are generally consistent with dry density after compaction ρ{sub core} by the core cutter method. This method differs depending on the boundary conditions of the bottom and the side of the construction place. Then, the reliability of the method was evaluated by applying it to two boundary conditions and comparing it with the core cutter method. The proposed estimation method showed good agreement with the core cutter method. Therefore, it was confirmed that it is possible to estimate dry density with high reliability. Although this method needs to be used in set with core cutter method, the measurement frequency of the core cutter method can be reduced, and the measurement of a vast area can be done in a short time without disturbance. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Spatial DIC Errors due to Pattern-Induced Bias and Grey Level Discretization

Digital image correlation (DIC) is an optical metrology method widely used in experimental mechanics for full-field shape, displacement and strain measurements. The required strain resolution for engineering applications of interest mandates DIC to have a high image displacement matching accuracy, on the order of 1/100th of a pixel, which necessitates an understanding of DIC errors. In this paper, we examine two spatial bias terms that have been almost completely overlooked. They cause a persistent offset in the matching of image intensities and thus corrupt DIC results. We name them pattern-induced bias (PIB), and intensity discretization bias (IDB). We show that the PIB error occurs in the presence of an undermatched shape function and is primarily dictated by the underlying intensity pattern for a fixed displacement field and DIC settings. The IDB error is due to the quantization of the gray level intensity values in the digital camera. In this paper we demonstrate these errors and quantify their magnitudes both experimentally and with synthetic images.

47 OTHER INSTRUMENTATION↗