Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Finite-frequency modeling of regional tropospheric infrasound using realistic atmospheres and terrain

Infrasonic waves have been observed to propagate to regional (greater than 15 km) distances through the troposphere. Infrasound propagation in the geometric acoustics approximation has shown that realistic terrain can scatter acoustic energy from tropospheric ducts; however, ray methods cannot intrinsically capture finite-frequency behavior such as diffraction. A two-dimensional finite-difference time-domain (FDTD) method has been developed to solve linearized equations for infrasound propagation with realistic terrain. Acoustic wave propagation over 100 km with both flat terrain and a Gaussian hill was first simulated in order to compare finite-frequency propagation with ray predictions. The effects of realistic terrain and atmospheres on infrasound signals from a 2012 surface explosion at the Utah Testing and Training Range are then investigated. Propagation through the troposphere is suggested by array processing results, but eigenrays are not predicted due to weak to nonexistent ducting conditions. FDTD modeling suggests that the inclusion of terrain and finite frequency effects helps explain much of the observed signal in a realistic scenario. Furthermore, these results suggest that geometric acoustics may underestimate propagation through the troposphere, and that recorded waveforms at regional distances may be noticeably affected by terrain.

58 GEOSCIENCES↗
NASA NTRS1987

A comparison of VLSI architecture of finite field multipliers using dual, normal or standard basis

Three different finite field multipliers are presented: (1) a dual basis multiplier due to Berlekamp; (2) a Massy-Omura normal basis multiplier; and (3) the Scott-Tavares-Peppard standard basis multiplier. These algorithms are chosen because each has its own distinct features which apply most suitably in different areas. Finally, they are implemented on silicon chips with nitride metal oxide semiconductor technology so that the multiplier most desirable for very large scale integration implementations can readily be ascertained.

Hsu, I. S.↗
DOE OSTI2020

Coupling between plasmonic and photonic crystal modes in suspended three-dimensional meta-films

A complementary metal oxide semiconductor (CMOS) compatible fabrication method for creating three-dimensional (3D) meta-films is presented. In contrast to metasurfaces, meta-films possess structural variation throughout the thickness of the film and can possess a sub-wavelength scale structure in all three dimensions. Here we use this approach to create 2D arrays of cubic silicon nitride unit cells with plasmonic inclusions of elliptical metallic disks in horizontal and vertical orientations with lateral array-dimensions on the order of millimeters. Fourier transform infrared (FTIR) spectroscopy is used to measure the infrared transmission of meta-films with either horizontally or vertically oriented ellipses with varying eccentricity. Shape effects due to the ellipse eccentricity, as well as localized surface plasmon resonance (LSPR) effects due to the effective plasmonic wavelength are observed in the scattering response. The structures were modeled using rigorous coupled wave analysis (RCWA), finite difference time domain (Lumerical), and frequency domain finite element (COMSOL). The silicon nitride support structure possesses a complex in-plane photonic crystal slab band structure due to the periodicity of the unit cells. We show that adjustments to the physical dimensions of the ellipses can be used to control the coupling to this band structure. The horizontally oriented ellipses show narrow, distinct plasmonic resonances while the vertically oriented ellipses possess broader resonances, with lower overall transmission amplitude for a given ellipse geometry. We attribute this difference in resonance behavior to retardation effects. The ability to couple photonic slab modes with plasmonic inclusions enables a richer space of optical functionality for design of metamaterial-inspired optical components.

42 ENGINEERING↗
NASA NTRS1996

Spectral Ambiguity of Allan Variance

We study the extent to which knowledge of Allan variance and other finite-difference variances determines the spectrum of a random process. The variance of first differences is known to determine the spectrum. We show that, in general, the Allan variance does not. A complete description of the ambiguity is given.

Allan Variance variances finite-difference varianc↗
DOE OSTI2020

Parametric study of optical transmission through plasmonic hole arrays modulated by the phase transition of vanadium dioxide

We have performed comprehensive electromagnetic simulations and preliminary experiments to explore the effects of geometrical and material parameters on the extraordinary optical transmission (EOT) through periodic arrays of subwavelength holes in a bilayer stack consisting of a gold or silver film atop a vanadium dioxide film (Au/Ag + VO 2 ), where the latter undergoes a semiconductor-to-metal phase transition. Using the finite-difference time-domain (FDTD) and finite-element methods (FEM), we vary iteratively the array periodicity, VO 2 film thickness and hole diameters, as well as the refractive index inside the VO 2 -layer holes and the VO 2 optical constants. For each variation, we compare the metallic-to-semiconducting ratios of the zero-order transmission ( T 00 ) peaks and find sharp maxima in these ratios within narrow parameter ranges. The maxima arise from Fabry-Perot and Fano-type resonances that minimize T 00 in the semiconducting phase of the perforated bilayers. At a fixed array period, the primary factors controlling the VO 2 -enabled EOT modulation are the VO 2 thickness, diameter of the VO 2 -layer holes, and absorption in the two VO 2 phases. Besides uncovering the origins of the higher metallic-phase T 00 , this study provides a protocol for optimizing the performance of the bilayer hole arrays for potential uses as dynamically tunable nano-optical devices.

Donev, Eugenii U. (ORCID:000000030555219X)↗
DOE OSTI2020

ATDM AD projects to ECP Report (Q3 FY2020)

The SPARC team has completed much of the workflow development needed for the coupled and uncoupled analyses needed to support our flight test validation work. The SPARC team has continued to make progress on development activities to support unsteady, full reentry vehicle aero analysis, which have focused on turbulence modeling, uniform mesh refinement, in-situ visualization (with the Paraview/Catalyst team), and performance analysis. The SPARC team is making solid progress towards our Q4 goal of demonstrating an improved Wall Modeled Large Eddy Simulation (WMLES) capability for standard geometries using high-order finite difference, discontinuous Galerkin, and low-dissipation finite volume methods. The SPARC team has met the minimum completion criteria for parts of the L1, including documenting performance and scaling of SPARC on Trinity, Sierra and Astra, and performing runs for the flight test validation case.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Study of hypervelocity meteoroid impact on orbital space stations

Structural damage resulting in hypervelocity impact of a meteorite on a spacecraft is discussed. Of particular interest is the backside spallation caused by such a collision. To treat this phenomenon two numerical schemes were developed in the course of this study to compute the elastic-plastic flow fracture of a solid. The numerical schemes are a five-point finite difference scheme and a four-node finite element scheme. The four-node finite element scheme proved to be less sensitive to the type of boundary conditions and loadings. Although further development work is needed to improve the program versatility (generalization of the network topology, secondary storage for large systems, improving of the coding to reduce the run time, etc.), the basic framework is provided for a utilitarian computer program which may be used in a wide variety of situations. Analytic results showing the program output are given for several test cases.

Leimbach, K. R.↗
NASA NTRS1976

Multidimensional difference schemes with fourth-order accuracy

An explicit finite-difference algorithm is presented for the solution of quasilinear divergence free multidimensional hyperbolic systems. The method consists of four steps per time level. The resulting scheme is fourth-order accurate in both space and time, though the intermediate steps are only first-order accurate. The family of schemes introduced is dissipative, and hence, suitable for both smooth flows and flows containing shocks. This method is compared, in several numerical examples, with both second-order schemes and others that are fourth order in space, but second order in time.

Turkel, E.↗
NASA NTRS1987

SINDA-NASTRAN interfacing program theoretical description and user's manual

The task of converting SINDA finite difference thermal model temperature results into NASTRAN finite element model thermal loads can be very labor intensive if there is not one node-to-one element, or systematic node-to-element. correlation between models. This paper describes the SINDA-NASTRAN Interfacing Program (SNIP), a FORTRAN computer code that generates NASTRAN structural model thermal load cards given by SINDA (or similar thermal model) temperature results and thermal model geometric data. SNIP generates NASTRAN thermal load cards for NASTRAN plate, shell, bar, and beam elements. The paper describes the interfacing procedures used by SNIP, and discusses set-up and operation of the program. Sample cases are included to demonstrate use of the program and show its performance under a variety of conditions. SNIP can provide structural model thermal loads that accurately reflect thermal model results while reducing the time required to interface thermal and structural models when compared to other methods.

Winegar, Steven R.↗
NASA NTRS1990

Combining Thermal And Structural Analyses

Computer code makes programs compatible so stresses and deformations calculated. Paper describes computer code combining thermal analysis with structural analysis. Called SNIP (for SINDA-NASTRAN Interfacing Program), code provides interface between finite-difference thermal model of system and finite-element structural model when no node-to-element correlation between models. Eliminates much manual work in converting temperature results of SINDA (Systems Improved Numerical Differencing Analyzer) program into thermal loads for NASTRAN (NASA Structural Analysis) program. Used to analyze concentrating reflectors for solar generation of electric power. Large thermal and structural models needed to predict distortion of surface shapes, and SNIP saves considerable time and effort in combining models.

Winegar, Steven R.↗
NASA NTRS1989

A high-order Lagrangian-decoupling method for the incompressible Navier-Stokes equations

A high-order Lagrangian-decoupling method is presented for the unsteady convection-diffusion and incompressible Navier-Stokes equations. The method is based upon: (1) Lagrangian variational forms that reduce the convection-diffusion equation to a symmetric initial value problem; (2) implicit high-order backward-differentiation finite-difference schemes for integration along characteristics; (3) finite element or spectral element spatial discretizations; and (4) mesh-invariance procedures and high-order explicit time-stepping schemes for deducing function values at convected space-time points. The method improves upon previous finite element characteristic methods through the systematic and efficient extension to high order accuracy, and the introduction of a simple structure-preserving characteristic-foot calculation procedure which is readily implemented on modern architectures. The new method is significantly more efficient than explicit-convection schemes for the Navier-Stokes equations due to the decoupling of the convection and Stokes operators and the attendant increase in temporal stability. Numerous numerical examples are given for the convection-diffusion and Navier-Stokes equations for the particular case of a spectral element spatial discretization.

Ho, Lee-Wing↗
NASA NTRS1993

Hybrid techniques for complex aerospace electromagnetics problems

Important aerospace electromagnetics problems include the evaluation of antenna performance on aircraft and the prediction and control of the aircraft's electromagnetic signature. Due to the ever increasing complexity and expense of aircraft design, aerospace engineers have become increasingly dependent on computer solutions. Traditionally, computational electromagnetics (CEM) has relied primarily on four disparate techniques: the method of moments (MoM), the finite-difference time-domain (FDTD) technique, the finite element method (FEM), and high frequency asymptotic techniques (HFAT) such as ray tracing. Each of these techniques has distinct advantages and disadvantages, and no single technique is capable of accurately solving all problems of interest on computers that are available now or will be available in the foreseeable future. As a result, new approaches that overcome the deficiencies of traditional techniques are beginning to attract a great deal of interest in the CEM community. Among these new approaches are hybrid methods which combine two or more of these techniques into a coherent model. During the ASEE Summer Faculty Fellowship Program a hybrid FEM/MoM computer code was developed and applied to a geometry containing features found on many modern aircraft.

Aberle, Jim↗
NASA NTRS1994

A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models

A benchmark calculation is proposed for evaluating the dynamical cores of atmospheric general circulation models (GCMs) independently of the physical parameterizations. The test focuses on the long-term statistical properties of a fully developed general circulation; thus, it is particularly appropriate for intercomparing the dynamics used in climate models. To illustrate the use of this benchmark, two very different atmospheric dynamical cores--one spectral, one finite difference--are compared. It is found that the long-term statistics produced by the two models are very similar. Selected results from these calculations are presented to initiate the intercomparison.

Held, Isaac M.↗
NASA NTRS1996

Drag reduction in turbulent MHD pipe flows

This is a preliminary study devoted to verifying whether or not direct simulations of turbulent Magneto-Hydro-Dynamic (MHD) flows in liquid metals reproduce experimental observations of drag reduction. Two different cases have been simulated by a finite difference scheme which is second order accurate in space and time. In the first case, an external azimuthal magnetic field is imposed. In this case, the magnetic field acts on the mean axial velocity and complete laminarization of the flow at N(sub a) = 30 has been achieved. In the second case, an axial magnetic field is imposed which affects only fluctuating velocities, and thus the action is less efficient. This second case is more practical, but comparison between numerical and experimental results is only qualitative.

Orlandi, P.↗
NASA NTRS1996

Wing-Body Aeroelasticity on Parallel Computers

This article presents a procedure for computing the aeroelasticity of wing-body configurations on multiple-instruction, multiple-data parallel computers. In this procedure, fluids are modeled using Euler equations discretized by a finite difference method, and structures are modeled using finite element equations. The procedure is designed in such a way that each discipline can be developed and maintained independently by using a domain decomposition approach. A parallel integration scheme is used to compute aeroelastic responses by solving the coupled fluid and structural equations concurrently while keeping modularity of each discipline. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Aeroelastic computations are illustrated for a high speed civil transport type wing-body configuration.

Guruswamy, Guru P.↗
NASA NTRS1997

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗
DOE OSTI2020

Direct E&M Field Solver for a Plane Wave without Sources

A three dimensional direct field solver for the electric and magnetic fields using the Finite Difference Time Domain (FDTD) method is described. Example simulations using periodic and absorbing boundary conditions in a vacuum medium for a plane wave propagating in three dimensions are explored and are compared to theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Enhancing two-photon spontaneous emission in rare earths using graphene and graphene nanoribbons

The enhancement of two-photon spontaneous emission (2PSE) from trivalent and divalent rare earth ions in proximity to graphene and graphene nanoribbons is calculated for achievable experimental conditions using a combination of finite difference time domain simulations and direct computation of transition rates between energy levels in rare earths. For Er 3+ , we find that the 2PSE rate is initially 8 orders lower than the single-photon spontaneous emission rate but that, with enhancement, 2PSE can reach 2.5% of the overall decay. When graphene nanoribbons are used, we also show that the emission of free-space photon pairs from Er 3+ at 3–3.2 μ⁢m via 2PSE can be increased by ~400. Our calculations show significantly less relative graphene-enhanced 2PSE than previous works, and we attribute this variation to differences in emitter size and assumed graphene mobility. Here, we also show that the internal energy structure of the ion can have an impact on the degree of 2PSE enhancement achievable and find that divalent rare earths are more favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗