Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite expression method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Multigrid techniques for the numerical solution of the diffusion equation

An accurate numerical solution of diffusion problems containing large local gradients can be obtained with a significant reduction in computational time by using a multigrid computational scheme. The spatial domain is covered with sets of uniform square grids of different sizes. The finer grid patterns overlap the coarse grid patterns. The finite-difference expressions for each grid pattern are solved independently by iterative techniques. Two interpolation methods were used to establish the values of the potential function on the fine grid boundaries with information obtained from the coarse grid solution. The accuracy and computational requirements for solving a test problem by a simple multigrid and a multilevel-multigrid method were compared. The multilevel-multigrid method combined with a Taylor series interpolation scheme was found to be best.

Phillips, R. E.↗
NASA NTRS1989

Determination of aerodynamic sensitivity coefficients in the transonic and supersonic regimes

The quasi-analytical approach is developed to compute airfoil aerodynamic sensitivity coefficients in the transonic and supersonic flight regimes. Initial investigation verifies the feasibility of this approach as applied to the transonic small perturbation residual expression. Results are compared to those obtained by the direct (finite difference) approach and both methods are evaluated to determine their computational accuracies and efficiencies. The quasi-analytical approach is shown to be superior and worth further investigation.

Elbanna, Hesham M.↗
NASA NTRS2001

Direct Calculation of the Scattering Amplitude Without Partial Wave Analysis

Two new developments in scattering theory are reported. We show, in a practical way, how one can calculate the full scattering amplitude without invoking a partial wave expansion. First, the integral expression for the scattering amplitude f(theta) is simplified by an analytic integration over the azimuthal angle. Second, the full scattering wavefunction which appears in the integral expression for f(theta) is obtained by solving the Schrodinger equation with the finite element method (FEM). As an example, we calculate electron scattering from the Hartree potential. With minimal computational effort, we obtain accurate and stable results for the scattering amplitude.

Shertzer, J.↗
NASA NTRS2002

Calculation of Scattering Amplitude Without Partial Analysis: Inclusion of Exchange - II

There was a method for calculating the whole scattering amplitude, f(Omega(sub k)), directly. The idea was to calculate the complete wave function Psi numerically, and use it in an integral expression for f, which can be reduced to a 2 dimensional quadrature. The original application was for e-H scattering without exchange. There the Schrodinger reduces a 2-d partial differential equation (pde), which was solved using the finite element method (FEM). Here we extend the method to the exchange approximation. The S.E. can be reduced to a pair of coupled pde's, which are again solved by the FEM. The formal expression for f(Omega(sub k)) consists two integrals, f+/- = f(sub d) +/- f(sub e); f(sub d) is formally the same integral as the no-exchange f. We have also succeeded in reducing f(sub e) to a 2-d integral. Results will be presented at the meeting.

Temkin, Aaron↗
NASA NTRS1992

Axially symmetric pulse propagation in semi-infinite hollow cylinders

The effects of dispersion and end conditions on axially symmetric pulse propagation in semiinfinite tubes are investigated. Numerical modeling of the dynamic response is accomplished using a stable and efficient finite element and wave propagation based method. This method yields the dispersive and modal data required to express displacements in a frequency domain modal expansion. Boundary conditions are applied to the tube end to obtain mode amplitudes. A fast Fourier transform is then used to get the time response. The specific end conditions considered are intended to give some insight into how the tube would behave as a member of a jointed truss structure.

Kohl, T.↗
DOE OSTI2023

Bayesian filtering for model predictive control of stochastic gene expression in single cells

This study describes a method for controlling the production of protein in individual cells using stochastic models of gene expression. By combining modern microscopy platforms with optogenetic gene expression, experimentalists are able to accurately apply light to individual cells, which can induce protein production. Here we use a finite state projection based stochastic model of gene expression, along with Bayesian state estimation to control protein copy numbers within individual cells. We compare this method to previous methods that use population based approaches. We also demonstrate the ability of this control strategy to ameliorate discrepancies between the predictions of a deterministic model and stochastic switching system.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1972

Asymptotic expansion of an impulse for an optimal finite burn

Closed form expressions are derived for the position and velocity of a spacecraft during a finite burn using the method involving the theory of asymptotic expansion of the optimal impulsive solution. The small parameter is given by the reciprocal of the mass flow rate. The expansion is given in terms of the impulsive solution and holds up through third-order for the velocity and fourth-order for the position.

Pines, S.↗
NASA NTRS1980

Grid generation using coarse, smooth finite elements

Grid quality solutions of a wide variety of equations can be generated to fit reasonable functional boundary conditions in two dimensions using very coarse rectangular finite elements. Some wavy operators (with no natural variational expression) were tried in order to demonstrate the method's versatility.

Dickson, L. J.↗
NASA NTRS1992

Sensitivity analysis of dynamic aeroelastic responses

This paper summarizes ongoing research on the sensitivity analysis of dynamic aeroelastic response of wings. Two approaches are being used to express the unsteady aerodynamic loads: (1) the frequency-domain approach, and (2) the state-space approach. The frequency-domain approach is demonstrated on a three-dimensional box wing and the state-space domain approach is demonstrated on a simple two-dimensional sectional model. Three different methods are used to find the sensitivities: (1) a purely finite difference approach, (2) a semi-analytical approach in which an analytical expression is used for calculating the sensitivity of an eigenvalue of the complex valued aeroelastic matrix, however, the derivatives of the components of the matrix are obtained using finite difference, and (3) a semi-analytic approach that differs from (2) in the sense that the sensitivity of the aerodynamic matrix is now obtained analytically. A good agreement is seen between the three sets of results. For the two-dimensional sectional model the results for the sensitivities of the flutter speed with respect to the various parameters are obtained.

Kapania, Rakesh K.↗
NASA NTRS1994

Using Finite-Element Analysis In Estimating Reliability

Method of estimating design survivability of structural component incorporates finite-element and probabilistic properties of materials. Involves evaluation of design parameters through direct comparisons of survivability of component expressed in terms of percentages of like components that survive at various lifetimes. Probabilistic properties of materials, given in terms of Weibull parameters, coupled with stress field computed by finite-element analysis to determine fatigue life based on initiation of cracks. Method applied to rotating disk containing bolt holes, representative of disks used in aerospace propulsion turbines. Also used in early stages of design process to optimize life-based designs, reducing testing of full-sized components needed to validate designs.

Zaretsky, Erwin V.↗
DOE OSTI2019

Diffraction of waves by a perfect electromagnetic conductor half-plane between isorefractive media

In this study, diffraction of waves by a perfectly electromagnetic conductor half-plane which is lying at the boundary of isorefractive media is investigated. The diffracted field expressions are constructed by using the relation between scattered geometrical optics fields and the transition boundaries. The high frequency asymptotic expressions of diffracted fields reach infinite magnitude values at the transition boundaries. The uniform expressions of fields are obtained by using the uniform theory of diffraction method in order to obtain the finite magnitude values at the transition boundaries and the results are analyzed numerically for different set of parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Mixed finite-difference scheme for free vibration analysis of noncircular cylinders

A mixed finite-difference scheme is presented for the free-vibration analysis of simply supported closed noncircular cylindrical shells. The problem is formulated in terms of eight first-order differential equations in the circumferential coordinate which possess a symmetric coefficient matrix and are free of the derivatives of the elastic and geometric characteristics of the shell. In the finite-difference discretization, two interlacing grids are used for the different fundamental unknowns in such a way as to avoid averaging in the difference-quotient expressions used for the first derivative. The resulting finite-difference equations are symmetric. The inverse-power method is used for obtaining the eigenvalues and eigenvectors.

Noor, A. K.↗
DOE OSTI2018

Exact excited states of nonintegrable models

We discuss a method of numerically identifying exact energy eigenstates for a finite system, whose form can then be obtained analytically. We demonstrate our method by identifying and deriving exact analytic expressions for several excited states, including an infinite tower, of the one-dimensional spin-1 Affleck-Kennedy-Lieb-Tasaki (AKLT) model, a celebrated nonintegrable model. The states thus obtained for the AKLT model can be interpreted as from one to an extensive number of quasiparticles on the ground state or on the highest excited state when written in terms of dimers. Included in these exact states is a tower of states spanning energies from the ground state to the highest excited state. Here, some of the states of the tower appear to be in the bulk of the energy spectrum, allowing us to make conjectures on the strong eigenstate thermalization hypothesis. We also generalize these exact states including the tower of states to the generalized integer spin AKLT models. Furthermore, we establish a correspondence between some of our states and those of the Majumdar-Ghosh model, yet another nonintegrable model, and extend our construction to the generalized integer spin AKLT models.

1-dimensional spin chains↗
DOE OSTI2019

Fluid–Solid Reaction Kinetics for Solids of Nonbasic Geometries and Determination of the Appropriate Shape Factors

The rate analysis involving a solid of nonbasic geometry in general requires a complex numerical solution. A closed-form expression of the reaction rate for such a solid would be very convenient and greatly simplify the analysis. A method for this purpose was developed in this work for solids with shapes of finite cylinders and finite rectangular parallelepipeds reacting under the rate control by chemical reaction or mass transfer, which may be extended to other irregular shapes. In this method, an appropriate shape factor is first obtained under the conditions of a topochemical reaction, which is then applied to the unified rate expressions previously obtained for solids of three basic shapes. Good agreements between the rates predicted with the appropriate shape factors and the exact or numerical solutions were obtained for both reaction-controlled and pore-diffusion-controlled reactions.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗
NASA NTRS1991

Incorporating finite element analysis into component life and reliability

A method for calculating a component's design survivability by incorporating finite element analysis and probabilistic material properties was developed. The method evaluates design parameters through direct comparisons of component survivability expressed in terms of Weibull parameters. The analysis was applied to a rotating disk with mounting bolt holes. The highest probability of failure occurred at, or near, the maximum shear stress region of the bolt holes. Distribution of failure as a function of Weibull slope affects the probability of survival. Where Weibull parameters are unknown for a rotating disk, it may be permissible to assume Weibull parameters, as well as the stress-life exponent, in order to determine the disk speed where the probability of survival is highest.

August, Richard↗
NASA NTRS1992

Fatigue criterion to system design, life and reliability: A primer

A method for estimating a component's design survivability by incorporating finite element analysis and probabilistic material properties was developed. The method evaluates design parameters through direct comparisons of component survivability expressed in terms of Weibull parameters. The analysis was applied to a rotating disk with mounting bolt holes. The highest probability of failure occurred at, or near, the maximum shear stress region of the bolt holes. Distribution of material failure as a function of Weibull slope affects the probability of survival. Where Weibull parameters are unknown for a rotating disk, it may be permissible to assume Weibull parameters, as well as the stress-life exponent, in order to determine the qualitative effect of disk speed on the probability of survival.

Zaretsky, Erwin V.↗