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,495 records · Page 83

Evaluation of Automated Model Calibration Techniques for Residential Building Energy Simulation

This simulation study adapts and applies the general framework described in BESTEST-EX (Judkoff et al 2010) for self-testing residential building energy model calibration methods. BEopt/DOE-2.2 is used to evaluate four mathematical calibration methods in the context of monthly, daily, and hourly synthetic utility data for a 1960's-era existing home in a cooling-dominated climate. The home's model inputs are assigned probability distributions representing uncertainty ranges, random selections are made from the uncertainty ranges to define "explicit" input values, and synthetic utility billing data are generated using the explicit input values. The four calibration methods evaluated in this study are: an ASHRAE 1051-RP-based approach (Reddy and Maor 2006), a simplified simulated annealing optimization approach, a regression metamodeling optimization approach, and a simple output ratio calibration approach. The calibration methods are evaluated for monthly, daily, and hourly cases; various retrofit measures are applied to the calibrated models and the methods are evaluated based on the accuracy of predicted savings, computational cost, repeatability, automation, and ease of implementation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Nonrad: Computing nonradiative capture coefficients from first principles

Point defects in semiconductor crystals provide a means for carriers to recombine nonradiatively. This recombination process impacts the performance of devices. We present the Nonrad code that implements the first-principles approach of Alkauskas et al. (2014) [8] for the evaluation of nonradiative capture coefficients based on a quantum-mechanical description of the capture process. An approach for evaluating electron-phonon coupling within the projector augmented wave formalism is presented. We also show that the common procedure of replacing Dirac delta functions with Gaussians can introduce errors into the resulting capture rate, and implement an alternative scheme to properly account for vibrational broadening. Finally, we assess the accuracy of using an analytic approximation to the Sommerfeld parameter by comparing with direct numerical evaluation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine processing of remotely acquired data

The two major branches of remote sensing are based on image orientation and numerical orientation. Numerically oriented systems tend to involve computers for data analysis. In designing an information-gathering system, the sensor as well as the means of analysis must be well mated to the type of system orientation. Attention is given to sensor types as related to system types, the multispectral approach and pattern recognition, the multispectral scanner as data source, an illustrative example, procedural details in the use of pattern recognition, the speed and cost of data processing, the use of spatial information, and data preprocessing steps.

Landgrebe, D. A.↗

On the estimation of the vertical ozone distribution from satellite measurements

A simple relationship has been established between the vertical ozone profile above the ozone maximum and satellite measured radiances for wavelengths shorter than 3000 A. The method is based on the computation of the penetration depths for solar radiation entering the atmosphere at wavelengths employed by the BUV satellite. Suitably defined penetration depths are shown to be linear functions of the normalized intensities for typical conditions, eliminating the requirement for iterative solutions and considerably simplifying the error analysis for the system. For ozone distributions which depart significantly from the exponential form, a noniterative extension of the analysis procedure has been developed and tested through application to simulated intensities.

Thomas, R. W. L.↗

A new method for mapping multidimensional data to lower dimensions

A multispectral mapping method is proposed which is based on the new concept of BEND (Bidimensional Effective Normalised Difference). The method, which involves taking one sample point at a time and finding the interrelationships between its features, is found very economical from the point of view of storage and processing time. It has good dimensionality reduction and clustering properties, and is highly suitable for computer analysis of large amounts of data. The transformed values obtained by this procedure are suitable for either a planar 2-space mapping of geological sample points or for making grayscale and color images of geo-terrains. A few examples are given to justify the efficacy of the proposed procedure.

Gowda, K. C.↗

Finite element solution techniques for large-scale problems in computational fluid dynamics

Element-by-element approximate factorization, implicit-explicit and adaptive implicit-explicit approximation procedures are presented for the finite-element formulations of large-scale fluid dynamics problems. The element-by-element approximation scheme totally eliminates the need for formation, storage and inversion of large global matrices. Implicit-explicit schemes, which are approximations to implicit schemes, substantially reduce the computational burden associated with large global matrices. In the adaptive implicit-explicit scheme, the implicit elements are selected dynamically based on element level stability and accuracy considerations. This scheme provides implicit refinement where it is needed. The methods are applied to various problems governed by the convection-diffusion and incompressible Navier-Stokes equations. In all cases studied, the results obtained are indistinguishable from those obtained by the implicit formulations.

Liou, J.↗

Method for experimental determination of flutter speed by parameter identification

A method for flight flutter testing is proposed which enables one to determine the flutter dynamic pressure from flights flown far below the flutter dynamic pressure. The method is based on the identification of the coefficients of the equations of motion at low dynamic pressures, followed by the solution of these equations to compute the flutter dynamic pressure. The initial results of simulated data reported in the present work indicate that the method can accurately predict the flutter dynamic pressure, as described. If no insurmountable difficulties arise in the implementation of this method, it may significantly improve the procedures for flight flutter testing.

Nissim, E.↗

Method for experimental determination of flutter speed by parameter identification

A method for flight flutter testing is proposed which enables one to determine the flutter dynamic pressure from flights flown far below the flutter dynamic pressure. The method is based on the identification of the coefficients of the equations of motion at low dynamic pressures, followed by the solution of these equations to compute the flutter dynamic pressure. The initial results of simulated data reported in the present work indicate that the method can accurately predict the flutter dynamic pressure, as described. If no insurmountable difficulties arise in the implementation of this method, it may significantly improve the procedures for flight flutter testing.

Nissim, E.↗

Time variant analysis of large scale constrained rotorcraft systems dynamics - An exploitation of IBM-3090 vector-processor's pipe-lining feature

A generalized algorithmic procedure is presented for handling the constraints in transmissions, which are treated as a multibody system of interconnected rigid/flexible bodies. The type of constraints are classified based on the interconnection of the bodies, assuming one or more points of contact to exist between them. The method is explained through flow charts and configuration/interaction tables. A significant increase in speed of execution is achieved by vectorizing the developed code in computationally intensive areas. The study of an example consisting of two meshing disks rotating at high angular velocity is carried out. The dynamic behavior of the constraint forces associated with the generalized coordinates of the system are plotted by selecting various modes. Applications are intended for the study of dynamic and subsequent prediction of constraint forces at the gear teeth contacting points in helicopter transmissions with the aim of improving performance dependability.

Amirouche, F. M. L.↗

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗

On the Use of CAD-Native Predicates and Geometry in Surface Meshing

Several paradigms for accessing CAD geometry during surface meshing for CFD are discussed. File translation, inconsistent geometry engines and non-native point construction are all identified as sources of non-robustness. The paper argues in favor of accessing CAD parts and assemblies in their native format, without translation, and for the use of CAD-native predicates and constructors in surface mesh generation. The discussion also emphasizes the importance of examining the computational requirements for exact evaluation of triangulation predicates during surface meshing. The native approach is demonstrated through an algorithm for the generation of closed manifold surface triangulations from CAD geometry. CAD parts and assemblies are used in their native format, and a part's native geometry engine is accessed through a modeler-independent application programming interface (API). In seeking a robust and fully automated procedure, the algorithm is based on a new physical space manifold triangulation technique specially developed to avoid robustness issues associated with poorly conditioned mappings. In addition, this approach avoids the usual ambiguities associated with floating-point predicate evaluation on constructed coordinate geometry in a mapped space. The technique is incremental, so that each new site improves the triangulation by some well defined quality measure. The algorithm terminates after achieving a prespecified measure of mesh quality and produces a triangulation such that no angle is less than a given angle bound, a or greater than pi - 2alpha. This result also sets bounds on the maximum vertex degree, triangle aspect-ratio and maximum stretching rate for the triangulation. In addition to the output triangulations for a variety of CAD parts, the discussion presents related theoretical results which assert the existence of such an angle bound, and demonstrate that maximum bounds of between 25 deg and 30 deg may be achieved in practice.

Aftosmis, M. J.↗

Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion plus Selected Inversion

We describe how to apply the recently developed pole expansion plus selected inversion (PEpSI) technique to Kohn-Sham density function theory (DFT) electronic structure calculations that are based on atomic orbital discretization. We give analytic expressions for evaluating charge density, total energy, Helmholtz free energy and atomic forces without using the eigenvalues and eigenvectors of the Kohn-Sham Hamiltonian. We also show how to update the chemical potential without using Kohn-Sham eigenvalues. The advantage of using PEpSI is that it has a much lower computational complexity than that associated with the matrix diagonalization procedure. We demonstrate the performance gain by comparing the timing of PEpSI with that of diagonalization on insulating and metallic nanotubes. For these quasi-1D systems, the complexity of PEpSI is linear with respect to the number of atoms. This linear scaling can be observed in our computational experiments when the number of atoms in a nanotube is larger than a few hundreds. Both the wall clock time and the memory requirement of PEpSI is modest. This makes it even possible to perform Kohn-Sham DFT calculations for 10,000-atom nanotubes on a single processor. We also show that the use of PEpSI does not lead to loss of accuracy required in a practical DFT calculation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Practical gust load alleviation and flutter suppression control laws based on a LQG methodology

A modified linear quadratic Gaussian (LQG) synthesis procedure has been used to design low-order robust multiloop controllers for a flexible airplane. The introduction of properly constructed fictitious Gauss-Markov processes in the control loops allowed meeting classical frequency-domain stability criteria using the direct synthesis procedures of modern time-domain control theory. Model reduction was used to simplify the control laws to the point where they could be easily implemented on onboard flight computers. These control laws provided excellent gust load and flutter mode control with good stability margins and compared very favorably to other control laws synthesized by the classical root-locus technique.

Gangsaas, D.↗

Parabolized reduced Navier-Stokes computational techniques

A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Rubin, Stanley G.↗

F-16XL and F-18 High Speed Acoustic Flight Test Databases

This report presents the recorded acoustic data and the computed narrow-band and 1/3-octave band spectra produced by F-18 and F-16XL aircraft in subsonic flight over an acoustic array. Both broadband-shock noise and turbulent mixing noise are observed in the spectra. Radar and c-band tracking systems provided the aircraft position which enabled directivity and smear angles from the aircraft to each microphone to be computed. These angles are based on source emission time and thus give some idea about the directivity of the radiated sound field due to jet noise. A follow-on static test was also conducted where acoustic and engine data were obtained. The acoustic data described in the report has application to community noise analysis, noise source characterization and validation of prediction models. A detailed description of the signal processing procedures is provided. Follow-on static tests of each aircraft were also conducted for which engine data and far-field acoustic data are presented.

Kelly, J. J.↗

Safety Aboard the International Space Station

As with any task that NASA takes on, safety is of utmost importaqce. There are pages of safety codes and procedures that must be followed before any idea can be brought to life. Unfortunately, the International Space Station s (ISS) safety regulations and procedures are based on lg standards rather than on Og. To aide in making this space age home away from home a less hazardous environment, I worked on several projects revolving around the dangers of flammable items in microgravity. The first task I was assigned was to track flames. This involves turning eight millimeter video recordings, of tests run in the five second drop tower, into avi format on the computer. The footage is then compressed and altered so that the flame can be seen more clearly. Using another program called Spotlight, line profiles were used to collect data describing the luminescence of the flame at different points. These raw data are saved as text files and run trough a macro so that a Matlab program can analyze it. By fitting the data to a curve and determining the areas of brightest luminescence, the behavior of the flame can be recorded numerically. After entering the data into a database, researchers can come back later and easily get information on flames resulting from different gas and liquid mixtures in microgravity. I also worked on phase two of the FATE project, which deals with safety aboard the ISS. This phase involves igniting projected droplets and determining how they react with secondary materials. Such simulations represent, on a small scale, the spread of onboard fires due to the effervescence of burning primary materials. I set up existing hardware to operate these experiments and ran tests with it, photographing the results. I also made CAD drawings of the apparatus and the area available on the (SF)2 rig for it to fit into. The experiment will later be performed on the KC-135, and the results gathered will be used to reanalyze current safety standards for the ISS, including the distance of required separation for flammable materials. Additional information is included in the original extended abstract.

Mintz, Shauna M.↗

Expansion of Tabulated Scattering Matrices in Generalized Spherical Functions

An efficient way to solve the vector radiative transfer equation for plane-parallel turbid media is to Fourier-decompose it in azimuth. This methodology is typically based on the analytical computation of the Fourier components of the phase matrix and is predicated on the knowledge of the coefficients appearing in the expansion of the normalized scattering matrix in generalized spherical functions. Quite often the expansion coefficients have to be determined from tabulated values of the scattering matrix obtained from measurements or calculated by solving the Maxwell equations. In such cases one needs an efficient and accurate computer procedure converting a tabulated scattering matrix into the corresponding set of expansion coefficients. This short communication summarizes the theoretical basis of this procedure and serves as the user guide to a simple public-domain FORTRAN program.

scattering matrix↗

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↗