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,117 records · Page 62

Computer aided modelling/compensator design for a flexible space antenna

Controller design algorithms are developed to produce simultaneously both a model of the plant and a compensator. The size of the model and properties of the compensator are driven by the performance requirements, the disturbance environment, and the location, number and type of sensors and actuators. The procedure is based on linear optimal control theory for distributed systems, and balanced realization theory is used to guide the development of the model and reduce the order of the compensator.

Gibson, J. S.↗

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models

Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.

59 BASIC BIOLOGICAL SCIENCES↗

Lung motion estimation using dynamic point shifting: An innovative model based on a robust point matching algorithm

Purpose: Image-guided radiotherapy is an advanced 4D radiotherapy technique that has been developed in recent years. However, respiratory motion causes significant uncertainties in image-guided radiotherapy procedures. To address these issues, an innovative lung motion estimation model based on a robust point matching is proposed in this paper. Methods: An innovative robust point matching algorithm using dynamic point shifting is proposed to estimate patient-specific lung motion during free breathing from 4D computed tomography data. The correspondence of the landmark points is determined from the Euclidean distance between the landmark points and the similarity between the local images that are centered at points at the same time. To ensure that the points in the source image correspond to the points in the target image during other phases, the virtual target points are first created and shifted based on the similarity between the local image centered at the source point and the local image centered at the virtual target point. Second, the target points are shifted by the constrained inverse function mapping the target points to the virtual target points. The source point set and shifted target point set are used to estimate the transformation function between the source image and target image. Results: The performances of the authors’ method are evaluated on two publicly available DIR-lab and POPI-model lung datasets. For computing target registration errors on 750 landmark points in six phases of the DIR-lab dataset and 37 landmark points in ten phases of the POPI-model dataset, the mean and standard deviation by the authors’ method are 1.11 and 1.11 mm, but they are 2.33 and 2.32 mm without considering image intensity, and 1.17 and 1.19 mm with sliding conditions. For the two phases of maximum inhalation and maximum exhalation in the DIR-lab dataset with 300 landmark points of each case, the mean and standard deviation of target registration errors on the 3000 landmark points of ten cases by the authors’ method are 1.21 and 1.04 mm. In the EMPIRE10 lung registration challenge, the authors’ method ranks 24 of 39. According to the index of the maximum shear stretch, the authors’ method is also efficient to describe the discontinuous motion at the lung boundaries. Conclusions: By establishing the correspondence of the landmark points in the source phase and the other target phases combining shape matching and image intensity matching together, the mismatching issue in the robust point matching algorithm is adequately addressed. The target registration errors are statistically reduced by shifting the virtual target points and target points. The authors’ method with consideration of sliding conditions can effectively estimate the discontinuous motion, and the estimated motion is natural. The primary limitation of the proposed method is that the temporal constraints of the trajectories of voxels are not introduced into the motion model. However, the proposed method provides satisfactory motion information, which results in precise tumor coverage by the radiation dose during radiotherapy.

60 APPLIED LIFE SCIENCES↗

A simulation of air pollution model parameter estimation using data from a ground-based LIDAR remote sensor

One way to obtain estimates of the unknown parameters in a pollution dispersion model is to compare the model predictions with remotely sensed air quality data. A ground-based LIDAR sensor provides relative pollution concentration measurements as a function of space and time. The measured sensor data are compared with the dispersion model output through a numerical estimation procedure to yield parameter estimates which best fit the data. This overall process is tested in a computer simulation to study the effects of various measurement strategies. Such a simulation is useful prior to a field measurement exercise to maximize the information content in the collected data. Parametric studies of simulated data matched to a Gaussian plume dispersion model indicate the trade offs available between estimation accuracy and data acquisition strategy.

Kibler, J. F.↗

Progressive fracture of fiber composites

Refined models and procedures are described for determining progressive composite fracture in graphite/epoxy angleplied laminates. Lewis Research Center capabilities are utilized including the Real Time Ultrasonic C Scan (RUSCAN) experimental facility and the Composite Durability Structural Analysis (CODSTRAN) computer code. The CODSTRAN computer code is used to predict the fracture progression based on composite mechanics, finite element stress analysis, and fracture criteria modules. The RUSCAN facility, CODSTRAN computer code, and scanning electron microscope are used to determine durability and identify failure mechanisms in graphite/epoxy composites.

Irvin, T. B.↗

Advanced feature extraction in remote sensing using artificial intelligence and geographic information systems

Traditional computer-assisted image-analysis techniques in remote sensing lag well behind human abilities in terms of both speed and accuracy. A fundamental limitation of computer-assisted techniques is their inability to assimilate a variety of different data types leading to an interpretation in a manner similar to human image interpretation. Expert systems and computer-vision techniques are proposed as a potential solution to these limitations. Some aspects of human expertise in image analysis may be codified into expert systems. Image understanding and symbolic reasoning provide a means of assimilating spatial information and spatial reasoning into the analysis procedure. Knowledge-based image-analysis systems incorporate many of these concepts and have been implemented for some well defined problem domains. Geographic information systems represent an excellent environment for this type of analysis, providing both analytic tools and contextual information to the analysis procedure.

Estes, John E.↗

A Rapid Aerodynamic Design Procedure Based on Artificial Neural Networks

An aerodynamic design procedure that uses neural networks to model the functional behavior of the objective function in design space has been developed. This method incorporates several improvements to an earlier method that employed a strategy called parameter-based partitioning of the design space in order to reduce the computational costs associated with design optimization. As with the earlier method, the current method uses a sequence of response surfaces to traverse the design space in search of the optimal solution. The new method yields significant reductions in computational costs by using composite response surfaces with better generalization capabilities and by exploiting synergies between the optimization method and the simulation codes used to generate the training data. These reductions in design optimization costs are demonstrated for a turbine airfoil design study where a generic shape is evolved into an optimal airfoil.

Rai, Man Mohan↗

Heats of Segregation of BCC Metals Using Ab Initio and Quantum Approximate Methods

Many multicomponent alloys exhibit surface segregation, in which the composition at or near a surface may be substantially different from that of the bulk. A number of phenomenological explanations for this tendency have been suggested, involving, among other things, differences among the components' surface energies, molar volumes, and heats of solution. From a theoretical standpoint, the complexity of the problem has precluded a simple, unified explanation, thus preventing the development of computational tools that would enable the identification of the driving mechanisms for segregation. In that context, we investigate the problem of surface segregation in a variety of bcc metal alloys by computing dilute-limit heats of segregation using both the quantum-approximate energy method of Bozzolo, Ferrante and Smith (BFS), and all-electron density functional theory. In addition, the composition dependence of the heats of segregation is investigated using a BFS-based Monte Carlo procedure, and, for selected cases of interest, density functional calculations. Results are discussed in the context of a simple picture that describes segregation behavior as the result of a competition between size mismatch and alloying effects

Good, Brian↗

Towards a Certified Lightweight Array Bound Checker for Java Bytecode

Dynamic array bound checks are crucial elements for the security of a Java Virtual Machines. These dynamic checks are however expensive and several static analysis techniques have been proposed to eliminate explicit bounds checks. Such analyses require advanced numerical and symbolic manipulations that 1) penalize bytecode loading or dynamic compilation, 2) complexify the trusted computing base. Following the Foundational Proof Carrying Code methodology, our goal is to provide a lightweight bytecode verifier for eliminating array bound checks that is both efficient and trustable. In this work, we define a generic relational program analysis for an imperative, stackoriented byte code language with procedures, arrays and global variables and instantiate it with a relational abstract domain as polyhedra. The analysis has automatic inference of loop invariants and method pre-/post-conditions, and efficient checking of analysis results by a simple checker. Invariants, which can be large, can be specialized for proving a safety policy using an automatic pruning technique which reduces their size. The result of the analysis can be checked efficiently by annotating the program with parts of the invariant together with certificates of polyhedral inclusions. The resulting checker is sufficiently simple to be entirely certified within the Coq proof assistant for a simple fragment of the Java bytecode language. During the talk, we will also report on our ongoing effort to scale this approach for the full sequential JVM.

Pichardie, David↗

A review of neutronics and thermal hydraulics–based screening methods applied to accelerated nuclear fuel qualification

This paper reviews the state-of-the-art engineering approach for using thermal hydraulic (TH) and neutronics modeling and simulation (M&S) tools to perform rapid screening studies of novel nuclear fuel concepts within the context of accelerated fuel qualification. Global research efforts have introduced nuclear fuel and material concepts that mark a significant departure from traditional reactor materials. The number of new technologies being considered for development for light water reactors and advanced reactor types has created the need for an accelerated fuel qualification procedure. A key component of this procedure is the rapid identification of the most promising fuel concepts using computational screening studies. Advanced TH and neutronic M&S tools should be leveraged to efficiently determine whether the reactor performance and safety characteristics of a given concept warrant additional studies or whether the concept requires modification or elimination. This paper reviews best practices for performing these TH and neutronics screening studies at various stages during a fuel concept's progression through the qualification procedure. The motivation behind standardizing this approach is to minimize time and resources spent on qualification activities for fuel concepts that could be quickly refined or eliminated from consideration based on their reactor physics and TH characteristics. Adoption of this screening procedure—which focuses primarily on nuclear fuels but may be applicable to other reactor materials—will also help accelerate new material qualification by generating boundary conditions crucial to fuel performance evaluations and highlighting needed areas of separate effects experimentation. This article reviews the motivation behind the introduction of novel nuclear fuel concepts, provides incentive for utilizing TH- and neutronics-based screening studies, describes the screening approach and methodology, and includes discussion on how to interpret screening results to provide recommendations for the continued development of a given concept.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An implicit marching procedure for the treatment of supersonic flow fields using the conservative full potential equation

A new aerodynamic prediction technique based on the conservative full potential equation is developed for the treatment of supersonic flow fields. This new technique bridges the gap between simplistic linear theory methods and complex Euler solvers. A novel local density linearization concept and a second order accurate retarded density scheme, both producing the correct artificial viscosity, are introduced in developing an implicit marching scheme for solving the scalar phi. The method produces results that compare well with Euler solvers and requires an order of magnitude less computer time and significantly less computer memory over existing nonlinear codes. The scalar phi formulation can be extended to handle subsonic pockets in the marching direction and also is suitable for developing inverse procedures where the shape corresponding to a prescribed loading is sought.

Shankar, V.↗

Algebraic Nonoverlapping Domain Decomposition Methods for Stabilized FEM and FV Discretizations

We consider preconditioning methods for convection dominated fluid flow problems based on a nonoverlapping Schur complement domain decomposition procedure for arbitrary triangulated domains. The triangulation is first partitioned into a number of subdomains and interfaces which induce a natural 2 x 2 partitioning of the p.d.e. discretization matrix. We view the Schur complement induced by this partitioning as an algebraically derived coarse space approximation. This avoids the known difficulties associated with the direct formation of an effective coarse discretization for advection dominated equations. By considering various approximations of the block factorization of the 2 x 2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 using MPI message passing protocol. A number of 2-D CFD calculations will be presented for both scalar advection-diffusion equations and the Euler equations discretized using stabilized finite element and finite volume methods. These results show very good scalability of the preconditioner for various discretizations as the number of processors is increased while the number of degrees of freedom per processor is fixed.

Barth, Timothy J.↗

SCALE Enhancements for Detailed Cask Dose Rate Analysis - Paper 37

The Nuclear Fuels Storage and Transportation Planning Project and the Used Fuel Disposition Campaign, both under the US Department of Energy (DOE) Office of Nuclear Energy, require the ability to predict k{sub eff}, heat generation and temperature distributions, and external dose rates from any combination of commercial spent nuclear fuel assemblies loaded into any type of licensed cask. Work in the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory has begun to develop the UNF-ST[and]DARDS system (Used Nuclear Fuel Storage, Transportation and Disposal Analysis Resource and Data System) that uses many of the nuclear analysis codes from the SCALE package, a comprehensive database of used nuclear fuel inventories from utilities, and a library of model template files to produce complete input files for the various types of UNF computational analyses. UNF-ST[and]DARDS uses the template files and input parameters that are applicable to a safety analysis to assemble and execute a complete input file, collect the output files and relevant information for subsequent calculations, and illustrate calculation results as a function of decay time and other relevant analysis parameters. This paper will focus on the shielding aspects and the ability to predict external dose rates from as-loaded casks. Realistic dose rate values based on actual canister fuel contents may be used in developing appropriate plans, procedures, and methods to keep personnel radiation exposures as low as reasonably achievable during cask repackaging and transportation operations, as well as in performing total system analyses and transportation risk evaluations. The MAVRIC shielding sequence of the SCALE package was created for the US Nuclear Regulatory Commission for analyzing dose rates outside of storage, transfer, or transportation casks loaded with spent nuclear fuel. MAVRIC contains automated variance reduction capabilities based on the CADIS and FW-CADIS methods that are effective for optimizing deep penetration Monte Carlo shielding problems for achieving low relative uncertainties in short computing times. MAVRIC uses the Denovo discrete ordinates code to compute importance maps and biased source terms for the Monaco fixed-source Monte Carlo code, which can use multi-group or continuous-energy cross section data. ORIGEN (Oak Ridge Isotope Generation and Depletion) is the SCALE computer code for calculating time-dependent isotopic concentrations during irradiation and decay for more than 2,200 nuclides that can be produced by nuclear fuel irradiation and activation. A series of SCALE enhancements was implemented to enable automatic decay heat and radiation source term generation and shielding safety analysis of as-loaded casks. The enhancements to SCALE made under this project will be a part of the SCALE 6.2 release expected in fall 2014. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SU-C-201-04: Quantification of Perfusion Heterogeneity Based On Texture Analysis for Fully Automatic Detection of Ischemic Deficits From Myocardial Perfusion Imaging

Purpose: Texture-based quantification of image heterogeneity has been a popular topic for imaging studies in recent years. As previous studies mainly focus on oncological applications, we report our recent efforts of applying such techniques on cardiac perfusion imaging. A fully automated procedure has been developed to perform texture analysis for measuring the image heterogeneity. Clinical data were used to evaluate the preliminary performance of such methods. Methods: Myocardial perfusion images of Thallium-201 scans were collected from 293 patients with suspected coronary artery disease. Each subject underwent a Tl-201 scan and a percutaneous coronary intervention (PCI) within three months. The PCI Result was used as the gold standard of coronary ischemia of more than 70% stenosis. Each Tl-201 scan was spatially normalized to an image template for fully automatic segmentation of the LV. The segmented voxel intensities were then carried into the texture analysis with our open-source software Chang Gung Image Texture Analysis toolbox (CGITA). To evaluate the clinical performance of the image heterogeneity for detecting the coronary stenosis, receiver operating characteristic (ROC) analysis was used to compute the overall accuracy, sensitivity and specificity as well as the area under curve (AUC). Those indices were compared to those obtained from the commercially available semi-automatic software QPS. Results: With the fully automatic procedure to quantify heterogeneity from Tl-201 scans, we were able to achieve a good discrimination with good accuracy (74%), sensitivity (73%), specificity (77%) and AUC of 0.82. Such performance is similar to those obtained from the semi-automatic QPS software that gives a sensitivity of 71% and specificity of 77%. Conclusion: Based on fully automatic procedures of data processing, our preliminary data indicate that the image heterogeneity of myocardial perfusion imaging can provide useful information for automatic determination of the myocardial ischemia.

60 APPLIED LIFE SCIENCES↗

Model selection and signal extraction using Gaussian Process regression

We present a novel computational approach for extracting localized signals from smooth background distributions. We focus on datasets that can be naturally presented as binned integer counts, demonstrating our procedure on the CERN open dataset with the Higgs boson signature, from the ATLAS collaboration at the Large Hadron Collider. Our approach is based on Gaussian Process (GP) regression — a powerful and flexible machine learning technique which has allowed us to model the background without specifying its functional form explicitly and separately measure the background and signal contributions in a robust and reproducible manner. Unlike functional fits, our GP-regression-based approach does not need to be constantly updated as more data becomes available. We discuss how to select the GP kernel type, considering trade-offs between kernel complexity and its ability to capture the features of the background distribution. We show that our GP framework can be used to detect the Higgs boson resonance in the data with more statistical significance than a polynomial fit specifically tailored to the dataset. Finally, we use Markov Chain Monte Carlo (MCMC) sampling to confirm the statistical significance of the extracted Higgs signature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fast computation algorithms for speckle pattern simulation

We present our development of a series of efficient computation algorithms, generally usable to calculate light diffraction and particularly for speckle pattern simulation. We use mainly the scalar diffraction theory in the form of Rayleigh-Sommerfeld diffraction formula and its Fresnel approximation. Our algorithms are based on a special form of the convolution theorem and the Fast Fourier Transform. They are able to evaluate the diffraction formula much faster than by direct computation and we have circumvented the restrictions regarding the relative sizes of the input and output domains, met on commonly used procedures. Moreover, the input and output planes can be tilted each to other and the output domain can be off-axis shifted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient algorithms for Hirshfeld-I charges

A new viewpoint on iterative Hirshfeld charges is presented, whereby the atomic populations obtained from such a scheme are interpreted as such populations which reproduce themselves. This viewpoint yields a self-consistent requirement for the Hirshfeld-I populations rather than being understood as the result of an iterative procedure. Based on this self-consistent requirement, much faster algorithms for Hirshfeld-I charges have been developed. In addition, new atomic reference densities for the Hirshfeld-I procedure are presented. The proposed reference densities are N-representable, display proper atomic shell structure and can be computed for any charged species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗