SEARCH · Search NASA

Results for “discretization 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Downdating a time-varying square root information filter

A new method to efficiently downdate an estimate and covariance generated by a discrete time Square Root Information Filter (SRIF) is presented. The method combines the QR factor downdating algorithm of Gill and the decentralized SRIF algorithm of Bierman. Efficient removal of either measurements or a priori information is possible without loss of numerical integrity. Moreover, the method includes features for detecting potential numerical degradation. Performance on a 300 parameter system with 5800 data points shows that the method can be used in real time and hence is a promising tool for interactive data analysis. Additionally, updating a time-varying SRIF filter with either additional measurements or a priori information proceeds analogously.

Muellerschoen, Ronald J.↗

NASA NTRS1975

A subvortex technique for the close approach to a discretized vortex sheet

The close-approach problem associated with flow calculation methods based on vortex-lattice theory was examined numerically using two-dimensional discretized vortex sheets. The analysis first yields a near-field radius of approximately the distance apart of the vortices in the lattice; only within this distance from the sheet are the errors arising from the discretization significant. Various modifications to the discrete vortices are then considered with the objective of reducing the errors. This leads to a near-field model in which a vortex splits into an increasing number of subvortices as it is approached. The subvortices, whose strengths vary linearly from the vortex position, are evenly distributed along an interpolated curve passing through the basic vortices. This subvortex technique can be extended to the three-dimensional case and is efficient because the number of vortices is effectively increased, but only where and when needed.

Maskew, B.↗

NASA NTRS1985

High resolution solutions of the Euler equations for vortex flows

Solutions of the Euler equations are presented for M = 1.5 flow past a 70-degree-swept delta wing. At an angle of attack of 10 degrees, strong leading-edge vortices are produced. Two computational approaches are taken, based upon fully three-dimensional and conical flow theory. Both methods utilize a finite-volume discretization solved by a pseudounsteady multistage scheme. Results from the two approaches are in good agreement. Computations have been done on a 16-million-word CYBER 205 using 196 x 56 x 96 and 128 x 128 cells for the two methods. A sizable data base is generated, and some of the practical aspects of manipulating it are mentioned. The results reveal many interesting physical features of the compressible vortical flow field and also suggest new areas needing research.

Murman, E. M.↗

NASA NTRS1988

Multigrid method for stability problems

The problem of calculating the stability of steady state solutions of differential equations is treated. Leading eigenvalues (i.e., having maximal real part) of large matrices that arise from discretization are to be calculated. An efficient multigrid method for solving these problems is presented. The method begins by obtaining an initial approximation for the dominant subspace on a coarse level using a damped Jacobi relaxation. This proceeds until enough accuracy for the dominant subspace has been obtained. The resulting grid functions are then used as an initial approximation for appropriate eigenvalue problems. These problems are being solved first on coarse levels, followed by refinement until a desired accuracy for the eigenvalues has been achieved. The method employs local relaxation on all levels together with a global change on the coarsest level only, which is designed to separate the different eigenfunctions as well as to update their corresponding eigenvalues. Coarsening is done using the FAS formulation in a non-standard way in which the right hand side of the coarse grid equations involves unknown parameters to be solved for on the coarse grid. This in particular leads to a new multigrid method for calculating the eigenvalues of symmetric problems. Numerical experiments with a model problem demonstrate the effectiveness of the method proposed. Using an FMG algorithm a solution to the level of discretization errors is obtained in just a few work units (less than 10), where a work unit is the work involved in one Jacobi relization on the finest level.

Taasan, Shlomo↗

NASA NTRS1987

Fast and stable recursive algorithms for continuous-time and discrete-time model conversions

Based on the Newton-Raphson method, this paper presents recursive algorithms that are rapidly convergent and more stable for modeling the equivalent continuous-time (discrete-time) model from the available discrete-time (continuous-time) model for a fixed sampling period. The newly developed recursive algorithms relax the constraints imposed upon the existing model conversion algorithms, and, thus, enhance the applications of microprocessors and associated microelectronics to digital control systems. A practical example is presented to demonstrate the effectiveness of the proposed procedures.

Shieh, L. S.↗

NASA NTRS2005

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗

NASA NTRS1977

A finite element analysis for acoustic transmission in nonuniform ducts

A multimodal finite-element analysis for two-dimensional and axisymmetric ducts with flow has been adopted to predict acoustic transmission in nonuniform ducts carrying a nonuniform flow (i.e., turbofan aircraft engine conditions). The analysis is based on a Galerkin method with a finite-element discretization, and multimodal matching at the duct ends. The element mesh and the number of modes required at each end are found to have a strong influence on the analysis. The numerical results compare well with results obtained from other techniques, including the method of weighted residuals.

Astley, R. J.↗

NASA NTRS1978

Global analysis of a buck regulator

Sufficient conditions for global stability of a buck regulator using a discrete control law are found. The method of paired systems and Liapunov functions are used to establish global stability and to study the convergence of the regulator. A heuristic argument is given that the optimal switching curves associated with the paired continuous systems approximate the optimal switching curves of the discrete systems.

Edwards, D. B.↗

NASA NTRS2018

Geometry and Joint Systems for Lattice-Based Reconfigurable Space Structures

We describe analytical methods for the design of the discrete elements of ultralight lattice structures. This modular, building block strategy allows for relatively simple element manufacturing, as well as relatively simple robotic assembly of low mass density structures on orbit, with potential for disassembly and reassembly into highly varying and large structures. This method also results in a structure that is easily navigable by relatively small mobile robots. The geometry of the cell can allow for high packing efficiency to minimize wasted payload volume while maximizing structural performance and constructability. We describe the effect of geometry choices on the final system mechanical properties, manufacturability of the components, and automated robotic constructability of a final system. Geometry choices considered include building block complexity, symmetry of the unit cell, and effects of vertex, edge, and face connectivity of the unit cell. Mechanical properties considered include strength scaling, modulus scaling, and structural performance of the joint, including proof load, shear load, mass, and loading area; as well as validation and verification opportunities. Manufacturability metrics include cost and time, manufacturing method (COTS versus custom), and tolerances required. Automated constructability metrics include local effects of loads imparted to the structure by the robot and assembly complexity, encompassing the ability of the robot to clamp and number of placement motions needed for assembly.

Ochalek, Megan E.↗

NASA NTRS2019

Geometry Systems for Lattice-Based Reconfigurable Space Structures

We describe analytical methods for the design of the discrete elements of ultralight lattice structures. This modular, building block strategy allows for relatively simple element manufacturing, as well as relatively simple robotic assembly of low mass density structures on orbit, with potential for disassembly and reassembly into highly varying and large structures. This method also results in a structure that is easily navigable by relatively small mobile robots. The geometry of the cell can allow for high packing efficiency to minimize wasted payload volume while maximizing structural performance and constructability. We describe the effect of geometry choices on the final system mechanical properties and automated robotic constructability of a final system. Geometric properties considered include number of attachments per voxel, number of attachments per coefficient of volume, and effects of vertex, edge, and face connectivity of the unit cell. Mechanical properties considered include strength scaling, modulus scaling, and packing efficiency of the lattice. Automated constructibility metrics include volume allowance for an end-effector, strut clearance angle for an end-effector, and packing efficiency. These metrics were applied to six lattice unit cell geometries: cube, cuboctahedron, octahedron, octet, rhombic dodecahedron, and truncated octahedron. A case study is presented to determine the most suitable lattice system for a specific set of strength and modulus scaling requirements while optimizing for ease of robotic assembly.

lattice geometry↗

DOE OSTI2013

A Variational Method for Interpolation Between Zone Centers in 2-D Geometry

Radiation-hydrodynamics codes are often based on a zone centered discretization strategy for the hydrodynamics, while transport requires a higher order interpolated treatment of the material state variable to produce a correct solution in optically thick media. In this paper, we use a variational method to generalize a previously developed 1-D interpolation method to 2-D geometries. While our interpolation method was developed to handle the gradient source term for thermal photons in the Difference Formulation, it may also be useful for the generation of an accurate zone centered discretization for diffusion methods. We provide a detailed description of the 2-D interpolation method and how it relates to the 1-D method that was developed earlier. We also discuss how this interpolation strategy can be extended to 3-D, but do not work through the details.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

42 ENGINEERING↗

NASA NTRS1984

On the solution of problems involving impact type loading

The response of beam and plate members to pulse and impact loading is investigated, and numerical methods are analyzed and compared. The equations of motion, the finite difference and finite element methods, the time integration/Runge-Kutta techniques, and material modeling are discussed in detail. It is found that both the finite difference and the finite element methods could be accurately employed to discretize the spatial variation in the displacements. Central differences or a fourth order Runge-Kutta algorithm could be used for the time integration. The total energy of the system would give the stability and accuracy of the solution. Results showed that the finite element method provided better efficiency in obtaining accurate solution than the finite difference method when a scalar processor is used. However, the finite difference method was more efficient on the vector processor. Therefore, the most efficient method of solution depends on the type of computer present for the analysis.

Moyer, E. T., Jr.↗

NASA NTRSDate not supplied

Geometry Systems for Lattice-Based Reconfigurable Space Structures

We describe analytic methods for the design of the discrete elements of ultralight lattice structures. This modular building block strategy allows for relatively simple element manufacturing as well as relatively simple robotic assembly of low mass-density structures on orbit, with potential for disassembly and reassembly into highly varying and large structures. This method also results in a structure that is easily navigable by relatively small, mobile robots. The geometry of the cell can allow for high packing efficiency to minimize wasted payload volume while maximizing structural performance and constructability. We describe the effect of geometry choices on the mechanical properties and automated robotic constructability of a final system. Geometric properties considered include number of attachments per voxel, number of attachments per coefficient of volume, and effects of vertex, edge, and face connectivity of the unit cell. Mechanical properties considered include strength scaling, modulus scaling, and packing efficiency of the lattice. Automated constructibility metrics include volume allowance for an end-effector, strut clearance angle for an end-effector, and packing efficiency. These metrics were applied to six lattice unit cell geometries: cube, cuboctahedron, octahedron, octet, rhombic dodecahedron, and truncated octahedron. A case study is presented to determine the most suitable lattice system for a specific set of strength and modulus scaling requirements while optimizing for ease of robotic assembly.

digital materials↗

DOE OSTI2013

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics

Multigrid methods can solve some classes of elliptic and parabolic equations to accuracy below the truncation error with a work-cost equivalent to a few residual calculations – so-called “textbook” multigrid efficiency. We investigate methods to solve the system of equations that arise in time dependent magnetohydrodynamics (MHD) simulations with textbook multigrid efficiency. We apply multigrid techniques such as geometric interpolation, full approximate storage, Gauss-Seidel smoothers, and defect correction for fully implicit, nonlinear, second-order finite volume discretizations of MHD. We apply these methods to a standard resistive MHD benchmark problem, the GEM reconnection problem, and add a strong magnetic guide field, which is a critical characteristic of magnetically confined fusion plasmas. We show that our multigrid methods can achieve near textbook efficiency on fully implicit resistive MHD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2022

Finite Difference Methods for Turbulence Simulations

The optimal finite difference discretization used in simulations of turbulent flows is influenced by both, the type of the scale resolving simulation (DNS or LES), as well as the flow-physics (hydrodynamic instabilities, shocks, acoustics, etc.) one expects to resolve. Insight into dispersion and dissipation error requirements for some common scale-resolving simulation scenarios help to highlight the issues faced in selecting a scheme.

Finite Difference Methods↗

DOE OSTI2022

A method for generating moving, orthogonal, area preserving polygonal meshes

A new method for generating locally orthogonal polygonal meshes from a set of generator points is presented in which polygon areas are a constraint. The area constraint property is particularly useful for particle methods where moving polygons track a discrete portion of material. Because Voronoi polygon meshes have some very attractive mathematical and numerical properties for numerical computation, a generalization of Voronoi polygon meshes was formulated that enforces a polygon area constraint. Area constrained moving polygonal meshes allow one to develop hybrid particle-mesh numerical methods that display some of the most attractive features of each approach. It is shown that this mesh construction method can continuously reconnect a moving, unstructured polygonal mesh in a pseudo-Lagrangian fashion without change in cell area/volume, and the method's ability to simulate various physical scenarios is shown. Overall, the advantages are identified for incompressible fluid flow calculations, with demonstration cases that include material discontinuities of all three phases of matter and large density jumps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1986

Wake shape and its effects on aerodynamic characteristics

The wake shape under symmetrical flight conditions and its effects on aerodynamic characteristics are examined. In addition, the effect of wake shape in sideslip and discrete vortices such as strake or forebody vortex on lateral characteristics is presented. The present numerical method for airplane configurations, which is based on discretization of the vortex sheet into vortex segments, verified the symmetrical and asymmetrical roll-up process of the trailing vortices. Also, the effect of wing wake on tail planes is calculated. It is concluded that at high lift the assumption of flat wake for longitudinal and lateral-directional characteristics should be reexamined.

Emdad, H.↗