Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic error”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Robust state estimation of feeding–blending systems in continuous pharmaceutical manufacturing

State estimation is a fundamental part of monitoring, control, and real-time optimization in continuous pharmaceutical manufacturing. For nonlinear dynamic systems with hard constraints, moving horizon estimation (MHE) can estimate the current state by solving a well-defined optimization problem where process complexities are explicitly considered as constraints. Traditional MHE techniques assume random measurement noise governed by some normal distributions. However, state estimates can be unreliable if noise is not normally distributed or measurements are contaminated with gross or systematic errors. In this paper, to improve the accuracy and robustness of state estimation, we incorporate robust estimators within the standard MHE skeleton, leading to an extended MHE framework. The proposed MHE approach is implemented on two pharmaceutical continuous feeding–blending system (FBS) configurations which include loss-in-weight (LIW) feeders and continuous blenders. Numerical results show that our MHE approach is robust to gross errors and can provide reliable state estimates when measurements are contaminated with outliers and drifts. Finally and moreover, the efficient solution of the MHE realized in this work, suggests feasible application of on-line state estimation on more complex continuous pharmaceutical processes.

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

A method for characterization of multiple dynamic constitutive parameters of FRCs

We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Dynamic load balancing for petascale quantum Monte Carlo applications: The Alias method

Diffusion Monte Carlo is the most accurate widely used Quantum Monte Carlo method for the electronic structure of materials, but it requires frequent load balancing or population redistribution steps to maintain efficiency and avoid accumulation of systematic errors on parallel machines. The load balancing step can be a significant factor affecting performance, and will become more important as the number of processing elements increases. We propose a new dynamic load balancing algorithm, the Alias Method, and evaluate it theoretically and empirically. An important feature of the new algorithm is that the load can be perfectly balanced with each process receiving at most one message. It is also optimal in the maximum size of messages received by any process. We also optimize its implementation to reduce network contention, a process facilitated by the low messaging requirement of the algorithm. Empirical results on the petaflop Cray XT Jaguar supercomputer at ORNL showing up to 30% improvement in performance on 120,000 cores. The load balancing algorithm may be straightforwardly implemented in existing codes. The algorithm may also be employed by any method with many near identical computational tasks that requires load balancing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Retrieval of water vapor column abundance and aerosol properties from ChemCam passive sky spectroscopy

In this work, we derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. This paper covers the methodology and initial results for water vapor and also provides preliminary results for aerosols. The data set presented here includes the results of 113 observations spanning from Mars Year 31 L s = 291° (March 30, 2013) to Mars Year 33 L s = 127° (March 24, 2016). Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO 2 absorptions and the known CO 2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. Our water vapor results give water vapor column abundance with a precision of ±0.6 precipitable microns and systematic errors no larger than ±0.3 precipitable microns, assuming uniform vertical mixing. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising interannual increase in water–ice cloud opacities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Analysis and improvement of the hot disk transient plane source method for low thermal conductivity materials

The hot disk transient plane source (TPS) method is a widely used standard technique (ISO 22007-2) for the characterization of thermal properties of materials, especially the thermal conductivity, k. Despite its well-established reliability for a wide variety of common materials, the hot disk TPS method is also known to suffer from a substantial systematic errors when applied to low-k thermal insulation materials. Here, we present a combined numerical and experimental study on the influence of the geometry of hot disk sensor on measured value of low-k materials. We demonstrate that the error is strongly affected by the finite thickness and thermal mass of the sensor's insulation layer was well as the corresponding increase of the effective heater size beyond the radius of the embedded metal heater itself. We also numerically investigate the dependence of the error on the sample thermal properties, confirming that the errors are worse in low-k samples. A simple correction function is also provided, which converts the apparent (erroneous) result from a standard hot disk TPS measurement to a more accurate value. A standard polyimide sensor was also optimized using both wet and dry etching to provide more accurate measurement directly. Experimentally corrected value of k for Airloy x56 aerogel and a commercial silica aerogel using the numerical correction factor derived based on the standard TPS sensor is in excellent agreement with the directly measured value from the TPS sensor using the optimized polyimide sensor. Finally, both of these methods can reduce the errors to less than 4% as compared to around 40% error of overestimation from raw values measured with the pristine sensor. Such results show that both the numerical correction to a pristine senor or an optimized sensor are capable of providing highly accurate value of thermal conductivity for such materials.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Deformable Image Registration for Adaptive Radiation Therapy of Head and Neck Cancer: Accuracy and Precision in the Presence of Tumor Changes

Purpose: To compare deformable image registration (DIR) accuracy and precision for normal and tumor tissues in head and neck cancer patients during the course of radiation therapy (RT). Methods and Materials: Thirteen patients with oropharyngeal tumors, who underwent submucosal implantation of small gold markers (average 6, range 4-10) around the tumor and were treated with RT were retrospectively selected. Two observers identified 15 anatomical features (landmarks) representative of normal tissues in the planning computed tomography (pCT) scan and in weekly cone beam CTs (CBCTs). Gold markers were digitally removed after semiautomatic identification in pCTs and CBCTs. Subsequently, landmarks and gold markers on pCT were propagated to CBCTs, using a b-spline-based DIR and, for comparison, rigid registration (RR). To account for observer variability, the pair-wise difference analysis of variance method was applied. DIR accuracy (systematic error) and precision (random error) for landmarks and gold markers were quantified. Time trend of the precisions for RR and DIR over the weekly CBCTs were evaluated. Results: DIR accuracies were submillimeter and similar for normal and tumor tissue. DIR precision (1 SD) on the other hand was significantly different (P<.01), with 2.2 mm vector length in normal tissue versus 3.3 mm in tumor tissue. No significant time trend in DIR precision was found for normal tissue, whereas in tumor, DIR precision was significantly (P<.009) degraded during the course of treatment by 0.21 mm/week. Conclusions: DIR for tumor registration proved to be less precise than that for normal tissues due to limited contrast and complex non-elastic tumor response. Caution should therefore be exercised when applying DIR for tumor changes in adaptive procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Head and Neck Margin Reduction With Adaptive Radiation Therapy: Robustness of Treatment Plans Against Anatomy Changes

Purpose: We set out to investigate loss of target coverage from anatomy changes in head and neck cancer patients as a function of applied safety margins and to verify a cone beam computed tomography (CBCT)–based adaptive strategy with an average patient anatomy to overcome possible target underdosage. Methods and Materials: For 19 oropharyngeal cancer patients, volumetric modulated arc therapy treatment plans (2 arcs; simultaneous integrated boost, 70 and 54.25 Gy; 35 fractions) were automatically optimized with uniform clinical target volume (CTV)–to–planning target volume margins of 5, 3, and 0 mm. We applied b-spline CBCT–to–computed tomography (CT) deformable registration to allow recalculation of the dose on modified CT scans (planning CT deformed to daily CBCT following online positioning) and dose accumulation in the planning CT scan. Patients with deviations in primary or elective CTV coverage >2 Gy were identified as candidates for adaptive replanning. For these patients, a single adaptive intervention was simulated with an average anatomy from the first 10 fractions. Results: Margin reduction from 5 mm to 3 mm to 0 mm generally led to an organ-at-risk (OAR) mean dose (D{sub mean}) sparing of approximately 1 Gy/mm. CTV shrinkage was mainly seen in the elective volumes (up to 10%), likely related to weight loss. Despite online repositioning, substantial systematic errors were present (>3 mm) in lymph node CTV, the parotid glands, and the larynx. Nevertheless, the average increase in OAR dose was small: maximum of 1.2 Gy (parotid glands, D{sub mean}) for all applied margins. Loss of CTV coverage >2 Gy was found in 1, 3, and 7 of 73 CTVs, respectively. Adaptive intervention in 0-mm plans substantially improved coverage: in 5 of 7 CTVs (in 6 patients) to <2 Gy of initially planned. Conclusions: Volumetric modulated arc therapy head and neck cancer treatment plans with 5-mm margins are robust for anatomy changes and show a modest increase in OAR dose. Margin reduction improves OAR sparing with approximately 1 Gy/mm at the expense of target coverage in a subgroup of patients. Patients at risk of CTV underdosage >2 Gy in 0-mm plans may be identified early in treatment using dose accumulation. A single intervention with an average anatomy derived from CBCT effectively mitigates discrepancies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Atomic Layer Thermopiles: Comprehensive static calibration, comparison and application in subsonic and supersonic flows

Here, this paper compares the calibration of different heat flux sensors in radiation- and convection-based sub- and supersonic operation. First, four heat flux sensors based on different principles: ALTP (Transverse Seebeck Effect), HFM-8E (differential-layer device), coaxial Thermocouple and a TG-2000 (circular-foil gage) are calibrated in a laser-based radiation setup. In a second step, all heat flux sensors are compared with a slug-calorimeter within a subsonic convection-dominated facility based on stagnation-point measurements of an impinging hot air jet. The obtained results indicate that a sensitivity transfer between a radiative calibrated sensor used in a mainly convective environment is not always possible and can lead to significant, systematic errors. Calibration in a subsonic shear flow with thermocouple readings demonstrate small divergence from the manufacturer provided sensitivity. Finally, supersonic testing with high frequency shock-boundary layer interactions highlight the need for ALTP rather than the conventional use of thermocouple arrays to resolve the high frequency phenomena associated with shock boundary layer interactions.

42 ENGINEERING↗
DOE OSTI2020

Effects of natural zeolites on field-scale geologic noble gas transport

Improving predictive models for noble gas transport through natural materials at the field-scale is an essential component of improving US nuclear monitoring capabilities. Several field-scale experiments with a gas transport component have been conducted at the Nevada National Security Site (Non-Proliferation Experiment, Underground Nuclear Explosion Signatures Experiment). However, the models associated with these experiments have not treated zeolite minerals as gas adsorbing phases. This is significant as zeolites are a common alteration mineral with a high abundance at these field sites and are shown here to significantly fractionate noble gases during field-scale transport. This fractionation and associated retardation can complicate gas transport predictions by reducing the signal-to-noise ratio to the detector (e.g. mass spectrometers or radiation detectors) enough to mask the signal or make the data difficult to interpret. Omitting adsorption-related retardation data of noble gases in predictive gas transport models therefore results in systematic errors in model predictions where zeolites are present.Herein is presented noble gas adsorption data collected on zeolitized and non-zeolitized tuff. Experimental results were obtained using a unique piezometric adsorption system designed and built for this study. Data collected were then related to pure-phase mineral analyses conducted on clinoptilolite, mordenite, and quartz. These results quantify the adsorption capacity of materials present in field-scale systems, enabling the modeling of low-permeability rocks as significant sorption reservoirs vital to bulk transport predictions.

58 GEOSCIENCES↗
DOE OSTI2017

Analyte-induced spectral filtering in femtosecond transient absorption spectroscopy

Here, we discuss the influence of spectral filtering by samples in femtosecond transient absorption measurements. Commercial instruments for transient absorption spectroscopy (TA) have become increasingly available to scientists in recent years and TA is becoming an established technique to measure the dynamics of photoexcited systems. Furthermore, we show that absorption of the excitation pulse by the sample can severely alter the spectrum and consequently the temporal pulse shape. This “spectral self-filtering” effect can lead to systematic errors and misinterpretation of data, most notably in concentration dependent measurements. Finally, the combination of narrow absorption peaks in the sample with ultrafast broadband excitation pulses is especially prone to this effect.

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

Multispectrum analysis of the oxygen A-band

Retrievals of atmospheric composition from near-infrared measurements require measurements of airmass to better than the desired precision of the composition. The oxygen bands are obvious choices to quantify airmass since the mixing ratio of oxygen is fixed over the full range of atmospheric conditions. The OCO-2 mission is currently retrieving carbon dioxide concentration using the oxygen A-band for airmass normalization. The 0.25% accuracy desired for the carbon dioxide concentration has pushed the required state-of-the-art for oxygen spectroscopy. To measure O2 A-band cross-sections with such accuracy through the full range of atmospheric pressure requires a sophisticated line-shape model (Rautian or Speed-Dependent Voigt) with line mixing (LM) and collision induced absorption (CIA). Models of each of these phenomena exist, however, this work presents an integrated self-consistent model developed to ensure the best accuracy. It is also important to consider multiple sources of spectroscopic data for such a study in order to improve the dynamic range of the model and to minimize effects of instrumentation and associated systematic errors. The techniques of Fourier Transform Spectroscopy (FTS) and Cavity Ring-Down Spectroscopy (CRDS) allow complimentary information for such an analysis. We utilize multispectrum fitting software to generate a comprehensive new database with improved accuracy based on these datasets. As a result, the extensive information will be made available as a multi-dimensional cross-section (ABSCO) table and the parameterization will be offered for inclusion in the HITRANonline database.

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

Couch height–based patient setup for abdominal radiation therapy

There are 2 methods commonly used for patient positioning in the anterior-posterior (A-P) direction: one is the skin mark patient setup method (SMPS) and the other is the couch height–based patient setup method (CHPS). This study compared the setup accuracy of these 2 methods for abdominal radiation therapy. The enrollment for this study comprised 23 patients with pancreatic cancer. For treatments (539 sessions), patients were set up by using isocenter skin marks and thereafter treatment couch was shifted so that the distance between the isocenter and the upper side of the treatment couch was equal to that indicated on the computed tomographic (CT) image. Setup deviation in the A-P direction for CHPS was measured by matching the spine of the digitally reconstructed radiograph (DRR) of a lateral beam at simulation with that of the corresponding time-integrated electronic portal image. For SMPS with no correction (SMPS/NC), setup deviation was calculated based on the couch-level difference between SMPS and CHPS. SMPS/NC was corrected using 2 off-line correction protocols: no action level (SMPS/NAL) and extended NAL (SMPS/eNAL) protocols. Margins to compensate for deviations were calculated using the Stroom formula. A-P deviation > 5 mm was observed in 17% of SMPS/NC, 4% of SMPS/NAL, and 4% of SMPS/eNAL sessions but only in one CHPS session. For SMPS/NC, 7 patients (30%) showed deviations at an increasing rate of > 0.1 mm/fraction, but for CHPS, no such trend was observed. The standard deviations (SDs) of systematic error (Σ) were 2.6, 1.4, 0.6, and 0.8 mm and the root mean squares of random error (σ) were 2.1, 2.6, 2.7, and 0.9 mm for SMPS/NC, SMPS/NAL, SMPS/eNAL, and CHPS, respectively. Margins to compensate for the deviations were wide for SMPS/NC (6.7 mm), smaller for SMPS/NAL (4.6 mm) and SMPS/eNAL (3.1 mm), and smallest for CHPS (2.2 mm). Achieving better setup with smaller margins, CHPS appears to be a reproducible method for abdominal patient setup.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Multiple-scattering Corrections to Measurements of the Prompt Fission Neutron Spectrum

The Chi-Nu project, conducted jointly by LANL and LLNL, aims to measure the shape of the prompt fission neutron spectrum (PFNS) for fission of {sup 239}Pu induced by neutrons from 50 keV to 15 MeV with accuracies of 3–5% in the outgoing energy from 50 keV to 9 MeV and 15% from 9 to 15 MeV. In order to meet this goal, detailed Monte Carlo simulations are being used to assess the importance and effect of every component in the experimental configuration. As part of this effort, we have also simulated some past PFNS measurements to identify possible sources of systematic error. We find that multiple scattering plays an important role in the target geometry, collimators, and detector response and that past experiments probably underestimated the extent of this effect.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Refining instrument response functions with 3-D Monte Carlo simulations of differential hard x-ray spectrometers

Three x-ray spectrometers – the hard x-ray detector (HXRD), the bremsstrahlung MeV x-ray spectrometer (BMXS), and the hard x-ray image plate (HXIP) – areused here to estimate hot-electron bremsstrahlung spectrum parameters for inertial confinement fusion experiments on OMEGA. Using 1-D deterministically calculated response functions, they show significantly different estimates of spectrum slope temperature. For all three detectors, 3-D Geant4 Monte Carlo simulations are developed and used to compute instrument response functions. Simulations explain features in raw image-plate data, but only significantly improve systematic errors in predicted detector channel signals for BMXS and do not entirely eliminate the discrepancy between estimates of slope temperature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

A modular apparatus for use in high-precision measurements of parity violation in polarized eV neutron transmission

We describe a modular apparatus for use in parity-violation measurements in epithermal neutron–nucleus resonances with high instantaneous neutron fluxes at the Manuel Lujan Jr. Neutron Scattering Center at Los Alamos National Laboratory. This apparatus is designed to conduct high-precision measurements of the parity-odd transmission asymmetry of longitudinally polarized neutrons through targets containing nuclei with p-wave neutron–nucleus resonances in the 0.1–10 eV energy regime and to accommodate a future search for time reversal violation in polarized neutron transmission through polarized nuclear targets. The apparatus consists of an adjustable neutron and gamma collimation system, a 3 He- 4 He ion chamber neutron flux monitor, two identical cryostats for target cooling, an adiabatic eV-neutron spin flipper, a near-unit efficiency 6 Li- 7 Li scintillation detector operated in current mode, a flexible CAEN data acquisition system, and a neutron spin filter based on spin-exchange optical pumping of 3 He gas. We describe the features of the apparatus design devoted to the suppression of systematic errors in parity-odd asymmetry measurements. We describe the configuration of the apparatus used to conduct a precision measurement of parity violation at the 0.7 eV p-wave resonance in 139 La which employs two identical 139 La targets, one to polarize the beam on the p-wave resonance using the weak interaction and one to analyze the polarization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Asymmetric errors

We present a procedure for handling asymmetric errors. Many results in particle physics are presented as values with different positive and negative errors, and there is no consistent procedure for handling them. We consider the difference between errors quoted, using pdfs and using likelihoods, and the difference between the rms spread of a measurement and the 68% central confidence region. We provide a comprehensive analysis of the possibilities, and software tools to enable their use.

Asymmetric↗
DOE OSTI2022

Harmonizing solar induced fluorescence across spatial scales, instruments, and extraction methods using proximal and airborne remote sensing: A multi-scale study in a soybean field

Solar-induced chlorophyll fluorescence (SIF) has been widely used to track vegetation photosynthesis at different scales ranging from in-situ measurements to satellite products. Airborne platforms sample SIF data at a spatial scale intermediate between in-situ and satellite, matching that of ground measurement (e.g. flux tower footprints and other field sampling), enabling us to explore causes of SIF variation and validate satellite-based SIF products. However, harmonizing SIF across sensors and platforms (correcting for systematic errors to yield a consistent, comparable SIF product) is challenging because SIF can be retrieved in different absorption windows, with different instruments and methods complicating the comparison between different observational levels (i.e., ground, airborne, satellites) and between sites equipped with different instruments with varying optical properties (spectral resolution and sampling intervals, spatial resolution). Additionally, the spatial and temporal variability of atmospheric properties can influence the retrieval of the weak SIF signal. Because of these complications, direct comparisons of airborne and ground SIF across scales are rarely attempted. Here, in this study, we combined airborne SIF data with simultaneous ‘ground truth’ data collected by stationary and mobile platforms in a soybean field in Nebraska, USA. In this effort, we tested several SIF extraction methods, including Fraunhofer Line Discrimination (FLD), improved Fraunhofer Line Discrimination (iFLD), Spectral Fitting Method (SFM), SpecFit, and a Singular Vector Decomposition (SVD) method. The SpecFit method was sensitive to the 715–740 nm water bands and removing the water bands in the fitting process yielded better agreement between the airborne and ground SIF spectra. Accurate estimation of the ground level downwelling irradiance obtained by ground measurements over a calibration target improved agreement between airborne and ground SIF retrievals at the O 2 A band, and allowed us to derive a SIF dataset with improved agreement across platforms and sampling scales. This experimental approach provided a method for generating comparable SIF signals across instruments, methods and platforms, which is critical to understanding the SIF-GPP relationship at different scales and to cross-validate the diversity of platforms used for satellite products calibration and validation.

54 ENVIRONMENTAL SCIENCES↗