Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational frameworks”

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 163 records · Page 9

An Efficient Numerical Algorithm for Solving Data Driven Feedback Control Problems

The goal of this paper is to solve a class of stochastic optimal control problems numerically, in which the state process is governed by an Itô type stochastic differential equation with control process entering both in the drift and the diffusion, and is observed partially. The optimal control of feedback form is determined based on the available observational data. In this work, we call this type of control problems the data driven feedback control. The computational framework that we introduce to solve such type of problems aims to find the best estimate for the optimal control as a conditional expectation given the observational information. To make our method feasible in providing timely feedback to the controlled system from data, we develop an efficient stochastic optimization algorithm to implement our computational framework.

97 MATHEMATICS AND COMPUTING↗

Computational Biorheology of Human Blood Flow in Health and Disease

Hematologic disorders arising from infectious diseases, hereditary factors and environmental influences can lead to, and can be influenced by, significant changes in the shape, mechanical and physical properties of red blood cells (RBCs), and the biorheology of blood flow. Hence, modeling of hematologic disorders should take into account the multiphase nature of blood flow, especially in arterioles and capillaries. We present here an overview of a general computational framework based on dissipative particle dynamics (DPD) which has broad applicability in cell biophysics with implications for diagnostics, therapeutics and drug efficacy assessments for a wide variety of human diseases. This computational approach, validated by independent experimental results, is capable of modeling the biorheology of whole blood and its individual components during blood flow so as to investigate cell mechanistic processes in health and disease. DPD is a Lagrangian method that can be derived from systematic coarse-graining of molecular dynamics but can scale efficiently up to arterioles and can also be used to model RBCs down to the spectrin level. We start from experimental measurements of a single RBC to extract the relevant biophysical parameters, using single-cell measurements involving such methods as optical tweezers, atomic force microscopy and micropipette aspiration, and cell-population experiments involving microfluidic devices. We then use these validated RBC models to predict the biorheological behavior of whole blood in healthy or pathological states, and compare the simulations with experimental results involving apparent viscosity and other relevant parameters. Finally, while the approach discussed here is sufficiently general to address a broad spectrum of hematologic disorders including certain types of cancer, this paper specifically deals with results obtained using this computational framework for blood flow in malaria and sickle cell anemia.

60 APPLIED LIFE SCIENCES↗

High-throughput identification of off-targets for the mechanistic study of severe adverse drug reactions induced by analgesics

Drugs may induce adverse drug reactions (ADRs) when they unexpectedly bind to proteins other than their therapeutic targets. Identification of these undesired protein binding partners, called off-targets, can facilitate toxicity assessment in the early stages of drug development. In this study, a computational framework was introduced for the exploration of idiosyncratic mechanisms underlying analgesic-induced severe adverse drug reactions (SADRs). The putative analgesic-target interactions were predicted by performing reverse docking of analgesics or their active metabolites against human/mammal protein structures in a high-throughput manner. Subsequently, bioinformatics analyses were undertaken to identify ADR-associated proteins (ADRAPs) and pathways. Using the pathways and ADRAPs that this analysis identified, the mechanisms of SADRs such as cardiac disorders were explored. For instance, 53 putative ADRAPs and 24 pathways were linked with cardiac disorders, of which 10 ADRAPs were confirmed by previous experiments. Moreover, it was inferred that pathways such as base excision repair, glycolysis/glyconeogenesis, ErbB signaling, calcium signaling, and phosphatidyl inositol signaling likely play pivotal roles in drug-induced cardiac disorders. In conclusion, our framework offers an opportunity to globally understand SADRs at the molecular level, which has been difficult to realize through experiments. It also provides some valuable clues for drug repurposing. - Highlights: • A novel computational framework was developed for mechanistic study of SADRs. • Off-targets of drugs were identified in large scale and in a high-throughput manner. • SADRs like cardiac disorders were systematically explored in molecular networks. • A number of ADR-associated proteins were identified.

60 APPLIED LIFE SCIENCES↗

A comparative study of shear band tracking strategies in three-dimensional finite elements with embedded weak discontinuities

We present a computational framework for the treatment of shear localization in metallic materials under dynamic loading, based on the integration of a shear band tracking strategy into an explicit 3D finite element formulation with embedded weak discontinuities. Within this computational framework, an embedded shear band’s mid-surface is represented by an iso-surface of a level set function, which is obtained by solving a heatconduction type boundary value problem (BVP). The solution of this BVP is carried out either globally over the entire problem domain, or locally at the level of individual elements in the mesh. In this paper, we present a detailed comparison of these global and local algorithmic implementations of the shear band tracking strategy. Numerical results obtained using these two implementations, and using a simplified formulation without shear band tracking, are presented and compared. Moreover, we compare the computational efficiency and parallel scaling performance of the local and global implementations. This comparative study shows that both implementations can simulate the initiation of two independent shear bands and their propagation past each other without merging, whereas only the global implementation can successfully simulate the merging of two branches of a single shear band. This study also confirms that the global implementation is more computationally intensive, since it requires a global system of linear equations associated with the level set BVP to be solved at each time step. Both implementations exhibit very good scalability in domain decomposition-based parallel simulations.

36 MATERIALS SCIENCE↗

Thermo-chemical Fluid Flow Simulation in Hot-Dip Galvanizing: The Evaluation of Dross Build-Up Formation

This paper presents a new computational framework to investigate the driving force for the formation of intermetallic dross particles in the zinc bath and dross build-up on the bath hardware. This is a major problem in continuous hot-dip galvanizing lines. The Computational Fluid Dynamics (CFD) model calculates the turbulent thermo-chemical flow conditions within the liquid melt. A detailed modeling of the steel strip–liquid interface enhances this approach, by means of a conjugated heat transfer calculation and a spatially resolved, temperature- and concentration-dependent iron dissolution and aluminum uptake. The heart of the computational framework is a thermodynamic model, which assesses the driving force for the formation or dissolution of dross particles in the bath and dross build-up on stationary and rotating equipment. The simulation results are validated with temperature and species depth profile measurements. The applicability of the CFD model is shown by investigating the locally resolved aluminum uptake and iron dissolution at the steel-strip surface, and the multi-physics conditions in the region near the roll. The novel approach of evaluating the thermodynamic driving force enables the assessment of the formation of dross build-up at the roll surface.

36 MATERIALS SCIENCE↗

Real-time x-ray fluoroscopy-based catheter detection and tracking for cardiac electrophysiology interventions

Purpose: X-ray fluoroscopically guided cardiac electrophysiology (EP) procedures are commonly carried out to treat patients with arrhythmias. X-ray images have poor soft tissue contrast and, for this reason, overlay of a three-dimensional (3D) roadmap derived from preprocedural volumetric images can be used to add anatomical information. It is useful to know the position of the catheter electrodes relative to the cardiac anatomy, for example, to record ablation therapy locations during atrial fibrillation therapy. Also, the electrode positions of the coronary sinus (CS) catheter or lasso catheter can be used for road map motion correction.Methods: In this paper, the authors present a novel unified computational framework for image-based catheter detection and tracking without any user interaction. The proposed framework includes fast blob detection, shape-constrained searching and model-based detection. In addition, catheter tracking methods were designed based on the customized catheter models input from the detection method. Three real-time detection and tracking methods are derived from the computational framework to detect or track the three most common types of catheters in EP procedures: the ablation catheter, the CS catheter, and the lasso catheter. Since the proposed methods use the same blob detection method to extract key information from x-ray images, the ablation, CS, and lasso catheters can be detected and tracked simultaneously in real-time.Results: The catheter detection methods were tested on 105 different clinical fluoroscopy sequences taken from 31 clinical procedures. Two-dimensional (2D) detection errors of 0.50 {+-} 0.29, 0.92 {+-} 0.61, and 0.63 {+-} 0.45 mm as well as success rates of 99.4%, 97.2%, and 88.9% were achieved for the CS catheter, ablation catheter, and lasso catheter, respectively. With the tracking method, accuracies were increased to 0.45 {+-} 0.28, 0.64 {+-} 0.37, and 0.53 {+-} 0.38 mm and success rates increased to 100%, 99.2%, and 96.5% for the CS, ablation, and lasso catheters, respectively. Subjective clinical evaluation by three experienced electrophysiologists showed that the detection and tracking results were clinically acceptable.Conclusions: The proposed detection and tracking methods are automatic and can detect and track CS, ablation, and lasso catheters simultaneously and in real-time. The accuracy of the proposed methods is sub-mm and the methods are robust toward low-dose x-ray fluoroscopic images, which are mainly used during EP procedures to maintain low radiation dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

IDAC-Dose 2.1, an internal dosimetry program for diagnostic nuclear medicine based on the ICRP adult reference voxel phantoms

To date, the estimated radiation-absorbed dose to organs and tissues in patients undergoing diagnostic examinations in nuclear medicine is derived via calculations based on models of the human body and the biokinetic behaviour of the radiopharmaceutical. An internal dosimetry computer program, IDAC-Dose2.1, was developed based on the International Commission on Radiological Protection (ICRP)-specific absorbed fractions and computational framework of internal dose assessment given for reference adults in ICRP Publication 133. The program uses the radionuclide decay database of ICRP Publication 107 and considers 83 different source regions irradiating 47 target tissues, defining the effective dose as presented in ICRP Publications 60 and 103. The computer program was validated against another ICRP dosimetry program, Dose and Risk Calculation (DCAL), that employs the same computational framework in evaluation of occupational and environmental intakes of radionuclides. IDAC-Dose2.1 has a sub-module for absorbed dose calculations in spherical structures of different volumes and composition; this sub-module is intended for absorbed dose estimates in radiopharmaceutical therapy. For nine specific alpha emitters, the absorbed dose contribution from their decay products is also included in the committed absorbed dose calculations. The absorbed doses and effective dose of 131 I-iodide determined by IDAC-Dose2.1 were validated against the dosimetry program DCAL, showing identical results. IDAC-Dose2.1 was used to calculate absorbed doses for intravenously administered 18 F-FDG and orally administered 99m Tc-pertechnetate and 131 I-iodide, three frequently used radiopharmaceuticals. Using the tissue weighting factors from ICRP Publication 103, the effective dose per administered activity was estimated to be 0.016 mSv/MBq for 18 F-FDG, 0.014 mSv/MBq for 99m Tc-pertechnetate, and 16 mSv/MBq for 131 I-iodide. The internal dosimetry program IDAC-Dose2.1 was developed and applied to three radiopharmaceuticals for validation against DCAL and to generate improved absorbed dose estimations for diagnostic nuclear medicine using specific absorbed fraction values of the ICRP computational voxel phantoms. The sub-module for absorbed dose calculations in spherical structures 1 mm to 9 cm in diameter and different tissue composition was included to broaden the clinical usefulness of the program. The IDAC-Dose2.1 program is free software for research and available for download at http://www.idac-dose.org.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

End-to-End Differentiable Modeling and Management of the Environment

Focal Area: (2) Data acquisition and assimilation enabled by machine learning, AI, and advanced methods including experimental/network design/optimization. We emphasize the importance of leveraging optimization techniques from AI/machine learning (ML) to solve challenging problems in Earth system modeling. Science Challenge: Automatic differentiation has had a transformative effect on ML by allowing the calculation of gradients of arbitrary functions in an incredibly large class of models. We can potentially realize similar improvements in parameter estimation and control for Earth system models (ESMs) by reimplementing them in computational frameworks from ML. Practitioners working with large (>10 7 parameters) models in ML and AI can obtain good predictive performance in a range of spatiotemporal tasks by making use of optimization via stochastic gradient descent and incorporating prior knowledge at multiple levels. We propose writing ESMs in open-source computational frameworks such as Torch, Tensorflow, and JAX to greatly expand the scope of environmental forecasting and management challenges, which can be addressed by leveraging automatic differentiation and gradient descent-like algorithms. We do not call for a wholesale replacement of physical models with data-driven surrogates, but rather advocate for interleaving physical and empirical equations in a manner that is most faithful to the extent of our scientific knowledge and observational data. Central to this topic is the merging of differentiable physical simulations with differentiable optimization layers, which are now both beginning to come to the forefront.

54 ENVIRONMENTAL SCIENCES↗

Passenger Aggregation Network with Very Efficient Listing (PANVEL) Ridesharing Model for Urban Air Mobility

This paper introduces a ridesharing model useful for examining the role of ridesharing in advancing Urban Air Mobility (UAM) operations. Although the success of UAM will likely require a host of new technologies in aircraft design, airspace management, autonomy, and more, operational innovations like ridesharing may also be key, especially to lower operating costs and attract a broader market. The ridesharing algorithm introduced, termed the Passenger Aggregation Network with Very Efficient Listing (PANVEL) model, estimates trip numbers and passenger occupancy aboard particular aircraft within a UAM network in a metropolitan area. This algorithm is applied within a transportation computational framework to determine the impact of ridesharing on mode choice, passenger travel patterns, and transportation costs. The algorithm aggregates passengers traveling between identical origin-destination pairs, factoring in individual-specific metrics, such as value of time and alternative mode options, to generate high load factors on UAM aircraft. The computational framework simulates passenger mode-choice selections using an effective-cost approach to determine the number of ridesharing UAM-preferred trips, and results with ridesharing are representative of a rough upper bound on UAM ridership potential for a given metro area. A case study in the Cleveland metro area is presented to demonstrate the algorithm and assess ridesharing potential.

advanced air mobility↗

Passenger Aggregation Network with Very Efficient Listing (PANVEL) Ridesharing Model for Urban Air Mobility

This presentation introduces a ridesharing model useful for examining the role of ridesharing in advancing Urban Air Mobility (UAM) operations. Although the success of UAM will likely require a host of new technologies in aircraft design, airspace management, autonomy, and more, operational innovations like ridesharing may also be key, especially to lower operating costs and attract a broader market. The ridesharing algorithm introduced, termed the Passenger Aggregation Network with Very Efficient Listing (PANVEL) model, estimates trip numbers and passenger occupancy aboard particular aircraft within a UAM network in a metropolitan area. This algorithm is applied within a transportation computational framework to determine the impact of ridesharing on mode choice, passenger travel patterns, and transportation costs. The algorithm aggregates passengers traveling between identical origin-destination pairs, factoring in individual-specific metrics, such as value of time and alternative mode options, to generate high load factors on UAM aircraft. The computational framework simulates passenger mode-choice selections using an effective-cost approach to determine the number of ridesharing UAM-preferred trips, and results with ridesharing are representative of a rough upper bound on UAM ridership potential for a given metro area. A case study in the Cleveland metro area is presented to demonstrate the algorithm and assess ridesharing potential.

advanced air mobility↗

Test and Validate Distributed Coaxial Cable Sensors for in situ Condition Monitoring of Coal-Fired Boiler Tubes

This project aims to test, validate, and advance the technology readiness level (from TRL5 to TRL7) of a novel low-cost distributed stainless-steel/ceramic coaxial cable sensing (SSC-CCS) technology for in situ monitoring of the boiler tube temperature in existing coal-fired power plants. The novel SSC-CCS sensing technology and associated condition-based monitoring (CBM) software to be demonstrated in this project will lead to an improved understanding of the boiler tube failure mechanisms and a prognostic system to improve the overall performance, reliability, and flexibility of the nation’s coal-fired power plant fleet. A boiler tube monitoring system with distributed coaxial cable temperature sensors and a sensor acquisition system was constructed. The high-temperature coaxial cable sensor with a length of 1.3m was made by using a quartz tube (1mm inner diameter (ID) and 6mm outer diameter (OD)) to concentrically separate a 304 stainless-steel (SS) rod (1mm OD) and SS tube (7.94mm OD and 6.16mm ID). The sensor acquisition system includes a vector network analyzer (VNA), a radio frequency (RF) power amplifier, multiple switches and a USB hub. The distributed stainless-steel quartz coaxial cable sensor (SSQ-CCS) had a linear response to temperature with a resolution uncertainty of σ = 0.77℃. To withstand the harsh conditions of 3,300 steam pressures and 800℃ high temperatures, the sensor was shielded by a protective tube made of the same material as the boiler tube. The protection tube had an OD of 1.5 inches and a thickness of 0.25 inches. In the laboratory tests, the sensor showed good sensitivity and fast response. The drift was bounded between +0.33% and -0.67% during a test at 600℃ for 350 hours, indicating good stability of the sensor. A field test was conducted where four sensors were welded on four superheat tubes (SH-Ts) at a coal-fired power station over 400 days. Conventional thermocouples were welded to the superheater tubes alongside the coaxial cable sensors for the purpose of comparison. Two sensors were capable of distributed sensing, with three multiplexed sensing sections. The other two sensors were single section. During the 400-day test period, the power plant experienced startups and shutdowns. At the steady state operations, the temperature of the boiler tube is about 600℃ (1112°F). The sensors recorded the entire coal-firing processes (start-up, steady state, and shut-down) and the glitch event. A GSM modem and a Watchdog were added to the system to ensure reliable data recording. The GSM modem sent daily messages to plant managers and Clemson team to inform the status of the sensor system. If the system was not normally working, the Watchdog would reboot the system automatically. The new coaxial cable based distributed sensing technology has been proven to be successful in both laboratory and field tests. A comprehensive four-stage multi-physics computational framework has been developed to assist the design, optimization, installation, and operation of SSQ-CCS. With the consideration of various operation conditions, we predict the distributions of flue gas temperatures within coal-fired boilers, the temperature correlation between the boiler tube and SSQ-CCS, and the safety of SSQ-CCS. A conditional-based monitoring system is implemented as well. The computational framework developed in this work can guide the future operation of coal-fired plants and other power plants for the safety prediction of boiler operations.

01 COAL, LIGNITE, AND PEAT↗

Numerical Simulation of a Flap-Edge Flowfield

In this paper we develop an approximate computational framework for simulation of the fluctuating flowfield associated with the complex vortex system seen at the side edge of a flap in a multielement high-lift airfoil system. The eventual goal of these simulations is to provide an estimate of the spectral content of these fluctuations, in order that the spectrum of the noise generated by such flowfields may be estimated. Results from simulations utilizing this computational framework are shown.

Streett, C. L.↗

Numerical Simulation of Fluctuations Leading to Noise in a Flap-Edge Flowfield

In this paper we develop an approximate computational framework for simulation of the fluctuating flowfield associated with the complex vortex system seen at the side edge of a flap in a multi-element high-lift airfoil system. The eventual goal of these simulations is to provide an estimate of the spectral content of these fluctuations, in order that the spectrum of the noise generated by such flowfields may be estimated. Results from simulations utilizing this computational framework are shown.

Streett, C. L.↗

A Component-based Programming Model for Composite, Distributed Applications

The nature of scientific programming is evolving to larger, composite applications that are composed of smaller element applications. These composite applications are more frequently being targeted for distributed, heterogeneous networks of computers. They are most likely programmed by a group of developers. Software component technology and computational frameworks are being proposed and developed to meet the programming requirements of these new applications. Historically, programming systems have had a hard time being accepted by the scientific programming community. In this paper, a programming model is outlined that attempts to organize the software component concepts and fundamental programming entities into programming abstractions that will be better understood by the application developers. The programming model is designed to support computational frameworks that manage many of the tedious programming details, but also that allow sufficient programmer control to design an accurate, high-performance application.

Eidson, Thomas M.↗

Four-Stage Multi-Physics Simulations to Assist Temperature Sensor Design for Industrial-Scale Coal-Fired Boiler

The growth of renewable energy sources presents a pressing challenge to the operation and maintenance of existing fossil fuel power plants, given that fossil fuel remains the predominant fuel source, responsible for over 60% of electricity generation in the United States. One of the main concerns within these fossil fuel power plants is the unpredictable failure of boiler tubes, resulting in emergency maintenance with significant economic and societal consequences. A reliable high-temperature sensor is necessary for in situ monitoring of boiler tubes and the safety of fossil fuel power plants. In this study, a comprehensive four-stage multi-physics computational framework is developed to assist the design, optimization installation, and operation of the high-temperature stainless-steel and quartz coaxial cable sensor (SSQ-CCS) for coal-fired boiler applications. With the consideration of various operation conditions, we predict the distributions of flue gas temperatures within coal-fired boilers, the temperature correlation between the boiler tube and SSQ-CCS, and the safety of SSQ-CCS. With the simulation-guided sensor installation plan, the newly designed SSQ-CCSs have been employed for field testing for more than 430 days. The computational framework developed in this work can guide the future operation of coal-fired plants and other power plants for the safety prediction of boiler operations.

42 ENGINEERING↗

Demonstrating GPU code portability and scalability for radiative heat transfer computations

High performance computing frameworks utilizing CPUs, Nvidia GPUs, and/or Intel Xeon Phis necessitate portable and scalable solutions for application developers. Nvidia GPUs in particular present numerous portability challenges with a different programming model, additional memory hierarchies, and partitioned execution units among streaming multiprocessors. This work presents modifications to the Uintah asynchronous many-task runtime and the Kokkos portability library to enable one single codebase for complex multiphysics applications to run across different architectures. Scalability and performance results are shown on multiple architectures for a globally coupled radiation heat transfer simulation, ranging from a single node to 16,384 Titan compute nodes.

Computer Science↗

Reaction–Transport Coupling in a Nanostructured Porous Electrode

The function of nanoporous electrodes involves chemical reactions coupled to transport of reactants and products into and out of the solid matrix. To understand how this coupling can depend on electrode structure and chemistry, a general multiscale, molecular-level computational framework is introduced and applied to a model system containing a nanoporous photoanode, the dye-sensitized solar cell. The pore structure and electrolyte composition are systematically varied to probe this coupling, revealing that in most cases the photocurrent density depends on dye coverage within the anode matrix, the specific architecture of an individual pore, and electrolyte gradients in between the anode and cathode. The large electrolyte gradients that develop in the anode are not kinetically controlling. In conclusion, the computational framework demonstrated is general and can be applied to other electrochemical, photochemical, and thermally driven nanoporous systems to evaluate the fundamental characteristics of their function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗