Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice calculation”

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

Monte Carlo calculations of elementary particle properties

The object of this project is to calculate the masses of the elementary particles. This ambitious goal apparently is not possible using analytic methods or known approximation methods. However, it is probable that the power of a modern super computer will make at least part of the low lying mass spectrum accessible through direct numerical computation. Initial attempts by several groups at calculating this spectrum on small lattices of space time points have been very promising. Using new methods and super computers considerable progress has been made towards evaluating the mass spectrum on comparatively large lattices. Only more time and faster machines with increased storage will allow calculations of systems with guaranteed minimal boundary effects. The ideas that currently go into this calculation are outlined.

Guralnik, G. S.↗
NASA NTRS1983

High angle-of-attack calculations of the subsonic vortex flow on slender bodies

A nonlinear vortex-lattice method is utilized for the calculation of separated vortex flows over slender bodies at high incidence. Symmetric and asymmetric vortex flow cases are calculated showing good agreement with recently obtained experimental data for an ogive-cyclinder body at high Reynolds numbers. The only input needed for these calculations is the positions of the separation lines on the body. Studies of some of the numerical aspects of this method are described and conclusions are implemented to improve the calculations.

Almosnino, D.↗
NASA NTRS1975

Aluminum under high pressure: 1 equation of state

A curve of applied pressure P versus lattice constant a is calculated for single crystal aluminum. It results from an application of the method of structural expansions for deriving the energies of simple metals, a method known to give reasonable results for the elastic constants even at second order in the effective electron-ion interaction. The method is taken from Fermi surface analysis . This essentially experimental information is used to verify that the extant face-centered cubic structure remains the preferred crystalline phase up to the highest pressures considered. Arguments are given to suggest that the P versus a curve should have reasonable a priori accuracy, which may be improved if experimental data in the intermediate pressure region can be provided to refine the energy dependent pseudo-potential. At 3 megabars the lattice constant is reduced by only 22%; the ion cores at this pressure are still very well separated.

Friedli, C.↗
NASA NTRS1991

A unified viscous theory of lift and drag of 2-D thin airfoils and 3-D thin wings

A unified viscous theory of 2-D thin airfoils and 3-D thin wings is developed with numerical examples. The viscous theory of the load distribution is unique and tends to the classical inviscid result with Kutta condition in the high Reynolds number limit. A new theory of 2-D section induced drag is introduced with specific applications to three cases of interest: (1) constant angle of attack; (2) parabolic camber; and (3) a flapped airfoil. The first case is also extended to a profiled leading edge foil. The well-known drag due to absence of leading edge suction is derived from the viscous theory. It is independent of Reynolds number for zero thickness and varies inversely with the square root of the Reynolds number based on the leading edge radius for profiled sections. The role of turbulence in the section induced drag problem is discussed. A theory of minimum section induced drag is derived and applied. For low Reynolds number the minimum drag load tends to the constant angle of attack solution and for high Reynolds number to an approximation of the parabolic camber solution. The parabolic camber section induced drag is about 4 percent greater than the ideal minimum at high Reynolds number. Two new concepts, the viscous induced drag angle and the viscous induced separation potential are introduced. The separation potential is calculated for three 2-D cases and for a 3-D rectangular wing. The potential is calculated with input from a standard doublet lattice wing code without recourse to any boundary layer calculations. Separation is indicated in regions where it is observed experimentally. The classical induced drag is recovered in the 3-D high Reynolds number limit with an additional contribution that is Reynold number dependent. The 3-D viscous theory of minimum induced drag yields an equation for the optimal spanwise and chordwise load distribution. The design of optimal wing tip planforms and camber distributions is possible with the viscous 3-D wing theory.

Yates, John E.↗
NASA NTRS1989

Aeroelastic modeling for the FIT team F/A-18 simulation

Some details of the aeroelastic modeling of the F/A-18 aircraft done for the Functional Integration Technology (FIT) team's research in integrated dynamics modeling and how these are combined with the FIT team's integrated dynamics model are described. Also described are mean axis corrections to elastic modes, the addition of nonlinear inertial coupling terms into the equations of motion, and the calculation of internal loads time histories using the integrated dynamics model in a batch simulation program. A video tape made of a loads time history animation was included as a part of the oral presentation. Also discussed is work done in one of the areas of unsteady aerodynamic modeling identified as needing improvement, specifically, in correction factor methodologies for improving the accuracy of stability derivatives calculated with a doublet lattice code.

Zeiler, Thomas A.↗
NASA NTRS1985

Deconfining phase transition and the continuum limit of lattice quantum chromodynamics

A large-scale Monte Carlo calculation is presented of the deconfining phase-transition temperature in lattice quantum chromodynamics without fermions. By using the Wilson action, it is found that the transition temperature as a function of the lattice coupling g is consistent with scaling behavior dictated by the perturbative beta function for 6/g-squared greater than 6.15.

Gottlieb, S. A.↗
NASA NTRS2001

Efficient Cache use for Stencil Operations on Structured Discretization Grids

We derive tight bounds on the cache misses for evaluation of explicit stencil operators on structured grids. Our lower bound is based on the isoperimetrical property of the discrete octahedron. Our upper bound is based on a good surface to volume ratio of a parallelepiped spanned by a reduced basis of the interference lattice of a grid. Measurements show that our algorithm typically reduces the number of cache misses by a factor of three, relative to a compiler optimized code. We show that stencil calculations on grids whose interference lattice have a short vector feature abnormally high numbers of cache misses. We call such grids unfavorable and suggest to avoid these in computations by appropriate padding. By direct measurements on a MIPS R10000 processor we show a good correlation between abnormally high numbers of cache misses and unfavorable three-dimensional grids.

Frumkin, Michael↗
NASA NTRS1984

Wind tunnel wall interference in closed, ventilated and adaptive test sections

A wall interference correction method for closed rectangular test sections was developed which uses measured wall pressures. Measurements with circular discs for blockage and a rectangular wing as a lift generator in a square closed test section validate this method. These measurements are intended to be a basis of comparison for measurements in the same tunnel using ventilated (in these case, slotted) walls. Using the vortex lattice method and homogeneous boundary conditions, calculations were performed which show sufficiently high pressure levels at the walls for correction purposes in test sections with porous walls. In Gottingen, an adaptive test section (which is a deformable rubber tube of 800 mm diameter) was built and a computer program was developed which is able to find the necessary wall adaptation for interference-free measurements in a single step. To check the program prior to the first run, the vortex lattice method was used to calculate wall pressure distributions in the nonadapted test section as input data for the one-step method. Comparison of the pressure distribution in the adapted test section with free-flight data shows nearly perfect agreement. An extension of the computer program can be made to evaluate the remaining interference corrections.

Holst, H.↗
NASA NTRS1990

The effect of frictional stress on the calculation of critical thickness in epitaxy

A model is developed for lattice-mismatched epitaxial layers, which makes it possible to calculate both the equilibrium critical thickness as well as the kinetic rate of misfit accommodations. The technique for calculating the equilibrium critical thickness is similar to that of Matthews (1966) in that it uses the force balance approach and the Volterra description of the dislocation line. For the kinetic component of misfit accommodation, the model uses a dynamic frictional force. A simple criterion is proposed for determining whether the generation of misfit dislocations will be controlled by the equilibrium or the kinetic factors.

Fox, Bradley A.↗
NASA NTRS1975

A combined representation method for use in band structure calculations. 1: Method

A representation was described whose basis levels combine the important physical aspects of a finite set of plane waves with those of a set of Bloch tight-binding levels. The chosen combination has a particularly simple dependence on the wave vector within the Brillouin Zone, and its use in reducing the standard one-electron band structure problem to the usual secular equation has the advantage that the lattice sums involved in the calculation of the matrix elements are actually independent of the wave vector. For systems with complicated crystal structures, for which the Korringa-Kohn-Rostoker (KKR), Augmented-Plane Wave (APW) and Orthogonalized-Plane Wave (OPW) methods are difficult to apply, the present method leads to results with satisfactory accuracy and convergence.

Friedli, C.↗
NASA NTRS1985

VORSTAB: A computer program for calculating lateral-directional stability derivatives with vortex flow effect

A computer program based on the Quasi-Vortex-Lattice Method of Lan is presented for calculating longitudinal and lateral-directional aerodynamic characteristics of nonplanar wing-body combination. The method is based on the assumption of inviscid subsonic flow. Both attached and vortex-separated flows are treated. For the vortex-separated flow, the calculation is based on the method of suction analogy. The effect of vortex breakdown is accounted for by an empirical method. A summary of the theoretical method, program capabilities, input format, output variables and program job control set-up are described. Three test cases are presented as guides for potential users of the code.

Lan, C. Edward↗
NASA NTRS1980

Higher order farfield drag minimization for a subcritical wing design code

A higher order Trefftz plane model of the undistorted interacting wing wakes of two symmetric subsonic planforms has been developed and integrated with an existing vortex lattice wing design code. The Trefftz plane calculation provides continuous, piecewise quadratically varying bound circulations, which are then interpolated in the nearfield vortex lattice representation of the wing. Integration of the surface slopes for a specified chord loading function provides the design wing camber surfaces. This paper outlines the theoretical development and presents comparisons between results obtained using the present theory and (1) previous exact solutions for minimum induced drag and (2) previous wing designs obtained using a discrete vortex wake moded.

Kuhlman, J. M.↗
NASA NTRS1982

Transonic flutter study of a wind-tunnel model of a supercritical wing with/without winglet

The scaled flutter model was a 1/6.5-size, semispan version of a supercritical wing (SCW) proposed for an executive-jet-transport airplane. The model was tested cantilever-mounted with a normal wingtip, a wingtip with winglet, and a normal wingtip ballasted to simulate the winglet mass properties. Flutter and aerodynamic data were acquired at Mach numbers from 0.6 to 0.95. The measured transonic flutter speed boundary for each wingtip configuration had roughly the same shape with a minimum flutter speed near M = 0.82. The winglet addition and wingtip mass ballast decreased the wing flutter speed by about 7 and 5%, respectively; thus, the winglet effect on flutter was more a mass effect than an aerodynamic effect. Flutter characteristics calculated using a doublet-lattice analysis (which included interference effects) were in good agreement with the experimental results up to M = 0.82. Comparisons of measured static aerodynamic data with predicted data indicated that the model was aerodynamically representative of the airplane SCW.

Ruhlin, C. L.↗
NASA NTRS1983

Low frequency dielectric relaxation in boracites

In order to elucidate the nature of the imperfections which adversely affect pyroelectric processes in boracites, the thermal and dispersive characteristics of the low frequency dielectric response in iron-iodide boracite (Fe3B7O13I) and copper-chloride boracite (Cu3B7O13Cl) have been investigated. These characteristics were measured as a function of crystallographic orientation and applied field in both the ferroelectric and paraelectric states. The low frequency dielectric relaxation of 100 line oriented multi-domain copper-chloride boracite clearly indicates the dipole nature of the lattice imperfections. The activation energies calculated from a noninteracting Debye model, are 0.53 eV in the ferroelectric phase and 0.10 eV in the paraelectric phase.

Kim, Q.↗
NASA NTRS1993

Anisotropic induced Compton scattering - Constraints on models of active galactic nuclei

A model for nonlinear radiative transfer on a lattice is used to explore observational signatures of anisotropic induced Compton scattering. Ways in which an intense beam of radiation may be altered in passing through a slab of plasma that is optically thin to spontaneous Thomson scattering are explained. A numerical approach to induced Compton scattering which involves solving the equations of nonlinear transfer on a lattice is described. Results of calculations for some simple geometrical configurations (spheres, disks, cylindrical 'jets', etc.) are presented. Some possible applications of these results to the interpretation of actual data on compact radio sources are outlined, and the FIR spectra of AGN are considered.

Coppi, P.↗
NASA NTRS1985

A vortex-lattice method for general, unsteady aerodynamics

A general method of calculating unsteady, incompressible, inviscid, three-dimensional flows around arbitrary planforms has been developed. The method is an extension of the vortex-lattice technique. It is not limited by aspect ratio, camber, or angle of attack, as long as vortex breakdown does not occur above the surface of the wing and separation occurs only along sharp edges. As the wing performs arbitrary maneuvers, the position of the wake and the distribution of circulation on the wing and in the wake are obtained as functions of time. One desirable feature of the present method is its ability to treat steady lifting flows very efficiently. Several examples of steady and unsteady flows are presented. These include rectangular wings, with and without flaps, delta, and cropped delta wings.

Konstadinopoulos, P.↗
NASA NTRS1992

Stiffnesses by (TtN) ensemble molecular dynamics

Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.

Tibbits, Patrick↗