Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian process”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stimulated emission tomography for efficient characterization of spatial entanglement

Stimulated-emission tomography (SET) is an excellent tool for characterizing the process of spontaneous parametric down conversion (SPDC), which is commonly used to create pairs of entangled photons for use in quantum information protocols. The use of stimulated emission increases the average number of detected photons by several orders of magnitude compared to the spontaneous process. In a SET measurement, the parametric down conversion is seeded by an intense signal field prepared with specified mode properties rather than by broadband multimodal vacuum fluctuations, as is the case for the spontaneous process. The SET process generates an intense idler field in a mode that is the complex conjugate to the signal mode. In this work we use SET to estimate the joint spatial mode distribution (JSMD) in the Laguerre-Gaussian (LG) basis of the two photons of an entangled photon pair. The pair is produced by parametric down conversion in a beta barium borate (BBO) crystal with type-II phase matching pumped at a wavelength of 405 nm along with a 780 nm seed signal beam prepared in a variety of LG modes to generate an 842 nm idler beam of which the spatial mode distribution is measured. We observe strong idler production and good agreement with the theoretical prediction of its spatial mode distribution. Our experimental procedure should enable the efficient determination of the photon-pair wavefunctions produced by low-brightness SPDC sources and the characterization of high-dimensional entangled-photon pairs. Published by the American Physical Society 2024

Xu, Yang (ORCID:0009000954547320)↗
DOE OSTI2016

TH-CD-209-08: Quantification of the Interplay Effect in Proton Pencil Beam Scanning Treatment of Lung

Purpose: To quantify the dose degradation caused by the interplay effect based on a beam specific motion analysis in proton pencil beam scanning (PBS) treatment of lung tumors Methods: PBS plans were optimized on average CT using a beam-specific PTV method for 10 consecutive patients with locally advanced non-small-cell-lung-cancer (NSCLC) treated with proton therapy to 6660/180 cGy. End inhalation (CT0) and end exhalation (CT50) were selected as the two extreme scenarios to acquire the relative stopping power ratio difference (Δrsp) for a respiration cycle. The water equivalent difference (ΔWET) per radiological path was calculated from the surface of patient to the iCTV by integrating the Δrsp of each voxel. The motion magnitude of each voxel within the target follows a quasi-Gaussian distribution. A motion index (MI (>5mm WET)), defined as the percentage of target voxels with an absolute integral ΔWET larger than 5 mm, was adopted as a metric to characterize interplay. To simulate the treatment process, 4D dose was calculated by accumulating the spot dose on the corresponding respiration phase to the reference phase CT50 by deformable image registration based on spot timing and patient breathing phase. Results: The study indicated that the magnitude of target underdose in a single fraction plan is proportional to the MI (p<0.001), with larger motion equating to greater dose degradation and standard deviations. The target homogeneity, minimum, maximum and mean dose in the 4D dose accumulations of 37 fractions varied as a function of MI. Conclusion: The MI quantification metric can predict the level of dose degradation in PBS lung cancer treatment, which potentially serves as a clinical decision tool to assess whether patients are suitable to receive PBS treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Resonant Slow Extraction Simulation using Bmad

Simulations of slow extraction and transport of charged particle beams from synchrotrons requires careful modeling and experimentation. In this paper, we outline how the Bmad modeling software was adapted and developed to run third-integer resonant extraction simulations of beams at Brookhaven National Laboratory’s Booster synchrotron. Further, we show experimental comparisons of the slowly extracted beams transferred to the NASA Space Radiation Laboratory (NSRL) transport line. In this process, beam passes through a stripping foil element at the extraction point, which, along with stripping remaining electrons from the beam ions, acts as a scatterer to modify the phase space, moving to a more Gaussian-like distribution. This modification to the beam helps generate a uniform beam at the beam line’s target location; which is necessary for the variety of experiments performed at the facility. During this work, beam energy loss and multiple scattering by foil routines were built in conjunction with the Bmad code developers and are now integrated into the software.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

NEW SUNS IN THE COSMOS. III. MULTIFRACTAL SIGNATURE ANALYSIS

In the present paper, we investigate the multifractality signatures in hourly time series extracted from the CoRoT spacecraft database. Our analysis is intended to highlight the possibility that astrophysical time series can be members of a particular class of complex and dynamic processes, which require several photometric variability diagnostics to characterize their structural and topological properties. To achieve this goal, we search for contributions due to a nonlinear temporal correlation and effects caused by heavier tails than the Gaussian distribution, using a detrending moving average algorithm for one-dimensional multifractal signals (MFDMA). We observe that the correlation structure is the main source of multifractality, while heavy-tailed distribution plays a minor role in generating the multifractal effects. Our work also reveals that the rotation period of stars is inherently scaled by the degree of multifractality. As a result, analyzing the multifractal degree of the referred series, we uncover an evolution of multifractality from shorter to larger periods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Gaussian mixture models as automated particle classifiers for fast neutron detectors

Pulse shape discrimination (PSD) is the task of classifying electronic pulse shapes for different particle types such as gamma rays and fast neutrons interacting in scintillators and read out by photo sensitive detectors. This field has been limited in its adoption of techniques found in the statistical learning community. Methods initially employed in the 1960s for analog electronic circuitry persist in the current PSD literature describing operations performed on digitized pulses, which are amenable to statistical rigor. Despite vast amounts of data collected at low energy levels, traditional PSD methods are unable to discriminate particles below a certain threshold. As such, in this work, Gaussian mixture models (GMMs) are used as a clustering technique for fast neutron detection in the absence of labeled data. GMMs yield improvements spanning the energy spectrum in a desirably efficient, unsupervised fashion. An extension, the Dirichlet Process GMM, provides further flexibility and classification improvement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Reliable noise measure in time-gated NMR data

Time gating is a commonly used approach in the pre-processing of nuclear magnetic resonance data before inversion. Gating suppresses spurious signals that can degrade recovered decay time distributions and therefore often stabilizes inversion. However, care must be taken in applying this technique when assumptions of uncorrelated Gaussianity break down and reduce stacking efficiency. If not properly accounted for, unreliable noise estimates introduce inversion artefacts through over- or underfitting of the data. Block-bootstrap resampling of noise realization proxies obtained through data phasing can be used to generate reliable noise estimates for the windows. Benefits of the approach are demonstrated through inversion of synthetic and borehole data. Finally, analysis confirms that bootstrapped noise metrics are more reliable under variable noise conditions and result in more stable inversion results.

58 GEOSCIENCES↗
NASA NTRS1974

Noncoherent detection of periodic signals

The optimal Bayes detector for a general periodic waveform having uniform delay and additive white Gaussian noise is examined. It is shown that the detector is much more complex than that for the well known cases of pure sine waves (i.e. classical noncoherent detection) and narrowband signals. An interpretation of the optimal processing is presented, and several implementations are discussed. The results have application to the noncoherent detection of optical square waves.

Gagliardi, R. M.↗
DOE OSTI2015

SU-E-J-130: Automating Liver Segmentation Via Combined Global and Local Optimization

Purpose: The aim of this work is to develop a robust algorithm for accurate segmentation of liver with special attention paid to the problems with fuzzy edges and tumor. Methods: 200 CT images were collected from radiotherapy treatment planning system. 150 datasets are selected as the panel data for shape dictionary and parameters estimation. The remaining 50 datasets were used as test images. In our study liver segmentation was formulated as optimization process of implicit function. The liver region was optimized via local and global optimization during iterations. Our method consists five steps: 1)The livers from the panel data were segmented manually by physicians, and then We estimated the parameters of GMM (Gaussian mixture model) and MRF (Markov random field). Shape dictionary was built by utilizing the 3D liver shapes. 2)The outlines of chest and abdomen were located according to rib structure in the input images, and the liver region was initialized based on GMM. 3)The liver shape for each 2D slice was adjusted using MRF within the neighborhood of liver edge for local optimization. 4)The 3D liver shape was corrected by employing SSR (sparse shape representation) based on liver shape dictionary for global optimization. Furthermore, H-PSO(Hybrid Particle Swarm Optimization) was employed to solve the SSR equation. 5)The corrected 3D liver was divided into 2D slices as input data of the third step. The iteration was repeated within the local optimization and global optimization until it satisfied the suspension conditions (maximum iterations and changing rate). Results: The experiments indicated that our method performed well even for the CT images with fuzzy edge and tumors. Comparing with physician delineated results, the segmentation accuracy with the 50 test datasets (VOE, volume overlap percentage) was on average 91%–95%. Conclusion: The proposed automatic segmentation method provides a sensible technique for segmentation of CT images. This work is supported by NIH/NIBIB (1R01-EB016777), National Natural Science Foundation of China (No.61471226 and No.61201441), Research funding from Shandong Province (No.BS2012DX038 and No.J12LN23), and Research funding from Jinan City (No.201401221 and No.20120109)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

A Collimated CZT Detector for Quantitative Gamma Assays - 17250

Cadmium-Zinc-Telluride [CZT] detectors for gamma spectroscopic measurements and room-temperature operation have had a long history of promising good performance. Very small detectors [0.05-0.5 cm{sup 3}] have been commercially available for many years, but their small size creates problems for spectroscopy because of very low photopeak efficiency, especially at medium and high energies. Large CZT [4-6 cm{sup 3}] detectors have been shown to work nicely in a research environment when complicated electronic and data processing algorithms are applied. We report here test results from a commercially available 1 cubic centimeter CZT detector housed in a portable custom Tungsten shield with a set of collimators with different apertures. The detector has 2-2.5% FWHM resolution at Cs-137 energies, and near-Gaussian shaped peaks. Efficiency calibrations can be quickly made for a wide variety of geometries using the ISOCS mathematical efficiency software. We show successful validation tests comparing accuracy of measurements on reference calibration point sources, line sources, and 200 liter drums. The low weight of the device [9 kg including the shield and collimator] and low power [powered via USB from the PC] makes these devices particularly suitable for human-deployed in-situ measurements of large objects containing medium levels of activity - e.g. characterization of objects encountered during decommissioning, emergency response, truck alarm investigation, or waste being shipped to a disposal facility. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Spatiotemporally dynamic implicit large eddy simulation using machine learning classifiers

In this article, we utilize machine learning to dynamically determine if a point on the computational grid requires implicit numerical dissipation for large eddy simulation (LES). The decision making process is learnt through a priori training on quantities derived from direct numerical simulation (DNS) data. In particular, we compute eddy-viscosities obtained through the coarse-graining of DNS quantities and utilize their projection onto a Gaussian distribution to categorize areas that may require dissipation. If our learning determines that closure is necessary, an upwinded scheme is utilized for computing the non-linear Jacobian. In contrast, if it is determined that closure is unnecessary, a symmetric and second-order accurate energy and enstrophy preserving Arakawa scheme is utilized instead. This results in a closure framework that precludes the specification of any model-form for the small scale contributions of turbulence but deploys an appropriate numerical dissipation from explicit closure driven hypotheses. Here, this methodology is deployed for the Kraichnan turbulence test-case and assessed through various statistical quantities such as angle-averaged kinetic energy spectra and vorticity structure functions. Our framework thus establishes a link between the use of explicit LES ideologies for closure and numerical dissipation-based modeling of turbulence leading to improved statistical fidelity of a posteriori simulations.

42 ENGINEERING↗
DOE OSTI2017

A backward Monte-Carlo method for time-dependent runaway electron simulations

Kinetic descriptions of runaway electrons (REs) are usually based on Fokker-Planck models that determine the probability distribution function of REs in 2-dimensional momentum space. Despite the simplification involved, the Fokker-Planck equation can rarely be solved analytically and direct numerical approaches [e.g., continuum and particle-based Monte Carlo (MC)] can be time consuming, especially in the computation of asymptotic-type observables including the runaway probability, the slowing-down and runaway mean times, and the energy limit probability. In this study, we present a novel backward MC approach to these problems based on backward stochastic differential equations that describe the dynamics of the runaway probability by means of the Feynman-Kac theory. The key ingredient of the backward MC algorithm is to place all the particles in a runaway state and simulate them backward from the terminal time to the initial time. As such, our approach can provide much faster convergence than direct MC methods (by significantly reducing the number of particles required to achieve a prescribed accuracy) while at the same time maintaining the advantages of particle-based methods (compared to continuum approaches). The proposed algorithm is unconditionally stable and can be parallelized as easy as the direct MC method, and its extension to dimensions higher than two is straightforward, thus paving the way for conducting large-scale RE simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Relative permeability as a stationary process: Energy fluctuations in immiscible displacement

Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work, we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ([Formula: see text]), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2004

Antenna Linear-Quadratic-Gaussian (LQG) Controllers: Properties, Limits of Performance, and Tuning Procedure

Wind gusts are the main disturbances that depreciate tracking precision of microwave antennas and radiotelescopes. The linear-quadratic-Gaussian (LQG) controllers - as compared with the proportional-and-integral (PI) controllers significantly improve the tracking precision in wind disturbances. However, their properties have not been satisfactorily understood; consequently, their tuning is a trial-and-error process. A control engineer has two tools to tune an LQG controller: the choice of coordinate system of the controller model and the selection of weights of the LQG performance index. This article analyzes properties of an open- and closed-loop antenna. It shows that the proper choice of coordinates of the open-loop model simplifies the shaping of the closed-loop performance. The closed-loop properties are influenced by the LQG weights. The article shows the impact of the weights on the antenna closed-loop bandwidth, disturbance rejection properties, and antenna acceleration. The bandwidth and the disturbance rejection characterize the antenna performance, while the acceleration represents the performance limit set by the antenna hardware (motors). The article presents the controller tuning procedure, based on the coordinate selection and the weight properties. The procedure rationally shapes the closed-loop performance, as an alternative to the trial-and-error approach.

Gawronski, W.↗
NASA NTRS2018

Realistic Covariance Generation for the GPM Spacecraft

We present several different methods for generating realistic predictive covariance, including Monte-Carlo simulations and more direct linear methods which require the addition of process noise. The Monte-Carlo simulation starts with an epoch uncertainty sample basis and propagates each trial to a time of interest in the future. The variance-covariance of the state elements as well as other higher order sample statistics can be readily computed from the propagated sample. While this method preserves the nonlinear effects on the propagated uncertainty, it is computationally intensive as a statistically significant sample size must be considered in the propagation process. Moreover, if the epoch covariance is optimistic, this can result in an underestimation of the prediction error. Another method is to propagate a state sensitivity matrix simultaneously with the satellite state, which allows the epoch covariances to be propagated forward in a linear fashion. This method does not preserve the non-linearity of the satellite state uncertainty but is much less computationally intensive. The propagated state covariance is the scaled to represent the realistic level of GPM state uncertainty via a "e-tuning process." The tuning process generates an inflation factor based on the observed error statistics of the predictive satellite trajectories when compared to the definitive ones. Difference tuning strategies are considered and compared via Goodness-of-Fit method testing for the Gaussian properties of the scaled covariance.

prediction error↗
DOE OSTI2026

Silicon-On-Sapphire Metasurfaces Generate Arrays of Dark and Bright Traps for Neutral Atoms

We demonstrated crystalline silicon-on-sapphire (c-SOS) metasurfaces that convert a Gaussian beam into arrays of complex optical traps, including arrays of optical bottle beams that trap atoms in dark regions interleaved with bright tweezer arrays. The high refractive index and indirect band gap of crystalline silicon make it possible to design high-resolution near-infrared (λ > 700 nm) metasurfaces that can be manufactured at scale using CMOS-compatible processes. Compared with active components like spatial light modulators (SLMs) that have become widely used to generate trap arrays, metasurfaces provide an indefinitely scalable number of pixels, enabling large arrays of complex traps in a very small form factor, as well as reduced dynamic noise. To design metasurfaces that can generate three-dimensional bottle beams to serve as dark traps, we modified the Gerchberg-Saxton algorithm to enforce complex-amplitude profiles at the focal plane of the metasurface and to optimize the uniformity of the traps across the array. We fabricated and measured c-SOS metasurfaces that convert a Gaussian laser beam into arrays of bright traps, dark traps, and interleaved bright/dark traps.

Gerchberg-Saxton↗
DOE OSTI2016

Parameterization of photon beam dosimetry for a linear accelerator

Purpose: In radiation therapy, accurate data acquisition of photon beam dosimetric quantities is important for (1) beam modeling data input into a treatment planning system (TPS), (2) comparing measured and TPS modeled data, (3) the quality assurance process of a linear accelerator’s (Linac) beam characteristics, (4) the establishment of a standard data set for comparison with other data, etcetera. Parameterization of the photon beam dosimetry creates a data set that is portable and easy to implement for different applications such as those previously mentioned. The aim of this study is to develop methods to parameterize photon beam dosimetric quantities, including percentage depth doses (PDDs), profiles, and total scatter output factors (S{sub cp}). Methods: S{sub cp}, PDDs, and profiles for different field sizes, depths, and energies were measured for a Linac using a cylindrical 3D water scanning system. All data were smoothed for the analysis and profile data were also centered, symmetrized, and geometrically scaled. The S{sub cp} data were analyzed using an exponential function. The inverse square factor was removed from the PDD data before modeling and the data were subsequently analyzed using exponential functions. For profile modeling, one halfside of the profile was divided into three regions described by exponential, sigmoid, and Gaussian equations. All of the analytical functions are field size, energy, depth, and, in the case of profiles, scan direction specific. The model’s parameters were determined using the minimal amount of measured data necessary. The model’s accuracy was evaluated via the calculation of absolute differences between the measured (processed) and calculated data in low gradient regions and distance-to-agreement analysis in high gradient regions. Finally, the results of dosimetric quantities obtained by the fitted models for a different machine were also assessed. Results: All of the differences in the PDDs’ buildup and the profiles’ penumbra regions were less than 2 and 0.5 mm, respectively. The differences in the low gradient regions were 0.20% ± 0.20% (<1% for all) and 0.50% ± 0.35% (<1% for all) for PDDs and profiles, respectively. For S{sub cp} data, all of the absolute differences were less than 0.5%. Conclusions: This novel analytical model with minimum measurement requirements was proved to accurately calculate PDDs, profiles, and S{sub cp} for different field sizes, depths, and energies.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

The predominance of Kohn–Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn–Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

High-energy sub-phonon lifetime pulse compression by stimulated Brillouin scattering in liquids

Multi-Joule level stimulated Brillouin scattering (SBS) pulse compression below the acoustic phonon lifetime is demonstrated with a energy-scalable generator-amplifier setup. Single-pass compression of pulses longer than 20τ B (τ B as phonon lifetime) to as short as 0.5τ B with ~100 mJ pulse energy is realized from the generator, by choosing the focusing length to match approximately with the full length at half maximum of the input Gaussian pulses. The interaction length is identified, both experimentally and numerically, as the key parameter in achieving sub-phonon lifetime pulse compression, with its main mechanism being the steepening of the Stokes pulse trailing edge via energy exchange process. After combining the generator with an amplifier that involves only collimated beams and serves as energy booster, the compression of 9 ns, 2 J pulses at 532 nm into 170 ps, 1.3 J per pulse is achieved in water, with very good stability in both pulse energy and duration. Furthermore, this work demonstrates for the first time the efficient high-energy SBS sub-phonon lifetime pulse compression, and paves a way to the reliable generation of sub-200 ps laser pulses with Joule-level energy.

74 ATOMIC AND MOLECULAR PHYSICS↗