Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis functions”

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 91 records · Page 5

Higher Order Modeling in Hybrid Approaches to the Computation of Electromagnetic Fields

Higher order geometry representations and interpolatory basis functions for computational electromagnetics are reviewed. Two types of vector-valued basis functions are described: curl-conforming bases, used primarily in finite element solutions, and divergence-conforming bases used primarily in integral equation formulations. Both sets satisfy Nedelec constraints, which optimally reduce the number of degrees of freedom required for a given order. Results are presented illustrating the improved accuracy and convergence properties of higher order representations for hybrid integral equation and finite element methods.

Wilton, Donald R.↗

Test functions for the moment method which yield the minimum mean square error

The method of moments has found extensive applications in the solution of a wide class of electromagnetic problems. Additionally, there are a number of basis functions which have found favor and with them different types of testing methods. This work examines the selection of test functions, given a selected set of basis functions, which serve to minimize the mean square error. These test functions will guarantee equal or lesser error compared with any other test function. Additionally, these test functions provide a monotonic improvement of the approximation as the increase in the order of the system. Hence while these test functions require the application of the differential or integral operator one additional time, they provide an absolute lower bound to the mean square error in the approximation and a means to systematically improve accuracy with each increase in the order of the unknowns.

Shamansky, H.↗

User's manual for CBS3DS, version 1.0

CBS3DS is a computer code written in FORTRAN 77 to compute the backscattering radar cross section of cavity backed apertures in infinite ground plane and slots in thick infinite ground plane. CBS3DS implements the hybrid Finite Element Method (FEM) and Method of Moments (MoM) techniques. This code uses the tetrahedral elements, with vector edge basis functions for FEM in the volume of the cavity/slot and the triangular elements with the basis functions for MoM at the apertures. By virtue of FEM, this code can handle any arbitrarily shaped three-dimensional cavities filled with inhomogeneous lossy materials; due to MoM, the apertures can be of any arbitrary shape. The User's Manual is written to make the user acquainted with the operation of the code. The user is assumed to be familiar with the FORTRAN 77 language and the operating environment of the computer the code is intended to run.

Reddy, C. J.↗

User's Manual for FEMOM3DS

FEMOM3DS is a computer code written in FORTRAN 77 to compute electromagnetic(EM) scattering characteristics of a three dimensional object with complex materials using combined Finite Element Method (FEM)/Method of Moments (MoM) technique. This code uses the tetrahedral elements, with vector edge basis functions for FEM in the volume of the cavity and the triangular elements with the basis functions similar to that described for MoM at the outer boundary. By virtue of FEM, this code can handle any arbitrarily shaped three-dimensional cavities filled with inhomogeneous lossy materials. The User's Manual is written to make the user acquainted with the operation of the code. The user is assumed to be familiar with the FORTRAN 77 language and the operating environment of the computers on which the code is intended to run.

Reddy, C.J.↗

A function approximation approach to anomaly detection in propulsion system test data

Ground test data from propulsion systems such as the Space Shuttle Main Engine (SSME) can be automatically screened for anomalies by a neural network. The neural network screens data after being trained with nominal data only. Given the values of 14 measurements reflecting external influences on the SSME at a given time, the neural network predicts the expected nominal value of a desired engine parameter at that time. We compared the ability of three different function-approximation techniques to perform this nominal value prediction: a novel neural network architecture based on Gaussian bar basis functions, a conventional back propagation neural network, and linear regression. These three techniques were tested with real data from six SSME ground tests containing two anomalies. The basis function network trained more rapidly than back propagation. It yielded nominal predictions with, a tight enough confidence interval to distinguish anomalous deviations from the nominal fluctuations in an engine parameter. Since the function-approximation approach requires nominal training data only, it is capable of detecting unknown classes of anomalies for which training data is not available.

Whitehead, Bruce A.↗

A diagnostic analysis of the VVP single-doppler retrieval technique

A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Boccippio, Dennis J.↗

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗

A multi-modal solution for the transmission of sound in nonuniform ducts

The method of weighted residuals in the form of a modified Galerkin method with boundary residuals is developed for the study of the transmission of sound in nonuniform ducts carrying a steady compressible flow. In this formulation the steady flow is modeled as essentially one-dimensional but with a kinematic modification to force tangency of the flow at the duct walls. Good agreement has been obtained with known results for transmission and reflection coefficients in hard walled ducts up to near sonic velocities. In ducts with acoustically compliant boundaries good comparisons have been more difficult to achieve except at low Mach numbers. The problem of transmission in a straight, acoustically treated duct with a uniform flow has been formulated and the Galerkin method used with basis functions derived from the case when flow is absent. Results indicate that favorable comparisons with exact computations can be obtained if care is taken in choosing the basis functions.

Eversman, W.↗

Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids

This paper is concerned with two important elements in the high-order accurate spatial discretization of finite volume equations over arbitrary grids. One element is the integration of basis functions over arbitrary domains, which is used in expressing various spatial integrals in terms of discrete unknowns. The other consists of quadrature approximations to those integrals. Only polynomial basis functions applied to polyhedral and polygonal grids are treated here. Non-triangular polygonal faces are subdivided into a union of planar triangular facets, and the resulting triangulated polyhedron is subdivided into a union of tetrahedra. The straight line segment, triangle, and tetrahedron are thus the fundamental shapes that are the building blocks for all integrations and quadrature approximations. Integrals of products up to the fifth order are derived in a unified manner for the three fundamental shapes in terms of the position vectors of vertices. Results are given both in terms of tensor products and products of Cartesian coordinates. The exact polynomial integrals are used to obtain symmetric quadrature approximations of any degree of precision up to five for arbitrary integrals over the three fundamental domains. Using a coordinate-free formulation, simple and rational procedures are developed to derive virtually all quadrature formulas, including some previously unpublished. Four symmetry groups of quadrature points are introduced to derive Gauss formulas, while their limiting forms are used to derive Lobatto formulas. Representative Gauss and Lobatto formulas are tabulated. The relative efficiency of their application to polyhedral and polygonal grids is detailed. The extension to higher degrees of precision is discussed.

Liu, Yen↗

Empirical Modeling of the Storm Time Innermost Magnetosphere Using Van Allen Probes and THEMIS Data: Eastward and Banana Currents

The structure of storm time currents in the inner magnetosphere, including its innermost region inside 4R(sub E), is studied for the first time using a modification of the empirical geomagnetic field model TS07D and new data from Van Allen Probes and Time History of Events and Macroscale Interactions during Substorms missions. It is shown that the model, which uses basis-function expansions instead of ad hoc current modules to approximate the magnetic field, consistently improves its resolution and magnetic field reconstruction with the increase of the number of basis functions and resolves the spatial structure and evolution of the innermost eastward current. This includes a connection between the westward ring current flowing largely at R > or approx. 3R(sub E) and the eastward ring current concentrated at R < or approx. 3R(sub E) resulting in a vortex current pattern. A similar pattern coined 'banana current' was previously inferred from the pressure distributions based on the energetic neutral atom imaging and first-principles ring current simulations. The morphology of the equatorial currents is dependent on storm phase. During the main phase, it is complex, with several asymmetries forming banana currents. Near SYM-H minimum, the banana current is strongest, is localized in the evening-midnight sector, and is more structured compared to the main phase. It then weakens during the recovery phase resulting in the equatorial currents to become mostly azimuthally symmetric.

Stephens, G. K.↗

Proper Orthogonal Decomposition in Optimal Control of Fluids

In this article, we present a reduced order modeling approach suitable for active control of fluid dynamical systems based on proper orthogonal decomposition (POD). The rationale behind the reduced order modeling is that numerical simulation of Navier-Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. We examine the possibility of obtaining reduced order models that reduce computational complexity associated with the Navier-Stokes equations while capturing the essential dynamics by using the POD. The POD allows extraction of certain optimal set of basis functions, perhaps few, from a computational or experimental data-base through an eigenvalue analysis. The solution is then obtained as a linear combination of these optimal set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations. We here use it in active control of fluid flows governed by the Navier-Stokes equations. We show that the resulting reduced order model can be very efficient for the computations of optimization and control problems in unsteady flows. Finally, implementational issues and numerical experiments are presented for simulations and optimal control of fluid flow through channels.

Ravindran, S. S.↗

Fault detection and diagnosis using neural network approaches

Neural networks can be used to detect and identify abnormalities in real-time process data. Two basic approaches can be used, the first based on training networks using data representing both normal and abnormal modes of process behavior, and the second based on statistical characterization of the normal mode only. Given data representative of process faults, radial basis function networks can effectively identify failures. This approach is often limited by the lack of fault data, but can be facilitated by process simulation. The second approach employs elliptical and radial basis function neural networks and other models to learn the statistical distributions of process observables under normal conditions. Analytical models of failure modes can then be applied in combination with the neural network models to identify faults. Special methods can be applied to compensate for sensor failures, to produce real-time estimation of missing or failed sensors based on the correlations codified in the neural network.

Kramer, Mark A.↗

Jig-Shape Optimization of a Low-Boom Supersonic Aircraft

A simple approach for optimizing the jig-shape is proposed in this study. This simple approach is based on an unconstrained optimization problem and applied to a low-boom supersonic aircraft. In this study, the jig-shape optimization is performed using the two-step approach. First, starting design variables are computed using the least-squares surface fitting technique. Next, the jig-shape is further tuned using a numerical optimization procedure based on an in-house object-oriented optimization tool. During the numerical optimization procedure, a design jig-shape is determined by the baseline jig-shape and basis functions. A total of 12 symmetric mode shapes of the cruise-weight configuration, rigid pitch shape, rigid left and right stabilator rotation shapes, and a residual shape are selected as sixteen basis functions. After three optimization runs, the trim shape error distribution is improved, and the maximum trim shape error of 0.9844 inches of the starting configuration becomes 0.00367 inch by the end of the third optimization run.

numerical analysi↗

Reduced-Order Modeling of a Heaving Airfoil

A reduced-order model of a flapping airfoil is developed using Proper Orthogonal Decomposition (POD). The proper basis functions, developed from snapshots of full Navier-Stokes simulations, are used for a Galerkin projection of the governing equations. The resulting coupled, nonlinear ordinary di.erential equations have a low dimension because the first few basis members capture most of the energy of the flow. The reduced-order model is used to simulate heaving motions that are both similar to and different from the motion(s) used to generate the basis functions, and the errors in the model are quantified. Several methods are used to generate mode sets that can be used over a range of heaving parameters, including snapshots from one, two, and multiple Navier-Stokes simulations. As snapshots from additional simulations are added to the decomposition, the mode sets become richer and can simulate a wider range of parameter space, at some computational cost. Whereas the POD method is fully applicable in three dimensions, the simulation technique based on a body-fixed and body-fitted grid suffers large overhead when extended to three dimensions. To reduce the overhead, an embedding technique is discussed which embeds the solid wing into a fixed Cartesian grid. The wing, which can now have multiple pieces and also be flexible, is represented by a distribution of body forces. This distribution is determined to give exactly the flow around a flapping wing.

Haj-Hariri, H.↗

Urn Models and Beta-splines

Some insight into the properties of beta-splines is gained by applying the techniques of urn models. Urn models are used to construct beta-spline basis functions and to derive the basic properties of these blending functions and the corresponding beta-spline curves. Only the simple notion of linear geometric continuity and with the most elementary beta parameter are outlined. Non-linear geometric continuity leads to additional beta parameters and to more complicated basis functions. Whether urn models can give us any insight into these higher order concepts still remains to be investigated.

Goldman, R. N.↗

Using EIGER for Antenna Design and Analysis

EIGER (Electromagnetic Interactions GenERalized) is a frequency-domain electromagnetics software package that is built upon a flexible framework, designed using object-oriented techniques. The analysis methods used include moment method solutions of integral equations, finite element solutions of partial differential equations, and combinations thereof. The framework design permits new analysis techniques (boundary conditions, Green#s functions, etc.) to be added to the software suite with a sensible effort. The code has been designed to execute (in serial or parallel) on a wide variety of platforms from Intel-based PCs and Unix-based workstations. Recently, new potential integration scheme s that avoid singularity extraction techniques have been added for integral equation analysis. These new integration schemes are required for facilitating the use of higher-order elements and basis functions. Higher-order elements are better able to model geometrical curvature using fewer elements than when using linear elements. Higher-order basis functions are beneficial for simulating structures with rapidly varying fields or currents. Results presented here will demonstrate curren t and future capabilities of EIGER with respect to analysis of installed antenna system performance in support of NASA#s mission of exploration. Examples include antenna coupling within an enclosed environment and antenna analysis on electrically large manned space vehicles.

Champagne, Nathan J.↗

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗

Reduced Order Methods for Prediction of Thermal-Acoustic Fatigue

The goal of this investigation is to assess the quality of high-cycle-fatigue life estimation via a reduced order method, for structures undergoing random nonlinear vibrations in a presence of thermal loading. Modal reduction is performed with several different suites of basis functions. After numerically solving the reduced order system equations of motion, the physical displacement time history is obtained by an inverse transformation and stresses are recovered. Stress ranges obtained through the rainflow counting procedure are used in a linear damage accumulation method to yield fatigue estimates. Fatigue life estimates obtained using various basis functions in the reduced order method are compared with those obtained from numerical simulation in physical degrees-of-freedom.

Przekop, A.↗