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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

TPSAS-NF1676L-13122-DND

An approach for assessing the delamination propagation and growth capabilities in commercial finite element codes was developed and demonstrated for the Virtual Crack Closure Technique (VCCT) implementations in ABAQUS?. The Double Cantilever Beam (DCB) specimen was chosen as an example. First, benchmark results to assess delamination propagation capabilities under static loading were created using models simulating specimens with different delamination lengths. Second, starting from an initially straight front, the delamination was allowed to propagate based on the algorithms implemented in the commercial finite element software. The load-displacement relationship obtained from the propagation analysis results and the benchmark results were compared. Good agreements could be achieved by selecting the appropriate input parameters, which were determined in an iterative procedure. A benchmark case to assess the fatigue growth prediction capabilities was also created. Input control parameters were varied to study the effect on the computed delamination growth. The benchmarks enabled the selection of the appropriate input parameters that yielded good agreements between the results obtained from the automated analysis and the benchmark results. Overall, the results are encouraging but further assessment for mixed-mode delamination is required.

Ronald Krueger↗
NASA NTRS1995

An Improved Ionospheric Correction Method for Wide-Area Augmentation Systems

A powerful approach for generating ionospheric corrections in wide area differential GPS applications has been developed, that can be applied to the Federal Aviation Administration's Wide Area Augmentation System (WAAS). This approach has been used to support ionospheric calibration for NASA's Deep Space Network and will be supporting real-time operations for global ionospheric specification and, possibly, single-frequency satellite altimeter calibration. It is a real-time, grid-based technique relying on a computationally efficient Kalman-type filter to produce accurate, smoothly varying ionospheric correction maps over the coverage area. Formal error maps are also computed, providing vertical delay errors over the WAAS grid, which are useful in integrity monitoring. This solution to the ionospheric correction problem is self-calibrating, since GPS transmitter inter-frequency biases are obtained as a by-product of the mapping procedure. This contrasts with other techniques in which bias values must be provided from some additional source.

Wide-Area↗
DOE OSTI2013

VALIDATION AND RESULTS OF A PSEUDO-MULTI-ZONE COMBUSTION TRAJECTORY PREDICTION MODEL FOR CAPTURING SOOT AND NOX FORMATION ON A MEDIUM DUTY DIESEL ENGINE

A pseudo-multi-zone phenomenological model has been created with the ultimate goal of supporting efforts to enable broader commercialization of low temperature combustion modes in diesel engines. The benefits of low temperature combustion are the simultaneous reduction in soot and nitric oxide emissions and increased engine efficiency if combustion is properly controlled. Determining what qualifies as low temperature combustion for any given engine can be difficult without expensive emissions analysis equipment. This determination can be made off-line using computer models or through factory calibration procedures. This process could potentially be simplified if a real-time prediction model could be implemented to run for any engine platform this is the motivation for this study. The major benefit of this model is the ability for it to predict the combustion trajectory, i.e. local temperature and equivalence ratio in the burning zones. The model successfully captures all the expected trends based on the experimental data and even highlights an opportunity for simply using the average reaction temperature and equivalence ratio as an indicator of emissions levels alone - without solving formation sub-models. This general type of modeling effort is not new, but a major effort was made to minimize the calculation duration to enable implementation as an input to real-time next-cycle engine controller Instead of simply using the predicted engine out soot and NOx levels, control decisions could be made based on the trajectory. This has the potential to save large amounts of calibration time because with minor tuning (the model has only one automatically determined constant) it is hoped that the control algorithm would be generally applicable.

Bittle, Joshua A.↗
DOE OSTI2016

SU-G-IeP3-05: Effects of Image Receptor Technology and Dose Reduction Software On Radiation Dose Estimates for Fluoroscopically-Guided Interventional (FGI) Procedures

Purpose: To investigate the effects of image receptor technology and dose reduction software on radiation dose estimates for most frequently performed fluoroscopically-guided interventional (FGI) procedures at a tertiary health care center. Methods: IRB approval was obtained for retrospective analysis of FGI procedures performed in the interventional radiology suites between January-2011 and December-2015. This included procedures performed using image-intensifier (II) based systems which were subsequently replaced, flat-panel-detector (FPD) based systems which were later upgraded with ClarityIQ dose reduction software (Philips Healthcare) and relatively new FPD system already equipped with ClarityIQ. Post procedure, technologists entered system-reported cumulative air kerma (CAK) and kerma-area product (KAP; only KAP for II based systems) in RIS; these values were analyzed. Data pre-processing included correcting typographical errors and cross-verifying CAK and KAP. The most frequent high and low dose FGI procedures were identified and corresponding CAK and KAP values were compared. Results: Out of 27,251 procedures within this time period, most frequent high and low dose procedures were chemo/immuno-embolization (n=1967) and abscess drainage (n=1821). Mean KAP for embolization and abscess drainage procedures were 260,657, 310,304 and 94,908 mGycm{sup 2}, and 14,497, 15,040 and 6307 mGycm{sup 2} using II-, FPD- and FPD with ClarityIQ- based systems, respectively. Statistically significant differences were observed in KAP values for embolization procedures with respect to different systems but for abscess drainage procedures significant differences were only noted between systems with FPD and FPD with ClarityIQ (p<0.05). Mean CAK reduced significantly from 823 to 308 mGy and from 43 to 21 mGy for embolization and abscess drainage procedures, respectively, in transitioning to FPD systems with ClarityIQ (p<0.05). Conclusion: While transitioning from II- to FPD- based systems was not associated with dose reduction for the most frequently performed FGI procedures, substantial dose reduction was noted with relatively newer systems and dose reduction software.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Computational Modeling of Fracture Toughness of Al-Si and Al-Zn-Mg-Cu Alloys with Detected Porosity

This article presents the results of an experimental and numerical analysis of the fracture resistance of cast aluminum alloys subjected to variable loads over time based on the fracture mechanics method. The subjects of the analysis were the AlSi9 and AlZn6MgCu alloys (usually used for plastic processing) in which the presence of porosities in the form of separate or dispersed voids was revealed experimentally. The material for the experimental procedures was obtained directly from a casting or from casting samples. Depending on the type of the defect detected, different numerical models were proposed. The results of the numerical calculations were verified experimentally. Static tensile tests, fatigue low-cycle tests, fracture mechanics tests, fatigue crack propagation tests, a microscopic examination of the sample fracture area, and an x-ray tomographic analysis were all performed. A correlation of the results between the numerical analysis and the experimental work has been found. On the basis of the proposed calculation models, the fatigue life of castings of the analyzed alloys can be predicted in the case of concentrated and dispersed defects.

36 MATERIALS SCIENCE↗
DOE OSTI2015

3D delivered dose assessment using a 4DCT-based motion model

Purpose: The purpose of this work is to develop a clinically feasible method of calculating actual delivered dose distributions for patients who have significant respiratory motion during the course of stereotactic body radiation therapy (SBRT). Methods: A novel approach was proposed to calculate the actual delivered dose distribution for SBRT lung treatment. This approach can be specified in three steps. (1) At the treatment planning stage, a patient-specific motion model is created from planning 4DCT data. This model assumes that the displacement vector field (DVF) of any respiratory motion deformation can be described as a linear combination of some basis DVFs. (2) During the treatment procedure, 2D time-varying projection images (either kV or MV projections) are acquired, from which time-varying “fluoroscopic” 3D images of the patient are reconstructed using the motion model. The DVF of each timepoint in the time-varying reconstruction is an optimized linear combination of basis DVFs such that the 2D projection of the 3D volume at this timepoint matches the projection image. (3) 3D dose distribution is computed for each timepoint in the set of 3D reconstructed fluoroscopic images, from which the total effective 3D delivered dose is calculated by accumulating deformed dose distributions. This approach was first validated using two modified digital extended cardio-torso (XCAT) phantoms with lung tumors and different respiratory motions. The estimated doses were compared to the dose that would be calculated for routine 4DCT-based planning and to the actual delivered dose that was calculated using “ground truth” XCAT phantoms at all timepoints. The approach was also tested using one set of patient data, which demonstrated the application of our method in a clinical scenario. Results: For the first XCAT phantom that has a mostly regular breathing pattern, the errors in 95% volume dose (D95) are 0.11% and 0.83%, respectively for 3D fluoroscopic images reconstructed from kV and MV projections compared to the ground truth, which is clinically comparable to 4DCT (0.093%). For the second XCAT phantom that has an irregular breathing pattern, the errors are 0.81% and 1.75% for kV and MV reconstructions, both of which are better than that of 4DCT (4.01%). In the case of real patient, although it is impossible to obtain the actual delivered dose, the dose estimation is clinically reasonable and demonstrates differences between 4DCT and MV reconstruction-based dose estimates. Conclusions: With the availability of kV or MV projection images, the proposed approach is able to assess delivered doses for all respiratory phases during treatment. Compared to the planning dose based on 4DCT, the dose estimation using reconstructed 3D fluoroscopic images was as good as 4DCT for regular respiratory pattern and was a better dose estimation for the irregular respiratory pattern.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Effective Ion Mobility Calculations for Macromolecules by Scattering off Electron Clouds

Broad commercialization and increasing resolving power of ion mobility spectrometry/mass spectrometry (IMS/MS) platforms has engendered an explosion of IMS applications to structural characterization of gas-phase biomolecules. That has renewed interest in more accurate and rapid ion mobility calculations needed to elicit ion geometries from the measurements. An approach based on scattering off electron density isosurfaces (SEDI) that mirrors the physics of molecular collisions was proven superior to the common methods involving atomic coordinates a decade ago, but has remained impractical for large ions because of extreme computational demands. Here, we accelerate SEDI by up to ~500 times using the fragment molecular orbital (FMO) approach for surface generation and the multiplexed scattering algorithm in conjunction with the new grid extrapolation procedure for cross section evaluations. Parallelization of the code on a supercomputer has produced major further speed gains, allowing SEDI calculations for proteins (defined by over a million surface points) with the precision of <0.1% in one minute. Initial tests reveal the anticipated dependences of mobility on the ion charge state and lower cross sections in view of reduced surface roughness. Present developments are expected to lead to broad application of SEDI in IMS studies of macromolecules, enabling more accurate and reliable structural assignments.

Alexeev, Yuri↗
DOE OSTI2016

Effects of anisotropic interaction-induced properties of hydrogen-rare gas compounds on rototranslational Raman scattering spectra: Comprehensive theoretical and numerical analysis

A comprehensive study is presented of many aspects of the depolarized anisotropic collision induced (CI) component of light scattered by weakly bound compounds composed of a dihydrogen molecule and a rare gas (Rg) atom, H{sub 2}–Rg. The work continues a series of earlier projects marking the revival of interest in linear light scattering following the development of new highly advanced tools of quantum chemistry and other theoretical, computational, and experimental means of spectral analyses. Sophisticated ab initio computing procedures are applied in order to obtain the anisotropic polarizability component’s dependence on the H{sub 2}–Rg geometry. These data are then used to evaluate the CI spectral lines for all types of Rg atoms ranging from He to Xe (Rn excluded). Evolution of the properties of CI spectra with growing polarizability/masses of the complexes studied is observed. Special attention is given to the heaviest, Kr and Xe based, scatterers. The influence of specific factors shaping the spectral lines (e.g., bound and metastable contribution, potential anisotropy) is discussed. Also the share of pressure broadened allowed rotational transitions in the overall spectral profile is taken into account and the extent to which it is separable from the pure CI contribution is discussed. We finish with a brief comparison between the obtained results and available experimental data.

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

Analysis of high-aspect-ratio jet-flap wings of arbitrary geometry

An analytical technique to compute the performance of an arbitrary jet-flapped wing is developed. The solution technique is based on the method of Maskell and Spence in which the well-known lifting-line approach is coupled with an auxiliary equation providing the extra function needed in jet-flap theory. The present method is generalized to handle straight, uncambered wings of arbitrary planform, twist, and blowing (including unsymmetrical cases). An analytical procedure is developed for continuous variations in the above geometric data with special functions to exactly treat discontinuities in any of the geometric and blowing data. A rational theory for the effect of finite wing thickness is introduced as well as simplified concepts of effective aspect ratio for rapid estimation of performance.

Lissaman, P. B. S.↗
NASA NTRS2006

Analytical Derivation of the Vegetation Optical Depth from the Microwave Polarization Difference Index

A numerical solution for the canopy optical depth in an existing microwave-based land surface parameter retrieval model is presented. The optical depth is derived from the microwave polarization difference index and the dielectric constant of the soil. The original procedure used an approximation in the form of a logarithmic decay function to define this relationship, and was derived through a series of lengthy polynomials. These polynomials had to be recalculated when the scattering albedo or antenna incidence angle changes. The new procedure is computationally more efficient and accurate.

Meesters, Antoon G. C. A.↗
NASA NTRSDate not supplied

F-ANG+: A 3-D Augmented-Stencil Face-Averaged Nodal-Gradient Cell-Centered Finite-Volume Method for Hypersonic Flows

We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.

Weighted Least-Squares↗
DOE OSTI2019

Correlations of the Elemental Compositions of Primary Coal Tar and Char

We report that simulations of large-scale coal combustors rely on accurate submodels to describe the chemical and physical changes in coal during reaction. Typically, researchers use simplified empirical submodels tuned to experimental data to reduce the computational complexity. When data are not readily available, researchers use simplifying assumptions, which can create inaccuracies and biases in a large simulation. One such simplifying assumption in coal research is how to describe the elemental composition of primary pyrolysis products. This paper explores several different empirical model forms to predict the dry, ash-free fractions of C, H, O, N, and S in both the char and the tar, using variables such as parent coal composition, reaction conditions (temperatures and particle residence times), and key coal structural parameters derived from NMR measurements. These model forms were correlated to existing data using a wide range of experimental data using a cross-validation procedure. Since coal structural values can be expensive to measure, several correlations from the literature were used to estimate these values based on information from the proximate and ultimate analyses of the parent coal, including a new correlation for the coal aromaticity. These model forms were tested against a set of measured elemental compositions of tar and char to find the best fit to use in the cross-validation process. Lastly, the best empirical models are presented that predict elemental composition of the coal char and tar after devolatilization.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2015

Variance-based sensitivity indices for stochastic models with correlated inputs

The goal of this article is the formulation of the principles of one of the possible strategies in implementing correlation between input random variables so as to be usable for algorithm development and the evaluation of Sobol’s sensitivity analysis. With regard to the types of stochastic computational models, which are commonly found in structural mechanics, an algorithm was designed for effective use in conjunction with Monte Carlo methods. Sensitivity indices are evaluated for all possible permutations of the decorrelation procedures for input parameters. The evaluation of Sobol’s sensitivity coefficients is illustrated on an example in which a computational model was used for the analysis of the resistance of a steel bar in tension with statistically dependent input geometric characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Vector tomography for reconstructing electric fields with non-zero divergence in bounded domains

In vector tomography (VT), the aim is to reconstruct an unknown multi-dimensional vector field using line integral data. In the case of a 2-dimensional VT, two types of line integral data are usually required. These data correspond to integration of the parallel and perpendicular projection of the vector field along the integration lines and are called the longitudinal and transverse measurements, respectively. In most cases, however, the transverse measurements cannot be physically acquired. Therefore, the VT methods are typically used to reconstruct divergence-free (or source-free) velocity and flow fields that can be reconstructed solely from the longitudinal measurements. In this paper, we show how vector fields with non-zero divergence in a bounded domain can also be reconstructed from the longitudinal measurements without the need of explicitly evaluating the transverse measurements. To the best of our knowledge, VT has not previously been used for this purpose. In particular, we study low-frequency, time-harmonic electric fields generated by dipole sources in convex bounded domains which arise, for example, in electroencephalography (EEG) source imaging. We explain in detail the theoretical background, the derivation of the electric field inverse problem and the numerical approximation of the line integrals. We show that fields with non-zero divergence can be reconstructed from the longitudinal measurements with the help of two sparsity constraints that are constructed from the transverse measurements and the vector Laplace operator. As a comparison to EEG source imaging, we note that VT does not require mathematical modeling of the sources. By numerical simulations, we show that the pattern of the electric field can be correctly estimated using VT and the location of the source activity can be determined accurately from the reconstructed magnitudes of the field. - Highlights: • Vector tomography is used to reconstruct electric fields generated by dipole sources. • Inverse solutions are based on longitudinal and transverse line integral measurements. • Transverse line integral measurements are used as a sparsity constraint. • Numerical procedure to approximate the line integrals is described in detail. • Patterns of the studied electric fields are correctly estimated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1963

Lunar Flight Study Series: Volume 7. Earth-Moon Transit Studies Based on Ephemeris Data and Using Best Available Computer Program: Principles for Reducing Earth-Moon Trajectory Analysis to Fundamentals - 2

This report presents an empirical investigation of earth-to-moon trajectories. The primary objective is to uncover relationships that result in the reduction of mission analysis or trajectory analysis problems to fundamentals. It is considered imperative that the results be accurate to the extent that any relationships that are uncovered are not brought about by simplifying or unrealistic assumptions. Consequently, the model used is as accurate as is presently available for use and the computations are performed under strict error controls. Procedures are reasonable from the engineering viewpoint. Some very helpful principles are uncovered. These may be briefly stated as follows: 1) Trajectories of constant flight time to the moon (in the 66 hour to 90 hour range), and arriving within a few minutes of the same time pass over a common point on the moon regardless of the arrival altitude and inclination, This common point of passage for such a family of trajectories is arbitrarily referred to as a VERTEX point; 2) These VERTEX points lie within a bounded region on the lunar surface as arrival time varies throughout the years, for a given flight time or flight time spread; 3) Increasing flight time results primarily in a longitude shift of the VERTEX point for a given arrival time; 4) Launch azimuth should be treated as a strong parameter regarding its influence on the vertex. Launch azimuth variations up to +/- 20 degrees about 90 degrees can shift the vertex as much as 3 degrees in latitude even when one is free to pick launch time, coast time, and S-IVB burn time appropriately with launch azimuth; and 5) The locus of periselenum, as arrival inclination takes on all possible values, is nearly circular about the vertex point for a given arrival altitude. A publication is now in preparation that applies these principles to much advantage in the solution of various trajectory analysis problems.

LUNAR FLIGHT↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Discrete Hahn polynomials for numerical solution of two-dimensional variable-order fractional Rayleigh–Stokes problem

The main aim of this paper is to find the numerical solutions of 2D Rayleigh–Stokes problem with the variable-order fractional derivatives in the Riemann–Liouville sense. The presented method is based on collocation procedure in combination with the new operational matrix of the variable-order fractional derivatives, in the Caputo sense, for the discrete Hahn polynomials. The main advantage of the proposed method is obtaining a global approximation for spatial and temporal discretizations, and it reduced the problem to an algebraic system, which is easier to solve. Also, the profit of approximating a continuous function by Hahn polynomials is that for computing the coefficients of the expansion, we only have to compute a summation and the calculation of coefficients is exact. The error bound for the approximate solution is estimated. Finally, we evaluate results of the presented method with other numerical methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

The Mount Wilson solar magnetograph - Scanning and data system

The paper describes a computer-operated image-scanning and data-collection system for the magnetograph at the Mt. Wilson 150-foot Tower telescope. The system is based on a minicomputer with a 32K word core memory and a generalized interface unit for controlling image motion, a keyboard, and an associated television screen. Operation of the solar image guider and the data-collection assembly is outlined along with the observation and data-reduction procedures. Advantages of the system include the ability to move the image in almost any conceivable fashion, a wide choice of integration times, and increased accuracy in magnetic and Doppler calibrations as well as in setting of the magnetic zero level.

Howard, R.↗