Search NASA⌕ Search

SEARCH · Search NASA

Results for “Problem structuring methods”

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,261 records · Page 70

A recurrence matrix solution for the dynamic response of aircraft in gusts

A systematic procedure developed for the calculation of the structural response of aircraft flying through a gust by use of difference equations in the solution of dynamic problems is first illustrated by means of a simple-damped-oscillator example. A detailed analysis is then given which leads to a recurrence matrix equation for the determination of the response of an airplane in a gust. The method takes into account wing bending and twisting deformations, fuselage deflection, vertical and pitching motion of the airplane, and some tail forces. The method is based on aerodynamic strip theory, but compressibility and three-dimensional aerodynamic effects can be taken into account approximately by means of over-all corrections. Either a sharp-edge gust or a gust of arbitrary shape in the spanwise or flight directions may be treated. In order to aid in the application of the method to any specific case, a suggested computational procedure is included. The possibilities of applying the method to a variety of transient aircraft problems, such as landing, are brought out. A brief review of matrix algebra, covering the extent to which it is used in the analysis, is also included. (author)

Houbolt, John C↗

A multidisciplinary approach to optimization of controlled space structures

A fundamental problem facing controls-structures analysts is a means of determining the trade-offs between structural design parameters and control design parameters in meeting some particular performance criteria. Developing a general optimization-based design methodology integrating the disciplines of structural dynamics and controls is a logical approach. The objective of this study is to develop such a method. Classical design methodology involves three phases. The first is structural optimization, wherein structural member sizes are varied to minimize structural mass, subject to open-loop frequency constraints. The next phase integrates control and structure design with control gains as additional design variables. The final phase is analysis of the 'optimal' integrated design phase considering 'real' actuators and 'standard' member sizes. The control gains could be further optimized for fixed structure, and actuator saturation constraints could be imposed. However, such an approach does not take full advantage of opportunities to tailor the structure and control system design as one system.

Woodard, Stanley E.↗

Computer Modeling of the Earliest Cellular Structures and Functions

In the absence of extinct or extant record of protocells (the earliest ancestors of contemporary cells). the most direct way to test our understanding of the origin of cellular life is to construct laboratory models of protocells. Such efforts are currently underway in the NASA Astrobiology Program. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs for molecules that perform proto-cellular functions. Many of these functions, such as import of nutrients, capture and storage of energy. and response to changes in the environment are carried out by proteins bound to membrane< We will discuss a series of large-scale, molecular-level computer simulations which demonstrate (a) how small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes, (b) how these peptides aggregate to form membrane-spanning structures (eg. channels), and (c) by what mechanisms such aggregates perform essential proto-cellular functions, such as proton transport of protons across cell walls, a key step in cellular bioenergetics. The simulations were performed using the molecular dynamics method, in which Newton's equations of motion for each item in the system are solved iteratively. The problems of interest required simulations on multi-nanosecond time scales, which corresponded to 10(exp 6)-10(exp 8) time steps.

Pohorille, Andrew↗

Overview of the NASA Glenn Flux Reconstruction Based High-Order Unstructured Grid Code

A computational fluid dynamics code based on the flux reconstruction (FR) method is currently being developed at NASA Glenn Research Center to ultimately provide a large- eddy simulation capability that is both accurate and efficient for complex aeropropulsion flows. The FR approach offers a simple and efficient method that is easy to implement and accurate to an arbitrary order on common grid cell geometries. The governing compressible Navier-Stokes equations are discretized in time using various explicit Runge-Kutta schemes, with the default being the 3-stage/3rd-order strong stability preserving scheme. The code is written in modern Fortran (i.e., Fortran 2008) and parallelization is attained through MPI for execution on distributed-memory high-performance computing systems. An h- refinement study of the isentropic Euler vortex problem is able to empirically demonstrate the capability of the FR method to achieve super-accuracy for inviscid flows. Additionally, the code is applied to the Taylor-Green vortex problem, performing numerous implicit large-eddy simulations across a range of grid resolutions and solution orders. The solution found by a pseudo-spectral code is commonly used as a reference solution to this problem, and the FR code is able to reproduce this solution using approximately the same grid resolution. Finally, an examination of the code's performance demonstrates good parallel scaling, as well as an implementation of the FR method with a computational cost/degree- of-freedom/time-step that is essentially independent of the solution order of accuracy for structured geometries.

High-Order Methods↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

Corbetta, M.↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

stochastic calculus↗

Composite Progressive Delamination Analysis of Thin Ply Buckling

The NASA Advanced Composites Project is a multi-year consortium of government, industry and academia to develop and validate technology that will reduce development and certification of advanced composite materials for aircraft structures. One of its tasks is to simulate and predict composite damage progression accurately through high fidelity finite element models. It included a subtask to verify an approach based on the Virtual Crack Closure Technology (VCCT) as well as Floating Node Method (FNM). This paper discusses predicted damage onset and growth in a thin-ply buckling problem. Sensitivity of the simulations to various initial crack length is studied. The compressive load-displacement curves are plotted and the crack growth critical load is predicted. Overall, the crack progression follows simulation benchmarks well. However, the critical-load overshoot phenomenon is observed in both VCCT and FNM, suggesting further investigation may be worth pursuing.

damage progression↗

Soot Formation in Freely-Propagating Laminar Premixed Flames

Soot formation within hydrocarbon-fueled flames is an important unresolved problem of combustion science. Thus, the present study is considering soot formation in freely-propagating laminar premixed flames, exploiting the microgravity environment to simplify measurements at the high-pressure conditions of interest for many practical applications. The findings of the investigation are relevant to reducing emissions of soot and continuum radiation from combustion processes, to improving terrestrial and spacecraft fire safety, and to developing methods of computational combustion, among others. Laminar premixed flames are attractive for studying soot formation because they are simple one-dimensional flows that are computationally tractable for detailed numerical simulations. Nevertheless, studying soot-containing burner-stabilized laminar premixed flames is problematical: spatial resolution and residence times are limited at the pressures of interest for practical applications, flame structure is sensitive to minor burner construction details so that experimental reproducibility is not very good, consistent burner behavior over the lengthy test programs needed to measure soot formation properties is hard to achieve, and burners have poor durability. Fortunately, many of these problems are mitigated for soot-containing, freely-propagating laminar premixed flames. The present investigation seeks to extend work in this laboratory for various soot processes in flames by observing soot formation in freely-propagating laminar premixed flames. Measurements are being made at both Normal Gravity (NG) and MicroGravity (MG), using a short-drop free-fall facility to provide MG conditions.

Lin, K.-C.↗

Structured background grids for generation of unstructured grids by advancing front method

A new method of background grid construction is introduced for generation of unstructured tetrahedral grids using the advancing-front technique. Unlike the conventional triangular/tetrahedral background grids which are difficult to construct and usually inadequate in performance, the new method exploits the simplicity of uniform Cartesian meshes and provides grids of better quality. The approach is analogous to solving a steady-state heat conduction problem with discrete heat sources. The spacing parameters of grid points are distributed over the nodes of a Cartesian background grid by interpolating from a few prescribed sources and solving a Poisson equation. To increase the control over the grid point distribution, a directional clustering approach is used. The new method is convenient to use and provides better grid quality and flexibility. Sample results are presented to demonstrate the power of the method.

Pirzadeh, Shahyar↗

Shock diffraction computations over complex structures

This work contains the results of a study aimed at the development of two- and three-dimensional numerical procedures for computing the flowfield generated by the interaction of a blast wave and a rigid body. A number of numerical procedures were applied to two-dimensional problems including both implicit and explicit algorithms. Each was tried on the blast wave-cylinder interaction problem. MacCormack's (1969) method with added fourth-order dissipation yielded the best results and was then applied to the blast wave-truck interaction problems in two dimensions. MacCormack's method was also used in three dimensions to determine the flowfield that results when a blast wave strikes a rectangular parallelepiped at an arbitrary angle. Both the twoand three-dimensional computations were compared with experiments in a number of ways. Two dimensional density contours show qualitative agreement for shock front location and Mach stem formation with spark shadowgraphs taken in a shock tube. Pressure-time histories indicate good quantitative agreement between theory and experiment both in two- and three-dimensions.

Mark, A.↗

Applications of Parallel Computation in Micro-Mechanics and Finite Element Method

This project discusses the application of parallel computations related with respect to material analyses. Briefly speaking, we analyze some kind of material by elements computations. We call an element a cell here. A cell is divided into a number of subelements called subcells and all subcells in a cell have the identical structure. The detailed structure will be given later in this paper. It is obvious that the problem is "well-structured". SIMD machine would be a better choice. In this paper we try to look into the potentials of SIMD machine in dealing with finite element computation by developing appropriate algorithms on MasPar, a SIMD parallel machine. In section 2, the architecture of MasPar will be discussed. A brief review of the parallel programming language MPL also is given in that section. In section 3, some general parallel algorithms which might be useful to the project will be proposed. And, combining with the algorithms, some features of MPL will be discussed in more detail. In section 4, the computational structure of cell/subcell model will be given. The idea of designing the parallel algorithm for the model will be demonstrated. Finally in section 5, a summary will be given.

Tan, Hui-Qian↗

Study of the modifications needed for efficient operation of NASTRAN on the Control Data Corporation STAR-100 computer

NASA structural analysis (NASTRAN) computer program is operational on three series of third generation computers. The problem and difficulties involved in adapting NASTRAN to a fourth generation computer, namely, the Control Data STAR-100, are discussed. The salient features which distinguish Control Data STAR-100 from third generation computers are hardware vector processing capability and virtual memory. A feasible method is presented for transferring NASTRAN to Control Data STAR-100 system while retaining much of the machine-independent code. Basic matrix operations are noted for optimization for vector processing.

Source record↗

Alternating direction adaptive grid generation

The present investigation is concerned with the development of an alternating direction method which can adaptively resolve numerical solutions involving physical problems by moving the points of a coordinate grid. The alternating direction movement algorithm is applied to adapt a grid to surface data and weights which are given with respect to the grid itself. Attention is given to details concerning the alternating direction method, aspects of algorithmic construction, the choice of weights which controls the coordinate movement and structure, and questions regarding the applications of the alternating direction movement algorithm. On an analytical basis, the method is seen to move points in a nonsingular fashion.

Eiseman, P. R.↗

The application of contraction theory to an iterative formulation of electromagnetic scattering

Contraction theory is applied to an iterative formulation of electromagnetic scattering from periodic structures and a computational method for insuring convergence is developed. A short history of spectral (or k-space) formulation is presented with an emphasis on application to periodic surfaces. To insure a convergent solution of the iterative equation, a process called the contraction corrector method is developed. Convergence properties of previously presented iterative solutions to one-dimensional problems are examined utilizing contraction theory and the general conditions for achieving a convergent solution are explored. The contraction corrector method is then applied to several scattering problems including an infinite grating of thin wires with the solution data compared to previous works.

Brand, J. C.↗

Solvability conditions for dendritic growth in the boundary-layer model with capillary anisotropy

This paper is concerned primarily with the development of an analytic approach to the theory of steady-state velocity selection in the boundary-layer model of dendritic solidification. The two-dimensional version of this model with a fourfold crystalline anisotropy alpha in the surface tension is considered. By extending a WKB method introduced in an earlier paper, the alpha dependence of the selected growth rate is determined in the limit of small alpha; and this rate is studied for large alphas in the limit in which the dimensionless undercooling approaches unity. Portions of the paper are devoted to a reinterpretation of the mathematical structure of the solvability condition in problems of this kind.

Langer, J. S.↗

Improving Search Properties in Genetic Programming

With the advancing computer processing capabilities, practical computer applications are mostly limited by the amount of human programming required to accomplish a specific task. This necessary human participation creates many problems, such as dramatically increased cost. To alleviate the problem, computers must become more autonomous. In other words, computers must be capable to program/reprogram themselves to adapt to changing environments/tasks/demands/domains. Evolutionary computation offers potential means, but it must be advanced beyond its current practical limitations. Evolutionary algorithms model nature. They maintain a population of structures representing potential solutions to the problem at hand. These structures undergo a simulated evolution by means of mutation, crossover, and a Darwinian selective pressure. Genetic programming (GP) is the most promising example of an evolutionary algorithm. In GP, the structures that evolve are trees, which is a dramatic departure from previously used representations such as strings in genetic algorithms. The space of potential trees is defined by means of their elements: functions, which label internal nodes, and terminals, which label leaves. By attaching semantic interpretation to those elements, trees can be interpreted as computer programs (given an interpreter), evolved architectures, etc. JSC has begun exploring GP as a potential tool for its long-term project on evolving dextrous robotic capabilities. Last year we identified representation redundancies as the primary source of inefficiency in GP. Subsequently, we proposed a method to use problem constraints to reduce those redundancies, effectively reducing GP complexity. This method was implemented afterwards at the University of Missouri. This summer, we have evaluated the payoff from using problem constraints to reduce search complexity on two classes of problems: learning boolean functions and solving the forward kinematics problem. We have also developed and implemented methods to use additional problem heuristics to fine-tune the searchable space, and to use typing information to further reduce the search space. Additional improvements have been proposed, but they are yet to be explored and implemented.

Janikow, Cezary Z.↗

Limb-darkening and the structure of the Jovian atmosphere

By observing the transit of various cloud features across the Jovian disk, limb-darkening curves were constructed for three regions in the 4.6 to 5.1 mu cm band. Several models currently employed in describing the radiative or dynamical properties of planetary atmospheres are here examined to understand their implications for limb-darkening. The statistical problem of fitting these models to the observed data is reviewed and methods for applying multiple regression analysis are discussed. Analysis of variance techniques are introduced to test the viability of a given physical process as a cause of the observed limb-darkening.

Newman, W. I.↗

An improved panel method for the solution of three-dimensional leading edge vortex flows Volume 2: User's guide and programmer's document

A computer program developed for solving the subsonic, three dimensional flow over wing-body configurations with leading edge vortex separation is presented. Instructions are given for the proper set up and input of a problem into the computer code. Program input formats and output are described, as well as the overlay structure of the program. The program is written in FORTRAN.

Tinoco, E. N.↗