Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High spatial resolution”

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 1,639 records · Page 91

Effects of head modeling errors on the spatial frequency representation of MEG

Abstract Objectives. We aim to investigate the effects of head model inaccuracies on signal and source reconstruction accuracies for various sensor array distances to the head. This allows for the assessment of the importance of head modeling for next-generation magnetoencephalography (MEG) sensors, optically-pumped magnetometers (OPM). Approach. A 1-shell boundary element method (BEM) spherical head model with 642 vertices of radius 9 cm and conductivity of 0.33 S m −1 was defined. The vertices were then randomly perturbed radially up to 2%, 4%, 6%, 8% and 10% of the radius. For each head perturbation case, the forward signal was calculated for dipolar sources located at 2 cm, 4 cm, 6 cm and 8 cm from the origin (center of the sphere), and for a 324 sensor array located at 10 cm to 15 cm from the origin. Equivalent current dipole (ECD) source localization was performed for each of these forward signals. The signal for each perturbed spherical head case was then analyzed in the spatial frequency domain, and the signal and ECD errors were quantified relative to the unperturbed case. Main results. In the noiseless and high signal-to-noise ratio (SNR) case of approximately ≥6 dB, inaccuracies in our spherical BEM head conductor models lead to increased signal and ECD inaccuracies when sensor arrays are placed closer to the head. This is true especially in the case of deep and superficial sources. In the noisy case however, the higher SNR for closer sensor arrays allows for an improved ECD fit and outweighs the effects of head geometry inaccuracies. Significance. OPMs may be placed directly on the head, as opposed to the more commonly used superconducting quantum interference device sensors which must be placed a few centimeters away from the head. OPMs thus allow for signals of higher spatial resolution to be captured, resulting in potentially more accurate source localizations. Our results suggest that an increased emphasis on accurate head modeling for OPMs may be necessary to fully realize its improved source localization potential.

Yeo, Wan-Jin (ORCID:0000000241417123)↗

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗

Short-lead seasonal precipitation forecast in northeastern Brazil using an ensemble of artificial neural networks

This study assesses the deterministic and probabilistic forecasting skill of a 1-month-lead ensemble of Artificial Neural Networks (EANN) based on low-frequency climate oscillation indices. The predictand is the February-April (FMA) rainfall in the Brazilian state of Ceará, which is a prominent subject in climate forecasting studies due to its high seasonal predictability. Additionally, the study proposes combining the EANN with dynamical models into a hybrid multi-model ensemble (MME). The forecast verification is carried out through a leave-one-out cross-validation based on 40 years of data. The EANN forecasting skill is compared with traditional statistical models and the dynamical models that compose Ceará’s operational seasonal forecasting system. A spatial comparison showed that the EANN was among the models with the smallest Root Mean Squared Error (RMSE) and Ranked Probability Score (RPS) in most regions. Moreover, the analysis of the area-aggregated reliability showed that the EANN is better calibrated than the individual dynamical models and has better resolution than Multinomial Logistic Regression for above-normal (AN) and below-normal (BN) categories. It is also shown that combining the EANN and dynamical models into a hybrid MME reduces the overconfidence of the extreme categories observed in a dynamically-based MME, improving the reliability of the forecasting system.

54 ENVIRONMENTAL SCIENCES↗

Deciphering phase evolution in complex metal oxide thin films via high-throughput materials synthesis and characterization

Discovery of structure-property relationships in thin film alloys of complex metal oxides enabled by high-throughput materials synthesis and characterization facilities is demonstrated here with a case-study. In this study, thin films of binary transition metal oxides (Ti–Zn) are prepared by pulsed laser deposition with continuously varying Ti:Zn ratio, creating combinatorial samples for exploration of the properties of this material family. The atomic structure and electronic properties are probed by spatially resolved techniques including x-ray absorption near edge structures (XANES) and x-ray fluorescence (XRF) at the Ti and Zn K-edge, x-ray diffraction, and spectroscopic ellipsometry. The observed properties as a function of Ti:Zn ratio are resolved into mixtures of five distinguishable phases by deploying multivariate curve resolution analysis on the XANES spectral series, under constraints set by results from the other characterization techniques. First-principles computations based on density function theory connect the observed properties of each distinct phase with structural and spectral characteristics of crystalline polymorphs of Ti–Zn oxide. Continuous tuning of the optical absorption edge as a function of Ti:Zn ratio, including the unusual observation of negative optical bowing, exemplifies a functional property of the film correlated to the phase evolution.

36 MATERIALS SCIENCE↗

ASCENT: a balloon-borne hard x-ray imaging spectroscopy telescope using transition edge sensor microcalorimeter detectors

Core collapse supernovae are thought to be one of the main sources in the galaxy of elements heavier than iron. Understanding the origin of the elements is thus tightly linked to our understanding of the explosion mechanism of supernovae and supernova nucleosynthesis. X-ray and gamma-ray observations of young supernova remnants, combined with improved theoretical modeling, have resulted in enormous improvements in our knowledge of these events. Here, the isotope Ti44 is one of the most sensitive probes of the innermost regions of the core collapse engine, and its spatial and velocity distribution are key observables. Hard x-ray imaging spectroscopy with the Nuclear Spectroscopic Telescope Array (NuSTAR) has provided new insights into the structure of the supernova remnant Cassiopeia A (Cas A), establishing the convective nature of the supernova engine. However, many questions about the details of this engine remain. We present here the concept for a balloon-borne follow-up mission called A SuperConducting ENergetic x-ray Telescope (ASCENT). ASCENT uses transition edge sensor gamma-ray microcalorimeter detectors with a demonstrated 55-eV full-width half maximum energy resolution at 97 keV. This 8- to 16-fold improvement in energy resolution over NuSTAR will allow for high-resolution imaging and spectroscopy of the Ti44 emission. This will allow for a detailed reconstruction of gamma-ray line redshifts, widths, and shapes, allowing us to address questions such as, What is the source of the neutron star kicks? What is the dominant production pathway for Ti44? Is the engine of Cas A unique?

79 ASTRONOMY AND ASTROPHYSICS↗

Global Dust Cycle and Direct Radiative Effect in E3SM Version 1: Impact of Increasing Model Resolution

Quantification of dust aerosols in Earth System Models (ESMs) has important implications for water cycle and biogeochemistry studies. This study examines the global life cycle and direct radiative effects (DREs) of dust in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1), and the impact of increasing model resolution both horizontally and vertically. The default 1° E3SMv1 captures the spatial and temporal variability in the observed dust aerosol optical depth (DAOD) reasonably well, but overpredicts dust absorption in the shortwave (SW). Simulations underestimate the dust vertical and long-range transport, compared with the satellite dust extinction profiles. After updating dust refractive indices and correcting for a bias in partitioning size-segregated emissions, both SW cooling and longwave (LW) warming of dust simulated by E3SMv1 are increased and agree better with other recent studies. The estimated net dust DRE of –0.42 Wm –2 represents a stronger cooling effect than the observationally based estimate –0.2 Wm –2 (–0.48 to +0.2), due to a smaller LW warming. Constrained by a global mean DAOD, model sensitivity studies of increasing horizontal and vertical resolution show strong influences on the simulated global dust burden and lifetime primarily through the change of dust dry deposition rate; there are also remarkable differences in simulated spatial distributions of DAOD, DRE, and deposition fluxes. Thus, constraining the global DAOD is insufficient for accurate representation of dust climate effects, especially in transitioning to higher- or variable-resolution ESMs. Better observational constraints of dust vertical profiles, dry deposition, size, and LW properties are needed.

54 ENVIRONMENTAL SCIENCES↗

Estimation of soil classes and their relationship to grapevine vigor in a Bordeaux vineyard: advancing the practical joint use of electromagnetic induction (EMI) and NDVI datasets for precision viticulture

Working within a vineyard in the Pessac Léognan Appellation of Bordeaux, France, this study documents the potential of using simple statistical methods with spatially-resolved and increasingly available electromagnetic induction (EMI) geophysical and normalized difference vegetation index (NDVI) datasets to accurately estimate Bordeaux vineyard soil classes and to quantitatively explore the relationship between vineyard soil types and grapevine vigor. First, co-located electrical tomographic tomography (ERT) and EMI datasets were compared to gain confidence about how the EMI method averaged soil properties over the grapevine rooting depth. Then, EMI data were used with core soil texture and soil-pit based interpretations of Bordeaux soil types (Brunisol, Redoxisol, Colluviosol and Calcosol) to estimate the spatial distribution of geophysically-identified Bordeaux soil classes. A strong relationship (r = 0.75, p < 0.01) was revealed between the geophysically-identified Bordeaux soil classes and NDVI (both 2 m resolution), showing that the highest grapevine vigor was associated with the Bordeaux soil classes having the largest clay fraction. The results suggest that within-block variability of grapevine vigor was largely controlled by variability in soil classes, and that carefully collected EMI and NDVI datasets can be exceedingly helpful for providing quantitative estimates of vineyard soil and vigor variability, as well as their covariation. The method is expected to be transferable to other viticultural regions, providing an approach to use easy-to-acquire, high resolution datasets to guide viticultural practices, including routine management and replanting.

54 ENVIRONMENTAL SCIENCES↗

Emu v1.1

Emu is a particle-in-cell code for solving the neutrino quantum kinetic equations in 1, 2, or 3 spatial dimensions with arbitrary angular resolution in order to simulate neutrino flavor transformation in neutron star merger and core collapse supernova environments. Emu represents the neutrino distribution function as a set of particles, each of which represent a collection of neutrinos and antineutrinos with unique position and momentum. Each particle carries two density matrices to define the flavor state of the neutrinos and antineutrinos it represents. Emu includes the vacuum, matter, and neutrino self-interaction potentials. Emu calculates the self-interaction potential using PIC deposition and interpolation algorithms that efficiently compute the local number density and flux at particle locations. Emu is implemented in C++ and is based on the AMReX library for high-performance, block-structured adaptive mesh refinement. Emu is parallelized with MPI + OpenMP for CPUs and MPI + CUDA for GPUs.

Willcox, Donald↗

Forward silicon vertex/tracking detector design and R$\&$D for the future Electron-Ion Collider

The proposed high-luminosity high-energy Electron-Ion Collider (EIC) will provide a clean environment to precisely study several fundamental questions in the fields of high-energy and nuclear physics . A low material budget and high granularity silicon vertex/tracking detector is critical to carry out a series of hadron and jet measurements at the future EIC especially for the heavy flavor product reconstruction or tagging. The conceptual design of a proposed forward silicon tracking detector with the pseudorapidity coverage from 1.2 to 3.5 has been developed in integration with different magnet options and the other EIC detector sub-systems. The tracking performance of this detector enables precise heavy flavor hadron and jet measurements in the hadron beam going direction. The detector R&D for the proposed silicon technology candidates: Low Gain Avalanche Diode (LGAD) and radiation hard depleted Monolithic Active Pixel Sensor (MALTA), which can provide good spatial and timing resolutions, is underway. Bench test results of the LGAD and MALTA prototype sensors will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Free space Thomson scattering to study high energy density shocks

A free space collective Thomson scattering system has been developed to study pulsed power produced plasmas. While most Thomson scattering diagnostics on pulsed power machines use a bundle of fibers to couple scattered light from the plasma to the spectrometer, this system used free space coupling of the light, which enabled a spatially continuous image of the plasma. Initial experiments with this diagnostic were performed on an inverse wire array generated by a 200 kA, 1100 ns rise time pulse power generator. The capabilities of this diagnostic were demonstrated by using the low frequency ion acoustic wave feature of the Thomson scattering spectra to measure the plasma flow velocity. The diagnostic was demonstrated to measure velocities between 20 and 40 km/s with an error of 4.7 km/s when fitting with a 600 μm spatial resolution or 8.9 km/s when fitting with a 150 μm spatial resolution. In some experiments, the diagnostic observed a bow shock in the plasma flow as the scattering intensity increased and flow velocity decreased.

47 OTHER INSTRUMENTATION↗

In-situ TEM investigation of void swelling in nickel under irradiation with analysis aided by computer vision

Understanding the stability of irradiation-induced voids in materials is important for engineering material's swelling behavior under irradiation. In-situ TEM offers a spatial and temporal resolution that is suitable for investigating the evolution of voids under irradiation. However, the in-situ videos have often been too large to be analyzed manually, leaving the valuable data underutilized. For this work, we developed a deep learning-based semantic segmentation model to consistently study the growth and shrinkage of voids in nickel under 1 MeV krypton ion irradiation at various temperatures from 525 °C to 650 °C. With a foil thickness near 100 nm and ion flux of 6.3 × 10 11 ions∙ cm –2 ∙s –1 , the pre-existing voids, which were created beforehand by irradiation at 600 °C to 0.5 dpa, shrank at low temperatures and grew at high temperatures under irradiation, where the transition occurred at 575 °C (~0.5 T M ). The observed stability transition provided new insight for the shrinkage mechanism of voids under irradiation. In addition, an annealing experiment on nickel, previously irradiated at 600 °C to 3 dpa, was performed sequentially at 650 °C to 720 °C to reveal the shrinkage rate of void as a function of temperature and void size. The advantage of combining computer vision and in-situ TEM to obtain comprehensive void evolution was demonstrated.

36 MATERIALS SCIENCE↗

Custom surface reflectance, shade mask, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study (2025)

This dataset contains land surface reflectance estimates and additional derived products generated from NEON Imaging Spectrometer (NIS) data collected in the Upper Gunnison river basin during June and July of 2025. Data was collected over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). These products were derived from radiance and LiDAR data collected by the NEON Airborne Observation Platform (AOP) campaign funded by the Colorado Headwaters Ecological Spectroscopy Study (CHESS) (doi:10.15485/3017965). Products include per-pixel surface reflectance (rfl) and reflectance uncertainty (rfl_unc), observational data (obs), canopy equivalent water thickness (ewt), and shade masks. Atmospheric correction was performed per flightline using the ISOFIT (Imaging Spectrometer Optimal FITting) optimal estimation framework to estimate surface reflectance and the associated per-band reflectance uncertainty. Reflectance retrievals achieved a mean absolute error of 1.5% across diverse validation surfaces (see validation report.pdf). Equivalent water thickness was calculated from surface reflectance using the Beer–Lambert absorption of liquid water. Shade masks were generated based on the geometry between the sun angle, ground surface, and sensor at the time of flight. Data products are provided per-flightline and as mosaics for each domain. Flightline data products are provided as ENVI-formatted binary files (rfl, rfl_unc, ewt) and GeoTIFFs (shade). Reflectance and uncertainty mosaics are provided as tiled NetCDFs, while all other mosaicked products are provided as cloud-optimized GeoTIFFs. These formats are supported by common geospatial software (e.g., QGIS, ArcGIS, ENVI) and programmatic libraries in Python (e.g., rasterio, xarray, spectral, netCDF4) and R (e.g., terra, ncdf4). Processing workflows were designed to be equivalent to those used to generate the 2018 CHESS campaign airborne imaging spectroscopy data products (doi:10.15485/3013527). All outputs were co-registered to a common spatial grid to support time series analyses. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgment: Data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). Computational research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004) and was funded by EMIT Extended Mission Phase E Science.

2018 NEON and 2025 CHESS Campaigns↗

The Blue Straggler Population of the Open Clusters Trumpler 5, Trumpler 20, and NGC 2477

We present a study based on Gaia DR2 of the population of blue straggler stars in the open clusters Trumpler 5, Trumpler 20, and NGC 2477. All candidates were selected according to their position in the color–magnitude diagram, their proper motion components, and their parallax. We also looked for yellow stragglers, i.e., possible evolved blue stragglers. We found that Trumpler 5 hosts a large blue straggler star population, which allowed us to analyze their radial distribution as a probe of the cluster’s dynamical status. The blue straggler star distribution was compared with that of red giant branch stars to evaluate mass segregation. Our results indicate that blue straggler stars are not more centrally concentrated than red giant branch stars stars in any of the clusters. The radial distribution of blue straggler stars in Trumpler 5 is flat. Additionally, using a multi-epoch radial velocity survey conducted with the high-resolution spectrograph FLAMES/GIRAFFE at the Very Large Telescope, we measured the radial velocities of a sample of stragglers to compare with the mean radial velocity and velocity dispersion of the clusters. Based on the radial velocity variations for different epochs, we roughly classified these stars as possible close or long-period binaries.

47 OTHER INSTRUMENTATION↗

Muon tracking in a LiquidO opaque scintillator detector

LiquidO is an innovative radiation detector concept. The core idea is to exploit stochastic light confinement in a highly scattering medium to self-segment the detector volume. In this paper, we demonstrate event-by-event muon tracking in a LiquidO opaque scintillator detector prototype. The detector consists of a 30 mm cubic scintillator volume instrumented with 64 wavelength-shifting fibres arranged in an 8 × 8 grid with a 3.2 mm pitch and read out by silicon photomultipliers. A wax-based opaque scintillator with a scattering length of approximately 0.5 mm is used. The tracking performance of this LiquidO detector is characterised with cosmic-ray muons and the position resolution is demonstrated to be 450 μm per row of fibres. These results highlight the potential of LiquidO opaque scintillator detectors to achieve fine spatial resolution, enabling precise particle tracking and imaging.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutrino Interactions observed with Large Area Picosecond Photodetectors in ANNIE

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton water-based neutrino detector located at Fermilab, approximately 110 m downstream of the Booster Neutrino Beam (BNB). ANNIE utilizes both photomultiplier tubes (PMTs) and advanced photodetectors, specifically Large Area Picosecond Photodetectors (LAPPDs), to detect Cherenkov light emitted by leptons produced in neutrino interactions within ANNIE. LAPPDs are a novel technology designed to detect photons with picosecond-level temporal resolution and sub-millimetre spatial precision. Multiple LAPPDs have been deployed in the ANNIE detector. This is the first use of this technology in a running particle physics experiment and has yielded the first detection of light from neutrino interactions in water with LAPPDs. In this poster, I will showcase the operational performance and functionality of LAPPDs in the ANNIE experiment. Neutrino beam data from the BNB is used to evaluate the timing precision, hit reconstruction performance, and beam response of deployed LAPPDs, demonstrating how this novel picosecond-resolution technology performs in a running neutrino water Cherenkov detector. This work highlights the successful integration of LAPPDs in ANNIE and provides quantitative benchmarks that inform their application in future neutrino experiments requiring high-resolution photon detection.

Aman, Mohammad Adil [Florida State U.] (ORCID:0009↗

Development of limited-view tomography for measurement of Spray G plume direction and liquid volume fraction

Abstract The method for direct injection of fuel in the cylinder of an IC engines is important to high-efficiency and low-emission performance. Optical spray diagnostics plays an important role in understanding plume movement and interaction for multi-hole injectors, and providing baseline understanding used for computational optimization of fuel delivery. Traditional planar or line-of-sight diagnostics fail to capture the liquid distribution because of optical thickness concerns. This work proposes a high-speed (67 kHz) extinction imaging technique at various injector rotations coupled to computed tomography (CT) for time-resolved reconstruction of liquid volume fraction in three dimensions. The number of views selected and processing were based on synthetic (modeled) liquid volume fraction data where extinction and CT adequately reconstructed each plume. The exercise showed that for an 8-hole, symmetric-design injector (ECN Spray G), only three different views are enough to reproduce the direction of each plume, and particularly the mean plume direction. Therefore, the number of views was minimized for experiments to save expense. Measurements applying this limited-view technique confirm plume–plume variations also detected with mechanical patternation, while providing better spatial and temporal resolution than achieved previously. Uncertainties due to the limited view within pressurized spray chambers, the droplet size, and optically thick regions are discussed. Graphic abstract

42 ENGINEERING↗

Characterization, comparison, and optimization of lattice light sheets

Lattice light sheet microscopy excels at the noninvasive imaging of three-dimensional (3D) dynamic processes at high spatiotemporal resolution within cells and developing embryos. Recently, several papers have called into question the performance of lattice light sheets relative to the Gaussian sheets most common in light sheet microscopy. Here, we undertake a theoretical and experimental analysis of various forms of light sheet microscopy, which demonstrates and explains why lattice light sheets provide substantial improvements in resolution and photobleaching reduction. The analysis provides a procedure to select the correct light sheet for a desired experiment and specifies the processing that maximizes the use of all fluorescence generated within the light sheet excitation envelope for optimal resolution while minimizing image artifacts and photodamage. We also introduce a new type of “harmonic balanced” lattice light sheet that improves performance at all spatial frequencies within its 3D resolution limits and maintains this performance over lengthened propagation distances allowing for expanded fields of view.

59 BASIC BIOLOGICAL SCIENCES↗

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗