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 451 records · Page 25

Chemical vapor deposition fluid flow simulation modelling tool

Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.

Bullister, Edward T.↗

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗

Thermal behavior characterization of a kilowatt-power-level cryogenically cooled Yb:YAG active mirror laser amplifier

We present an analysis of the thermal behavior of a high-energy kilowatt-average-power diode-pumped cryogenically cooled Yb:YAG active mirror laser amplifier based on measurements and simulations. Maps of the temperature distribution of the laser material at pump powers up to 1 kW were obtained for the first time by spatially and spectrally resolving the fluorescence induced by a scanning probe beam. The wavefront distortion resulting from the front surface deformation and the overall deformation of the gain medium assembly were measured using a Mach–Zehnder interferometer. The measured deformations agree well with the results of thermomechanical modeling using finite element method simulations, and with the results of focal length shift measurements. The relative contributions to the optical path difference (OPD) of the mechanical deformations, refractive index changes, and electronic contribution are discussed. We propose that the Cr 4+ :YAG cladding plays a significant role in both the temperature distribution and the overall OPD changes. Moreover, the pump-induced mechanical deformations of the assembly dominate the OPD changes in this kilowatt-average-pump-power cryogenically cooled Yb:YAG active mirror laser.

36 MATERIALS SCIENCE↗

Enstrophy budget in decaying compressible turbulence

The interaction of vorticity and shock waves in decaying compressible turbulence is investigated by direct numerical simulation of the Navier-Stokes equation using a spectral method with 64-cubed collocation points. The Taylor microscale Reynolds number decreases from 36 to 17 and the rms Mach number from 0.8 to 0.4 during the decay. It is observed that vorticity is created inside curved shocks through baroclinic interaction and intensified by the compression of fluid elements. Vortex stretching is, on the other hand, dominant outside shocks and is responsible for enhancing the strength of vorticity. A dipole of vorticity is typically created at the vertex of curved shocks, the structure of which is explained in terms of jump relations of velocity derivatives across a curved shock.

Kida, Shigeo↗

Large tuning of narrow-beam terahertz plasmonic lasers operating at 78 K

A new tuning mechanism is demonstrated for single-mode metal-clad plasmonic lasers, in which the refractive-index of the laser’s surrounding medium affects the resonant-cavity mode in the same vein as the refractive-index of gain medium inside the cavity. Reversible, continuous, and mode-hop-free tuning of ~57 GHz is realized for single-mode narrow-beam terahertz plasmonic quantum-cascade lasers (QCLs), which is demonstrated at a much more practical temperature of 78 K. The tuning is based on post-process deposition/etching of a dielectric (silicon-dioxide) on a QCL chip that has already been soldered and wire-bonded onto a copper mount. This is a considerably larger tuning range compared to previously reported results for terahertz QCLs with directional far-field radiation patterns. The key enabling mechanism for tuning is a recently developed antenna-feedback scheme for plasmonic lasers, which leads to the generation of hybrid surface-plasmon-polaritons propagating outside the cavity of the laser with a large spatial extent. The effect of dielectric deposition on QCL’s characteristics is investigated in detail including that on maximum operating temperature, peak output power, and far-field radiation patterns. Single-lobed beam with low divergence (<7°) is maintained through the tuning range. The antenna-feedback scheme is ideally suited for modulation of plasmonic lasers and their sensing applications due to the sensitive dependence of spectral and radiative properties of the laser on its surrounding medium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NCAR Release of CAM‐SE in CESM2.0: A Reformulation of the Spectral Element Dynamical Core in Dry‐Mass Vertical Coordinates With Comprehensive Treatment of Condensates and Energy

Abstract It is the purpose of this paper to provide a comprehensive documentation of the new NCAR (National Center for Atmospheric Research) version of the spectral element (SE) dynamical core as part of the Community Earth System Model (CESM2.0) release. This version differs from previous releases of the SE dynamical core in several ways. Most notably the hybrid sigma vertical coordinate is based on dry air mass, the condensates are dynamically active in the thermodynamic and momentum equations (also referred to as condensate loading), and the continuous equations of motion conserve a more comprehensive total energy that includes condensates. Not related to the vertical coordinate change, the hyperviscosity operators and the vertical remapping algorithms have been modified. The code base has been significantly reduced, sped up, and cleaned up as part of integrating SE as a dynamical core in the CAM (Community Atmosphere Model) repository rather than importing the SE dynamical core from High‐Order Methods Modeling environment as an external code.

54 ENVIRONMENTAL SCIENCES↗

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Enhancement of terahertz fields in LiTaO 3 waveguides using a conical pulse front

The development of methods for the generation of strong ultrafast electromagnetic pulses in the terahertz (THz) spectral range has led to a surge of progress in nonlinear THz spectroscopy and THz control of molecular and collective responses. For spectroscopy in the 1-THz range, the submillimeter wavelengths and associated large spot sizes, large optical elements, and short distances between final focusing elements and samples can lead to cumbersome experimental setups that are incompatible with some sample environments. Here, we introduce a novel terahertz ring excitation (TREx) optical pumping geometry to generate superposing, focusing fields in planar THz waveguides made out of the electro-optic material lithium tantalate. High THz fields, >175 kV/cm, are generated and measured optically with no free-space THz propagation. The field level achieved by pumping with a sequence of concentric rings of excitation light exceeds by about 20× the result of a single cylindrically focused line of pump light that has been used routinely in previous work. The technique opens new prospects for compact waveguide-based linear and nonlinear THz spectroscopy and signal processing.

47 OTHER INSTRUMENTATION↗

Temperature Measurements Using the Dielectronic Satellite Lines of Fe xvii

The coronal temperature of Capella and Procyon has been derived from the spectral emission of Fe XVII observed in high resolution with Chandra's High Energy Transmission Grating Spectrometer. The measurement makes use of well resolved dielectronic satellite features with a spectator electron in the n = 4 shell. This method does not require a priori knowledge of the ionization equilibrium or of the elemental abundances. We find the temperature of Fe XVII emission is ${264}_{-51}^{+142}$ eV for Procyon and ${514}_{-25}^{+27}$ eV for Capella, i.e., 3.1 and 6.0 MK, respectively. As a result, the Procyon measurement is based on three X-rays observed in the dielectronic satellite feature, which illustrates the power of the method even for low-signal observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High spectral resolution remote sensing of the land

Spectral remote sensing has been practiced on a large scale since the launch of Landsat 1 in 1972. The limited information contained in this spectrally undersampled data set has led to the development of sophisticated statistical-inferential methods for data analysis. The results are usually limited by the availability of ground truth information. Recent technological developments have made it feasible to create narrow-band, contiguous, spectral image data sets that make possible the identification of surface cover materials based on the complete reflectance spectrum for each picture element. This capability will revolutionize the use of remote sensing data and require new deterministic image processing techniques to extract the full information content from the data. Sensors, based on the concept of imaging spectrometry and the new technology of area array infrared detectors, have been constructed and are candidates for Shuttle and space platform flights.

Goetz, A. F. H.↗

Using Directional Diffusion Coefficients for Nonlinear Diffusion Acceleration of the First Order S{sub N} Equations in Near-Void Regions

The common definition of the diffusion coefficient as the inverse of three times the transport cross section is not compatible with voids. Morel introduced a non-local tensor diffusion coefficient that remains finite in voids. It can be obtained by solving an auxiliary transport problem without scattering or fission. Larsen and Trahan successfully applied this diffusion coefficient for enhancing the accuracy of diffusion solutions of very high temperature reactor (VHTR) problems that feature large, optically thin channels in the z-direction. It was demonstrated that a significant reduction of error can be achieved in particular in the optically thin region. Along the same line of thought, non-local diffusion tensors have been applied to modeling the TREAT reactor confirming the findings of Larsen and Trahan. Previous work of the authors have introduced a flexible Nonlinear-Diffusion Acceleration (NDA) method for the first order S{sub N} equations discretized with the discontinuous finite element method (DFEM). This NDA method uses a scalar diffusion coefficient in the low-order system that is obtained as the flux weighted average of the inverse transport cross section. Hence, it suffers from very large and potentially unbounded diffusion coefficients in the low order problem. However, it was noted that the choice of the diffusion coefficient does not influence consistency of the method at convergence and hence the diffusion coefficient is essentially a free parameter. The choice of the diffusion coefficient does, however, affect the convergence behavior of the nonlinear diffusion iterations. Within this work we use Morel's non-local diffusion coefficient in the aforementioned NDA formulation in lieu of the flux weighted inverse of three times the transport cross section. The goal of this paper is to demonstrate that significant enhancement of the spectral properties of NDA can be achieved in near void regions. For testing the spectral properties of the NDA with non-local diffusion coefficients, the periodical horizontal interface problem is used. This problem consists of alternating stripes of optically thin and thick materials both of which feature scattering ratios close to unity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Application of distance correction to ChemCam laser-induced breakdown spectroscopy measurements

Laser-induced breakdown spectroscopy (LIBS) provides chemical information from atomic, ionic, and molecular emissions from which geochemical composition can be deciphered. Analysis of LIBS spectra in cases where targets are observed at different distances, as is the case for the ChemCam instrument on the Mars rover Curiosity, which performs analyses at distances between 2 and 7.4 m is not a simple task. Previously, we showed that spectral distance correction based on a proxy spectroscopic standard created from first-shot dust observations on Mars targets ameliorates the distance bias in multivariate-based elemental-composition predictions of laboratory data. In this work, we correct an expanded set of neutral and ionic spectral emissions for distance bias in the ChemCam data set. By using and testing different selection criteria to generate multiple proxy standards, we find a correction that minimizes the difference in spectral intensity measured at two different distances and increases spectral reproducibility. When the quantitative performance of distance correction is assessed, there is improvement for SiO 2 , Al 2 O 3 , CaO, FeOT, Na 2 O, K 2 O, that is, for most of the major rock forming elements, and for the total major-element weight percent predicted. But, for MgO the method does not provide improvements while for TiO 2 , it yields inconsistent results. Additionally, we observed that many emission lines do not behave consistently with distance, evidenced from laboratory analogue measurements and ChemCam data. This limits the effectiveness of the method.

36 MATERIALS SCIENCE↗

Imaging system and method for enhanced visualization of near surface vascular structures

A sub-dermal structure visualization system and method is disclosed. The system may use an illumination module including: a near-infrared (NIR) light source adapted to substantially uniformly illuminate an imaged area including sub-dermal regions thereof with NIR light; and a first optical system including at least one optical element for controlling at least one of spectral and polarization properties of the NIR light prior to illuminating the imaged area. An imaging module rejects unwanted signals from an imaged area while passing desired optical signals that are to be received by an image acquisition module. The desired optical signals can comprise a vein visualization signal to assist in visualizing a vascular structure below a skin layer of a patient.

Demos, Stavros G.↗

Spectral Element Applications in Complex Nuclear Reactor Geometries: Tet-to-Hex Meshing

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a let-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYSICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

ANSYS ICEM↗

Achromatic registration of quadrature components of the optical spectrum in spectral domain optical coherence tomography

We have thoroughly investigated the method of simultaneous reception of spectral components with the achromatised quadrature phase shift between two portions of a reference wave, designed for the effective suppression of the 'mirror' artefact in the resulting image obtained by means of spectral domain optical coherence tomography (SD OCT). We have developed and experimentally tested a phase-shifting element consisting of a beam divider, which splits the reference optical beam into the two beams, and of delay lines being individual for each beam, which create a mutual phase difference of π/2 in the double pass of the reference beam. The phase shift achromatism over a wide spectral range is achieved by using in the delay lines the individual elements with different dispersion characteristics. The ranges of admissible adjustment parameters of the achromatised delay line are estimated for exact and inexact conformity of the geometric characteristics of its components to those calculated. A possibility of simultaneous recording of the close-to-quadrature spectral components with a single linear photodetector element is experimentally confirmed. The suppression of the artefact mirror peak in the OCT-signal by an additional 9 dB relative to the level of its suppression is experimentally achieved when the air delay line is used. Two-dimensional images of the surface positioned at an angle to the axis of the probe beam are obtained with the correction of the 'mirror' artefact while maintaining the dynamic range of the image. (laser biophotonics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

WE-G-18C-07: Accelerated Water/fat Separation in MRI for Radiotherapy Planning Using Multi-Band Imaging Techniques

Purpose: Dixon sequences are used to characterize disease processes, obtain good fat or water separation in cases where fat suppression fails and to obtain pseudo-CT datasets. Dixon's method uses at least two images acquired with different echo times and thus requires prolonged acquisition times. To overcome associated problems (e.g., for DCE/cine-MRI), we propose to use a method for water/fat separation based on spectrally selective RF pulses. Methods: Two alternating RF pulses were used, that imposes a fat selective phase cycling over the phase encoding lines, which results in a spatial shift for fat in the reconstructed image, identical to that in CAIPIRINHA. Associated aliasing artefacts were resolved using the encoding power of a multi-element receiver array, analogous to SENSE. In vivo measurements were performed on a 1.5T clinical MR-scanner in a healthy volunteer's legs, using a four channel receiver coil. Gradient echo images were acquired with TE/TR = 2.3/4.7ms, flip angle 20°, FOV 45×22.5cm{sup 2}, matrix 480×216, slice thickness 5mm. Dixon images were acquired with TE,1/TE,2/TR=2.2/4.6/7ms. All image reconstructions were done in Matlab using the ReconFrame toolbox (Gyrotools, Zurich, CH). Results: RF pulse alternation yields a fat image offset from the water image. Hence the water and fat images fold over, which is resolved using in-plane SENSE reconstruction. Using the proposed technique, we achieved excellent water/fat separation comparable to Dixon images, while acquiring images at only one echo time. Conclusion: The proposed technique yields both inphase water and fat images at arbitrary echo times and requires only one measurement, thereby shortening the acquisition time by a factor 2. In future work the technique may be extended to a multi-band water/fat separation sequence that is able to achieve single point water/fat separation in multiple slices at once and hence yields higher speed-up factors.

60 APPLIED LIFE SCIENCES↗

SOLAR OPACITY CALCULATIONS USING THE SUPER-TRANSITION-ARRAY METHOD

A new opacity model has been developed based on the Super-Transition-Array (STA) method for the calculation of monochromatic opacities of plasmas in local thermodynamic equilibrium. The atomic code, named STAR (STA-Revised), is described and used to calculate spectral opacities for a solar model implementing the recent AGSS09 composition. Calculations are carried out throughout the solar radiative zone. The relative contributions of different chemical elements and atomic processes to the total Rosseland mean opacity are analyzed in detail. Monochromatic opacities and charge-state distributions are compared with the widely used Opacity Project (OP) code, for several elements near the radiation–convection interface. STAR Rosseland opacities for the solar mixture show a very good agreement with OP and the OPAL opacity code throughout the radiation zone. Finally, an explicit STA calculation was performed of the full AGSS09 photospheric mixture, including all heavy metals. It was shown that, due to their extremely low abundance, and despite being very good photon absorbers, the heavy elements do not affect the Rosseland opacity.

79 ASTRONOMY AND ASTROPHYSICS↗