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 1,783 records · Page 99

Ablation Volume Measurement After Percutaneous Cryoablation Using a Two-cryo-probe Technique for Small Hepatocellular Carcinomas

Purpose: To calculate the ablation volume of percutaneous cryoablation (PCA) with a two-cryo-probe technique for small hepatocellular carcinomas (HCCs) and to assess risk factors for local tumor recurrence (LTR). Materials and Methods: Between January 2013 and June 2017, 96 patients (mean age, 60.7 years; range, 37–83 years) with 106 small HCCs were retrospectively analyzed. The ablation volume, technical success, LTR rates, and complications were evaluated. Ablation volume was measured after the first freezing and immediately after completing PCA twice via computed tomography imaging. Prognostic factors associated with LTR were analyzed. Results: Technical success was achieved in all patients. The mean final ablation volume was 19.1 ± 4.8 cm{sup 3}, which was significantly higher than the first freezing ablation volume (14.7 ± 4.3 cm{sup 3}, p < 0.001). During the median follow-up period of 16.7 months (range, 5–52 months), LTR-free survival rates based on Kaplan–Meier analyses at 6 months, 1 year, and 2 years were 87.7%, 84.0%, and 80.2%, respectively. Only one major complication of post-procedural arterial bleeding occurred. A final ablation volume/tumor volume < 10 was a significant risk factor for LTR (p = 0.044). Conclusions: A ratio of final ablation volume to tumor volume < 10 was a significant prognostic factor for LTR. Therefore, the measurement and prediction of the final ablation volume are important to reduce LTR.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Implementation of a hybrid particle code with a PIC description in r–z and a gridless description in ϕ into OSIRIS

For many plasma physics problems, three-dimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry and the dominant three-dimensional physics is captured by the inclusion of only a few azimuthal harmonics. Recently, it was proposed [1] to model one such problem, laser wakefield acceleration, by expanding the fields and currents in azimuthal harmonics and truncating the expansion. The complex amplitudes of the fundamental and first harmonic for the fields were solved on an r–z grid and a procedure for calculating the complex current amplitudes for each particle based on its motion in Cartesian geometry was presented using a Marder's correction to maintain the validity of Gauss's law. In this paper, we describe an implementation of this algorithm into OSIRIS using a rigorous charge conserving current deposition method to maintain the validity of Gauss's law. We show that this algorithm is a hybrid method which uses a particles-in-cell description in r–z and a gridless description in ϕ. We include the ability to keep an arbitrary number of harmonics and higher order particle shapes. Examples for laser wakefield acceleration, plasma wakefield acceleration, and beam loading are also presented and directions for future work are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reliable and efficient solution of genome-scale models of Metabolism and macromolecular Expression

Currently, Constraint-Based Reconstruction and Analysis (COBRA) is the only methodology that permits integrated modeling of Metabolism and macromolecular Expression (ME) at genome-scale. Linear optimization computes steady-state flux solutions to ME models, but flux values are spread over many orders of magnitude. Data values also have greatly varying magnitudes. Furthermore, standard double-precision solvers may return inaccurate solutions or report that no solution exists. Exact simplex solvers based on rational arithmetic require a near-optimal warm start to be practical on large problems (current ME models have 70,000 constraints and variables and will grow larger). We also developed a quadrupleprecision version of our linear and nonlinear optimizer MINOS, and a solution procedure (DQQ) involving Double and Quad MINOS that achieves reliability and efficiency for ME models and other challenging problems tested here. DQQ will enable extensive use of large linear and nonlinear models in systems biology and other applications involving multiscale data.

59 BASIC BIOLOGICAL SCIENCES↗

Development of precision laser goniometer systems

The concept of constructing precision laser goniometer systems, based on integrating a ring laser and an optical angle sensor with the holographic principle of angular scale recording is considered. The concept implies the application of the cross-calibration procedure, aimed at determination of systematic components of the errors of angle sensors, used in the system. The results of the presented system studies demonstrate the error of angular measurements amounting to {approx}0.01''. The results of implementing the proposed concept in the creation of a standard system of the plane angle unit of rigid rotation and the measuring and computing complex for automated control of digital angle transducers with high digit capacity are briefly presented. (laser optics 2012)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-T-04: 3D Dose Based Patient Compensator QA Procedure for Proton Radiotherapy

Purpose: In proton double-scattering radiotherapy, compensators are the essential patient specific devices to contour the distal dose distribution to the tumor target. Traditional compensator QA is limited to checking the drilled surface profiles against the plan. In our work, a compensator QA process was established that assess the entire compensator including its internal structure for patient 3D dose verification. Methods: The fabricated patient compensators were CT scanned. Through mathematical image processing and geometric transformations, the CT images of the proton compensator were combined with the patient simulation CT images into a new series of CT images, in which the imaged compensator is placed at the planned location along the corresponding beam line. The new CT images were input into the Eclipse treatment planning system. The original plan was calculated to the combined CT image series without the plan compensator. The newly computed patient 3D dose from the combined patientcompensator images was verified against the original plan dose. Test plans include the compensators with defects intentionally created inside the fabricated compensators. Results: The calculated 3D dose with the combined compensator and patient CT images reflects the impact of the fabricated compensator to the patient. For the test cases in which no defects were created, the dose distributions were in agreement between our method and the corresponding original plans. For the compensator with the defects, the purposely changed material and a purposely created internal defect were successfully detected while not possible with just the traditional compensator profiles detection methods. Conclusion: We present here a 3D dose verification process to qualify the fabricated proton double-scattering compensator. Such compensator detection process assesses the patient 3D impact of the fabricated compensator surface profile as well as the compensator internal material and structure changes. This research receives funding support from CURA Medical Technologies.

60 APPLIED LIFE SCIENCES↗

Stochastic characterization of phase detection algorithms in phase-shifting interferometry

Phase-shifting interferometry (PSI) is the preferred non-contact method for profiling sub-nanometer surfaces. Based on monochromatic light interference, the method computes the surface profile from a set of interferograms collected at separate stepping positions. Errors in the estimated profile are introduced when these positions are not located correctly. In order to cope with this problem, various algorithms that minimize the effects of certain types of stepping errors (linear, sinusoidal, etc.) have been developed. Despite the relatively large number of algorithms suggested in the literature, there is no unified way of characterizing their performance when additional unaccounted random errors are present. Here, we suggest a procedure for quantifying the expected behavior of each algorithm in the presence of independent and identically distributed (i.i.d.) random stepping errors, which can occur in addition to the systematic errors for which the algorithm has been designed. As a result, the usefulness of this method derives from the fact that it can guide the selection of the best algorithm for specific measurement situations.

42 ENGINEERING↗

Hot, cold, and annual reference atmospheres for Edwards Air Force Base, California (1975 version)

Reference atmospheres pertaining to summer (hot), winter (cold), and mean annual conditions for Edwards Air Force Base, California, are presented from surface to 90 km altitude (700 km for the annual model). Computed values of pressure, kinetic temperature, virtual temperature, and density and relative differences percentage departure from the Edwards reference atmospheres, 1975 (ERA-75) of the atmospheric parameters versus altitude are tabulated in 250 m increments. Hydrostatic and gas law equations were used in conjunction with radiosonde and rocketsonde thermodynamic data in determining the vertical structure of these atmospheric models. The thermodynamic parameters were all subjected to a fifth degree least-squares curve-fit procedure, and the resulting coefficients were incorporated into Univac 1108 computer subroutines so that any quantity may be recomputed at any desired altitude using these subroutines.

Johnson, D. L.↗

An efficient algorithm for solution of the unsteady transonic small-disturbance equation

A time accurate approximate factorization (AF) algorithm is formulated for solution of the three dimensional unsteady transonic small-disturbance equation. The AF algorithm consists of a time linearization procedure coupled with a Newton iteration technique. Superior stability characteristics of the new algorithm are demonstrated through applications to steady and oscillatory flows at subsonic and supersonic freestream conditions for an F-5 fighter wing. For steady flow calculations, the size of the time step is cycled to achieve rapid convergence. For unsteady flow calculations, the AF algorithm is sufficiently robust to allow the step size to be selected based on accuracy rather than on stability considerations. Therefore, accurate solutions are obtained in only several hundred time steps yielding a significant computational cost savings when compared to alternative methods.

Batina, John T.↗

An efficient algorithm for solution of the unsteady transonic small-disturbance equation

A time accurate approximation factorization (AF) algorithm is formulated for solution of the three-dimensional unsteady transonic small-disturbance equation. The AF algorithm consists of a time linearization procedure coupled with a Newton iteration technique. Superior stability characteristics of the new algorithm are demonstrated through applications to steady and oscillatory flows at subsonic and supersonic freestream conditions for an F-5 fighter wing. For steady flow calculations, the size of the time step is cycled to achieve rapid convergence. For unsteady flow calculations, the AF algorithm is sufficiently robust to allow the step size to be selected based on accuracy rather than on stability considerations. Therefore, accurate solutions are obtained in only several hundred time steps yielding a significant computational cost savings when compared to alternative methods.

Batina, John T.↗

An experiment on the use of disposable plastics as a reinforcement in concrete beams

Illustrated here is the concept of reinforced concrete structures by the use of computer simulation and an inexpensive hands-on design experiment. The students in our construction management program use disposable plastic as a reinforcement to demonstrate their understanding of reinforced concrete and prestressed concrete beams. The plastics used for such an experiment vary from plastic bottles to steel reinforced auto tires. This experiment will show the extent to which plastic reinforcement increases the strength of a concrete beam. The procedure of using such throw-away plastics in an experiment to explain the interaction between the reinforcement material and concrete, and a comparison of the test results for using different types of waste plastics are discussed. A computer analysis to simulate the structural response is used to compare the test results and to understand the analytical background of reinforced concrete design. This interaction of using computers to analyze structures and to relate the output results with real experimentation is found to be a very useful method for teaching a math-based analytical subject to our non-engineering students.

Chowdhury, Mostafiz R.↗

Line relaxation methods for the solution of 2D and 3D compressible flows

An implicit finite element based algorithm for the compressible Navier-Stokes equations is outlined, and the solution of the resulting equation by a line relaxation on general meshes of triangles or tetrahedra is described. The problem of generating and adapting unstructured meshes for viscous flows is reexamined, and an approach for both 2D and 3D simulations is proposed. An efficient approach appears to be the use of an implicit/explicit procedure, with the implicit treatment being restricted to those regions of the mesh where viscous effects are known to be dominant. Numerical examples demonstrating the computational performance of the proposed techniques are given.

Hassan, O.↗

Fragmented Land Cover Types and Estimation of Area with Course Spatial Resolution Imagery

Imagery of coarse resolution, such weather satellite imagery with 1 sq km pixels, is increasingly used to monitor dynamic and fragmented types of land surface types, such as scars from recent fires and ponds in wetlands. Accurate estimates of these land cover types at regional to global scales are required to assess the roles of fires and wetlands in global warming, yet difficult to compute when much of the area is accounted for by fragments about the same size as the pixels. In previous research, we found that size distribution of the fragments in several example scenes fit simple two-parameter models and related effects of coarse resolution to errors in area estimates based on pixel counts. We report on progress to develop accurate area estimates based on modelling the size distribution of the fragments, including analysis of size distributions on an expanded set of maps developed from digital imagery and a test of a procedure to correct for effects of coarse spatial resolution.

Hlavka, Chris↗

New method for fitting coefficients in standard model effective theory

We present an alternative method for carrying out a principal-component analysis of Wilson coefficients in standard model effective field theory (SMEFT). The method is based on singular-value decomposition (SVD). The SVD method provides information about the sensitivity of experimental observables to physics beyond the standard model that is not accessible in the Fisher-information method. In principle, the SVD method can also have computational advantages over diagonalization of the Fisher information matrix. We demonstrate the SVD method by applying it to the dimension-6 coefficients for the process of top-quark decay to a b quark and a W boson and use this example to illustrate some pitfalls in widely used fitting procedures. We also outline an iterative procedure for applying the SVD method to dimension-8 SMEFT coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Efficient Scheduling of Recursive Control Flow on GPUs

Graphics processing units (GPUs) have rapidly emerged as a very significant player in high performance computing. Single instruction multiple thread (SIMT) pipelines are typically used in GPUs to exploit parallelism and maximize performance. Although support for unstructured control flow has been included in GPUs, efficiently managing thread divergence for arbitrary parallel programs remains a critical challenge. In this paper, we focus on the problem of supporting recursion in modern GPUs. We design and comparatively evaluate various algorithms to manage thread divergence encountered in recursive programs. The results improve upon traditional post-dominator based reconvergence mechanisms designed to handle thread divergence due to control flow within a procedure.

GPUs↗

Proactive Data Containers for Scientific Storage (Final Report)

Emerging HPC systems are expected to be deployed with an unprecedented level of complexity, due to a deep system memory/storage hierarchy and heterogeneity of the storage hardware. This hierarchy is expected to range from CPU cache through several levels of volatile memory to nonvolatile memory, traditional hard disks, and tape. Simple and efficient methods of data management and movement through this hierarchy is critical for scientific applications using exascale systems. Existing storage system and I/O (SSIO) technologies face severe challenges in dealing with these requirements. POSIX and MPI I/O standards that are the basis for existing I/O libraries and parallel file systems present fundamental challenges in the areas of scalable metadata operations, semantics-based data movement performance tuning, asynchronous operation, and support for scalable consistency of distributed operations. Moving toward new paradigms for SSIO in the extreme-scale era, we have proposed to investigate novel object-based data abstractions and storage mechanisms that take advantage of the deep storage hierarchy and enable proactive automated performance tuning. In order to achieve these overarching goals, we initiated an effort to develop a fundamental new data abstraction, called Proactive Data Containers (PDC). A PDC is a container within a locus of storage (memory, NVRAM, disk, etc.) that stores science data in an object-centric manner. Managing data as objects enables powerful optimization opportunities for data movement and transformations. The R&D focus of this project are: 1) formulation of object-oriented PDCs and their mapping in different levels of the exascale storage hierarchy; 2) efficient strategies for moving data in deep storage hierarchies using PDCs; 3) techniques for transforming and reorganizing data based on application requirements; and 4) novel analysis paradigms for enabling data transformations and user-defined analysis on data in PDCs. Toward achieving these overarching goals, we designed an object-centric application programing interface (API) for HPC, scalable metadata management for object-centric storage systems, and data movement optimizations such as Data Elevator for moving data between two levels of storage devices and TAPIOCA for efficient aggregation of data on compute nodes. We then implemented several components of the PDC system. They include metadata management, data placement services, remote procedure calls, data aggregation, etc. We have put them together into the overall PDC framework.

97 MATHEMATICS AND COMPUTING↗

Cyclic plasticity models and application in fatigue analysis

An analytical procedure for prediction of the cyclic plasticity effects on both the structural fatigue life to crack initiation and the rate of crack growth is presented. The crack initiation criterion is based on the Coffin-Manson formulae extended for multiaxial stress state and for inclusion of the mean stress effect. This criterion is also applied for the accumulated damage ahead of the existing crack tip which is assumed to be related to the crack growth rate. Three cyclic plasticity models, based on the concept of combination of several yield surfaces, are employed for computing the crack growth rate of a crack plane stress panel under several cyclic loading conditions.

Kalev, I.↗

Investigation of the transient fuel preburner manifold and combustor

A computational fluid dynamics (CFD) model with finite rate reactions, FDNS, was developed to study the start transient of the Space Shuttle Main Engine (SSME) fuel preburner (FPB). FDNS is a time accurate, pressure based CFD code. An upwind scheme was employed for spatial discretization. The upwind scheme was based on second and fourth order central differencing with adaptive artificial dissipation. A state of the art two-equation k-epsilon (T) turbulence model was employed for the turbulence calculation. A Pade' Rational Solution (PARASOL) chemistry algorithm was coupled with the point implicit procedure. FDNS was benchmarked with three well documented experiments: a confined swirling coaxial jet, a non-reactive ramjet dump combustor, and a reactive ramjet dump combustor. Excellent comparisons were obtained for the benchmark cases. The code was then used to study the start transient of an axisymmetric SSME fuel preburner. Predicted transient operation of the preburner agrees well with experiment. Furthermore, it was also found that an appreciable amount of unburned oxygen entered the turbine stages.

Wang, Ten-See↗

Method for calculating the rolling and yawing moments due to rolling for unswept wings with or without flaps or ailerons by use of nonlinear section lift data

The methods of NACA reports 865 and 1090 have been applied to the calculation of the rolling-moment and yawing-moment coefficients due to rolling for unswept wings with or without flaps or ailerons. The methods are based on lifting-line theory and allow the use of nonlinear section lift data. The method presented in this report permits calculations to be made somewhat beyond maximum lift for wings having no twist or continuous twist and employing airfoil sections which do not display large discontinuities in the lift curves. Calculations can be made up to maximum lift for wings with discontinuous twist such as that produced by partial-span flaps or ailerons, or both. Two calculated examples are presented in simplified computing forms in order to illustrate the procedures involve.

Martina, Albert P↗