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.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Computation of transonic vortex flows past delta wings Integral equation approach

The steady full-potential equation is written in the form of Poisson's equation, and the solution of the velocity field is expressed in terms of an integral equation. The solution consists of a surface integral of vorticity distribution on the wing and its free-vortex sheets and a volume integral of source distribution within a volume around the wing and its free-vortex sheets. The solution is obtained through successive iteration cycles. The source distribution is computed by using a mixed finite-difference scheme of the Murman-Cole type. The method is applied to delta wings. Numerical examples show that a conical shock is captured on the suction side of the wing. It is attached to the lower surface of the leading-edge vortex but does not necessarily reach to the wing surface.

Kandil, O. A.↗

NASA NTRS1993

Application of symbolic and algebraic manipulation software in solving applied mechanics problems

As its name implies, symbolic and algebraic manipulation is an operational tool which not only can retain symbols throughout computations but also can express results in terms of symbols. This report starts with a history of symbolic and algebraic manipulators and a review of the literatures. With the help of selected examples, the capabilities of symbolic and algebraic manipulators are demonstrated. These applications to problems of applied mechanics are then presented. They are the application of automatic formulation to applied mechanics problems, application to a materially nonlinear problem (rigid-plastic ring compression) by finite element method (FEM) and application to plate problems by FEM. The advantages and difficulties, contributions, education, and perspectives of symbolic and algebraic manipulation are discussed. It is well known that there exist some fundamental difficulties in symbolic and algebraic manipulation, such as internal swelling and mathematical limitation. A remedy for these difficulties is proposed, and the three applications mentioned are solved successfully. For example, the closed from solution of stiffness matrix of four-node isoparametrical quadrilateral element for 2-D elasticity problem was not available before. Due to the work presented, the automatic construction of it becomes feasible. In addition, a new advantage of the application of symbolic and algebraic manipulation found is believed to be crucial in improving the efficiency of program execution in the future. This will substantially shorten the response time of a system. It is very significant for certain systems, such as missile and high speed aircraft systems, in which time plays an important role.

Tsai, Wen-Lang↗

NASA NTRS1994

Parameter dimension of turbulence-induced phase errors and its effects on estimation in phase diversity

The method of phase diversity has been used in the context of incoherent imaging to estimate jointly an object that is being imaged and phase aberrations induced by atmospheric turbulence. The method requires a parametric model for the phase-aberration function. Typically, the parameters are coefficients to a finite set of basis functions. Care must be taken in selecting a parameterization that properly balances accuracy in the representation of the phase-aberration function with stability in the estimates. It is well known that over parameterization can result in unstable estimates. Thus a certain amount of model mismatch is often desirable. We derive expressions that quantify the bias and variance in object and aberration estimates as a function of parameter dimension.

Thelen, Brian J.↗

NASA NTRS1977

Over-relaxation applied to the MacCormack finite-difference scheme

An over-relaxation procedure is applied to the MacCormack finite-difference scheme in order to reduce the computation time required to obtain a steady-state solution. The implementation of this acceleration procedure to an existing computer program using the regular MacCormack method is extremely simple and does not require additional storage. The over-relaxation procedure does not alter the steady-state solution, which is second-order accurate. The method is first applied to Burgers' equation. A stability condition and an expression for the increase in the rate of convergence are derived. The method is then applied to the calculation of the hypersonic viscous flow over a flat plate, using the complete Navier-Stokes equations, and the inviscid flow over a wedge. Reductions in computing time by factors of 3 and 1.5, respectively, are obtained by over-relaxation.

Desideri, J.-A.↗

NASA NTRS1985

Axisymmetric solid elements by a rational hybrid stress method

Four-node axisymmetric solid elements are derived by a new version of hybrid method for which the assumed stresses are expressed in complete polynomials in natural coordinates. The stress equilibrium conditions are introduced through the use of additional displacements as Lagrange multipliers. A rational procedure is to choose the displacement terms such that the resulting strains are also of complete polynomials of the same order. Example problems all indicate that elements obtained by this procedure lead to better results in displacements and stresses than that by other finite elements.

Tian, Z.↗

NASA NTRS1989

Multi-grid for structures analysis

In structural analysis the amount of computational time necessary for a solution is proportional to the number of degrees of freedom times the bandwidth squared. In implicit time analysis, this must be done at each discrete point in time. If, in addition, the problem is nonlinear, then this solution must be iterated at each point in time. If the bandwidth is large, the size of the problem that can be analyzed is severely limited. The multi-grid method is a possible algorithm which can make this solution much more computationally efficient. This method has been used for years in computational fluid mechanics. It works on the fact that relaxation is very efficient on the high frequency components of the solution (nearest neighbor interactions) and not very good on low frequency components of the solution (far interactions). The multi-grid method is then to relax the solution on a particular model until the residual stops changing. This indicates that the solution contains the higher frequency components. A coarse model is then generated for the lower frequency components to the solution. The model is then relaxed for the lower frequency components of the solution. These lower frequency components are then interpolated to the fine model. In computational fluid mechanics the equations are usually expressed as finite differences.

Kascak, Albert F.↗

NASA NTRS1977

Impedance characteristics of coaxial and planar magnetoplasma capacitors

A theory has been developed for the impedance of a homogeneous magnetoplasma enclosed between two specular reflecting coaxial electrodes, with a static magnetic field parallel to the electrode axes. The parallel-plate magnetoplasma capacitor is treated as a sub-case. Starting with the Vlasov equation, an integral equation is derived for the electric field. Solving this equation, and integrating to obtain the voltage, gives the capacitor impedance. This includes a capacitive component, and a resistive component expressing the Landau damping associated with the open orbits of electrons reflected at the electrodes. A direct numerical solution of the field integral equation has been carried out for a range of values of magnetic field, plasma density, and signal frequency. The values of impedance so obtained are compared with the predictions of macroscopic theory, and of an approximate microscopic theory in which open orbits are ignored and solutions are obtained using finite Fourier transform methods. The mathematical relations between these theories are demonstrated.

Harker, K. J.↗

NASA NTRS1984

Analytical description of an imperfect Fabry-Perot etalon

A novel method is presented for the investigation of Fabry-Perot etalon defects, incorporating the effects of microsurface imperfection, bowing, and departure from parallelism, as well as the effect of finite aperture on the etalon, into the perturbation function. The tolerance imperfections of reflective layers can be incorporated into the transfer coefficients of a noncorrugated thin-film system, and any defect can be investigated either separately or in combination with others through the use of the appropriate analytical or numerical simulation of the corrugated function. In the context of numerical computation, the present method permits the simulation and evaluation of different specific situations. The expressions derived for the transfer coefficients of a corrugated, equidistant thin film system may be applicable to other branches of applied optics.

Prikryl, I.↗

NASA NTRS1990

Panel methods: An introduction

Panel methods are numerical schemes for solving (the Prandtl-Glauert equation) for linear, inviscid, irrotational flow about aircraft flying at subsonic or supersonic speeds. The tools at the panel-method user's disposal are (1) surface panels of source-doublet-vorticity distributions that can represent nearly arbitrary geometry, and (2) extremely versatile boundary condition capabilities that can frequently be used for creative modeling. Panel-method capabilities and limitations, basic concepts common to all panel-method codes, different choices that were made in the implementation of these concepts into working computer programs, and various modeling techniques involving boundary conditions, jump properties, and trailing wakes are discussed. An approach for extending the method to nonlinear transonic flow is also presented. Three appendices supplement the main test. In appendix 1, additional detail is provided on how the basic concepts are implemented into a specific computer program (PANAIR). In appendix 2, it is shown how to evaluate analytically the fundamental surface integral that arises in the expressions for influence-coefficients, and evaluate its jump property. In appendix 3, a simple example is used to illustrate the so-called finite part of the improper integrals.

Erickson, Larry L.↗

NASA NTRS1987

Velocity selection at large undercooling in a two-dimensional nonlocal model of solidification

The formation of needle-crystal dendrites from an undercooled melt is investigated analytically, applying the method of Caroli et al. (1986) to Langer's (1980) symmetric two-dimensional nonlocal model of solidification with finite anisotropy in the limit of large undercooling. A solution based on the WKB approximation is obtained, and a saddle-point evaluation is performed. It is shown that needle-crystal solutions exist only if the capillary anisotropy is nonzero, in which case a particular value of the growth velocity can be selected. This finding and the expression for the dependence of the selected velocity on the singular perturbation parameter and the strength of the anisotropy are found to be in complete agreement with the results of a boundary-layer model (Langer and Hong, 1986).

Barbieri, Angelo↗

NASA NTRS2002

Dynamic Stability of Uncertain Laminated Beams Under Subtangential Loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. In this research, we focus on the deterministic and probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. Thus a shear-deformable laminated beam element, including warping effects, is derived to study the deterministic and probabilistic response of laminated beams. This twenty-one degrees of freedom element can be used for solving both static and dynamic problems. In the first-order shear deformable model used here we have employed a more accurate method to obtain the transverse shear correction factor. The dynamic version of the principle of virtual work for laminated composites is expressed in its nondimensional form and the element tangent stiffness and mass matrices are obtained using analytical integration The stability is studied by giving the structure a small disturbance about an equilibrium configuration, and observing if the resulting response remains small. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: Exact Monte Carlo Simulation, Sensitivity Based Monte Carlo Simulation, and Probabilistic FEA. These methods were integrated into the developed finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems, as well as to study the stability analysis, using the dynamic criterion.

Goyal, Vijay K.↗

NASA NTRS1989

Radiative transfer theory for active remote sensing of a forested canopy

A canopy is modeled as a two-layer medium above a rough interface. The upper layer stands for the forest crown, with the leaves modeled as randomly oriented and distributed disks and needles and the branches modeled as randomly oriented finite dielectric cylinders. The lower layer contains the tree trunks, modeled as randomly positioned vertical cylinders above the rough soil. Radiative-transfer theory is applied to calculate EM scattering from such a canopy, is expressed in terms of the scattering-amplitude tensors (SATs). For leaves, the generalized Rayleigh-Gans approximation is applied, whereas the branch and trunk SATs are obtained by estimating the inner field by fields inside a similar cylinder of infinite length. The Kirchhoff method is used to calculate the soil SAT. For a plane wave exciting the canopy, the radiative-transfer equations are solved by iteration to the first order in albedo of the leaves and the branches. Numerical results are illustrated as a function of the incidence angle.

Karam, M. A.↗

NASA NTRS1987

The effect of wave breaking on wave spectrum in water of finite depth

An approximate method is devised to compute the energy-containing portion of the spectrum of waves in water of finite depth, taking into account the effect of wave breaking. It is assumed that there exists a linear and Gaussian ideal wave train whose spectrum is first calculated using the wave energy flux balance equation without considering wave breaking. The Miche wave-breaking criterion for waves in water of finite depth is then applied to limit the wave elevation and establish an expression for the breaking wave elevation in terms of the elevation and elevation's second time derivative of the ideal waves. Simple expressions for the mean value, the mean square value, and the spectrum of the breaking waves are then obtained, and numerical results are presented graphically.

Tung, C. C.↗

NASA NTRS1982

A theoretical study of a two-wavelength lidar technique for the measurement of atmospheric temperature profiles

The theory of differential absorption lidar measurements for lines with a Voigt profile is given and applied to a two-wavelength technique for measuring the atmospheric temperature profile using a high J line in the oxygen A band. Explicit expressions for the temperature and pressure dependence of the absorption coefficient are developed for lines with a Voigt profile. An iteration procedure for calculating the temperature for narrow laser bandwidths is described which has an accuracy better than 0.2 K for bandwidths less than 0.01/cm. To reduce the errors in lidar measurements due to uncertainties in pressure, a method for estimating the pressure from the temperature profile is described. A procedure for extending the differential absorption technique to the case of finite laser bandwidth with good accuracy is also described. Simulation results show that a knowledge of the laser frequency is needed to 0.005/cm for accurate temperature measurements. Evaluation of the sensitivity for both ground- and Shuttle-based measurements shows accuracies generally better than 1 K. This technique allows up to an order of magnitude improvement in sensitivity compared to other differential absorption lidar techniques.

Korb, C. L.↗

NASA NTRS1987

Exact image theory for the problem of dielectric/magnetic slab

Exact image method, recently introduced for the exact solution of electromagnetic field problems involving homogeneous half spaces and microstrip-like geometries, is developed for the problem of homogeneous slab of dielectric and/or magnetic material in free space. Expressions for image sources, creating the exact reflected and transmitted fields, are given and their numerical evaluation is demonstrated. Nonradiating modes, guided by the slab and responsible for the loss of convergence of the image functions, are considered and extracted. The theory allows, for example, an analysis of finite ground planes in microstrip antenna structures.

Lindell, I. V.↗

NASA NTRS1987

Effect of current on spectrum of breaking waves in water of finite depth

This paper presents an approximate method to compute the mean value, the mean square value and the spectrum of waves in water of finite depth taking into account the effect of wave breaking with or without the presence of current. It is assumed that there exists a linear and Gaussian ideal wave train whose spectrum is first obtained using the wave energy flux balance equation without considering wave breaking. The Miche wave breaking criterion for waves in finite water depth is used to limit the wave elevation and establish an expression for the breaking wave elevation in terms of the elevation and its second time derivative of the ideal waves. Simple expressions for the mean value, the mean square value and the spectrum are obtained. These results are applied to the case in which a deep water unidirectional wave train, propagating normally towards a straight shoreline over gently varying sea bottom of parallel and straight contours, encounters an adverse steady current whose velocity is assumed to be uniformly distributed with depth. Numerical results are obtained and presented in graphical form.

Tung, C. C.↗

NASA NTRS2010

Repeatability and Reproducibility of Compression Strength Measurements Conducted According to ASTM E9

Ten commercial laboratories participated in an interlaboratory study to establish the repeatability and reproducibility of compression strength tests conducted according to ASTM International Standard Test Method E9. The test employed a cylindrical aluminum AA2024-T351 test specimen. Participants measured elastic modulus and 0.2 % offset yield strength, YS(0.2 % offset), using an extensometer attached to the specimen. The repeatability and reproducibility of the yield strength measurement, expressed as coefficient of variations were cv(sub r)= 0.011 and cv(sub R)= 0.020 The reproducibility of the test across the laboratories was among the best that has been reported for uniaxial tests. The reported data indicated that using diametrically opposed extensometers, instead of a single extensometer doubled the precision of the test method. Laboratories that did not lubricate the ends of the specimen measured yield stresses and elastic moduli that were smaller than those measured in laboratories that lubricated the specimen ends. A finite element analysis of the test specimen deformation for frictionless and perfect friction could not explain the discrepancy, however. The modulus measured from stress-strain data were reanalyzed using a technique that finds the optimal fit range, and applies several quality checks to the data. The error in modulus measurements from stress-strain curves generally increased as the fit range decreased to less than 40 % of the stress range.

Luecke, William E.↗

NASA NTRS1982

Nonlinear Green's function method for unsteady transonic flows

Advantages to employing Green's function in describing unsteady three-dimensional transonic flows are explored. The development of the function for application to linear subsonic and supersonic unsteady aerodynamics is reviewed. It is shown that unique solutions are possible for external flows, with all functional expressions being defined in Prandtl-Glauert space. The development of methods of using the Green's function for transonic flows is traced, noting the necessity of including the effects of significant nonlinear terms. The steady-state problem is considered to demonstrate the shock-capturing ability of the method and the usefulness of the function in the incompressible, subsonic, transonic, and supersonic areas of potential unsteady three-dimensional flows around complex configurations. Computational time is asserted to be an order of magnitude less than with finite difference methods.

Tseng, K.↗