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 397 records · Page 22

Flux Cube Reconstruction from Slitless Spectroscopy

Slitless spectroscopy enables efficient, large-area surveys without target preselection, yet it faces challenges from source blending, higher noise, and lost spatial–spectral information. We present an advanced, nonparametric, data-driven algorithm that leverages multiple dispersion angles to reconstruct three-dimensional flux distributions, providing low-resolution integral field unit capabilities from slitless data. By treating each pixel as an independent element, our method naturally handles source confusion without requiring prior assumptions regarding redshifts, templates, or model libraries. We validate the algorithm using simulated Roman Space Telescope wide-field slitless spectroscopy images that are equivalent to what is expected from the High-Latitude Time-Domain Survey. First, we demonstrate that a host-galaxy model reconstructed from multiple dispersion angles can be used to accurately subtract host light from a transient, recovering a Type Ia supernova spectrum with minimal bias. Second, we showcase a high-fidelity flux-cube reconstruction of a complex galaxy, successfully measuring the redshift and recovering continuum, emission, and absorption features. This approach highlights the potential of multi-dispersion-angle slitless data to provide spatially resolved spectral information in a nonparametric way, which is traditionally accessible only with integral field spectroscopy, opening a new window into large, unbiased, and spatially resolved studies of galaxy evolution.

Griggio, M. [Space Telescope Science Institute, Ba↗

Heat Transfer Measurements on Surfaces with Natural Ice Castings and Modeled Roughness

An experimental method is described to measure and compare the convective heat transfer coefficient of natural and simulated ice accretion roughness and to provide a rational means for determining accretion-related enhanced heat transfer coefficients. The natural ice accretion roughness was a sample casting made from accretions at the NASA Lewis Icing Research Tunnel (IRT). One of these castings was modeled using a Spectral Estimation Technique (SET) to produce three roughness elements patterns that simulate the actual accretion. All four samples were tested in a flat-plate boundary layer at angle of attack in a "dry" wind tunnel test. The convective heat transfer coefficient was measured using infrared thermography. It is shown that, dispite some problems in the current data set, the method does show considerable promise in determining roughness-induced heat transfer coefficients, and that, in addition to the roughness height and spacing in the flow direction, the concentration and spacing of elements in the spanwise direction are important parameters.

Breuer, Kenneth S.↗

Changes in SSL Device Efficiency and Optical Performance Under Accelerated Aging Conditions

Lighting application efficiency (LAE) describes the efficient delivery of light from the light source to the lighted task and is viewed as a new frontier—increasing energy savings with solid-state lighting (SSL) technologies. The framework for LAE that is proposed by the U.S. Department of Energy (DOE) consists of four major elements: light source efficiency, optical delivery efficiency, spectral efficiency, and intensity effectiveness. This report focuses on a sampling of the available SSL products that can be broadly defined as having modified spectral output because the method of spectra modification has a significant impact on light source efficiency and long-term optical delivery and spectral efficiencies. This report focused on the changes in light source, optical delivery, and spectral efficiencies that occur during aging of SSL devices. An accelerated stress test (AST) regiment was developed for the devices under test (DUTs) examined in this report to provide insights into how the performance of SSL devices change with aging. The AST protocols demonstrate that the SSL products discussed in this report often reduce light source efficiency to achieve different spectral characteristics. In addition, aging of the optical components (e.g., lenses, solder masks) in SSL devices can produce increased light absorption, which negatively impacts optical delivery efficiency and light source efficiency. The modified spectral outputs of the selected SSL products discussed in this report come in a variety of form factors and achieve enhanced optical performance by using a variety of methods. All of the products examined in this study use mid-power light-emitting diodes (MP-LEDs), although the number, configuration, phosphor content, and light-emitting diode (LED) pump of the MP-LEDs differ. Product MS-1 is a 60-watt (W) replacement A19 lamp with a hermetically sealed glass globe, which contains an embedded optical filter to absorb green and yellow emissions, thereby creating a “sunlike” modified spectrum. Product MS-2 is a 60-W replacement A19 lamp with 30 MP-LEDs, and it uses a violet LED pump, along with green and red phosphor emissions, to produce a “healthy” spectrum. This spectrum omits blue emissions in an effort to reduce melanopic lux. Product MS-3 is an LED module consisting of 21 MP-LEDs that use a violet LED pump, along with blue, green, and red phosphors, to produce a “sunlike” spectrum. Products MS-4 and MS-5 are both 6-inch (in) downlights that use a manual switching mechanism so that users can select application-specific correlated color temperatures before installation. Products MS-4 and MS-5 both contain two LED primaries (2,700 Kelvin [K] and 5,000 K) for spectral tuning. Product MS-4 contains 12 MP-LEDs for each LED primary, and Product MS-5 contains 10 MP-LEDs for each LED primary. This report summarizes the overall findings from up to 8,000 hours (hrs) of AST on the lamp DUTs (Products MS-1 and MS-2); up to 5,000 hrs of AST on the LED light engine DUTs (Product MS-3); and up to 7,000 hrs of AST on the downlight DUTs (Products MS-4 and MS-5). The AST procedures used in this study included a room temperature operational life (RTOL) test, an operational life test conducted at 45 degrees Celsius (°C; 45OL) test, an operational life test conducted at 75°C (75OL), a wet high-temperature operational life test performed at 65°C and 90% relative humidity (6590), and a wet high-temperature operational life test performed at 75°C and 75% relative humidity (7575). The AST procedures used for Products MS-1 and MS-2 were RTOL, 45OL, and 6590. The AST procedures used for Products MS-3, MS-4, and MS-5 were RTOL, 75OL, and 7575. During the ASTs described herein, separate populations of each product (three DUTs in each population for Products MS-1, MS-2, MS-4, and MS-5; four DUTs in each population for Product MS-3) were subjected to power cycling of 1 hr on and 1 hr off. The key findings from this study include the following. Enhanced optical performance came at the cost of reduced light source efficiency for the lamps and light engines examined in this study. The optical filter used to produce a “sunlike” spectrum for Product MS-1 reduced light source efficiency by 26%, from 113 lumens per watt (lm/W) to 85 lm/W. Products MS-2 and MS-3 that used a violet-pumped LED to achieve “healthy” and “sunlike” enhanced optical performance, respectively, suffered the largest reduction in initial light source efficiency (49 lm/W for Product MS-2 and 68 lm/W for Products MS-3) perhaps because of the use of violet LED as the optical pump. Product MS-2 also had the poorest color fidelity because of undersaturation of blues and oversaturation of greens and yellows. Chromaticity maintenance of the products in this study was generally good, with parametric failure only occurring for one product at one AST condition (i.e., Product MS-3 at 7575). The chromaticity shift for Product MS-3 in 7575 test conditions resulted from reduced emissions of the broad green and red phosphors used to mimic sunlight in the 500–750 nanometer (nm) range. As a result of these phosphor emission losses, chromaticity shifted toward the more stable violet-pumped LED and blue phosphor. Although chromaticity maintenance was acceptable for most products tested, different chromaticity shift mechanisms were observed for temperature and humidity tests compared with temperature alone for Products MS-1, MS-2, and MS-3. For Product MS-1, a relative increase in green emissions was observed with humidity, which may indicate humidity-accelerated degradation of the optical filter at green wavelengths (optical delivery efficiency reduction) or degradation of the emitters in the red region (spectral efficiency reduction). The violet-pumped products (i.e., MS-2 and MS-3) exhibited different chromaticity shift behavior in the temperature-humidity environments, with Product MS-2 shifting generally yellow because of photo-induced oxidation of the plastic globe (i.e., a reduction of optical delivery efficiency) and Product MS-3 shifting toward violet and blue emitters because of faster loss of emission from green and red phosphors (i.e., a reduction in spectral efficiency). Because of an initial drop in power consumption, the light source efficiency of downlights (Products MS-4 and MS-5) increased at RTOL throughout the test duration. The light source efficiency initially increased during 75OL until a decay in luminous flux maintenance (LFM) dominated light source efficiency, whereas the light source efficiency decreased for the entire 7575 test duration. It was found that increasing the ambient environment of an SSL device from 25°C to 75°C decreases the LFM by 2.7 to 5.5, depending on design of the SSL device. Adding humidity to the system (75OL to 7575) was found to decrease LFM by another factor of approximately 3.0. A common location of failure was identified for the downlights (Products MS-4 and MS-5) operated in the 7575 environment. This failure abruptly occurred between 3,000 to 5,000 hrs of testing on the film capacitor of the electromagnetic interference (EMI) filter near the diode bridge for 11 of the 12 downlights. It is likely that this failure was caused by a combination of the high electrical voltage across the capacitor, the location of the capacitor near the heat sources (e.g., transformers, diode bridges, power resistors), and the high stress environment of 7575. Longer test times are needed to fully understand the LFM, luminous efficacy changes, and chromaticity shifts for some of the products. Because of the high reliability, light source efficiency, and spectral tuning capabilities of LEDs, the expectations of LAE for SSL products are greatly increased over traditional lighting products. The data gathered in this report begin to provide an understanding of the tradeoffs that current SSL products undergo when optimizing the different efficiency elements of LAE. As discussed in this report, the findings from the tests conducted show that the optimization of spectral efficiency can come at the cost of initial light source efficiency. Furthermore, the introduction of violet LEDs can promote long-term optical delivery efficiency degradation, and the introduction of optical filters or new phosphors can lead to unwanted spectral efficiency changes because the ratio of phosphor emitters changes with aging. The data presented in this report also identified a common failure location in 6-in downlights. The results regarding long-term behavior of the modified spectra devices studied provide valuable information about changes to the light source, spectral, and optical delivery efficiencies as the devices age. This information can be used to improve future SSL designs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Expansion dynamics and chemistry evolution in ultrafast laser filament produced plasmas

Laser ablation in conjunction with optical emission spectroscopy is a potential non-contact, standoff detection method for all elements in the periodic table and certain isotopes such as radionuclides. Currently, significant development efforts are on-going to use ultrafast laser filaments for remote sensing. The application of filaments is of particular interest in extending the range of stand-off capability associated with elemental and isotopic detection via laser-induced breakdown spectroscopy. In this study, we characterize the expansion dynamics and chemical evolution of filament-produced uranium (U) plasmas. Laser filaments are generated in the laboratory by loosely focusing 35 femtosecond (fs), 6 milli Joule (mJ) pulses in air. Two-dimensional spectrally-integrated and time-resolved imaging was performed to study hydrodynamics and evolution of U atomic and UO molecular emission in filament-produced U plasmas. Our results highlight that filament ablation of U plasmas gives a cylindrical plasma plume morphology with an appearance of plume splitting into slow and fast moving components at later times of its evolution. Emission from the slow-moving region shows no distinct spectral features (i.e. broadband-like) and is contributed in part by nanoparticles generated during ultrafast laser ablation. Additionally, we find U atoms and U oxide molecules (i.e. UO, UxOy) co-exist in the filament produced plasma, which can be attributed to the generation of low-temperature plasma conditions during filament ablation

(140.3440) Laser-induced breakdown, (280.5395) Pla↗

Measuring interstellar dust grains from the haloes of binary X-ray sources

A method for determining the location, size distribution, and composition of interstellar dust grains from the haloes of binary X-ray sources is described. The method is based on the time dependence of the halo X-ray emission during an eclipse of the point source. The distance of the dust cloud is inferred directly from observations of the growth of the dark circle around the eclipsed point source. Grain size distribution in the interstellar dust clouds is derived from the halo brightness profile at a given photon energy. The composition of the grains is obtained by observing the spectral features in the total optical depth due to X-ray scattering near the photoionization thresholds of the elements in the grains. The haloes become dim just below the ionization threshold, and sharply increase in function brightness as the photon energy increases beyond the absorption edge of the relevant atom. The method makes it possible to observe and analyze the halo emission during eclipse even without imaging the source. It is suggested that haloes probably account for the soft X-rays observed during eclipses and dips observed in sources including 4U1915-05 and X1755-338.

Xu, Y.↗

Deciphering a removed text on a 16th c. centrepiece: The benefits of using AGLAEMap and AXSIA to process PIXE maps.

In the field of archaeometry, it is not uncommon to be presented with art objects that contain inscriptions, signatures and other writings that are nearly impossible to read. Scanned microbeam PIXE offers an attractive approach to attack this problem, but even then the distribution of characteristic X-rays of the element(s) used in these writings can remain illegible. We show in this paper that two methods were used to reveal the inscription: first the use of a GUPIXWin, TRAUPIXE and AGLAEMap software suite enables to make quantitative analysis of each pixel, to visualize the results and to select X-ray peaks that could enable to distinguish letters. Then, the Automated eXpert Spectral Image Analysis (AXSIA) program developed at Sandia, which analyzes the x-ray intensity vs. Energy and (X, Y) position “datacubes”, was used to factor the datacube into 1) principle component spectral shapes and 2) the weighting images of these components. The specimen selected for this study was a silver plaque representing a scroll from the so-called “MerkelscheTafelaufsatz,” a centrepiece made by the Nuremberg goldsmith Wenzel Jamnitzer in 1549. X-ray radiography of the plaque shows lines of different silver thicknesses, meaning that a text has been removed. The PIXE analysis used a 3-MeV proton beam focused to 50μm and scanned across the sample on different areas of interest of several cm². This analysis showed major elements of Cu and Ag, and minor elements such as Pb, Au, Hg. X-ray intensity maps were then made by setting windows on the various x-ray peaks but the writing on the centrepiece was not revealed even if the map of Cu after data treatment at AGLAE enabled to distinguish some letters. The AXSIA program enabled to factor two main spectral shapes from the datacube that were quite similar and involved virtually all of the X-rays being generated. Nevertheless, small differences between these factors were observed for the Cu K X-rays, Pb, Bi and Au L X-rays. The plot of the factor with the highest Au signal gave also information on the shape of some letters. The comparison of the results obtained by the two methods shows that they both drastically improve the resolution and contrast of such writings and that each of the method can also bring different information on the composition and thus the techniques used for the writing.

AGLAEMap↗

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗

Surface-plasmon enhanced photodetection at communication band based on hot electrons

Surface plasmons can squeeze light into a deep-subwavelength space and generate abundant hot electrons in the nearby metallic regions, enabling a new paradigm of photoconversion by the way of hot electron collection. Unlike the visible spectral range concerned in previous literatures, we focus on the communication band and design the infrared hot-electron photodetectors with plasmonic metal-insulator-metal configuration by using full-wave finite-element method. Titanium dioxide-silver Schottky interface is employed to boost the low-energy infrared photodetection. The photodetection sensitivity is strongly improved by enhancing the plasmonic excitation from a rationally engineered metallic grating, which enables a strong unidirectional photocurrent. With a five-step electrical simulation, the optimized device exhibits an unbiased responsivity of ∼0.1 mA/W and an ultra-narrow response band (FWHM = 4.66 meV), which promises to be a candidate as the compact photodetector operating in communication band.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Advances in scientific literature mining for interpreting materials characterization

Abstract Using synchrotron light sources, such as the National Synchrotron Light Source II at Brookhaven National Laboratory, scientists in fields as diverse as physics, biology, and materials science, identify the atomic structure, chemical composition, or other important properties of varied specimens. x-ray spectroscopy from light sources is particularly valuable for materials research with vast information available about reference spectra in the scientific literature. However, as the technique is applicable to many science domains, searching for information about select x-ray spectroscopy spectra is impeded by the sheer number of publications. Moreover, useful information about the context of an experiment or figures presented in papers can be buried among the details, which takes time to assess. This work presents a scientific literature mining system that supports data acquisition, information extraction, and user interaction for referencing x-ray spectra identification and spectral interpretation. The goal is to provide efficient access to useful spectral data to researchers who may spend only a few days at a synchrotron light source. With this system, users browse a classification tree for papers arranged according to x-ray spectroscopic methods, chemical elements, and x-ray absorption spectroscopy edges. Relevant figures are extracted with sentences from the paper that explain them, known as ‘figure explanatory text.’ Notably, this system focuses on semantic aspects (logical analysis) to find figure explanatory text using deep contextualized word embeddings techniques and contains an interface to obtain labeled data from domain experts that is used to evaluate and improve the model.

Park, Gilchan (ORCID:0000000201536646)↗

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗

Implementation of stacked ensemble machine learning for the detection of surrogate plutonium contamination in soil via LIBS

Supervised machine learning methods have demonstrated increased utility for the quantification of lanthanide and actinide elements in atomic spectroscopy applications. This study implements laser-induced breakdown spectroscopy (LIBS) for the identification of plutonium surrogate material (CeO 2 ) in soil matrices by training supervised machine learning methods on the recorded spectral data. A bagged ensemble using Random Forest yields the highest sensitivity predictions with a detection limit of 0.015 wt.% CeO 2 . However, high precision in Ce content prediction required the use of a stacked ensemble regression, which provided the superlative Ce quantification model with an error of 0.107% and a detection limit of 0.022 wt.%. Furthermore, the high performance of the stacked ensemble demonstrates its potential to enhance the accuracy and sensitivity of nuclear contaminant detection using field-deployable spectroscopic analyzers in real-world scenarios.

47 OTHER INSTRUMENTATION↗

Time-Resolved Optical Pump-Resonant X-ray Probe Spectroscopy of 4-Thiouracil: A Simulation Study

We theoretically monitor the photoinduced ππ* → nπ* internal conversion process in 4-thiouracil (4TU), triggered by an optical pump. The element-sensitive spectroscopic signatures are recorded by a resonant X-ray probe tuned to the sulfur, oxygen, or nitrogen K-edge. We employ high-level electronic structure methods optimized for core-excited electronic structure calculation combined with quantum nuclear wavepacket dynamics computed on two relevant nuclear modes, fully accounting for their quantum nature of nuclear motions. We critically discuss the capabilities and limitations of the resonant technique. For sulfur and nitrogen, we document a pre-edge spectral window free from ground-state background and rich with ππ* and nπ* absorption features. The lowest sulfur K-edge shows strong absorption for both ππ* and nπ*. In the lowest nitrogen K-edge window, we resolve a state-specific fingerprint of the ππ* and an approximate timing of the conical intersection via its depletion. A spectral signature of the nπ* transition, not accessible by UV–vis spectroscopy, is identified. The oxygen K-edge is not sensitive to molecular deformations and gives steady transient absorption features without spectral dynamics. Here, the ππ*/nπ* coherence information is masked by more intense contributions from populations. Altogether, element-specific time-resolved resonant X-ray spectroscopy provides a detailed picture of the electronic excited-state dynamics and therefore a sensitive window into the photophysics of thiobases.

97 MATHEMATICS AND COMPUTING↗

The virtual element method for eigenvalue problems with potential terms on polytopic meshes

In this work, we extend the conforming virtual element method (VEM) to the numerical resolution of eigenvalue problems with potential terms on a polytopic mesh. An important application is that of the Schrödinger equation with a pseudopotential term. This model is a fundamental element in the numerical resolution of more complex problems from the Density Functional Theory. The VEM is based on the construction of the discrete bilinear forms of the variational formulation through certain polynomial projection operators that are directly computable from the degrees of freedom. The method shows a great flexibility with respect to the meshes and provides a correct spectral approximation with optimal convergence rates. This point is discussed from both the theoretical and the numerical viewpoint. The performance of the method is numerically investigated by solving the quantum harmonic oscillator problem with the harmonic potential and a singular eigenvalue problem with zero potential for the first eigenvalues.

97 MATHEMATICS AND COMPUTING↗

Combining Astrometry and Elemental Abundances: The Case of the Candidate Pre-Gaia Halo Moving Groups G03-37, G18-39, and G21-22

While most moving groups are young and nearby, a small number have been identified in the Galactic halo. Understanding the origin and evolution of these groups is an important piece of reconstructing the formation history of the halo. Here we report on our analysis of three putative halo moving groups: G03-37, G18-39, and G21-22. Based on Gaia EDR3 data, the stars associated with each group show some scatter in velocity (e.g., Toomre diagram) and integrals of motion (energy, angular momentum) spaces, counter to expectations of moving-group stars. We choose the best candidate of the three groups, G21-22, for follow-up chemical analysis based on high-resolution spectroscopy of six presumptive members. Using a new Python code that uses a Bayesian method to self-consistently propagate uncertainties from stellar atmosphere solutions in calculating individual abundances and spectral synthesis, we derive the abundances of α- (Mg, Si, Ca, Ti), Fe-peak (Cr, Sc, Mn, Fe, Ni), odd-Z (Na, Al, V), and neutron-capture (Ba, Eu) elements for each star. We find that the G21-22 stars are not chemically homogeneous. Based on the kinematic analysis for all three groups and the chemical analysis for G21-22, we conclude the three are not genuine moving groups. The case for G21-22 demonstrates the benefit of combining kinematic and chemical information in identifying conatal populations when either alone may be insufficient. Comparing the integrals of motion and velocities of the six G21-22 stars with those of known structures in the halo, we tentatively associate them with the Gaia-Enceladus accretion event.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Physics–Dynamics Coupling with Element-Based High-Order Galerkin Methods: Quasi-Equal-Area Physics Grid

Atmospheric modeling with element-based high-order Galerkin methods presents a unique challenge to the conventional physics–dynamics coupling paradigm, due to the highly irregular distribution of nodes within an element and the distinct numerical characteristics of the Galerkin method. The conventional coupling procedure is to evaluate the physical parameterizations ( physics) on the dynamical core grid. Evaluating the physics at the nodal points exacerbates numerical noise from the Galerkin method, enabling and amplifying local extrema at element boundaries. Grid imprinting may be substantially reduced through the introduction of an entirely separate, approximately isotropic finite-volume grid for evaluating the physics forcing. Integration of the spectral basis over the control volumes provides an area-average state to the physics, which is more representative of the state in the vicinity of the nodal points rather than the nodal point itself and is more consistent with the notion of a “large-scale state” required by conventional physics packages. This study documents the implementation of a quasi-equal-area physics grid into NCAR’s Community Atmosphere Model Spectral Element and is shown to be effective at mitigating grid imprinting in the solution. The physics grid is also appropriate for coupling to other components within the Community Earth System Model, since the coupler requires component fluxes to be defined on a finite-volume grid, and one can be certain that the fluxes on the physics grid are, indeed, volume averaged.

54 ENVIRONMENTAL SCIENCES↗

Simplified dichromated gelatin hologram recording process

A simplified method for making dichromated gelatin (DCG) holographic optical elements (HOE) has been discovered. The method is much less tedious and it requires a period of processing time comparable with that for processing a silver halide hologram. HOE characteristics including diffraction efficiency (DE), linearity, and spectral sensitivity have been quantitatively investigated. The quality of the holographic grating is very high. Ninety percent or higher diffraction efficiency has been achieved in simple plane gratings made by this process.

Georgekutty, Tharayil G.↗

Robust Metric-Aligned Quad-Dominant Meshing Using L(sub p) Centroidal Voronoi Tessellation

We introduce a meshing algorithm that can be used to both generate and adapt meshes for bounded domains in an anisotropic manner. This is particularly beneficial when anisotropic flow features like shock waves or contact discontinuities are present in the computational domain. The algorithm presented in this paper is based upon meshing under the imposed Riemannian metric tensor, which controls the orientation and size of the mesh elements. In this way there is no need for user intervention to recognize these features. We demonstrate that the method indeed aligns the elements with the underlying metric and produces right-angled simplices that can be recombined into quadrilateral elements. The aim is to eventually incorporate this meshing strategy in the monolithic high-order spectral element solver that is currently being developed at NASA Ames. This paper has two main contributions: First, we demonstrate that we can generate quad-dominant metric-aligned meshes for bounded domains using a generalized form of L(sub p)-Centroidal Voronoi Tessellation (L(sub p)-CVT). Unlike previous works, we do not rely on a background mesh and discretize the bounded domain in a hierarchical way by first discretizing the boundaries and then the volume using the underlying metric. Second, we present an alternative for clipping the Voronoi cells on the boundary, which is common practice in CVT-based meshing algorithms, by reconstructing the Voronoi cells using the defined metric field. In this way we avoid the geometrical complexity of the clipping procedure and we show that we evaluate the energy and its gradient correctly. We show that the reconstruction of the computational domain is consistent with the Lloyds’ algorithm that is used to compute the L(sub p)-CVT.

Centroidal↗

Beyond Generalized Eigenvalues in Lattice Quantum Field Theory

Two analysis techniques, the generalized eigenvalue method (GEM) or Prony's (or related) method (PM), are commonly used to analyze statistical estimates of correlation functions produced in lattice quantum field theory calculations. GEM takes full advantage of the matrix structure of correlation functions but only considers individual pairs of time separations when much more data exists. PM can be applied to many time separations and many individual matrix elements simultaneously but does not fully exploit the matrix structure of the correlation function. We combine both these methods into a single framework based on matrix polynomials. As these algebraic methods are well known for producing extensive spectral information about statistically-noisy data, the method should be paired with some information criteria, like the recently proposed Bayesean model averaging.

Fleming, George T.↗