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 289 records · Page 16

Preliminary results of the analysis of CaII K spectroheliograms

The preliminary results of the photometry of CaII K spectroheliograms are presented. From the spectrograms for 1992, plages, the magnetic network, intranetwork elements and the chromospheric background were separated using the histogram method. The intensity and area of these separated features, as well as the full disk intensity, were derived. The spatial K index was compared to the spectral CaII K index derived from line profiles. It was found that the spatial K index and intensity of plages, the network elements and the intranetwork and background regions were highly correlated with the MgII h and k c/w ratio.

Kariyappa, R.↗

TU-AB-BRC-03: Accurate Tissue Characterization for Monte Carlo Dose Calculation Using Dual-and Multi-Energy CT Data

Purpose: To develop a general method for human tissue characterization with dual-and multi-energy CT and evaluate its performance in determining elemental compositions and the associated proton stopping power relative to water (SPR) and photon mass absorption coefficients (EAC). Methods: Principal component analysis is used to extract an optimal basis of virtual materials from a reference dataset of tissues. These principal components (PC) are used to perform two-material decomposition using simulated DECT data. The elemental mass fraction and the electron density in each tissue is retrieved by measuring the fraction of each PC. A stoichiometric calibration method is adapted to the technique to make it suitable for clinical use. The present approach is compared with two others: parametrization and three-material decomposition using the water-lipid-protein (WLP) triplet. Results: Monte Carlo simulations using TOPAS for four reference tissues shows that characterizing them with only two PC is enough to get a submillimetric precision on proton range prediction. Based on the simulated DECT data of 43 references tissues, the proposed method is in agreement with theoretical values of protons SPR and low-kV EAC with a RMS error of 0.11% and 0.35%, respectively. In comparison, parametrization and WLP respectively yield RMS errors of 0.13% and 0.29% on SPR, and 2.72% and 2.19% on EAC. Furthermore, the proposed approach shows potential applications for spectral CT. Using five PC and five energy bins reduces the SPR RMS error to 0.03%. Conclusion: The proposed method shows good performance in determining elemental compositions from DECT data and physical quantities relevant to radiotherapy dose calculation and generally shows better accuracy and unbiased results compared to reference methods. The proposed method is particularly suitable for Monte Carlo calculations and shows promise in using more than two energies to characterize human tissue with CT.

60 APPLIED LIFE SCIENCES↗

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

Efficient algorithms for use in probabilistic finite element analysis

This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.

Dias, J. B.↗

Enhancing Thomson scattering polychromator performance with multi-pass spectral filters

In photon-deficient, noncollective Thomson scattering diagnostics, filter polychromators are typically employed in the spectral analysis of Thomson-scattered signals to achieve acceptable signal-to-noise performance. Currently, the most common polychromator filter configuration employs a set of single-passband optical filters that define individual spectral channels. Here, we introduce a new spectral analysis method for Thomson scattering based on spectral filters with multiple passbands, referred to as Thomson scattering spectral multiplexing. Implementing multi-bandpass spectral filters on polychromators increases the achievable range of electron temperature measurement for a given number of filters employed. In addition, Thomson scattering spectral multiplexing reduces systematic measurement uncertainty, with fewer required spectral channels, thereby decreasing light loss from reduced optical element interactions. A multi-bandpass filter set, optimized by a genetic algorithm, has been successfully installed and tested on the Helically Symmetric eXperiment (HSX), demonstrating the benefits of the Thomson scattering spectral multiplexing method.

Instruments & Instrumentation↗

Preliminary Results of the Analysis

The preliminary results of the photometry of CaII K spectroheliograms taken at the NationalSolar Observatory at Sacramento peak are presented in this paper. We have digitizedspectroheliograms for 1980 (maximum of SC21), 1985 (minimum of SC21), 1987 (beginning of theascending phase of SC22), 1988 and 1989 (ascending phase and maximum of SC22), and 1992(declining phase of SC22). We have analyzed images for 1992 and separated the plages, the magneticnetwork, internetwork elements and the chromospheric background using histogram method. Wehave derived the intensity and area of these features as well as the full disk intensity (Spatial KIndex). The Spatial K Index has been compared to the spectral Ca K index derived from the lineprofiles and total solar and UV irradiance measured by the UARS and NOAA9 satellites. Thecontribution of plages, the magnetic network and internetwork element to the changes observed intotal solar and UV irradiances are also estimated.

CaII K Spectroheliograms↗

Quantifying spectral albedo effects on bifacial photovoltaic module measurements and system model predictions

We provide a comprehensive analysis of the effect of spectral albedo on photovoltaic (PV) module measurements and system model predictions. We demonstrate how to account for albedo in indoor bifacial device measurements by adjusting the applied irradiance using the scaled rear irradiance method, exemplified on fabricated silicon heterojunction (SHJ) modules. System model performance is studied using a detailed 3D finite-element model, DUET, for fixed-tilt and horizontal single-axis tracked (SAT) arrays between 15 and 75°N. Spectral effects cause variations in measured SHJ module short-circuit current up to 2% and efficiency variation up to 0.3% abs. We further demonstrate that rear-side spectral mismatch factors (SMMs) resulting from including or omitting spectral albedo in PV system modeling vary between ±13%, while total (front+rear) SMMs vary up to 3%, depending on the deployment configuration and latitude. SAT array SMMs are weakly correlated with latitude, while fixed-tilt array SMMs increase with latitude, driven by an increasing proportion of ground-reflected light on the front-side of modules. Ground-reflections can constitute between 2% and 32% of total incident module irradiance, with notably high (>10%) contributions for fixed-tilt arrays at high latitude. Effects of spectral albedo are most significant for: (1) fixed-tilt deployments at high latitudes, (2) wide bandgap technologies such as perovskite and cadmium telluride cells, (3) albedos which vary steeply over the technology's absorption range, and (4) high albedo ground covers. Overall, we demonstrate that omitting spectral albedo effects can result in PV measurement and system-level modeling uncertainties on the order of several percent in these cases.

14 SOLAR ENERGY↗

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗

Hierarchical model reduction driven by a proper orthogonal decomposition for parametrized advection-diffusion-reaction problems

This work combines the Hierarchical Model (HiMod) reduction technique with a standard Proper Orthogonal Decomposition (POD) to solve parametrized partial differential equations for the modeling of advection-diffusion-reaction phenomena in elongated domains (e.g., pipes). This combination leads to what we define as HiPOD model reduction, which merges the reliability of HiMod reduction with the computational efficiency of POD. Two HiPOD techniques are presented and assessed by an extensive numerical verification.

97 MATHEMATICS AND COMPUTING↗

Spectral-domain moment-method analysis of coplanar microstrip parasitic subarrays

Basic characteristics of several configurations of coplanar microstrip parasitic subarrays consisting of one fed patch and two or more parasitic patches were investigated by means of a spectral-domain full-wave analysis and the moment method analysis. Results are presented for radiating- and nonradiating edge-coupled three-element linear subarrays and for a five-patch cross. A comparison of the theoretical input impedance results obtained by the analysis of a three-element linear array showed a reasonable agreement between computed and measured R and X values.

Chen, Wei↗

Improved precision in As speciation analysis with HERFD-XANES at the As K -edge: the case of As speciation in mine waste

High-energy-resolution fluorescence-detected (HERFD) X-ray absorption near-edge spectroscopy (XANES) is a spectroscopic method that allows for increased spectral feature resolution, and greater selectivity to decrease complex matrix effects compared with conventional XANES. XANES is an ideal tool for speciation of elements in solid-phase environmental samples. Accurate speciation of As in mine waste materials is important for understanding the mobility and toxicity of As in near-surface environments. In this study, linear combination fitting (LCF) was performed on synthetic spectra generated from mixtures of eight measured reference compounds for both HERFD-XANES and transmission-detected XANES to evaluate the improvement in quantitative speciation with HERFD-XANES spectra. The reference compounds arsenolite (As 2 O 3 ), orpiment (As 2 S 3 ), getchellite (AsSbS 3 ), arsenopyrite (FeAsS), kaňkite (FeAsO 4 ·3.5H 2 O), scorodite (FeAsO 4 ·2H 2 O), sodium arsenate (Na 3 AsO 4 ), and realgar (As 4 S 4 ) were selected for their importance in mine waste systems. Statistical methods of principal component analysis and target transformation were employed to determine whether HERFD improves identification of the components in a dataset of mixtures of reference compounds. LCF was performed on HERFD- and total fluorescence yield (TFY)-XANES spectra collected from mine waste samples. Arsenopyrite, arsenolite, orpiment, and sodium arsenate were more accurately identified in the synthetic HERFD-XANES spectra compared with the transmission-XANES spectra. In mine waste samples containing arsenopyrite and either scorodite or kaňkite, LCF with HERFD-XANES measurements resulted in fits with smaller R -factors than concurrently collected TFY measurements. The improved accuracy of HERFD-XANES analysis may provide enhanced delineation of As phases controlling biogeochemical reactions in mine wastes, contaminated soils, and remediation systems.

58 GEOSCIENCES↗

Optimal Polynomial Smoothers and One‐Sided V‐Cycles for Poisson Problems

The solution to the Poisson equation arising from the spectral element discretization of the incompressible Navier‐Stokes equations needs robust preconditioning strategies. One such strategy is multigrid. To realize the potential of multigrid methods, effective smoothing strategies are needed. Chebyshev polynomial smoothers, in conjunction with pointwise Jacobi or additive Schwarz methods (ASMs), prove to be an effective smoother. Other polynomial smoothers, however, may provide superior convergence to the multigrid preconditioner. The authors compare the standard Chebyshev polynomial smoothers to both the novel fourth‐kind Chebyshev polynomial smoothers proposed by Lottes as well as smoothers based on the polynomial of best uniform approximation to as proposed by Kraus, Vassilevski, and Zikatanov. At the cost of symmetry, further improvements may be made. For example, a order polynomial smoother on both sides of the V‐cycle may be substituted with an order polynomial smoother on one side at no additional cost. The choice of omitting the postsmoother in favor of higher‐order polynomial presmoothing is advantageous in cases where the multigrid approximation property constant is large. The authors consider a 2D model problem based on finite differences to motivate the choice of polynomial smoother, order, and whether to apply postsmoothing for the target application of high‐order ‐geometric multigrid methods for GPU architectures. Results from both domains demonstrate the substantial improvement of these approaches over the standard Chebyshev polynomial smoother with a symmetric V‐cycle.

97 MATHEMATICS AND COMPUTING↗

Spectral methods for inviscid, compressible flows

Report developments in the application of spectral methods to two dimensional compressible flows are reviewed. A brief introduction to spectral methods -- their history and especially their implementation -- is provided. The stress is on those techniques relevant to transonic flow computation. The spectral multigrid iterative methods are discussed with application to the transonic full potential equation. Discontinuous solutions of the Euler equations are considered. The key element is the shock fitting technique which is briefly explained.

Hussaini, M. Y.↗

Spectral methods for inviscid, compressible flows

Report developments in the application of spectral methods to two dimensional compressible flows are reviewed. A brief introduction to spectral methods - their history and especially their implementation - is provided. The stress is on those techniques relevant to transonic flow computation. The spectral multigrid iterative methods are discussed with application to the transonic full potential equation. Discontinuous solutions of the Euler equations are considered. The key element is the shock fitting technique which is briefly explained.

Hussaini, M. Y.↗

X-Ray Fluorescence Solvent Detection at the Substrate-Adhesive Interface

With environmental regulations limiting the use of volatile organic compounds, low-vapor pressure solvents have replaced traditional degreasing solvents for bond substrate preparation. When used to clean and prepare porous bond substrates such as phenolic composites, low vapor pressure solvents can penetrate deep into substrate pore networks and remain there for extended periods. Trapped solvents can interact with applied adhesives either prior to or during cure, potentially compromising bond properties. Currently, methods for characterizing solvent time-depth profiles in bond substrates are limited to bulk gravimetric or sectioning techniques. While sectioning techniques such as microtome allow construction of solvent depth profiles, their depth resolution and reliability are limited by substrate type. Sectioning techniques are particularly limited near the adhesive-substrate interface where depth resolution is further limited by adhesive-substrate hardness and, in the case of a partially cured adhesive, mechanical properties differences. Additionally, sectioning techniques cannot provide information about lateral solvent diffusion. Cross-section component mapping is an alternative method for measuring solvent migration in porous substrates that eliminates the issues associated with sectioning techniques. With cross-section mapping, the solvent-wiped substrate is sectioned perpendicular rather than parallel to the wiped surface, and the sectioned surface is analyzed for the solvent or solvent components of interest using a two-dimensional mapping or imaging technique. Solvent mapping can be performed using either direct or indirect methods. With a direct method, one or more solvent components are mapped using red or Raman spectroscopy together with a moveable sample stage and/or focal plane array detector. With an indirect method, an elemental "tag" not present in the substrate is added to the solvent before the substrate is wiped. Following cross sectioning, the tag element can then be mapped by its characteristic x-ray emission using either x-ray fluorescence, or electron-beam energy-and wavelength-dispersive x-ray spectrometry. The direct mapping techniques avoid issues of different diffusion or migration rates of solvents and elemental tags, while the indirect techniques avoid spectral resolution issues in cases where solvents and substrates have adjacent or overlapping peaks. In this study, cross-section component indirect mapping is being evaluated as a method for measuring migration of d-limonene based solvents in glass-cloth phenolic composite (GCP) prior to and during subsequent bonding and epoxy adhesive cure.

Wurth, Laura↗

Preliminary design study. Shuttle modular scanning spectroradiometer

Fundamental concepts on which to base a detailed design for a Shuttle Modular Scanning Spectroradiometer were developed, and a preliminary design is presented. The recommended design features modularity and flexibility. It includes a 75-cm f/1.7-telescope assembly in an all-reflective Schmidt configuration, a solid state scan system (pushbroom) with high resolution over a 15 deg field of view, and ten detector channels covering the spectral range from 0.45 to 12.5 micrometers. It uses charge transfer device techniques to accommodate a large number of detector elements for earth observation measurements. Methods for in-flight radiometric calibration, for image motion compensation, and for data processing are described. Recommendations for ground support equipment are included, and interfaces with the shuttle orbiter vehicle are illustrated.

Source record↗

Stabilization of numerical interchange in spectral-element magnetohydrodynamics

In this study, auxiliary numerical projections of the divergence of flow velocity and vorticity parallel to magnetic field are developed and tested for the purpose of suppressing unphysical interchange instability in magnetohydrodynamic simulations. The numerical instability arises with equal-order C 0 finite- and spectral-element expansions of the flow velocity, magnetic field, and pressure and is sensitive to behavior at the limit of resolution. The auxiliary projections are motivated by physical field-line bending, and coercive responses to the projections are added to the flow-velocity equation. Their incomplete expansions are limited to the highest-order orthogonal polynomial in at least one coordinate of the spectral elements. Cylindrical eigenmode computations show that the projections induce convergence from the stable side with first-order ideal-MHD equations during h-refinement and p-refinement. Hyperbolic and parabolic projections and responses are compared, together with different methods for avoiding magnetic divergence error. Lastly, the projections are also shown to be effective in linear and nonlinear time-dependent computations with the NIMROD code [C. R. Sovinec, et al., J. Comput. Phys. 195 (2004) 355-386], provided that the projections introduce numerical dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Absolute nuclear charge radius by Na-like spectral line separation in high-Z elements

Abstract We describe a novel technique to determine absolute nuclear radii of high- Z nuclides. Utilizing accurate theoretical atomic structure calculations together with precise measurements of extreme ultraviolet transitions in highly charged ions this method allows for precise determinations of absolute nuclear charge radii based upon the well-known nuclear radii of their neighboring elements. This method can work for elements without stable isotopes, and its accuracy may be competitive with current methods (electron scattering and muonic x-ray spectroscopy).

Hosier, A. (ORCID:0000000328366605)↗