SEARCH · Engineering Papers

Results for “temporal computing”

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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Downscaling SSP-consistent global spatial urban land projections from 1/8-degree to 1-km resolution 2000–2100

Abstract Long-term, spatial urban land projections that simultaneously offer global coverage and local-scale empirical accuracy are rare. Recently a set of such projections was produced using data-science-based simulations and the Shared Socioeconomic Pathways (SSPs). These projections update at decadal time intervals from 2000 to 2100 with a spatial resolution of 1/8 degree, while many socio-environmental studies customarily run their analysis and modelling at finer spatial resolutions, e.g. 1-km. Here we develop and validate an algorithm to downscale the 1/8-degree spatial urban land projections to the 1-km resolution. The algorithm uses an iterative process to allocate the decadal amount of urban land expansion originally projected for each 1/8-degree grid to its constituent 1-km grids. The results are a set of global maps showing urban land fractions at the 1-km resolution, updated at decadal intervals from 2000 to 2100, under five different urban land expansion scenarios consistent with the SSPs. The data can support studies of potential interactions between future urbanization and environmental changes across spatial and temporal scales.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

DOME: Directional medical embedding vectors from Electronic Health Records

Motivation: The increasing availability of Electronic Health Record (EHR) systems has created enormous potential for translational research. Recent developments in representation learning techniques have led to effective large-scale representations of EHR concepts along with knowledge graphs that empower downstream EHR studies. However, most existing methods require training with patient-level data, limiting their abilities to expand the training with multi-institutional EHR data. On the other hand, scalable approaches that only require summary-level data do not incorporate temporal dependencies between concepts. Methods: We introduce a DirectiOnal Medical Embedding (DOME) algorithm to encode temporally directional relationships between medical concepts, using summary-level EHR data. Specifically, DOME first aggregates patient-level EHR data into an asymmetric co-occurrence matrix. Then it computes two Positive Pointwise Mutual Information (PPMI) matrices to correspondingly encode the pairwise prior and posterior dependencies between medical concepts. Following that, a joint matrix factorization is performed on the two PPMI matrices, which results in three vectors for each concept: a semantic embedding and two directional context embeddings. They collectively provide a comprehensive depiction of the temporal relationship between EHR concepts. Results: We highlight the advantages and translational potential of DOME through three sets of validation studies. First, DOME consistently improves existing direction-agnostic embedding vectors for disease risk prediction in several diseases, for example achieving a relative gain of 5.5% in the area under the receiver operating characteristic (AUROC) for lung cancer. Second, DOME excels in directional drug-disease relationship inference by successfully differentiating between drug side effects and indications, correspondingly achieving relative AUROC gain over the state-of-the-art methods by 10.8% and 6.6%. Finally, DOME effectively constructs directional knowledge graphs, which distinguish disease risk factors from comorbidities, thereby revealing disease progression trajectories. The source codes are provided at https://github.com/celehs/Directional-EHRembedding.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗

DOE OSTI2018

Evaluation of non-invasive biomonitoring of 2,4-Dichlorophenoxyacetic acid (2,4-D) in saliva

The objective of this study was to evaluate the potential for noninvasive biomonitoring of 2,4-dichlorophenoxyacetic acid (2,4-D) in saliva. Using an in vitro rat salivary gland epithelial cell (SGEC) system, a collection of experiments investigating chemical protein binding, temporal and directional transport, as well as competitive transport with para-aminohippuric acid (PAH), a substrate for renal organic anion transporters, was conducted to identify cellular transport parameters required to computationally model salivary transport of 2,4-D. Additionally, a physiological protein gradient was implemented to mimic physiologically relevant concentrations of protein in rat plasma and saliva, and under these conditions the transfer of 2,4-D was markedly slower, driven by increased protein binding (i.e. reduced free 2,4-D species available to cross salivary barrier). The rate of transfer was directly proportional to the amount of unbound 2,4-D and demonstrated no indication of active transport. An in vivo assessment of 2,4-D exposure in rats revealed a non-linear protein binding response in plasma, indicating saturated protein binding and increased levels of free 2,4-D species at higher doses. A strong correlation between concentrations in saliva and free unbound species in plasma was observed (Pearson correlation coefficient = 0.95). Saliva:plasma 2,4-D ratios measured in vivo (0.0079) were comparable over dose, but significantly different than ratios measured in vitro (physiological conditions) (0.034), possibly due to 2,4-D concentrations in saliva not being at equilibrium with 2,4-D concentrations in blood, as well as physiological features absent in in vitro settings (e.g. blood flow). We demonstrated that 2,4-D is consistently transported into saliva using both in vitro and in vivo models, making 2,4-D a potential candidate for human non-invasive salivary biomonitoring. Further work is needed to understand whether current sensor limits of detection are sufficient to measure occupationally relevant exposures.

Carver, Zana A.↗

DOE OSTI2017

Hydropower Optimization Using Artificial Neural Network Surrogate Models of a High-Fidelity Hydrodynamics and Water Quality Model

Hydropower operations optimization subject to environmental constraints is limited by challenges associated with dimensionality and spatial and temporal resolution. The need for high-fidelity hydrodynamic and water quality models within optimization schemes is driven by improved computational capabilities, increased requirements to meet specific points of compliance with greater resolution, and the need to optimize operations of not just single reservoirs but systems of reservoirs. This study describes an important advancement for computing hourly power generation schemes for a hydropower reservoir using high-fidelity models, surrogate modeling techniques, and optimization methods. The predictive power of the high-fidelity hydrodynamic and water quality model CE-QUAL-W2 is successfully emulated by an artificial neural network, then integrated into a genetic algorithm optimization approach to maximize hydropower generation subject to constraints on dam operations and water quality. This methodology is applied to a multipurpose reservoir near Nashville, Tennessee, USA. The model successfully reproduced high-fidelity reservoir information while enabling 6.8% and 6.6% increases in hydropower production value relative to actual operations for dissolved oxygen (DO) limits of 5 and 6 mg/L, respectively, while witnessing an expected decrease in power generation at more restrictive DO constraints. Exploration of simultaneous temperature and DO constraints revealed capability to address multiple water quality constraints at specified locations. Here, the reduced computational requirements of the new modeling approach demonstrated an ability to provide decision support for reservoir operations scheduling while maintaining high-fidelity hydrodynamic and water quality information as part of the optimization decision support routines.

13 HYDRO ENERGY↗

DOE OSTI2023

Using defects as a ‘fossil record’ to help interpret complex processes during additive manufacturing: as applied to raster-scanned electron beam powder bed additively manufactured Ti–6Al–4V

Abstract Defects in parts produced by additive manufacturing, instead of simply being perceived as deleterious, can act as important sources of information associated with the complex physical processes that occur during materials deposition and subsequent thermal cycles. Indeed, they act as materials-state ‘fossil’ records of the dynamic AM process. The approach of using defects as epoch-like records of prior history has been developed while studying additively manufactured Ti–6Al–4V and has given new insights into processes that may otherwise remain either obscured or unquantified. Analogous to ‘epochs,’ the evolution of these defects often is characterized by physics that span across a temporal length scale. To demonstrate this approach, a broad range of analyses including optical and electron microscopy, X-ray computed tomography, energy-dispersive spectroscopy, and electron backscatter diffraction have been used to characterize a raster-scanned electron beam Ti–6Al–4V sample. These analysis techniques provide key characteristics of defects such as their morphology, location within the part, complex compositional fields interacting with the defects, and structures on the free surfaces of defects. Observed defects have been classified as banding, spherical porosity, and lack of fusion. Banding is directly related to preferential evaporation of Al, which has an influence on mechanical properties. Lack-of-fusion defects can be used to understand columnar grain growth, fluid flow of melt pools, humping, and spattering events. Graphical abstract

36 MATERIALS SCIENCE↗

DOE OSTI2022

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Sub-100-fs energy transfer in coenzyme NADH is a coherent process assisted by a charge-transfer state

Excitation energy transfer (EET) is a key photoinduced process in biological chromophoric assemblies. Here we investigate the factors which can drive EET into efficient ultrafast sub-ps regimes. We demonstrate how a coherent transport of electronic population could facilitate this in water solvated NADH coenzyme and uncover the role of an intermediate dark charge-transfer state. High temporal resolution ultrafast optical spectroscopy gives a 54 ± 11 fs time constant for the EET process. Nonadiabatic quantum dynamical simulations computed through the time-evolution of multidimensional wavepackets suggest that the population transfer is mediated by photoexcited molecular vibrations due to strong coupling between the electronic states. The polar aqueous solvent environment leads to the active participation of a dark charge transfer state, accelerating the vibronically coherent EET process in favorably stacked conformers and solvent cavities. Our work demonstrates how the interplay of structural and environmental factors leads to diverse pathways for the EET process in flexible heterodimers and provides general insights relevant for coherent EET processes in stacked multichromophoric aggregates like DNA strands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Time-dependent, multimode interaction analysis of the gyroklystron amplifier

In this paper, a time-dependent multimode nonlinear analysis for the gyroklystron amplifier has been developed by extending the analysis of gyrotron oscillators by employing the self-consistent approach. The nonlinear analysis developed here has been validated by taking into account the reported experimental results for a 32.3 GHz, three cavity, second harmonic gyroklystron operating in the TE{sub 02} mode. The analysis has been used to estimate the temporal RF growth in the operating mode as well as the nearby competing modes. Device gain and bandwidth have been computed for different drive powers and frequencies. The effect of various beam parameters, such as beam voltage, beam current, and pitch factor, has also been studied. The computational results have estimated the gyroklystron saturated RF power ∼319 kW at 32.3 GHz with efficiency ∼23% and gain ∼26.3 dB with device bandwidth ∼0.027% (8 MHz) for a 70 kV, 20 A electron beam. The computed results are found to be in agreement with the experimental values within 10%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Three-dimensional kinetic modeling of streamer propagation in a nitrogen/helium gas mixture

A fully resolved kinetic model (Particle-in-Cell and Direct Simulation Monte Carlo for particle/photon collisions) of a near atmospheric pressure ionization wave is presented here. Fully resolving the required numerical spatial (sub-um) and temporal scales (tens of fs) for atmospheric pressure discharges in three dimensions is still a challenging task on modern super computers. To keep the overall problem tractable, the total number of elements are reduced by only simulating a 10 degree wedge rather than a full 360 degree geometry. The ionization wave is generated in a needle-plane configuration with a gap size of 250 um and a background of nitrogen and helium gas. A voltage of 1500 V is applied to the anode and an initial electron and ion density of 10 9 cm 3 is seeded in a region near the anode electrode tip and extending towards the cathode. As these initial electrons are swept away, photoionization and photoemission create new electrons and allow the ionization front to propagate towards the cathode. Results from the 90% N 2 , 10% He discharge indicate that photoionization has minimal impact on plasma formation processes and cathode photoemission is the dominant mechanism for new electrons. In the 90% He, 10% N 2 discharge case, however, photoionization likely has an impact as the observed locations of photoionization occur far enough away from the ionization front to allow for sufficient avalanche processes that contribute to the propagation of the ionization wave. Additionally, the electron energy distribution functions in the 90% He, 10% N 2 case indicate that there is less energy loss to the low lying molecular N 2 electronic states as well as the vibrational and rotational modes. This leads to higher electron energies and faster plasma development times of ~0.4 ns for the 90% He, 10% N 2 case, and ~1.5 ns for the 90% N 2 , 10% He case. In addition to analysis of the ionization wave results, the overall challenges associated with simulations near atmospheric pressure discharges in three dimensions are discussed, including the limitations of the 10 degree wedge that produces, at least qualitatively, minimal 3D effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Disentangling the impact of quasiparticles and two-level systems on the statistics of superconducting-qubit lifetime

Temporal fluctuations in the superconducting qubit lifetime, T 1 , present additional challenges in the pursuit of fault-tolerant quantum computing. Although the exact mechanisms remain unclear, T 1 fluctuations are generally attributed to strong coupling between the qubit and a few near-resonant two-level systems (TLSs), which can exchange energy with an ensemble of thermally fluctuating two-level fluctuators (TLFs) at low frequencies. Here, we report T 1 measurements of qubits with varying geometrical footprints and surface dielectrics as a function of temperature. By analyzing the noise spectrum of the qubit depolarization rate, Γ 1 = 1 / T 1 , we disentangle the contributions of TLSs, nonequilibrium quasiparticles (QPs), and equilibrium (thermally excited) QPs to the variance in Γ 1 . We find that the Γ 1 variance in qubits with smaller footprints is more susceptible to QP and TLS fluctuations than that in larger-footprint qubits. Furthermore, the QP-induced variances in all qubits align with the theoretical framework of QP diffusion and fluctuation. These findings offer valuable insights for future qubit design and engineering optimization.

Zhu, Shaojiang [Fermilab] (ORCID:0000000293180092)↗

DOE OSTI2016

TH-AB-209-05: Validating Hemoglobin Saturation and Dissolved Oxygen in Tumors Using Photoacoustic Computed Tomographic Spectroscopic Imaging

Purpose: Photoacoustic computed tomographic spectroscopy (PCT-S) provides intra-tumor measurements of oxygenation with high spatial resolution (0.2mm) and temporal fidelity (1–2 minutes) without the need for exogenous agents or ionizing radiation, thus providing a unique in vivo assay to measure SaO{sub 2} and investigate acute and chronic forms of hypoxia. The goal of this study is to validate in vivo SaO{sub 2} levels within tail artery of mice and the relationship between SaO{sub 2} and pO{sub 2} within subcutaneous breast tumors using PCT-S imaging, pulse oximetry and an OxyLite probe. Methods: A closed circuit phantom was fabricated to control blood oxygenation levels, where SaO{sub 2} was measured using a co-oximeter and pO{sub 2} using an Oxylite probe. Next, SaO{sub 2} levels within the tail arteries of mice (n=3) were measured using PCT-S and pulse oximetry while breathing high-to-low oxygen levels (6-cycles). Finally, PCT-S was used to measure SaO{sub 2} levels in MCF-7, MCF-7-VEGF165, and MDA-MB-231 xenograft breast tumors and compared to Oxylite pO{sub 2} levels values. Results: SaO{sub 2} and pO{sub 2} data obtained from the calibration phantom was fit to Hill’s equation: aO{sub 2} levels between 88 and 52% demonstrated a linear relationship (r2=0.96) and a 3.2% uncertainty between PCT-S values relative to pulse oximetry. Scatter plots of localized PCT-S measured SaO2 and Oxylite pO{sub 2} levels in MCF-7/MCF-7-VEGF165 and MDA-MD-231 breast tumors were fit to Hill’s equation: P50=17.2 and 20.7mmHg, and n=1.76 and 1.63. These results are consistent with sigmoidal form of Hill’s equation, where the lower P{sub 50} value is indicative of an acidic tumor microenvironment. Conclusion: The results demonstrate photoacoustic imaging can be used to measure SaO{sub 2} cycling and intra-tumor oxygenation, and provides a powerful in vivo assay to investigate the role of hypoxia in radiation, anti-angiogenic, and immunotherapies.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Route Sanitizer: Connected Vehicle Trajectory De-Identification Tool

Route Sanitizer is ORNL's connected vehicle moving object database de-identification tool and a graphical user interface to ORNL's connected vehicle de-identification algorithm. It uses the Google Chrome (soon to be Electron) platform so it will run on different computing platforms. The basic de-identification strategy is record redaction: portions of a vehicle trajectory (e.g. sequences of precise temporal spatial records) are removed. It does not alter retained records. The algorithm uses custom techniques to find areas within trajectories that may be considered private, then it suppresses those in addition to enough of the trajectory surrounding those locations to protect against "inference attacks" in a mathematically sound way. Map data is integrated into the process to make this possible.

Carter, JasonM↗

DOE OSTI2023

In Situ Investigation of Upper Airway Occlusion in Sleep Disordered Breathing Using Ultrasonic Transducer Arrays

This work presents a novel application of ultrasound for the real-time, non-invasive investigation of occlusion of the upper airway during events of obstructive sleep apnea/hypopnea syndrome. It is hypothesized that ultrasonic pulses applied to the neck during apneic events produce spectral and temporal features that can detect apnea occurrence. Theoretical models of ultrasound propagation and an in vitro test were conducted to test this hypothesis in both transmission and reflection modes. Complete specifications and technical details of the system design and fabrication, which is mounted on each subject’s neck, are presented, including the methodology. Nine patients (seven male and two female, mean age of 42 years, with a range of 25 to 56 years, and body mass index 37.6 ± 6.6 kg/m2) were recruited for a full night study, which included simultaneous nocturnal polysomnography for the validation of the results. Nine temporal features and four spectral features were extracted from the envelope of the received pulse waveform. These were used to compute 26 metrics to quantify the changes in the ultrasonic waveforms between normal breathing and apneic events. The statistical analysis of the collected ultrasonic data showed that at least two or more of the proposed features could detect apneic events in all subjects. The findings establish the feasibility of the proposed method as a cost-effective and non-invasive OSAHS screening tool.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Spatiotemporal Downscaling Model for Solar Irradiance Forecast Using Nearest-Neighbor Random Forest and Gaussian Process

Accurate solar photovoltaic (PV) capacity estimation requires high-resolution, site-specific solar irradiance data to account for localized variability. However, global datasets, such as the National Solar Radiation Database (NSRDB), provide regional averages that fail to capture the fine-scale fluctuations critical for large-scale grid integration. This limitation is particularly relevant in the context of increasing distributed energy resources (DERs) penetration, such as rooftop PV. Additionally, it is critical to the implementation of the U.S. Federal Energy Regulatory Commission (FERC) Order 2222, which facilitates DER participation in U.S. bulk power markets. To address this challenge, this study evaluates Nearest-Neighbor Random Forest (NNRF) and Nearest-Neighbor Gaussian Process (NNGP) models for spatiotemporal downscaling of global solar irradiance data. By leveraging historical irradiance and meteorological data, these models incorporate spatial, temporal, and feature-based correlations to enhance local irradiance predictions. The NNRF model, a machine-learning approach, prioritizes computational efficiency and predictive accuracy, while the NNGP model offers a level of interpretability and prediction uncertainty by numerically quantifying correlations and dependencies in the data. Model validation was conducted using day-ahead predictions. The results showed that the average Goodness of Fit (GoF) of the NNRF model of 90.61% across all eight sites outperformed the GoF of the NNGP of 85.88%. Additionally, the computational speed of NNRF was 2.5 times faster than the NNGP. Finally, the NNGP displayed polynomial scaling while the NNRF scaled linearly with increasing number of nearest neighbors. Additional validation of the model on five sites in Puerto Rico further confirmed the superiority of the NNRF model over the NNGP model. These findings highlight the robustness and computational efficiency of NNRF for large-scale solar irradiance downscaling, making it a strong candidate for improving PV capacity estimation and real-time electricity market integration for DERs.

Asiedu, Shadrack (ORCID:0009000646004826)↗

DOE OSTI2017

Route Sanitizer: Connected Vehicle Trajectory De-Identification Tool

Carter, Jason M↗

DOE OSTI2015

DIATOM (Data Initialization and Modification) Library Version 7.0

DIATOM is a library that provides numerical simulation software with a computational geometry front end that can be used to build up complex problem geometries from collections of simpler shapes. The library provides a parser which allows for application-independent geometry descriptions to be embedded in simulation software input decks. Descriptions take the form of collections of primitive shapes and/or CAD input files and material properties that can be used to describe complex spatial and temporal distributions of numerical quantities (often called “database variables” or “fields”) to help define starting conditions for numerical simulations. The capability is designed to be general purpose, robust and computationally efficient. By using a combination of computational geometry and recursive divide-and-conquer approximation techniques, a wide range of primitive shapes are supported to arbitrary degrees of fidelity, controllable through user input and limited only by machine resources. Through the use of call-back functions, numerical simulation software can request the value of a field at any time or location in the problem domain. Typically, this is used only for defining initial conditions, but the capability is not limited to just that use. The most recent version of DIATOM provides the ability to import the solution field from one numerical solution as input for another.

Crawford, David A.↗

DOE OSTI2016

SU-D-BRA-04: Computerized Framework for Marker-Less Localization of Anatomical Feature Points in Range Images Based On Differential Geometry Features for Image-Guided Radiation Therapy

Purpose: To propose a computerized framework for localization of anatomical feature points on the patient surface in infrared-ray based range images by using differential geometry (curvature) features. Methods: The general concept was to reconstruct the patient surface by using a mathematical modeling technique for the computation of differential geometry features that characterize the local shapes of the patient surfaces. A region of interest (ROI) was firstly extracted based on a template matching technique applied on amplitude (grayscale) images. The extracted ROI was preprocessed for reducing temporal and spatial noises by using Kalman and bilateral filters, respectively. Next, a smooth patient surface was reconstructed by using a non-uniform rational basis spline (NURBS) model. Finally, differential geometry features, i.e. the shape index and curvedness features were computed for localizing the anatomical feature points. The proposed framework was trained for optimizing shape index and curvedness thresholds and tested on range images of an anthropomorphic head phantom. The range images were acquired by an infrared ray-based time-of-flight (TOF) camera. The localization accuracy was evaluated by measuring the mean of minimum Euclidean distances (MMED) between reference (ground truth) points and the feature points localized by the proposed framework. The evaluation was performed for points localized on convex regions (e.g. apex of nose) and concave regions (e.g. nasofacial sulcus). Results: The proposed framework has localized anatomical feature points on convex and concave anatomical landmarks with MMEDs of 1.91±0.50 mm and 3.70±0.92 mm, respectively. A statistically significant difference was obtained between the feature points on the convex and concave regions (P<0.001). Conclusion: Our study has shown the feasibility of differential geometry features for localization of anatomical feature points on the patient surface in range images. The proposed framework might be useful for tasks involving feature-based image registration in range-image guided radiation therapy.

60 APPLIED LIFE SCIENCES↗