Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A comprehensive simulation model of the performance of photochromic films in absorbance-modulation-optical-lithography

Optical lithography is the most prevalent method of fabricating micro-and nano-scale structures in the semiconductor industry due to the fact that patterning using photons is fast, accurate and provides high throughput. However, the resolution of this technique is inherently limited by the physical phenomenon of diffraction. Absorbance-Modulation-Optical Lithography (AMOL), a recently developed technique has been successfully demonstrated to be able to circumvent this diffraction limit. AMOL employs a dual-wavelength exposure system in conjunction with spectrally selective reversible photo-transitions in thin films of photochromic molecules to achieve patterning of features with sizes beyond the far-field diffraction limit. We have developed a finite-element-method based full-electromagnetic-wave solution model that simulates the photo-chemical processes that occur within the thin film of the photochromic molecules under illumination by the exposure and confining wavelengths in AMOL. This model allows us to understand how the material characteristics influence the confinement to sub-diffraction dimensions, of the transmitted point spread function (PSF) of the exposure wavelength inside the recording medium. The model reported here provides the most comprehensive analysis of the AMOL process to-date, and the results show that the most important factors that govern the process, are the polarization of the two beams, the ratio of the intensities of the two wavelengths, the relative absorption coefficients and the concentration of the photochromic species, the thickness of the photochromic layer and the quantum yields of the photoreactions at the two wavelengths. The aim of this work is to elucidate the requirements of AMOL in successfully circumventing the far-field diffraction limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Modeling Resonances of Dielectric Spheres in the Terahertz Range

We have studied the resonance properties of dielectric spheres with subwavelength dimensions in the optical and terahertz frequency range. Distribution of the electric field vector in a sphere irradiated by a plane polarized wave has been numerically simulated using the finite element method. The obtained results clearly demonstrate that it is in principle possible to use resonances of various orders in dielectric spheres for the amplification of current in a conducting antenna with dimensions several times smaller than the radiation wavelength. In the case of resonances within short spectral intervals, it is possible to provide a current density gain up to 25–30 dB as compared to that in the absence of a dielectric sphere surrounding the antenna.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1975

Optimization of structures undergoing harmonic or stochastic excitation

The optimal design was investigated of simple structures subjected to dynamic loads, with constraints on the structures' responses. Optimal designs were examined for one dimensional structures excited by harmonically oscillating loads, similar structures excited by white noise, and a wing in the presence of continuous atmospheric turbulence. The first has constraints on the maximum allowable stress while the last two place bounds on the probability of failure of the structure. Approximations were made to replace the time parameter with a frequency parameter. For the first problem, this involved the steady state response, and in the remaining cases, power spectral techniques were employed to find the root mean square values of the responses. Optimal solutions were found by using computer algorithms which combined finite elements methods with optimization techniques based on mathematical programming. It was found that the inertial loads for these dynamic problems result in optimal structures that are radically different from those obtained for structures loaded statically by forces of comparable magnitude.

Johnson, E. H.↗
NASA NTRS1976

Biplex, pipelined FFT

A method of implementing a pipelined fast Fourier transform (FFT) is presented that makes full use of the elements that comprise it. For a given bandwidth and spectral resolution, this approach requires 25 percent less memory, and needs only half the logic speed required for the conventional implementation. Since the spectrum is broken into two parts, two analog-to-digital converters, and attendant mixer-filters, are needed. These two input elements need only operate at half the speed of the one they replaced.

Emerson, R. F.↗
NASA NTRS1993

High-order cyclo-difference techniques: An alternative to finite differences

The summation-by-parts energy norm is used to establish a new class of high-order finite-difference techniques referred to here as 'cyclo-difference' techniques. These techniques are constructed cyclically from stable subelements, and require no numerical boundary conditions; when coupled with the simultaneous approximation term (SAT) boundary treatment, they are time asymptotically stable for an arbitrary hyperbolic system. These techniques are similar to spectral element techniques and are ideally suited for parallel implementation, but do not require special collocation points or orthogonal basis functions. The principal focus is on methods of sixth-order formal accuracy or less; however, these methods could be extended in principle to any arbitrary order of accuracy.

Carpenter, Mark H.↗
NASA NTRSDate not supplied

Modal Propellant Gauging: High-resolution and non-invasive gauging of both settled and unsettled liquids in reduced gravity.

The modal response of a liquid-filled tank to external acoustic excitation can be used to infer with high resolutionthe mass of contained liquid, the mass flow rate of liquids into and out of the tank, and changes in tank pressure.Both contained liquid mass and internal ullage pressure affect the modal response of the tank walls through fluidmass-loading of the tank walls and pressure-induced wall stiffening, respectively. Modal Propellant Gauging refersto the technology that exploits these shifts in modal frequencies to infer the mass of propellant in a tank. MPG is anon-invasive gauging technology that has demonstrated gauging resolutions of 1% for settled propellants and 2-3%for unsettled, sloshing propellants. Extensive parabolic flight testing of the MPG system on model tanks has beenconducted to validate the technology in reduced gravity. MPG testing on a qualification tank for the Orion Program'sEuropean Service Module has also been conducted and is reported here. Finite element modeling of the OrionESM "upper" tank is discussed and compared with measurement data. Three computational approaches to massdetermination, Peak Tracking, Point Sensor, and Spectral Density methods, are described here. Use cases are definedand analyzed in the context of the Orion ESM Qualification tank data, and an implementation scheme for continuousmass gauging on the Orion ESM is discussed.Keywords: Propellant Mass Gauging, Modal Analysis, Orion ESM

Propellant Mass Gauging↗
NASA NTRS2016

Extracting Damping Ratio from Dynamic Data and Numerical Solutions

There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.

Casiano, M. J.↗
DOE OSTI2024

Optical encoder devices and systems

Devices, systems and methods for encoding information using optical components are described. An example photonic filtered sampler includes a spectral shaper configured to receive an optical pulse train, a dispersive element positioned to receive an output of the spectral shaper and to expand spectral contents thereof in time, and a modulator configured to receive an output of the dispersive element and a radio frequency (RF) signal, and to produce a modulated output optical signal in accordance with the RF signal. In this configuration, one or more characteristics of the modulated output optical signal is determined based on a spectral shape provided by the spectral shaper and dispersive properties of the dispersive element.

Buckley, Brandon Walter↗
DOE OSTI2017

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Spatial-Spectral Representation for X-Ray Fluorescence Image Super-Resolution

X-ray fluorescence (XRF) scanning of works of art is becoming an increasing popular nondestructive analyticalmethod. The high quality XRF spectra is necessary to obtain significant information on both major and minor elements used for characterization and provenance analysis. However, there is a tradeoff between the spatial resolution of an XRF scan and the signalto- noise ratio (SNR) of each pixel’s spectrum, due to the limited scanning time. In this project, we propose an XRF image super-resolution method to address this tradeoff; thus, obtaining a high spatial resolution XRF scan with high SNR. We fuse a lowresolution XRF image and a conventional RGB high-resolution image into a product of both high spatial and high spectral resolution XRF image. There is no guarantee of a one to one mapping between XRF spectrum and RGB color since, for instance, paintings with hidden layers cannot be detected in visible but can in X-ray wavelengths. We separate the XRF image into the visible and nonvisible components. The spatial resolution of the visible component is increased utilizing the high-resolution RGB image, whereas the spatial resolution of the non-visible component is increased using a total variation super-resolution method. Finally, the visible and nonvisible components are combined to obtain the final result.

Dai, Qiqin↗
DOE OSTI2022

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗
NASA NTRS1972

Automatic classification of soils and vegetation with ERTS-1 data

Preliminary results of a test of a computerized analysis method using ERTS 1 data are presented. The method consisted of a four-spectral-band supervised, maximum likelihood, Gaussian classifier with training statistics derived through a combination of clustering and manual methods. The multivariate analysis method leads to the assignment of each resolution element of the data to one of a preselected set of discrete classes. The data frame was an area over the Texas-Oklahoma border including Lake Texoma. The study suggests that multispectral scanner data coupled with machine processing shows promise for earth surface cover surveys. Futhermore, the processing time is short and consequently the costs are low; a full frame can be analyzed completely within 48 hours.

Landgrebe, D. A.↗
NASA NTRS1991

Hierarchical classifier design in high-dimensional, numerous class cases

As progress in new sensor technology continues, increasingly high spectral resolution sensors are being developed. These sensors give more detailed and complex data for each picture element and greatly increase the dimensionality of data over past systems. Three methods for designing a decision tree classifier are discussed; a top down approach, a bottom up approach, and a hybrid approach. Three feature extraction techniques are implemented. Canonical and extended canonical techniques are mainly dependent on the mean difference between two classes. An autocorrelation technique is dependent on the correlation differences. The mathematical relationship among sample size, dimensionality, and risk value is derived.

Kim, Byungyong↗
NASA NTRS1993

A priori error estimates for an hp-version of the discontinuous Galerkin method for hyperbolic conservation laws

A priori error estimates are derived for hp-versions of the finite element method for discontinuous Galerkin approximations of a model class of linear, scalar, first-order hyperbolic conservation laws. These estimates are derived in a mesh dependent norm in which the coefficients depend upon both the local mesh size h(sub K) and a number p(sub k) which can be identified with the spectral order of the local approximations over each element.

Bey, Kim S.↗
DOE OSTI2019

Optical properties of laser-modified diamond: From visible to microwave range

We have measured the optical properties of a graphitised layer with dc conductivity of 305 Ω{sup −1} cm{sup −1} produced by excimer laser irradiation (KrF, τ = 20 ns, λ = 248 nm) of a polycrystalline diamond surface grown by chemical vapour deposition. Studies have been conducted in a wide range of wavelengths: from far-IR to visible spectral regions. The resulting constants of the Drude model are the basic parameters for calculating the electromagnetic response of diamond/graphite photon elements or metamaterials fabricated using the direct laser writing method. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Helium abundance of the sun.

Critical review of the four methods that have been used to estimate the ratio of helium to hydrogen in the sun - the solar neutrino flux, spectral intensity of helium lines in prominences and the chromosphere, elemental abundance of solar cosmic rays, and variations of solar wind He/H. At present the measured value for the neutrino flux cannot be explained by any of the theoretical models of the sun, and thus the neutrino measurements cannot yield an estimate of He/H. Helium line intensity measurements give the most accurate values for He/H now available - i.e., 5 to 8% helium by number. However, these results refer to the chromosphere and prominences only, and questions of spatial and temporal variations of the observed quantities and of the inferred helium abundances have not yet been dealt with adequately.

Hirshberg, J.↗
NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗
NASA NTRSDate not supplied

Quick and Correlative TOF-SIMS Analysis of Dispersoid Content in Powder Feedstock and Printed Oxide Dispersion Strengthened Alloys

Oxide dispersion strengthening (ODS) of metal alloys, when combined with additive manufacturing (AM), can be used to solve high-temperature material application problems. The performance of ODS materials is dependent on nano-scale oxide being distributed uniformly throughout the metal matrix; a homogenous distribution is critical to the creep and yield strength. Therefore, quantifying and analyzing these oxide nanoparticles is essential for ODS alloy development. Assessing the chemistry and morphology of dispersoids in a material can be conducted using a variety of X-ray and electron-based methods. Time of Flight Secondary Ion Mass Spectrometry (TOF-SIMS), which assesses chemical concentrations from ion masses, provides the necessary spatial resolution and avoids potential spectral overlap of all elements present. In this work, preparing ODS powder feedstock and Laser Powder Bed Fusion AM was conducted with a NiCoCr alloy and yttria dispersoid, mixed with a Resodyn resonance mixer and printed on an EOS M100 3D printer at NASA Glenn. TOF-SIMS was performed on a Physical Electronics NanoTOF TRIFT V TOF-SIMS at Case Western Reserve University. The TOF-SIMS results collected on the feedstock powder samples correspond to the qualitative interpretation of the homogeneity of the oxide coating on the metal powder from SEM images. This provides quantitative support to the qualitative microscopy images collected earlier in the process design phase. The variance in the yttrium distribution on the powder feedstock samples also correlated with the variance of the yttrium distribution detected in the printed samples. This correlation supported the hypothesis that an even coating of yttria on the metal powders contributes to an even distribution of yttria in the printed material.

Laura G. Wilson↗