SEARCH · Search NASA

Results for “generalized finite differences”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1990

A multiblock multigrid three-dimensional Euler equation solver

Current aerodynamic designs are often quite complex (geometrically). Flexible computational tools are needed for the analysis of a wide range of configurations with both internal and external flows. In the past, geometrically dissimilar configurations required different analysis codes with different grid topologies in each. The duplicity of codes can be avoided with the use of a general multiblock formulation which can handle any grid topology. Rather than hard wiring the grid topology into the program, it is instead dictated by input to the program. In this work, the compressible Euler equations, written in a body-fitted finite-volume formulation, are solved using a pseudo-time-marching approach. Two upwind methods (van Leer's flux-vector-splitting and Roe's flux-differencing) were investigated. Two types of explicit solvers (a two-step predictor-corrector and a modified multistage Runge-Kutta) were used with multigrid acceleration to enhance convergence. A multiblock strategy is used to allow greater geometric flexibility. A report on simple explicit upwind schemes for solving compressible flows is included.

Cannizzaro, Frank E.↗

NASA NTRS1992

An installed nacelle design code using a multiblock Euler solver. Volume 1: Theory document

An efficient multiblock Euler design code was developed for designing a nacelle installed on geometrically complex airplane configurations. This approach employed a design driver based on a direct iterative surface curvature method developed at LaRC. A general multiblock Euler flow solver was used for computing flow around complex geometries. The flow solver used a finite-volume formulation with explicit time-stepping to solve the Euler Equations. It used a multiblock version of the multigrid method to accelerate the convergence of the calculations. The design driver successively updated the surface geometry to reduce the difference between the computed and target pressure distributions. In the flow solver, the change in surface geometry was simulated by applying surface transpiration boundary conditions to avoid repeated grid generation during design iterations. Smoothness of the designed surface was ensured by alternate application of streamwise and circumferential smoothings. The capability and efficiency of the code was demonstrated through the design of both an isolated nacelle and an installed nacelle at various flow conditions. Information on the execution of the computer program is provided in volume 2.

Chen, H. C.↗

NASA NTRS2001

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗

DOE OSTI2021

Durability of plasma sprayed Thermal Barrier Coatings with controlled properties part Ⅱ: Effects of geometrical curvature

The present investigation elucidates the effects of substrate curvature on the durability of Air Plasma Sprayed (APS) Thermal Barrier Coatings (TBCs). Traditionally, planar disk specimens are utilized in Furnace Cycle Testing (FCT). In most cases, delamination is initiated at the disk's free edge and then propagates along the TBC-bond coat interface. However, in turbine components (e.g., blades/vanes), they lack significant free edges and the coatings are deposited on non-flat surfaces of varying radii of curvature. These geometrical discrepancies imply significant differences in the stress states as compared to planar disk specimens. Therefore, the part geometry inevitably affects the TBC failure mechanisms in the turbine system, and consequently influences TBC durability. Nevertheless, these effects have not been fully explored in the past, despite the anecdotal knowledge which suggests the leading edge of turbine blades to be one of the most common failure/spallation locations. Here in this study, representative TBC systems were deposited onto superalloy disks with flat substrate and rods with curved substrate. The experimental results from FCT suggest TBC durability on rods is significantly lower than when sprayed on disks. Furthermore, the durability of TBCs on rods appears to increase with higher porosity, which is contradictory to reported trends for TBCs on disks. In addition, the effects of various bond coats are not consistent with those observed in disks. To clarify the underlying effects of curved geometries, detailed stress analyses were performed. They revealed a state of thermal stresses which is unique to curved geometries. It was found that TBC on rods experience tensile radial stresses as well as tensile hoop stresses during the cooling. These stresses result in a complex interplay in failure processes in rods, which explains the different observed trends in the durability from those obtained with disks.

36 MATERIALS SCIENCE↗

DOE OSTI2016

2016 KIVA-hpFE Development: A Robust and Accurate Engine Modeling Software

Los Alamos National Laboratory and its collaborators are facilitating engine modeling by improving accuracy and robustness of the modeling, and improving the robustness of software. We also continue to improve the physical modeling methods. We are developing and implementing new mathematical algorithms, those that represent the physics within an engine. We provide software that others may use directly or that they may alter with various models e.g., sophisticated chemical kinetics, different turbulent closure methods or other fuel injection and spray systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Max-independent set and the quantum alternating operator ansatz

he maximum-independent set (MIS) problem of graph theory using the quantum alternating operator ansatz is studied. We perform simulations on the Rigetti Forest simulator for the square ring, K 2,3 , and K3,3 graphs and analyze the dependence of the algorithm on the depth of the circuit and initial states. The probability distribution of observation of the feasible states representing maximum-independent sets is observed to be asymmetric for the MIS problem, which is unlike the Max-Cut problem where the probability distribution of feasible states is symmetric. For asymmetric graphs, it is shown that the algorithm clearly favors the independent set with the larger number of elements even for finite circuit depth. Finally, we also compare the approximation ratios for the algorithm when we choose different initial states for the square ring graph and show that it is dependent on the choice of the initial state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Flexural phonons in supported graphene: from pinning to localization

We identify graphene layer on a disordered substrate as a system where localization of phonons can be observed. Generally, observation of localization for scattering waves is not simple, because the Rayleigh scattering is inversely proportional to a high power of wavelength. The situation is radically different for the out of plane vibrations, so-called flexural phonons, scattered by pinning centers induced by a substrate. In this case, the scattering time for vanishing wave vector tends to a finite limit. One may, therefore, expect that physics of the flexural phonons exhibits features characteristic for electron localization in two dimensions, albeit without complications caused by the electron-electron interactions. We confirm this idea by calculating statistical properties of the Anderson localization of flexural phonons for a model of elastic sheet in the presence of the pinning centers. Finally, we discuss possible manifestations of the flexural phonons, including the localized ones, in the electronic thermal conductance.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Review of the Mechanical Properties of Materials Used in Nb 3 Sn Magnets for Particle Accelerators

Superconducting magnets experience significant thermo-mechanical loads throughout their life cycle. These are introduced by the electro-magnetic forces during powering, but also by the prestress applied in many magnet designs. Further to this, the large thermal excursion that components of different materials experience can generate significant internal forces. The loads are also experienced by the superconducting coils, whose critical current can decrease as a consequence of the applied strain. It is then crucial to predict the overall mechanical behavior and conservatively design a magnet, avoiding failure of the mechanical components and of the superconducting coils. Finite Element Analysis (FEA) is generally used to perform these tasks, but its results rely heavily on the material properties and models used. This is in particular true for the coil composite, which is simplified to allow reasonable model sizes in full magnet models. Here we present the state-of-art knowledge of the mechanical properties of the materials mostly used in superconducting magnet construction. We review elastic and plastic properties at room and cryogenic temperature, thermal contraction, and summarize the state-of-art failure criteria for these materials. Finally, the paper summarizes the present understanding of the mechanical behavior and limits of Nb3Sn coils. For the first time, an orthotropic failure criteria is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1973

Structural expansions for the ground state energy of a simple metal

A structural expansion for the static ground state energy of a simple metal is derived. An approach based on single particle band structure which treats the electron gas as a non-linear dielectric is presented, along with a more general many particle analysis using finite temperature perturbation theory. The two methods are compared, and it is shown in detail how band-structure effects, Fermi surface distortions, and chemical potential shifts affect the total energy. These are of special interest in corrections to the total energy beyond third order in the electron ion interaction, and hence to systems where differences in energies for various crystal structures are exceptionally small. Preliminary calculations using these methods for the zero temperature thermodynamic functions of atomic hydrogen are reported.

Hammerberg, J.↗

NASA NTRS1981

On eigensolutions for discontinuous liners in a duct containing uniform mean flow

Sound attenuation in a rectangular acoustically lined duct containing uniform mean flow is analytically investigated using the generalized Wiener-Hopf technique. Uniqueness of the solution is enforced for lined sections of the finite axial extent by imposing edge conditions at the liner interface. Possible edge conditions are considered, including the Kutta condition, and the causal solution corresponding to edge conditions is considered the best choice. Solution methods such as the mode matching and singularity methods imply differing edge conditions, and results show that power attenuation is insensitive to the imposed edge conditions, although significant differences are observed for the reflection coefficient. The amplitude of the exponentially increasing instability mode in the lined section must be set to zero as a first approximation to the nonlinear situation, and results indicate that measurements of the reflection factor can be used to make a more definite decision about physically appropriate edge conditions.

Koch, W.↗

NASA NTRS1985

New stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = Au sub x + Bu + f in the quarter plane x greater than or equal to 0, t greater than or equal to 0. The criteria are used to easily establish stability for numerous combinations of well known basic schemes and boundary conditions, thus generalizing many special cases studied in the recent literature. A number of examples are examined, including the unitary unconditionally stable Crank-Nicholson scheme and an almost-dissipative unconditionally stable backward Euler scheme.

Goldberg, M.↗

NASA NTRS1989

Shear buckling of specially orthotropic plates with centrally located cutouts

The shear loading of a composite rectangular plate with a centrally located circular cutout was analyzed in order to predict the buckling load of the plate. The first step in this analysis is to calculare the in-plane stress distribution of the plate prior to buckling. This problem can be solved using finite element methods, but a quicker and more cost efficient method is desired. The method chosen to determine the pre-buckling stresses is that of boundary collocation using complex variables. The real valued force functions are written in terms of two complex valued functions, each of which is a function of a different complex variable. The solution of the generalized biharmonic function is a superposition of functions of these two complex variables. For this analysis, the force functions are represented with a Laurent series. The constants in these series are found using boundary collocation. Two force equations for complex variables are satisfied over the plate boundaries to obtain these constants. The stresses in the plate are then known for any particular locations on the plate given the applied shear on the edges of the plate.

Klang, Eric C.↗

NASA NTRS1992

Finite-difference algorithms for the time-domain Maxwell's equations - A numerical approach to RCS analysis

The applications of two CFD-based finite-difference methods to computational electromagnetics are investigated. In the first method, the time-domain Maxwell's equations are solved using the explicit Lax-Wendroff scheme and in the second method, the second-order wave equations satisfying the Maxwell's equations are solved using the implicit Crank-Nicolson scheme. The governing equations are transformed to a generalized curvilinear coordinate system and solved on a body-conforming mesh using the scattered-field formulation. The induced surface current and the bistatic radar cross section are computed and the results are validated for several two-dimensional test cases involving perfectly-conducting scatterers submerged in transverse-magnetic plane waves.

Vinh, Hoang↗

NASA NTRS1993

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗

NASA NTRS1996

Nonlinear Analysis of the Space Shuttle Super-Lightweight External Fuel Tank

The results of buckling and nonlinear analyses of the Space Shuttle External Tank super-lightweight liquid oxygen (LOX) tank are presented. Modeling details and results are presented for two prelaunch loading conditions and for two full-scale structural tests conducted on the original external tank. These results illustrate three distinctly different types of nonlinear responses for thin-walled shells subjected to combined mechanical and thermal loads. These nonlinear response phenomena consist of bifurcation-type buckling, short-wavelength nonlinear bending, and nonlinear collapse associated with a limit point. For each case, the results show that accurate predictions of nonlinear behavior generally require a large scale high-fidelity finite element model. Results are also presented that show that a fluid filled launch vehicle shell can be highly sensitive to initial geometric imperfections. In addition, results presented for two full scale structural tests of the original standard weight external tank suggest that the finite element modeling approach used in the present study is sufficient for representing the nonlinear behavior of the super lightweight LOX tank.

Nemeth, Michael P.↗

DOE OSTI2017

On entropic uncertainty relations in the presence of a minimal length

Entropic uncertainty relations for the position and momentum within the generalized uncertainty principle are examined. Studies of this principle are motivated by the existence of a minimal observable length. Then the position and momentum operators satisfy the modified commutation relation, for which more than one algebraic representation is known. One of them is described by auxiliary momentum so that the momentum and coordinate wave functions are connected by the Fourier transform. However, the probability density functions of the physically true and auxiliary momenta are different. As the corresponding entropies differ, known entropic uncertainty relations are changed. Using differential Shannon entropies, we give a state-dependent formulation with correction term. State-independent uncertainty relations are obtained in terms of the Rényi entropies and the Tsallis entropies with binning. Such relations allow one to take into account a finiteness of measurement resolution. - Highlights: • Entropic uncertainty bounds increase for position and physically true momentum. • Correction terms are always nonzero for wave packets of a finite width. • An increase of uncertainty bounds depends on acceptance functions of apparatuses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1989

An experimental and analytical investigation of the effect of spanwise curvature on wing flutter at Mach number of 0.7

An experimental and analytical study was conducted at Mach 0.7 to investigate the effects of spanwise curvature on flutter. Two series of rectangular planform wings of aspect ration 1.5 and curvature ranging from zero (uncurved) to 1.04/ft were flutter tested in the NASA Langley Transonic Dynamics Tunnel (TDT). One series consisted of models with a NACA 65 A010 airfoil section and the other of flat plate cross section models. Flutter analyses were conducted for correlation with the experimental results by using structural finite element methods to perform vibration analysis and two aerodynamic theories to obtain unsteady aerodynamic load calculations. The experimental results showed that for one series of models the flutter dynamic pressure increased significantly with curvature while for the other series of models the flutter dynamic pressure decreased with curvature. The flutter analyses, which generally predicted the experimental results, indicated that the difference in behavior of the two series of models was primarily due to differences in their structural properties.

Rivera, Jose A., Jr.↗