Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local ordering”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: A careful study of the boundary error

The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Carpenter, Mark H.↗
NASA NTRS2002

The Connection Between Local Icosahedral Order in Metallic Liquids and the Nucleation of Ordered Phases

Over fifty years ago, David Turnbull showed that the temperature of many metallic liquids could be decreased far below their equilibrium melting temperature before crystallization occurred. To explain those surprising results, Charles Frank hypothesized that the local structures of undercooled metallic liquids are different from those of crystal phases, containing a significant degree of icosahedral order that is incompatible with extended periodicity. Such structural differences must create a barrier to the formation crystal phases, explaining the observed undercooling behavior. If true, the nucleation from the liquid of phases with extended icosahedral order should be easier. Icosahedral order is often favored in small clusters, as observed recently in liquid-like clusters of pure Pb on the (111) surface of Si(3), for example. However, it has never been shown that an increasing preference for icosahedral phase formation can be directly linked with the development of icosahedral order in the undercooled liquid. Owing to the combination of very recent advances in levitation techniques and the availability of synchrotron X-ray and high flux neutron facilities.

Kelton, K. F.↗
NASA NTRS2003

The Connection Between Local Icosahedral Order in Metallic Liquids and the Nucleation of Ordered Phases

Over fifty years ago, David Turnbull showed that the temperature of many metallic liquids could be decreased far below their equilibrium melting temperature before crystallization occurred. To explain those surprising results, Charles Frank hypothesized that the local structures of undercooled metallic liquids are different from those of crystal phases, containing a significant degree of icosahedral order that is incompatible with extended periodicity. Such structural differences must create a barrier to the formation crystal phases, explaining the observed undercooling behavior. If true, the nucleation from the liquid of phases with extended icosahedral order should be easier. Icosahedral order is often favored in small clusters, as observed recently in liquid-like clusters of pure Pb on the (111) surface of Si, for example. However, it has never been shown that an increasing preference for icosahedral phase formation can be directly linked with the development of icosahedral order in the undercooled liquid. Owing to the combination of very recent advances in levitation techniques and the availability of synchrotron x-ray and high flux neutron facilities, this is shown here.

Curreri, Peter A.↗
NASA NTRS2003

The Connection Between Local Icosahedral Order in Metallic Liquids and the Nucleation Behavior of Ordered Phases

Over fifty years ago, David Turnbull showed that the temperature of many metallic liquids could be decreased far below their equilibrium melting temperature before crystallization occurred. To explain those surprising results, Charles Frank hypothesized that the local structures of undercooled metallic liquids are different from those of crystal phases, containing a significant degree of icosahedral order that is incompatible with extended periodicity. Such structural differences must create a barrier to the formation crystal phases, explaining the observed undercooling behavior. If true, the nucleation from the liquid of phases with extended icosahedral order should be easier. Icosahedral order is often favored in small clusters, as observed recently in liquid-like clusters of pure Pb on the (111) surface of Si[3], for example. However, it has never been shown that an increasing preference for icosahedral phase formation can be directly linked with the development of icosahedral order in the undercooled liquid. Owing to the combination of very recent advances in levitation techniques and the availability of synchrotron x-ray and high flux neutron facilities, this is shown here.

Kelton, K. F.↗
NASA NTRS1985

Texture classification by local rank correlation

A new approach to texture classification based on local rank correlation is proposed here. Its performance is compared with Laws' method which uses local convolution with feature masks. In the experiments, texture samples are classified based on their distribution of local statistics, either rank correlations or convolutions. The new method achieves generally optimal classification rates. It appears to be more robust because local order statistics are unaffected by local sample differences due to monotonic shifts of texture gray values and are less sensitive to noise.

Harwood, D.↗
NASA NTRS1981

Influence of boundary approximations and conditions on finite difference solutions

Numerical representations of boundary approximations and conditions for three problems are investigated to determine the resulting global accuracy of the steady state solution. Numerical accuracy with various boundary approximations is determined for quasi-one-dimensional inviscid flow in a duct with the interior grid points evaluated using the MacCormack scheme. When an extrapolation approximation with first order local truncation error is used, the global second order accuracy of the difference scheme can be destroyed. For one dimensional flow in a porous medium, an implicit midpoint difference scheme which is consistent with the boundary conditions is developed without the need of boundary approximations. A dissipative model problem is solved with the boundary conditions discretized with first and second order accuracy. The overall second order accuracy of the difference scheme is destroyed if first order numerical representation of one of the boundary conditions is used. With a boundary approximation, the second order global accuracy of the model problem is retained if either second order extrapolation or first order representation of the governing equation is used.

Blottner, F. G.↗
NASA NTRS1994

An HP Adaptive Discontinuous Galerkin Method for Hyperbolic Conservation Laws

This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.

Bey, Kim S.↗
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗
NASA NTRS2005

Multi-Dimensional Spectral Difference Method for Unstructured Grids

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. It combines the best features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. Universal reconstructions are obtained by distributing unknowns in a geometrically similar manner for all unstructured cells. Placements of the unknown and flux points with various order of accuracy are given for the line, triangular and tetrahedral elements. The data structure of the new method permits an optimum use of cache memory, resulting in further computational efficiency on modern computers. A new pointer system is developed that reduces memory requirements and simplifies programming for any order of accuracy. Numerical solutions are presented and compared with the exact solutions for wave propagation problems in both two and three dimensions to demonstrate the capability of the method. Excellent agreement has been found. The method is simpler and more efficient than previous discontinuous Galerkin and spectral volume methods for unstructured grids.

Liu, Yen↗
NASA NTRS1996

Defect Structure of Beta NiAl Using the BFS Method for Alloys

The semiempirical BFS method for alloys is generalized by replacing experimental input with first-principles results thus allowing for the study of complex systems. In order to examine trends and behavior of a system in the vicinity of a given point of the phase diagram a search procedure based on a sampling of selected configurations is employed. This new approach is applied to the study of the beta phase of the Ni-Al system, which exists over a range of composition from 45-60 at.% Ni. This methodology results in a straightforward and economical way of reproducing and understanding the basic features of this system. At the stoichiometric composition, NiAl should exist in a perfectly ordered B2 structure. Ni-rich alloys are characterized by antisite point defects (with Ni atoms in the Al sites) with a decrease in lattice parameters. On the Al-rich side of stoichiometry there is a steep decrease in lattice parameter and density with increasing Al content. The presence of vacancies in Ni sites would explain such behavior. Recent X-ray diffraction experiments suggest a richer structure: the evidence, while strongly favoring the presence of vacancies in Ni sites, also suggests the possibility of some vacancies in Al sites in a 3:1 ratio. Moreover, local ordering of vacant sites may be preferred over a random distribution of individual point defects.

Bozzolo, Guillermo↗
NASA NTRS1991

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary tasks performed are: (1) the development of a second order local thermodynamic nonequilibrium (LTNE) model for atoms; (2) the continued development of vibrational nonequilibrium models; and (3) the development of a new multicomponent diffusion model. In addition, studies comparing these new models with previous models and results were conducted and reported.

Carlson, Leland A.↗
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. A discussion on the Discontinuous Spectral Difference (SD) Method, locations of the unknowns and flux points and numerical results are also presented.

Liu, Yen↗
NASA NTRS1998

Electronics Manufacturer Provided With Testing and Evaluation Data Necessary to Obtain Additional Orders

A local electronics manufacturer, the Sterling Manufacturing Company, was presented with the opportunity to supply 30,000 automotive cellular antennas to a European subsidiary of a large U.S. auto manufacturer. Although the company built an antenna that they believed would meet the auto manufacturer's specifications, they were unable to conduct the necessary validation tests in-house. They decided to work with NASA Lewis Research Center's Space Electronics Division, which, as part of its technology development program, evaluates the performance of antennas in its Microwave Systems Lab to assess their capabilities for space communications applications. Data measured in Lewis' Microwave Systems Lab proved that Sterling's antenna performed better than specified by the auto manufacturer.

Source record↗
NASA NTRS1993

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗
NASA NTRS1986

Embedded mesh solution of the 2-D Euler equations - Evaluation of interface formulations

Solution of the steady 2-D Euler equations using mesh embedding, or local grid refinement, with a cell-centered finite volume scheme is investigated. Embedded regions which are topologically similar to the global grid are considered. An isoenergetic model for the governing equations is used in Jameson's finite volume multistage scheme with modifications to the boundary conditions and smoothing. A detailed study of the embedding interface flux and smoothing formulations is conducted. Taylor expansion analysis reveals that local second order spatial accuracy is not possible if a conservative interface flux formulation is used. The analysis also gives constraints for local first order accuracy. An energy stability analysis indicates that downwind weighting of interface fluxes causes local instabilities. Analysis shows that conservative interface smoothing formulations must have a locally convective component, but that correct interface formulations allow globally dissipative smoothing. Embedded mesh solutions obtained with this scheme are presented for a transonic airfoil. They show that if embedding interfaces are close to the shocks, then small modifications in the interface location can have large effects on converge and solution accuracy.

Allmaras, S. R.↗
NASA NTRS1994

Accurate finite difference methods for time-harmonic wave propagation

Finite difference methods for solving problems of time-harmonic acoustics are developed and analyzed. Multidimensional inhomogeneous problems with variable, possibly discontinuous, coefficients are considered, accounting for the effects of employing nonuniform grids. A weighted-average representation is less sensitive to transition in wave resolution (due to variable wave numbers or nonuniform grids) than the standard pointwise representation. Further enhancement in method performance is obtained by basing the stencils on generalizations of Pade approximation, or generalized definitions of the derivative, reducing spurious dispersion, anisotropy and reflection, and by improving the representation of source terms. The resulting schemes have fourth-order accurate local truncation error on uniform grids and third order in the nonuniform case. Guidelines for discretization pertaining to grid orientation and resolution are presented.

Harari, Isaac↗
NASA NTRS1990

A cost-effective strategy for nonoscillatory convection without clipping

Clipping of narrow extrema and distortion of smooth profiles is a well known problem associated with so-called high resolution nonoscillatory convection schemes. A strategy is presented for accurately simulating highly convective flows containing discontinuities such as density fronts or shock waves, without distorting smooth profiles or clipping narrow local extrema. The convection algorithm is based on non-artificially diffusive third-order upwinding in smooth regions, with automatic adaptive stencil expansion to (in principle, arbitrarily) higher order upwinding locally, in regions of rapidly changing gradients. This is highly cost effective because the wider stencil is used only where needed-in isolated narrow regions. A recently developed universal limiter assures sharp monotonic resolution of discontinuities without introducing artificial diffusion or numerical compression. An adaptive discriminator is constructed to distinguish between spurious overshoots and physical peaks; this automatically relaxes the limiter near local turning points, thereby avoiding loss of resolution in narrow extrema. Examples are given for one-dimensional pure convection of scalar profiles at constant velocity.

Leonard, B. P.↗