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 487 records · Page 27

A robust and efficient stepwise regression method for building sparse polynomial chaos expansions

Polynomial Chaos (PC) expansions are widely used in various engineering fields for quantifying uncertainties arising from uncertain parameters. The computational cost of classical PC solution schemes is unaffordable as the number of deterministic simulations to be calculated grows dramatically with the number of stochastic dimension. This considerably restricts the practical use of PC at the industrial level. A common approach to address such problems is to make use of sparse PC expansions. This paper presents a non-intrusive regression-based method for building sparse PC expansions. The most important PC contributions are detected sequentially through an automatic search procedure. The variable selection criterion is based on efficient tools relevant to probabilistic method. Two benchmark analytical functions are used to validate the proposed algorithm. The computational efficiency of the method is then illustrated by a more realistic CFD application, consisting of the non-deterministic flow around a transonic airfoil subject to geometrical uncertainties. To assess the performance of the developed methodology, a detailed comparison is made with the well established LAR-based selection technique. The results show that the developed sparse regression technique is able to identify the most significant PC contributions describing the problem. Moreover, the most important stochastic features are captured at a reduced computational cost compared to the LAR method. The results also demonstrate the superior robustness of the method by repeating the analyses using random experimental designs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multi-Objective Optimization of a Wrought Magnesium Alloy for High Strength and Ductility

An optimization technique is coupled with crystal plasticity based finite element (CPFE) computations to aid the microstructural design of a wrought magnesium alloy for improved strength and ductility. The initial microstructure consists of a collection of sub-micron sized grains containing deformation twins. The variables used in the simulations are crystallographic texture, and twin spacing within the grains. It is assumed that plastic deformation occurs mainly by dislocation slip on two sets of slip systems classified as hard and soft modes. The hard modes are those slip systems that are inclined to the twin planes and the soft mode consists of dislocation glide along the twin plane. The CPFE code calculates the stress-strain response of the microstructure as a function of the microstructural parameters and the length-scale of the features. A failure criterion based on a critical shear strain and a critical hydrostatic stress is used to define ductility. The optimization is based on the sequential generation of an initial population defined by the texture and twin spacing variables. The CPFE code and the optimizer are coupled in parallel so that new generations are created and analyzed dynamically. In each successive generation, microstructures that satisfy at least 90% of the mean strength and mean ductility in the current generation are retained. Multiple generation runs based on the above procedure are carried out in order to obtain maximum strength-ductility combinations. Finally, the implications of the computations for the design of a wrought magnesium alloy are discussed. Research sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U. S. Department of Energy.

36 MATERIALS SCIENCE↗

Velocity-split Navier-Stokes solution procedure for incompressible high Reynolds number external flows

A method for solving the Navier-Stokes equations based on splitting the velocity vector into its rotational and irrotational parts was sucessfully applied to internal flow computations. The applicability of the method to external flows is examined by studying several model problems. The model problems are those of laminar and turbulent incompressible flow past a semi-infinite flat plate and laminar incompressible flow past a finite flat plate. For these problems, the procedure accurately reproduces the known solutions and is computationally very efficient even at high Reynolds numbers. Computational aspects of the method are discussed along with the possibility of using the procedure to retrofit a viscous capability into existing potential flow codes.

Dwoyer, D. L.↗

Reducing the error of geoid undulation computations by modifying Stokes' function

The truncation theory as it pertains to the calculation of geoid undulations based on Stokes' integral, but from limited gravity data, is reexamined. Specifically, the improved procedures of Molodenskii et al. are shown through numerical investigations to yield substantially smaller errors than the conventional method that is often applied in practice. In this improved method, as well as in a simpler alternative to the conventional approach, the Stokes' kernel is suitably modified in order to accelerate the rate of convergence of the error series. These modified methods, however, effect a reduction in the error only if a set of low-degree potential harmonic coefficients is utilized in the computation. Consider, for example, the situation in which gravity anomalies are given in a cap of radius 10 deg and the GEM 9 (20,20) potential field is used. Then, typically, the error in the computed undulation (aside from the spherical approximation and errors in the gravity anomaly data) according to the conventional truncation theory is 1.09 m; with Meissl's modification it reduces to 0.41m, while Molodenskii's improved method gives 0.45 m. A further alteration of Molodenskii's method is developed and yields an RMS error of 0.33 m. These values reflect the effect of the truncation, as well as the errors in the GEM 9 harmonic coefficients. The considerable improvement, suggested by these results, of the modified methods over the conventional procedure is verified with actual gravity anomaly data in two oceanic regions, where the GEOS-3 altimeter geoid serves as the basis for comparison. The optimal method of truncation, investigated by Colombo, is extremely ill-conditioned. It is shown that with no corresponding regularization, this procedure is inapplicable.

Jekeli, C.↗

A Computed Tomography Automation Architecture Developed from a Comprehensive Literature Review

This paper presents a literature review on techniques related to the computed tomography procedure that incorporate automation elements in their research investigations or industrial applications. Computed tomography (CT) is a non-destructive testing (NDT) technique in that the imaging and inspection are performed without damaging the sample, allowing for additional or repeated analysis if necessary. The reviewed literature is organized based on the steps associated with a general NDT task in order to define an end-to-end computed tomography automation architecture. The process steps include activities prior to image collection, during the scan, and after the data are collected. It further reviews efforts related to repeating this process based on a previous scan result. By analyzing the multiple existing but disparate efforts found in the literature, we present a framework for fully automating NDT procedures and discuss the remaining technical gaps in the developed framework.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

HADY-I, a FORTRAN program for the compressible stability analysis of three-dimensional boundary layers

A computer program HADY-I for calculating the linear incompressible or compressible stability characteristics of the laminar boundary layer on swept and tapered wings is described. The eigenvalue problem and its adjoint arising from the linearized disturbance equations with the appropriate boundary conditions are solved numerically using a combination of Newton-Raphson interative scheme and a variable step size integrator based on the Runge-Kutta-Fehlburh fifth-order formulas. The integrator is used in conjunction with a modified Gram-Schmidt orthonormalization procedure. The computer program HADY-I calculates the growth rates of crossflow or streamwise Tollmien-Schlichting instabilities. It also calculates the group velocities of these disturbances. It is restricted to parallel stability calculations, where the boundary layer (meanflow) is assumed to be parallel. The meanflow solution is an input to the program.

El-Hady, N. M.↗

ENSAERO - A multidisciplinary program for fluid/structural interaction studies of aerospace vehicles

A computational procedure is developed that uses a moving zonal grid concept to model complex flexible aerospace vehicles. The Euler/Navier-Stokes equations are used to model the flow, and computations are made using efficient methods based on both central and upwind schemes. The structure is represented by a finite element method which can model general aerospace vehicles. Provisions are made to accommodate other disciplines such as controls and thermal loads. The code is capable of computing unsteady flows on flexible wings with vortical flows. Adaptation of this procedure for parallel processing and validation for complete aerospace configurations is in progress.

Guruswamy, G. P.↗

Preliminary design of composite wings for buckling, strength and displacement constraints

An unstiffened panel buckling constraint for balanced, symmetric laminated composites is included on the global design level in a mathematical programming structural optimization procedure for designing wing structures. Constraints are introduced by penalty functions, and Newton's method based on approximate second derivatives of the penalty terms is used as the search algorithm to obtain minimum-mass designs. Constraint approximations used during the optimization process contribute to the computational efficiency of the procedure. A criterion is developed that identifies the appropriate conservative form of the constraint approximations that are used with the optimization procedure. Minimum-mass design results are obtained for a multispar high-aspect-ratio wing subjected to material strength, minimum-gage, displacement, panel buckling and twist constraints. The material systems considered for the examples are all graphite-epoxy, graphite-epoxy with boron-epoxy spar caps, and all aluminum. The composite material designs are shown to have an advantage over the aluminum designs since they can often satisfy additional constraints with only small mass increases.

Starnes, J. H., Jr.↗

Nonlinear smoothing identification algorithm with application to data consistency checks

A parameter identification algorithm for nonlinear systems is presented. It is based on smoothing test data with successively improved sets of model parameters. The smoothing, which is iterative, provides all of the information needed to compute the gradients of the smoothing performance measure with respect to the parameters. The parameters are updated using a quasi-Newton procedure, until convergence is achieved. The advantage of this algorithm over standard maximum likelihood identification algorithms is the computational savings in calculating the gradient. This algorithm was used for flight-test data consistency checks based on a nonlinear model of aircraft kinematics. Measurement biases and scale factors were identified. The advantages of the presented algorithm and model are discussed.

Idan, M.↗

Enabling Computational Nanotechnology through JavaGenes in a Cycle Scavenging Environment

A genetic algorithm procedure is developed and implemented for fitting parameters for many-body inter-atomic force field functions for simulating nanotechnology atomistic applications using portable Java on cycle-scavenged heterogeneous workstations. Given a physics based analytic functional form for the force field, correlated parameters in a multi-dimensional environment are typically chosen to fit properties given either by experiments and/or by higher accuracy quantum mechanical simulations. The implementation automates this tedious procedure using an evolutionary computing algorithm operating on hundreds of cycle-scavenged computers. As a proof of concept, we demonstrate the procedure for evaluating the Stillinger-Weber (S-W) potential by (a) reproducing the published parameters for Si using S-W energies in the fitness function, and (b) evolving a "new" set of parameters using semi-empirical tightbinding energies in the fitness function. The "new" parameters are significantly better suited for Si cluster energies and forces as compared to even the published S-W potential.

Globus, Al↗

Static and Dynamic Model Update of an Inflatable/Rigidizable Torus Structure

The present work addresses the development of an experimental and computational procedure for validating finite element models. A torus structure, part of an inflatable/rigidizable Hexapod, is used to demonstrate the approach. Because of fabrication, materials, and geometric uncertainties, a statistical approach combined with optimization is used to modify key model parameters. Static test results are used to update stiffness parameters and dynamic test results are used to update the mass distribution. Updated parameters are computed using gradient and non-gradient based optimization algorithms. Results show significant improvements in model predictions after parameters are updated. Lessons learned in the areas of test procedures, modeling approaches, and uncertainties quantification are presented.

Horta, Lucas G.↗

Prediction of static aerodynamic characteristics for space-shuttle-like and other bodies at angles of attack from 0 deg to 180 deg

An engineering-type procedure is presented for computing normal-force, axial-force, and pitching-moment coefficients for bodies at angles of attack from 0 deg to 180 deg. The procedure is ideally suited for estimating the aerodynamic characteristics of space shuttle booster-like bodies because of the wide range of angles of attack, Mach numbers, and Reynolds numbers to be encountered. The analytical formulas, plots, and references given are also applicable for shuttle orbiter, missile, and aircraft-like bodies of both circular and noncircular cross section. The method for computing normal-force and pitching-moment coefficients is based upon the original proposal of Allen that the crossflow or lift distribution over a body can be expressed as the sum of slender-body potential term and an empirical viscous crossflow term. Although experimental data from which to verify the procedure at very high angles of attack are extremely limited, the comparisons made thus far of computed with experimental results are good. In this report the procedure has been shown to be capable of predicting reasonably well the experimental variation of C sub N, C sub A, C sub m, and x sub ac/l with angle of attack for nine bodies of revolution at a free-stream Mach number of 2.86.

Jorgensen, L. H.↗

Resilience-based Recovery Scheduling of Transportation Network in Mixed Traffic Environment: A Deep-Ensemble-Assisted Active Learning Approach

Devising effective post-hazard recovery strategies is critical in enhancing the resilience of transportation networks (TNs). However, existing work does not consider the multiclass users’ travel behavior in network functionality quantification and the metaheuristic solution procedures often suffer from extensive computational burden due to the exploration need in large solution space and the expensive functionality quantification. Here, this study develops a bilevel decision-making framework for the resilience-based recovery scheduling of the TN in a mixed traffic environment with connected and autonomous vehicles (CAVs) and human-driven vehicles (HDVs). The lower level quantifies the TN's functionality over time considering different travel behavior of CAV and HDV users arisen from their different levels of information perception. The upper level presents a novel deep-ensemble-assisted active learning approach to balance optimization performance and computational cost. This framework can help decision makers better quantify the TN's functionality to support effective recovery scheduling of TN with different mixed traffic scenarios ranging from HDV-only to future CAV-dominant traffic. The optimization approach bears the potential to be extended to solving general large-scale network recovery scheduling problems effectively and efficiently. The proposed methodology is demonstrated using a real-world traffic network in Southern California under earthquake considering deterministic and stochastic repair durations.

42 ENGINEERING↗

Simplified cyclic structural analyses of SSME turbine blades

Anisotropic high-temperature alloys are used to meet the safety and durability requirements of turbine blades for high-pressure turbopumps in reusable space propulsion systems. The applicability to anisotropic components of a simplified inelastic structural analysis procedure developed at the NASA Lewis Research Center is assessed. The procedure uses as input the history of the total strain at the critical crack initiation location computed from elastic finite-element analyses. Cyclic heat transfer and structural analyses are performed for the first stage high-pressure fuel turbopump blade of the space shuttle main engine. The blade alloy is directionally solidified MAR-M 246 (nickel base). The analyses are based on a typical test stand engine cycle. Stress-strain histories for the airfoil critical location are computed using both the MARC nonlinear finite-element computer code and the simplified procedure. Additional cases are analyzed in which the material yield strength is arbitrarily reduced to increase the plastic strains and, therefore, the severity of the problem. Good agreement is shown between the predicted stress-strain solutions from the two methods. The simplified analysis uses about 0.02 percent (5 percent with the required elastic finite-element analyses) of the CPU time used by the nonlinear finite element analysis.

Kaufman, A.↗

Design optimization of transonic airfoils

Numerical optimization procedures were considered for the design of airfoils in transonic flow based on the transonic small disturbance (TSD) and Euler equations. A sequential approximation optimization technique was implemented with an accurate approximation of the wave drag based on the Nixon's coordinate straining approach. A modification of the Euler surface boundary conditions was implemented in order to efficiently compute design sensitivities without remeshing the grid. Two effective design procedures producing converged designs in approximately 10 global iterations were developed: interchanging the role of the objective function and constraint and the direct lift maximization with move limits which were fixed absolute values of the design variables.

Joh, C.-Y.↗

Least-squares/parabolized Navier-Stokes procedure for optimizing hypersonic wind tunnel nozzles

A new procedure is demonstrated for optimizing hypersonic wind-tunnel-nozzle contours. The procedure couples a CFD computer code to an optimization algorithm, and is applied to both conical and contoured hypersonic nozzles for the purpose of determining an optimal set of parameters to describe the surface geometry. A design-objective function is specified based on the deviation from the desired test-section flow-field conditions. The objective function is minimized by optimizing the parameters used to describe the nozzle contour based on the solution to a nonlinear least-squares problem. The effect of the changes in the nozzle wall parameters are evaluated by computing the nozzle flow using the parabolized Navier-Stokes equations. The advantage of the new procedure is that it directly takes into account the displacement effect of the boundary layer on the wall contour. The new procedure provides a method for optimizing hypersonic nozzles of high Mach numbers which have been designed by classical procedures, but are shown to produce poor flow quality due to the large boundary layers present in the test section. The procedure is demonstrated by finding the optimum design parameters for a Mach 10 conical nozzle and a Mach 6 and a Mach 15 contoured nozzle.

Korte, John J.↗