Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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 523 records · Page 29

Multiscale Modeling of Hematologic Disorders

Parasitic infectious diseases and other hereditary hematologic disorders are often associated with major changes in the shape and viscoelastic properties of red blood cells (RBCs). Such changes can disrupt blood flow and even brain perfusion, as in the case of cerebral malaria. Modeling of these hematologic disorders requires a seamless multiscale approach, where blood cells and blood flow in the entire arterial tree are represented accurately using physiologically consistent parameters. In this chapter, we present a computational methodology based on dissipative particle dynamics (DPD) which models RBCs as well as whole blood 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 small arteries and can also be used to model RBCs down to spectrin level. To this end, we present two complementary mathematical models for RBCs and describe a systematic procedure on extracting the relevant input parameters from optical tweezers and microfluidic experiments for single RBCs. We then use these validated RBC models to predict the behavior of whole healthy blood and compare with experimental results. The same procedure is applied to modeling malaria, and results for infected single RBCs and whole blood are presented.

60 APPLIED LIFE SCIENCES↗

Use of C-Arm Cone Beam CT During Hepatic Radioembolization: Protocol Optimization for Extrahepatic Shunting and Parenchymal Enhancement

PurposeTo optimize a C-arm computed tomography (CT) protocol for radioembolization (RE), specifically for extrahepatic shunting and parenchymal enhancement.Materials and MethodsA prospective development study was performed per IDEAL recommendations. A literature-based protocol was applied in patients with unresectable and chemorefractory liver malignancies undergoing an angiography before radioembolization. Contrast and scan settings were adjusted stepwise and repeatedly reviewed in a consensus meeting. Afterwards, two independent raters analyzed all scans. A third rater evaluated the SPECT/CT scans as a reference standard for extrahepatic shunting and lack of target segment perfusion.ResultsFifty scans were obtained in 29 procedures. The first protocol, using a 6 s delay and 10 s scan, showed insufficient parenchymal enhancement. In the second protocol, the delay was determined by timing parenchymal enhancement on DSA power injection (median 8 s, range 4–10 s): enhancement improved, but breathing artifacts increased (from 0 to 27 %). Since the third protocol with a 5 s scan decremented subjective image quality, the second protocol was deemed optimal. Median CNR (range) was 1.7 (0.6–3.2), 2.2 (−1.4–4.0), and 2.1 (−0.3–3.0) for protocol 1, 2, and 3 (p = 0.80). Delineation of perfused segments was possible in 57, 73, and 44 % of scans (p = 0.13). In all C-arm CTs combined, the negative predictive value was 95 % for extrahepatic shunting and 83 % for lack of target segment perfusion.ConclusionAn optimized C-arm CT protocol was developed that can be used to detect extrahepatic shunts and non-perfusion of target segments during RE.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Quality control for community-based sea-ice model development

A new collaborative organization for sea-ice model development, the CICE Consortium, has devised quality control procedures to maintain the integrity of its numerical codes' physical representations, enabling broad participation from the scientific community in the Consortium's open software development environment. Using output from five coupled and uncoupled configurations of the Los Alamos Sea Ice Model, CICE, we formulate quality control methods that exploit common statistical properties of sea-ice thickness, and test for significant changes in model results in a computationally efficient manner. New additions and changes to CICE are graded into four categories, ranging from bit-for-bit amendments to significant, answer-changing upgrades. These modifications are assessed using criteria that account for the high level of autocorrelation in sea-ice time series, along with a quadratic skill metric that searches for hemispheric changes in model answers across an array of different CICE configurations. These metrics also provide objective guidance for assessing new physical representations and code functionality.

58 GEOSCIENCES↗

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable implementation of polynomial filtering for density functional theory calculation in PARSEC

In this work, we present an efficient implementation of polynomial filtering methods in PARSEC, a real-space pseudopotential based Kohn–Sham density functional theory solver. The implementation described here improves upon a Chebyshev-filtered subspace iteration algorithm used in the previous version of PARSEC. We present a hybrid polynomial filtering scheme that combines Chebyshev-filtered subspace iteration and a spectrum slicing method that partitions the spectrum into several spectral slices and uses bandpass-filtered subspace iteration to compute approximate eigenpairs within each interior slice simultaneously. We describe a procedure to partition a spectrum and construct polynomial filters. We also discuss a number of practical issues such as the use of appropriate data layouts for carrying out the computation on a two-dimensional process grid and how to achieve good load balance by allocating an appropriate number of process groups to each spectral slice. Numerical examples are presented to demonstrate the effectiveness of the hybrid polynomial filtering method as well as the superior parallel scalability of spectrum slicing in comparison to that of Chebyshev-filtered subspace iteration.

97 MATHEMATICS AND COMPUTING↗

CoRE MOF DB: A curated experimental metal-organic framework database with machine-learned properties for integrated material-process screening

Here, we present an updated version of the Computation-Ready, Experimental (CoRE) Metal-Organic Framework (MOF) database, which includes a curated set of computation-ready MOF crystal structures designed for high-throughput computational materials discovery. Data collection and curation procedures were improved from the previous version to enable more frequent updates in the future. Machine-learning-predicted properties, such as stability metrics and heat capacities, are included in the dataset to streamline screening activities. An updated version of MOFid was developed to provide detailed information on metal nodes, organic linkers, and topologies of an MOF structure. DDEC6 partial atomic charges of MOFs were assigned based on a machine-learning model. Gibbs ensemble Monte Carlo simulations were used to classify the hydrophobicity of MOFs. The finalized dataset was subsequently used to perform integrated material-process screening for various carbon-capture conditions using high-fidelity temperature-swing adsorption (TSA) simulations. Our workflow identified multiple MOF candidates that are predicted to outperform CALF-20 for these applications.

CoRE MOF database↗

Radiation forces and torque on a rigid elliptical cylinder in acoustical plane progressive and (quasi)standing waves with arbitrary incidence

This paper presents two key contributions; the first concerns the development of analytical expressions for the axial and transverse acoustic radiation forces exerted on a 2D rigid elliptical cylinder placed in the field of plane progressive, quasi-standing, or standing waves with arbitrary incidence. The second emphasis is on the acoustic radiation torque per length. The rigid elliptical cylinder case is important to be considered as a first-order approximation of the behavior of a cylindrical fluid column trapped in air because of the significant acoustic impedance mismatch at the particle boundary. Based on the rigorous partial-wave series expansion method in cylindrical coordinates, non-dimensional acoustic radiation force and torque functions are derived and defined in terms of the scattering coefficients of the elliptic cylinder. A coupled system of linear equations is obtained after applying the Neumann boundary condition for an immovable surface in a non-viscous fluid and solved numerically by matrix inversion after performing a single numerical integration procedure. Computational results for the non-dimensional force components and torque, showing the transition from the progressive to the (equi-amplitude) standing wave behavior, are performed with particular emphasis on the aspect ratio a/b, where a and b are the semi-axes of the ellipse, the dimensionless size parameter, as well as the angle of incidence ranging from end-on to broadside incidence. The results show that the elliptical geometry has a direct influence on the radiation force and torque, so that the standard theory for circular cylinders (at normal incidence) leads to significant miscalculations when the cylinder cross section becomes non-circular. Moreover, the elliptical cylinder experiences, in addition to the acoustic radiation force, a radiation torque that vanishes for the circular cylinder case. The application of the formalism presented here may be extended to other 2D surfaces of arbitrary shape, such as Chebyshev cylindrical particles with a small deformation, stadiums (with oval shape), or other non-circular geometries.

42 ENGINEERING↗

SU-D-BRF-07: Ultrasound and Fluoroscopy Based Intraoperative Image-Guidance System for Dynamic Dosimetry in Prostate Brachytherapy

Purpose: Prostate brachytherapy is a common treatment method for low-risk prostate cancer patients. Intraoperative treatment planning is known to improve the treatment procedure and the outcome. The current limitation of intraoperative treatment planning is the inability to localize the seeds in relation to the prostate. We developed an image-guidance system to fulfill this need to achieve intraoperative dynamic dosimetry in prostate brachytherapy. Methods: Our system is based on standard imaging equipments available in the operating room, including the transrectal ultrasound (TRUS) and the mobile C-arm. A simple fiducial is added to compute the C-arm pose. Three fluoroscopic images and an ultrasound volume of the seeds and the prostate are acquired and processed by four image processing algorithms: seed segmentation, fiducial detection with pose estimation, seed reconstruction, and seeds-to-TRUS registration. The updated seed positions allow the physician to assess the quality of implantation and dynamically adjust the treatment plan during the course of surgery to achieve improved exit dosimetry. Results: The system was tested on 10 phantoms and 37 patients. Seed segmentation resulted in a 1% false negative and 2% false positive rates. Fiducial detection with pose estimation resulted in a detection rate of 98%. Seed reconstruction had a mean reconstruction error of 0.4 mm. Seeds-to-TRUS registration had a mean registration error of 1.3 mm. The total processing time from image acquisition to registration was approximately 1 minute. Conclusion: We present an image-guidance system for intraoperative dynamic dosimetry in prostate brachytherapy. Using standard imaging equipments and a simple fiducial, our system can be easily adopted in any clinics. Robust image processing algorithms enable accurate and fast computation of the delivered dose. Especially, the system enables detection of possible hot/cold spots during the surgery, allowing the physician to address these in the operating room rather than requiring additional therapy afterwards. This research was supported in part by the National Institutes of Health/National Cancer Institute (NIH/NCI) under grant 2R44CA099374 and grant 1R01CA151395, and in part by the Department of Defense (DoD) under grant W81XWH-05-1-0407.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Transmission Constraint Filtering in Large-Scale Security-Constrained Unit Commitment

When solving the security-constrained unit commitment problem, one of the most complicating factors is handling the large number of transmission constraints, corresponding to both base case and N - 1 line contingency scenarios. Although it is well known that only a few of these constraints need to be enforced, identifying this critical subset of constraints efficiently remains a challenge. In this letter, we propose a novel and simple iterative contingency-screening procedure that is able to eliminate 99.4% of the constraints selected by existing iterative methods, allowing for the solution of much larger-scale problems. We report computational results in realistic instances with up to 6468 buses and 9000 transmission lines. Furthermore, the method was also independently implemented and evaluated at MISO, where it performed faster than alternative methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An Iterative Response-Surface-Based Approach for Chance-Constrained AC Optimal Power Flow Considering Dependent Uncertainty

A modern power system is characterized by a stochastic variation of the loads and an increasing penetration of renewable energy generation, which results in large uncertainties in its states. These uncertainties bring formidable challenges to the power system planning and operation process. To address these challenges, we propose a cost-effective, iterative response-surface-based approach for the chance-constrained AC optimal power-flow problem that aims to ensure the secure operation of the power systems considering dependent uncertainties. Starting from a stochastic-sampling-based framework, we first utilize the copula theory to simulate the dependence among multivariate uncertain inputs. Then, to reduce the prohibitive computational time required in the traditional Monte-Carlo method, we propose, instead of using the original complicated power-system model, to rely on a polynomial-chaos-based response surface. This response surface allows us to efficiently evaluate the time-consuming power-system model at arbitrary distributed sampled values with a negligible computational cost. This further enables us to efficiently conduct an online stochastic testing for the system states that not only screens out the statistical active constraints, but also assists in a better design of the tightened bounds without using any Gaussian or symmetric assumption. Finally, an iterative procedure is executed to fine-tune the optimal solution that better satisfies a predefined probability. The simulations conducted in multiple test systems demonstrate the excellent performance of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Two-link Staggered Quark Smearing in QUDA

Gauge covariant smearing based on the 3D lattice Laplacian can be used to create extended operators that have better overlap with hadronic ground states. For staggered quarks, we make use of two-link parallel transport to preserve taste properties. We have implemented the procedure in QUDA. We present the performance of this code on the NVIDIA A100 GPUs in Indiana University's Big Red 200 supercomputer and on the AMD MI250X GPUs in Oak Ridge Leadership Computer Facility's (OLCF's) Crusher and discuss its scalability. We also study the performance improvement from using NVSHMEM on OLCF's Summit. Reusing precomputed two-link products for all sources and sinks, it reduces the total smearing time for a baryon correlator measurement by a factor of 100-120 as compared with the original MILC code and reduces the overall time by 60-70%.

Gottlieb, Steven↗

Neutron source in the MCNPX shielding calculating for electron accelerator driven facility

Argonne National Laboratory (ANL) of USA and Kharkov Inst. of Physics and Technology (KIPT) of Ukraine have been collaborating on the design development of an experimental neutron source facility. It is an accelerator driven system (ADS) utilizing a subcritical assembly driven by electron accelerator. The facility will be utilized for performing basic and applied nuclear researches, producing medical isotopes, and training young nuclear specialists. Monte Carlo code MCNPX has been utilized as a design tool due to its capability to transport electrons, photons, and neutrons at high energies. However the facility shielding calculations with MCNPX need enormous computational resources and the small neutron yield per electron makes sampling difficulty for the Monte Carlo calculations. A method, based on generating and utilizing neutron source file, was proposed and tested. This method reduces significantly the required computer resources and improves the statistics of the calculated neutron dose outside the shield boundary. However the statistical errors introduced by generating the neutron source were not directly represented in the results, questioning the validity of this methodology, because an insufficiently sampled neutron source can cause error on the calculated neutron dose. This paper presents a procedure for the validation of the generated neutron source file. The impact of neutron source statistic on the neutron dose is examined by calculating the neutron dose as a function of the number of electron particles used for generating the neutron source files. When the value of the calculated neutron dose converges, it means the neutron source has scored sufficient records and statistic does not have apparent impact on the calculated neutron dose. In this way, the validity of neutron source and the shield analyses could be verified. (authors)

43 PARTICLE ACCELERATORS↗

Test-Based Calibration of Safety Factors for Capacity Models

In this paper, a simple procedure to calibrate the safety factor of capacity models is presented. The calibration can be carried out based on any available database of experimental tests, even of limited size. The procedure aims to assess the model capability of predicting the test results and to calibrate the safety factor so that the capacity equation meets the target reliability level required by the code or sought by the calibrator. After predicting each test of the database with the capacity equation under consideration, the test-prediction pairs are checked for the property of linearity, and the relative error for the properties of homoscedasticity and normality. Once these properties are fulfilled—which may require a nonlinear transformation of the test values and/or the predictions—the closed-form equation proposed in this paper is employed to compute a target design value. The model safety factor is finally obtained by comparing such target design value with the design value obtained from the code. The paper also proposes two approximate analytical equations to compute the tolerance factor, used to attain any given fractile, as a function of the (even small) number of tests, with any assigned confidence level. A fundamental outcome of the procedure is that it yields an objective indicator of the model accuracy, measured by the standard deviation of its error, which may be regarded as a parameter useful for selecting the most reliable model among different competing ones. In the long run, the application of the proposed procedure will allow achieving a uniform reliability level throughout all capacity models used in codes and guidelines. A further advantage is that the partial safety factors so derived can be straightforwardly updated when more experiments become available. Finally, as an example, the proposed procedure is herein applied to the ACI 318 shear design capacity equation for concrete members unreinforced in shear.

42 ENGINEERING↗

A review of privacy in energy applications

As the distribution system continues to experience an increase in distributed energy resource (DER) and electric vehicle (EV) penetration, so does the need for new solutions that can help grid operators manage and leverage their capabilities. This will undoubtedly lead to new operational schemes and business opportunities that will transform the traditional consumer into a prosumer who will be more actively engaged in grid operations. Although the field is still under active development, many of the potential use cases presented in literature or industry are built upon edge computing, two-way communications, and other innovative computational constructs to attain their goals. However, at their core, many use cases assume a great level of data access to aid with the decision-making process, an assumption that may need to be revised to ensure fair and equitable operational processes are maintained. This may be particularly true as edge resources are predicted to participate in retail-side, many-to-many, or peer-to-peer markets and thus may lead to financial impacts if data access considerations are ignored. The need to revise data access mechanisms can be further justified by the introduction of new participants into the operational process, who do not have the same level of trust, nor the incentives to focus on energy delivery as their primary objective. At the same time, more consumers are becoming aware of their own data, and the potential impacts of its abuse. To help solution developers better understand these risks, this report has been developed to offer an initial introduction to the topic of privacy. This is achieved by 1) Highlighting the need for privacy-aware solutions; 2) Encouraging system designers to be inquisitive about the status quo; 3) Documenting the existing threat space; 4) Presenting and evaluating tools that may be helpful towards enabling better privacy postures; and 5) Making recommendations to encourage the adoption of better practices. From a technical perspective, the report focuses on evaluating two potential techniques by applying them to the Transactive Energy Space. Based on the obtained results, it can be established that differential privacy (DP) methods may have limited applicability when highly correlated, time-series data records need to be protected. However, DP may be a powerful tool when it is used to aggregate and analyze mid-size and large-size data sets in a more traditional statistical environment. The second tool under evaluation is threshold cryptography, which can guarantee complete secrecy (and thus privacy) but requires the establishment of key management procedures and dedicated communication channels for key coordination. Therefore, due to its increased computational overhead, the use of threshold cryptography must be weighted using a cost/benefit analysis on a per-application basis.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Graph Theory Approach to High-Throughput Surface Adsorption Structure Generation

We present a methodology for graph based enumeration of surfaces and unique chemical adsorption structures bonded to those surfaces. Utilizing the graph produced from a bulk structure, we create a unique graph representation for any general slab cleave and further extend that representation to include a large variety of catalytically relevant adsorbed molecules. Here, we also demonstrate simple geometric procedures to generate 3D initial guesses of these enumerated structures. While generally useful for generating a wide variety of structures used in computational surface science and heterogeneous catalysis, these techniques are also key to facilitating an informatics approach to the high-throughput search for more effective catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Framework for Integrating Distributed Energy Resources into Distribution Systems

This paper focuses on developing a novel theoretical framework for effective coordination and control of a large number of distributed energy resources in distribution systems in order to more reliably manage the future U.S. electric power grid under the high penetration of renewable generation. The proposed framework provides a systematic view of the overall structure of the future distribution systems along with the underlying information flow, functional organization, and operational procedures. It is characterized by the features of being open, flexible and interoperable with the potential to support dynamic system configuration. Under the proposed framework, the energy consumption of various DERs is coordinated and controlled in a hierarchical way by using market-based approaches. The real-time voltage control is simultaneously considered to complement the real power control in order to keep nodal voltages stable within acceptable ranges during real time. In addition, computational challenges associated with the proposed framework are also discussed with recommended practices.

Lian, Jianming↗

Real-Time and Adaptive Reservoir Computing With Application to Profile Prediction in Fusion Plasma

Nuclear fusion is a promising alternative to address the problem of sustainable energy production. The tokamak is an approach to fusion based on magnetic plasma confinement, constituting a complex physical system with many control challenges. Here we study the characteristics and optimization of reservoir computing (RC) for real-time and adaptive prediction of plasma profiles in the DIII-D tokamak. Our experiments demonstrate that RC achieves comparable results to state-of-the-art (deep) convolutional neural networks (CNNs) and long short-term memory (LSTM) models, with a significantly easier and faster training procedure. This efficient approach allows for fast and frequent adaptation of the model to new situations, such as changing plasma conditions or different fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hyperspectral Image Classification Using Graph Clustering Methods

Hyperspectral imagery is a challenging modality due to the dimension of the pixels which can range from hundreds to over a thousand frequencies depending on the sensor. Most methods in the literature reduce the dimension of the data using a method such as principal component analysis, however this procedure can lose information. More recently methods have been developed to address classification of large datasets in high dimensions. This paper presents two classes of graph-based classification methods for hyperspectral imagery. Using the full dimensionality of the data, we consider a similarity graph based on pairwise comparisons of pixels. The graph is segmented using a pseudospectral algorithm for graph clustering that requires information about the eigenfunctions of the graph Laplacian but does not require computation of the full graph. We develop a parallel version of the Nyström extension method to randomly sample the graph to construct a low rank approximation of the graph Laplacian. With at most a few hundred eigenfunctions, we can implement the clustering method designed to solve a variational problem for a graph-cut-based semi-supervised or unsupervised classification problem. We implement OpenMP directive-based parallelism in our algorithms and show performance improvement and strong, almost ideal, scaling behavior. The method can handle very large datasets including a video sequence with over a million pixels, and the problem of segmenting a data set into a pre-determined number of classes.

47 OTHER INSTRUMENTATION↗