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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SE2-ANL Subchannel Analysis with New Heating Calculation Based on VARIANT Transport Code

The multi-assembly, steady-state sub-channel analysis code SE2-ANL has been widely used in thermal-hydraulic analysis of sodium cooled fast reactors (SFR). It is used to compute the core-wide temperature profiles for a given flow distribution and to determine the assembly flow rates yielding a desired (assembly-wise) peak temperature distribution. It is a modified version of SUPERENERGY-2 coupled with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In the current version of SE2-ANL, the neutron and gamma fluxes are calculated using the finite-difference diffusion theory option of the DIF3D code with six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. In order to determine the pin power distributions more accurately by eliminating these limitations, a new coupled heating calculation procedure has recently developed based on the neutron and gamma transport calculations with the VARIANT code. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at the center of each fuel pin and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The SE2-ANL code has been modified to interface with the new, coupled neutron and gamma heating calculation scheme. In addition, several improvements have been made on the computational methods and models of SE2-ANL. To eliminate the limit on the axial mesh size due to the numerical instability problem of the current explicit scheme, the numerical scheme for discretizing the energy equation in the axial direction has been changed to the {theta}-method of the SLTHEN code. An automated flow allocation scheme has also been implemented in SE2-ANL to determine assembly flow rates such that the peak cladding mid-wall temperatures of individual orifice zones are equalized over the burn cycle. The impacts of the new heating calculation scheme based on VARIANT transport solutions on the cladding and fuel temperatures have been tested against the existing method of SE2-ANL. Using the existing and new heating calculation schemes, full-core subchannel analyses were performed for the metal core design of the 1000 MWt Advanced Burner Reactor (ABR-1000), and the resulting coolant, cladding and fuel temperatures were compared. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1990

An adaptive refinement procedure for transient thermal analysis using nodeless variable finite elements

An adaptive mesh refinement procedure that uses nodeless variables and quadratic interpolation functions is presented for analyzing transient thermal problems. A temperature based finite element scheme with Crank-Nicolson time marching is used to obtain the thermal solution. The strategies used for mesh adaption, computing refinement indicators, and time marching are described. Examples in one and two dimensions are presented and comparisons are made with exact solutions. The effectiveness of this procedure for transient thermal analysis is reflected in good solution accuracy, reduction in number of elements used, and computational efficiency.

Ramakrishnan, R.↗
NASA NTRS1979

Computational optimization and wind tunnel test of transonic wing designs

A practical procedure for the optimum design of transonic wings is demonstrated. The procedure uses an optimization program based on the method of feasible directions coupled with an aerodynamic analysis program which solves the three-dimensional potential equation for subsonic through transonic flow. Two new wings for the A-7 aircraft were designed by using the optimization procedure to achieve specified surface pressure distributions. The new wings, along with the existing A-7 wing, were tested in the Ames 11-foot transonic wind tunnel. The experimental data show that all of the performance goals were met.

Haney, H. P.↗
NASA NTRS1987

Algebraic grid generation

Algebraic grid generation is the direct expression of a physical coordinate system as a function of a uniform grid in a rectangular computational coordinate system. Algebraic grid generation is based on mathematical interpolation and is presented in general terms of multivariate transfinite interpolation. The multisurface method and the two-boundary technique are described as univariate procedures that can be applied within the context of transfinite interpolation. A technique for grid clustering is described. Problems that are commonly encountered in three-dimensional grid generation are discussed and approaches for dealing with complex physical domains using multiple computational grid blocks are presented.

Smith, Robert E.↗
NASA NTRS1993

The effects of a dynamic graphical model during simulation-based training of console operation skill

LOADER is a Windows-based simulation of a complex procedural task. The task requires subjects to execute long sequences of console-operation actions (e.g., button presses, switch actuations, dial rotations) to accomplish specific goals. The LOADER interface is a graphical computer-simulated console which controls railroad cars, tracks, and cranes in a fictitious railroad yard. We hypothesized that acquisition of LOADER performance skill would be supported by the representation of a dynamic graphical model linking console actions to goal and goal states in the 'railroad yard'. Twenty-nine subjects were randomly assigned to one of two treatments (i.e., dynamic model or no model). During training, both groups received identical text-based instruction in an instructional-window above the LOADER interface. One group, however, additionally saw a dynamic version of the bird's-eye view of the railroad yard. After training, both groups were tested under identical conditions. They were asked to perform the complete procedure without guidance and without access to either type of railroad yard representation. Results indicate that rather than becoming dependent on the animated rail yard model, subjects in the dynamic model condition apparently internalized the model, as evidenced by their performance after the model was removed.

Farquhar, John D.↗
DOE OSTI2016

A model-based scatter artifacts correction for cone beam CT

Purpose: Due to the increased axial coverage of multislice computed tomography (CT) and the introduction of flat detectors, the size of x-ray illumination fields has grown dramatically, causing an increase in scatter radiation. For CT imaging, scatter is a significant issue that introduces shading artifact, streaks, as well as reduced contrast and Hounsfield Units (HU) accuracy. The purpose of this work is to provide a fast and accurate scatter artifacts correction algorithm for cone beam CT (CBCT) imaging. Methods: The method starts with an estimation of coarse scatter profiles for a set of CBCT data in either image domain or projection domain. A denoising algorithm designed specifically for Poisson signals is then applied to derive the final scatter distribution. Qualitative and quantitative evaluations using thorax and abdomen phantoms with Monte Carlo (MC) simulations, experimental Catphan phantom data, and in vivo human data acquired for a clinical image guided radiation therapy were performed. Scatter correction in both projection domain and image domain was conducted and the influences of segmentation method, mismatched attenuation coefficients, and spectrum model as well as parameter selection were also investigated. Results: Results show that the proposed algorithm can significantly reduce scatter artifacts and recover the correct HU in either projection domain or image domain. For the MC thorax phantom study, four-components segmentation yields the best results, while the results of three-components segmentation are still acceptable. The parameters (iteration number K and weight β) affect the accuracy of the scatter correction and the results get improved as K and β increase. It was found that variations in attenuation coefficient accuracies only slightly impact the performance of the proposed processing. For the Catphan phantom data, the mean value over all pixels in the residual image is reduced from −21.8 to −0.2 HU and 0.7 HU for projection domain and image domain, respectively. The contrast of the in vivo human images is greatly improved after correction. Conclusions: The software-based technique has a number of advantages, such as high computational efficiency and accuracy, and the capability of performing scatter correction without modifying the clinical workflow (i.e., no extra scan/measurement data are needed) or modifying the imaging hardware. When implemented practically, this should improve the accuracy of CBCT image quantitation and significantly impact CBCT-based interventional procedures and adaptive radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Orbital-dependent second-order scaled-opposite-spin correlation functionals in the optimized effective potential method

The performance of correlated optimized effective potential (OEP) functionals based on the spin-resolved second-order correlation energy is analysed. The relative importance of singly- and doubly- excited contributions as well as the effect of scaling the same- and opposite- spin components is investigated in detail comparing OEP results with Kohn–Sham (KS) quantities determined via an inversion procedure using accurate ab initio electronic densities. Special attention is dedicated in particular to the recently proposed scaled-opposite–spin OEP functional [I. Grabowski, E. Fabiano, and F. Della Sala, Phys. Rev. B 87, 075103 (2013)] which is the most advantageous from a computational point of view. We find that for high accuracy, a careful, system dependent, selection of the scaling coefficient is required. We analyse several size-extensive approaches for this selection. Finally, we find that a composite approach, named OEP2-SOSh, based on a post-SCF rescaling of the correlation energy can yield high accuracy for many properties, being comparable with the most accurate OEP procedures previously reported in the literature but at substantially reduced computational effort.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Feynman formulae and phase space Feynman path integrals for tau-quantization of some Lévy-Khintchine type Hamilton functions

Evolution semigroups generated by pseudo-differential operators are considered. These operators are obtained by different (parameterized by a number τ) procedures of quantization from a certain class of functions (or symbols) defined on the phase space. This class contains Hamilton functions of particles with variable mass in magnetic and potential fields and more general symbols given by the Lévy-Khintchine formula. The considered semigroups are represented as limits of n-fold iterated integrals when n tends to infinity. Such representations are called Feynman formulae. Some of these representations are constructed with the help of another pseudo-differential operator, obtained by the same procedure of quantization; such representations are called Hamiltonian Feynman formulae. Some representations are based on integral operators with elementary kernels; these are called Lagrangian Feynman formulae. Langrangian Feynman formulae provide approximations of evolution semigroups, suitable for direct computations and numerical modeling of the corresponding dynamics. Hamiltonian Feynman formulae allow to represent the considered semigroups by means of Feynman path integrals. In the article, a family of phase space Feynman pseudomeasures corresponding to different procedures of quantization is introduced. The considered evolution semigroups are represented as phase space Feynman path integrals with respect to these Feynman pseudomeasures, i.e., different quantizations correspond to Feynman path integrals with the same integrand but with respect to different pseudomeasures. This answers Berezin’s problem of distinguishing a procedure of quantization on the language of Feynman path integrals. Moreover, the obtained Lagrangian Feynman formulae allow also to calculate these phase space Feynman path integrals and to connect them with some functional integrals with respect to probability measures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

High temperature cyclic oxidation furnace testing at NASA Lewis Research Center

A standardized method of testing the cyclic oxidation resistance of various alloys in static air to 1200 C was developed and is routinely used at NASA Lewis Research Center. Test samples are automatically raised and lowered into a resistance wound furnace for a series of fixed interval heating and cooling cycles. Spall catchers collect the accumulated spall from each sample. The samples are weighed intermittently to generate specific weight change/time data. At various test times the samples and the accumulated spall are analyzed by X-ray diffraction. A computer program uses this gravimetric and X-ray data as input to print out the oxidation curves and specific weight change/time and X-ray results in a published format, organizes, and indexes the data. So far, several hundred Fe, Ni, and Co base alloys were tested using this same basic procedure and results form the basis of a series of cyclic oxidation handbooks to be published by NASA. Such specific weight change/time data were used to estimate the oxidative metal consumption by several computer modeling techniques to rank alloys and for use in life testing estimates.

Barrett, C. A.↗
DOE OSTI2015

Development of efficient time-evolution method based on three-term recurrence relation

The advantage of the real-time (RT) propagation method is a direct solution of the time-dependent Schrödinger equation which describes frequency properties as well as all dynamics of a molecular system composed of electrons and nuclei in quantum physics and chemistry. Its applications have been limited by computational feasibility, as the evaluation of the time-evolution operator is computationally demanding. In this article, a new efficient time-evolution method based on the three-term recurrence relation (3TRR) was proposed to reduce the time-consuming numerical procedure. The basic formula of this approach was derived by introducing a transformation of the operator using the arcsine function. Since this operator transformation causes transformation of time, we derived the relation between original and transformed time. The formula was adapted to assess the performance of the RT time-dependent Hartree-Fock (RT-TDHF) method and the time-dependent density functional theory. Compared to the commonly used fourth-order Runge-Kutta method, our new approach decreased computational time of the RT-TDHF calculation by about factor of four, showing the 3TRR formula to be an efficient time-evolution method for reducing computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Determination of solid/porous wall boundary conditions from wind tunnel data for computational fluid dynamics codes

A computational and experimental study has been undertaken to investigate methods of modeling solid and porous wall boundary conditions in CFD codes. The procedure utilizes experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution is then imposed as a boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the boundary conditions at the wind tunnel walls using only sparse experimental measurements has been investigated. Verification of the approach using computational studies has been carried out using an incompressible flow solver. The current work demonstrates this technique for low speed flows and compares the result with experimental data obtained from a heavily instrumented variable porosity test section. Position and refinement of experimental measurements required to describe porous wall boundary conditions has also been considered for application to other porous wall wind tunnels.

Beutner, Thomas J.↗
NASA NTRS1987

Non-linear analysis of symmetric structures with unsymmetric boundary conditions

Two computational techniques are developed analytically to generate global approximation vectors for use in the nonlinear large-deflection analysis of symmetric structures with asymmetric boundary conditions. One method is based on linear combinations of FEM-generated symmetric and antisymmetric vector components, while the other employs a preconditioned-conjugate-gradient technique. The derivation of the procedures is explained in detail, and their effectiveness is evaluated in sample computations for an elliptic toroid and a cylindrical panel; the results are presented in tables and graphs and characterized.

Noor, Ahmed K.↗
DOE OSTI2016

TU-FG-201-12: Designing a Risk-Based Quality Assurance Program for a Newly Implemented Y-90 Microspheres Procedure

Purpose: To create a quality assurance program based upon a risk-based assessment of a newly implemented SirSpheres Y-90 procedure. Methods: A process map was created for a newly implemented SirSpheres procedure at a community hospital. The process map documented each step of this collaborative procedure, as well as the roles and responsibilities of each member. From the process map, different potential failure modes were determined as well as any current controls in place. From this list, a full failure mode and effects analysis (FMEA) was performed by grading each failure mode’s likelihood of occurrence, likelihood of detection, and potential severity. These numbers were then multiplied to compute the risk priority number (RPN) for each potential failure mode. Failure modes were then ranked based on their RPN. Additional controls were then added, with failure modes corresponding to the highest RPNs taking priority. Results: A process map was created that succinctly outlined each step in the SirSpheres procedure in its current implementation. From this, 72 potential failure modes were identified and ranked according to their associated RPN. Quality assurance controls and safety barriers were then added for failure modes associated with the highest risk being addressed first. Conclusion: A quality assurance program was created from a risk-based assessment of the SirSpheres process. Process mapping and FMEA were effective in identifying potential high-risk failure modes for this new procedure, which were prioritized for new quality assurance controls. TG 100 recommends the fault tree analysis methodology to design a comprehensive and effective QC/QM program, yet we found that by simply introducing additional safety barriers to address high RPN failure modes makes the whole process simpler and safer.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1999

The Application of a Boundary Integral Equation Method to the Prediction of Ducted Fan Engine Noise

The prediction of ducted fan engine noise using a boundary integral equation method (BIEM) is considered. Governing equations for the BIEM are based on linearized acoustics and describe the scattering of incident sound by a thin, finite-length cylindrical duct in the presence of a uniform axial inflow. A classical boundary value problem (BVP) is derived that includes an axisymmetric, locally reacting liner on the duct interior. Using potential theory, the BVP is recast as a system of hypersingular boundary integral equations with subsidiary conditions. We describe the integral equation derivation and solution procedure in detail. The development of the computationally efficient ducted fan noise prediction program TBIEM3D, which implements the BIEM, and its utility in conducting parametric noise reduction studies are discussed. Unlike prediction methods based on spinning mode eigenfunction expansions, the BIEM does not require the decomposition of the interior acoustic field into its radial and axial components which, for the liner case, avoids the solution of a difficult complex eigenvalue problem. Numerical spectral studies are presented to illustrate the nexus between the eigenfunction expansion representation and BIEM results. We demonstrate BIEM liner capability by examining radiation patterns for several cases of practical interest.

Dunn, M. H.↗
NASA NTRS2001

Self-motion perception: assessment by real-time computer-generated animations

We report a new procedure for assessing complex self-motion perception. In three experiments, subjects manipulated a 6 degree-of-freedom magnetic-field tracker which controlled the motion of a virtual avatar so that its motion corresponded to the subjects' perceived self-motion. The real-time animation created by this procedure was stored using a virtual video recorder for subsequent analysis. Combined real and illusory self-motion and vestibulo-ocular reflex eye movements were evoked by cross-coupled angular accelerations produced by roll and pitch head movements during passive yaw rotation in a chair. Contrary to previous reports, illusory self-motion did not correspond to expectations based on semicircular canal stimulation. Illusory pitch head-motion directions were as predicted for only 37% of trials; whereas, slow-phase eye movements were in the predicted direction for 98% of the trials. The real-time computer-generated animations procedure permits use of naive, untrained subjects who lack a vocabulary for reporting motion perception and is applicable to basic self-motion perception studies, evaluation of motion simulators, assessment of balance disorders and so on.

Non-NASA Center↗
NASA NTRS2002

Evaluation of Operational Procedures for Using a Time-Based Airborne Inter-arrival Spacing Tool

An airborne tool has been developed based on the concept of an aircraft maintaining a time-based spacing interval from the preceding aircraft. The Advanced Terminal Area Approach Spacing (ATAAS) tool uses Automatic Dependent Surveillance-Broadcast (ADS-B) aircraft state data to compute a speed command for the ATAAS-equipped aircraft to obtain a required time interval behind another aircraft. The tool and candidate operational procedures were tested in a high-fidelity, full mission simulator with active airline subject pilots flying an arrival scenario using three different modes for speed control. The objectives of this study were to validate the results of a prior Monte Carlo analysis of the ATAAS algorithm and to evaluate the concept from the standpoint of pilot acceptability and workload. Results showed that the aircraft was able to consistently achieve the target spacing interval within one second (the equivalent of approximately 220 ft at a final approach speed of 130 kt) when the ATAAS speed guidance was autothrottle-coupled, and a slightly greater (4-5 seconds), but consistent interval with the pilot-controlled speed modes. The subject pilots generally rated the workload level with the ATAAS procedure as similar to that with standard procedures, and also rated most aspects of the procedure high in terms of acceptability. Although pilots indicated that the head-down time was higher with ATAAS, the acceptability of head-down time was rated high. Oculometer data indicated slight changes in instrument scan patterns, but no significant change in the amount of time spent looking out the window between the ATAAS procedure versus standard procedures.

Oseguera-Lohr, Rosa M.↗
NASA NTRS1987

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗
DOE OSTI2015

A mass conserving level set method for detailed numerical simulation of liquid atomization

An improved mass conserving level set method for detailed numerical simulations of liquid atomization is developed to address the issue of mass loss in the existing level set method. This method introduces a mass remedy procedure based on the local curvature at the interface, and in principle, can ensure the absolute mass conservation of the liquid phase in the computational domain. Three benchmark cases, including Zalesak's disk, a drop deforming in a vortex field, and the binary drop head-on collision, are simulated to validate the present method, and the excellent agreement with exact solutions or experimental results is achieved. It is shown that the present method is able to capture the complex interface with second-order accuracy and negligible additional computational cost. The present method is then applied to study more complex flows, such as a drop impacting on a liquid film and the swirling liquid sheet atomization, which again, demonstrates the advantages of mass conservation and the capability to represent the interface accurately.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗