SEARCH · Engineering Papers

Results for “spatial discretization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Evolution of Local Microstructures: Spatial Instabilities of Coarsening Clusters

This work examines the diffusional growth of discrete phase particles dispersed within a matrix. Engineering materials are often microstructurally heterogeneous, and the details of the microstructure determine how well that material performs in a given application. Critical to the development of designing multiphase microstructures with long-term stability is the process of Ostwald ripening. Ripening, or phase coarsening, is diffusion-limited and arises in polydisperse multiphase materials. Growth and dissolution occur because fluxes of solute, driven by chemical potential gradients at the interfaces of the dispersed phase material, depend on particle size. Competitive kinetics of these processes dictates that larger particles grow at the expense of smaller ones, overall leading to an increase of the average particle size. The classical treatment of phase coarsening was done by Todes, Lifshitz, and Slyozov, (TLS) in the limit of zero volume fraction, V(sub V)=0 of the dispersed phase. Since the publication of TLS theory, there have been numerous investigations, many of which sought to describe the kinetic scaling behavior over a range of volume fractions. Some studies in the literature report that the relative increase in coarsening rate at low (but not zero) volume fractions compared to that predicted by TLS is proportional to V(exp 1/2)(sub v) whereas others suggest V(exp 1/3)(sub v).This issue has been addressed recently by simulation studies at low volume fractions in three dimensions by members of the Rensselaer/MSFC team.

Glicksman, Martin E.↗

NASA NTRS2010

Contribution of Small-Scale Correlated Fluctuations of Microstructural Properties of a Spatially Extended Geophysical Target Under the Assessment of Radar Backscatter

The study of the collective effects of radar scattering from an aggregation of discrete scatterers randomly distributed in a space is important for better understanding the origin of the backscatter from spatially extended geophysical targets (SEGT). We consider the microstructure irregularities of a SEGT as the essential factor that affect radar backscatter. To evaluate their contribution this study uses the "slice" approach: particles close to the front of incident radar wave are considered to reflect incident electromagnetic wave coherently. The radar equation for a SEGT is derived. The equation includes contributions to the total backscatter from correlated small-scale fluctuations of the slice's reflectivity. The correlation contribution changes in accordance with an earlier proposed idea by Smith (1964) based on physical consideration. The slice approach applied allows parameterizing the features of the SEGT's inhomogeneities.

Yurchak, Boris S.↗

DOE OSTI2023

Imaging Nanoscale Energy Transport and Conversion with Ultrafast Electron Microscopy (Final Technical Report)

Light-matter interactions are ubiquitous in nature and reside at the heart of innumerable technologies. The cascade of processes that occur when a material absorbs a photon of light are exceedingly complex and are interwoven in both space and time, rendering precise determination of the atomic-scale and ultrafast mechanisms immensely challenging. The advent of methods for generating short pulses of light several decades ago led to major advances in understanding the initial moments of light absorption and the resultant effects, though directly interrogating the response of the atoms within the material continued to prove challenging. More recently, methods for generating ultrashort pulses of X-rays and fast electrons have opened the way to probing photoinduced structural dynamics of a wide range of matter in multiple phases. An especially promising laboratory-scale method is ultrafast electron microscopy (UEM), wherein the modalities of conventional transmission electron microscopes (imaging, diffraction, spectroscopy) are extended into the femtosecond temporal regime. Here, the ultrafast imaging and diffraction modalities of UEM were used to study the transient structural aspects of photoexcitation of semiconducting materials (e.g., spatially-resolved electron-phonon coupling, excitation and emission of acoustic phonons, and discrete nanoscale scattering processes). The project had three main objectives: (1) determination of the excitation mechanisms of dense, hypersonic charge-carrier waves and the spatially mediated means by which they couple to the lattice via coherent phonon emission, especially with UEM imaging, wherein effects of nanoscale structural and morphological features on the coupling and relaxation dynamics are expected to dictate nucleation sites and preferred wave vectors, (2) elucidation of photoinduced acoustic-phonon seeding, emergence, propagation, and decay over nanoscale crystal regions and especially with respect to local strain fields and atomic-scale disruptions in lattice order, and (3) interwoven with the first two were efforts aimed at realizing combined angstrom-femtosecond spatiotemporal imaging with UEM. The outcomes and impacts of this project were the generation of new knowledge with respect to fundamental light-matter interactions and, especially, the spatially-mediated excitation and evolution of the structural response of materials following coherent photoexcitation. Importantly, the spatial and temporal resolutions of the UEM imaging modalities used are well-suited for such studies and enable spatially-resolved mechanisms to be determined as a function of atomic order, structural features, and morphology. In addition, ultrafast crystallographic measurements were used to correlate real- and reciprocal-space dynamics in order to determine atomic-scale preferential wave vectors and ultrafast scattering mechanisms, especially as dictated by specimen boundary conditions. The obtained results, as detailed in peer-reviewed publications and presentations, illustrate the importance of ultrafast, angstrom-scale real-space imaging for developing a comprehensive understanding of energy transport and conversion in materials.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Gauging anomalous unitary operators

Boundary theories of static bulk topological phases of matter are obstructed in the sense that they cannot be realized on their own as isolated systems. The obstruction can be quantified/characterized by quantum anomalies, in particular when there is a global symmetry. Similarly, topological Floquet evolutions can realize obstructed unitary operators at their boundaries. In this paper, we discuss the characterization of such obstructions by using quantum anomalies. As a particular example, here we discuss time-reversal symmetric boundary unitary operators in one and two spatial dimensions, where the anomaly emerges as we gauge the so-called Kubo-Martin-Schwinger (KMS) symmetry. We also discuss mixed anomalies between particle number conserving U(1) symmetry and discrete symmetries, such as C and CP, for unitary operators in odd spatial dimensions that can be realized at the boundaries of topological Floquet systems in even spatial dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2000

High Spatial Resolution Airborne Multispectral Thermal Infrared Remote Sensing Data for Analysis of Urban Landscape Characteristics

We have used airborne multispectral thermal infrared (TIR) remote sensing data collected at a high spatial resolution (i.e., 10m) over several cities in the United States to study thermal energy characteristics of the urban landscape. These TIR data provide a unique opportunity to quantify thermal responses from discrete surfaces typical of the urban landscape and to identify both the spatial arrangement and patterns of thermal processes across the city. The information obtained from these data is critical to understanding how urban surfaces drive or force development of the Urban Heat Island (UHI) effect, which exists as a dome of elevated air temperatures that presides over cities in contrast to surrounding non-urbanized areas. The UHI is most pronounced in the summertime where urban surfaces, such as rooftops and pavement, store solar radiation throughout the day, and release this stored energy slowly after sunset creating air temperatures over the city that are in excess of 2-4'C warmer in contrast with non-urban or rural air temperatures. The UHI can also exist as a daytime phenomenon with surface temperatures in downtown areas of cities exceeding 38'C. The implications of the UHI are significant, particularly as an additive source of thermal energy input that exacerbates the overall production of ground level ozone over cities. We have used the Airborne Thermal and Land Applications Sensor (ATLAS), flown onboard a Lear 23 jet aircraft from the NASA Stennis Space Center, to acquire high spatial resolution multispectral TIR data (i.e., 6 bandwidths between 8.2-12.2 (um) over Huntsville, Alabama, Atlanta, Georgia, Baton Rouge, Louisiana, Salt Lake City, Utah, and Sacramento, California. These TIR data have been used to produce maps and other products, showing the spatial distribution of heating and cooling patterns over these cities to better understand how the morphology of the urban landscape affects development of the UHI. In turn, these data have been used by government officials, urban planners, and other decision-makers, to make more informed decisions on how to mitigate the UHI and its subsequent impacts.

Quattrochi, Dale A.↗

DOE OSTI2023

Improved treatment of multi-material cells in thermal radiation transport codes

High-energy-density physics simulations with non-conformal meshes of materials require a multi-material (MM) closure that affects the thermal radiation transport (TRT). We propose a set of novel closures that work for an arbitrary number of materials, both grey and multigroup energy discretizations, and any angular discretization (such as Sn, IMC, or diffusion). For each spatial cell, our closures let each species (ion and electron) of each material have its own temperature, density, and internal energy, but use a single radiation distribution that interacts with all materials within the cell. Our closures maintain energy conservation, do not incur increased computational cost in the TRT solve itself, are compatible with single-material TRT solvers, and do not make any temperature-equilibrium assumptions. Here we test our closures on a wide range of increasingly realistic problems and find them to be robust.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Capability Demonstration of a 3D CdZnTe Detector on a High-Altitude Balloon Flight

In collaboration between the University of Michigan and Los Alamos National Laboratory, a 3D position-sensing CdZnTe (CZT) detector prototype was built and integrated into a high-altitude balloon platform to evaluate its performance in a space-like mixed-radiation environment. The detector prototype, Orion Eagle, was designed to operate in near-vacuum environments without any temperature regulation. Orion Eagle was hand-launched from NASA’s Columbia Scientific Balloon Facility (CSBF) at Fort Sumner, NM on September 26, 2021, and successfully operated throughout a 9-hour flight, which reached 38.5 km in altitude. The flight met its objectives, successfully detecting atmospheric gamma rays and galactic cosmic rays, and raising the Technical Readiness Level from 4 to 6 for large-volume 3D CZT detector technology for space applications. Ionization tracks produced by charged particles create spatial signatures in the detector that are distinguishable from discrete gamma-ray interactions. Therefore, the 3D position-sensing capabilities using pixelated electrodes on a CZT detector can help enable discrimination of background charged particles from gamma-ray events without an anticoincidence shield. The potential for background rejection capability, ambient-temperature operation, gamma-ray coded-aperture and Compton imaging, and near High Purity Germanium (HPGe) energy resolution motivate the use of large-volume 3D CZT imaging spectrometers in future space missions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗

NASA NTRS1980

Aliasing and blurring in 2-D sampled imagery

The quality of image reconstructions from discrete data suffers not only from the blurring of spatial detail caused by limitations in the spatial frequency response of electrooptical systems, but also from the aliasing generated if spatial detail has been undersampled. Previous studies have observed that reasonable spot intensity profiles and photosensor aperture shapes of equivalent size result in about equal blurring but that some profiles and shapes suppress aliasing better than others. This paper presents quantitative results of the magnitude of aliasing and blurring as a function of random radiance fields typical for natural scenes and of spatial responses and sampling intervals typical for TV cameras and optical-mechanical scanners. These results indicate that aliasing may often be a larger source of degradation than either blurring or electronic noise.

Huck, F. O.↗

NASA NTRS1991

Do clusters of galaxies affect the spectrum of the microwave background?

The distortion, averaged over the sky, due to the Compton scattering of background photons with electrons in the hot gas in clusters of galaxies is calculated. Using an existing sample of X-ray clusters, various values of the density parameter Omega, and plausible models for cluster evolution, Monte Carlo realizations of the microwave sky are generated. The spatial structure of these simulations shows a network of discrete sources whose properties can be a strong function of both evolution and Omega. The amount of spectral distortion in the models is greatest for models characterized by self-similar cluster evolution in an open universe and is within an order of magnitude of the current upper limits. Thus, improved observational sensitivity must inevitably detect some deviation from a blackbody spectrum.

Markevitch, M.↗

NASA NTRS2013

Remote Observations of Reentering Spacecraft Including the Space Shuttle Orbiter

Flight measurement is a critical phase in development, validation and certification processes of technologies destined for future civilian and military operational capabilities. This paper focuses on several recent NASA-sponsored remote observations that have provided unique engineering and scientific insights of reentry vehicle flight phenomenology and performance that could not necessarily be obtained with more traditional instrumentation methods such as onboard discrete surface sensors. The missions highlighted include multiple spatially-resolved infrared observations of the NASA Space Shuttle Orbiter during hypersonic reentry from 2009 to 2011, and emission spectroscopy of comparatively small-sized sample return capsules returning from exploration missions. Emphasis has been placed upon identifying the challenges associated with these remote sensing missions with focus on end-to-end aspects that include the initial science objective, selection of the appropriate imaging platform and instrumentation suite, target flight path analysis and acquisition strategy, pre-mission simulations to optimize sensor configuration, logistics and communications during the actual observation. Explored are collaborative opportunities and technology investments required to develop a next-generation quantitative imaging system (i.e., an intelligent sensor and platform) with greater capability, which could more affordably support cross cutting civilian and military flight test needs.

Horvath, Thomas J.↗

NASA NTRSDate not supplied

Indigenous Organic Matter Associated with Unusual Low Temperature Aqueous Secondary Phases in the Mars Meteorite Nakhla

In freshly fractured chips of the Mars meteorite Nakhla, we found carbonaceous phases intimately associated with secondary aqueous alteration phases, both mineral and amorphous, interpreted to have formed through the low-temperature, mineral-water dissolution of the host Martian basalt.We have characterized the nature of carbonaceous matter and spatially associated secondary phases within Nakhla using a multidisciplinary suite of optical and electron beam instrumentation(Thomas-Keprta et al., 2022). Secondary Phases: The preserved secondary mineral phases are best described as ‘juvenile’ hydrated amorphous silica formed when water interacts with mafic rocks. Although amorphous silica on Earth rarely persists through geologic time due to alteration into other phases, e.g.,microcrystalline opal-CT and eventually quartz, the lack of maturation of primarily amorphous silica with embedded nano phase ferrihydrite in Nakhla infers these phases formed in a limited water environment by low temperature, circumneutral pH fluids stabilized by the presence of coexisting amorphous silica.Additionally, we report the presence of the ferrous hydroxycarbonate mineral chukanovite (Fe2(OH)2CO3)within some of the secondary mineral aggregates studied. Organic Matter:Indigenous Nakhla carbonaceous matter is present both in condensed phases and in a dispersed state spatially associated with secondary alteration phases. In the former, discrete refractory micron to submicron assemblages were identified that appear macromolecular in nature and, in several cases, associated with F and, in one case, significant N. While organically bound Fhas yet to be detected remotely on Mars, the identification of organic N within one of the carbonaceous phases is considered a prerequisite for Martian astrobiological potential. Summary: Neither the intimate association of carbonaceous matter with secondary phases nor the identification of chukanovite have been previously reported in any of the Martian meteorites.A variety of habitability-related metastable solid-phase, oxidation-reduction, and organic-molecular sample attributes in Nakhla formed hundreds of millions of years ago near Mars’ surface; these alteration products persisted on Mars until very recentlyand resemble the amorphous materials observed at Gale crater by Mars Science Laboratory CuriosityRover (Rampe et al., 2020).Similar phases may be expected to occur elsewhere on Mars (e.g., Jezerocrater). References: Rampe E. B. et al.(2020) Mineralogy and geochemistry of sedimentary rocks and eolian sediments in Gale crater, Mars: A review after six Earth years of exploration with Curiosity.Geochemistry80, 125605; Thomas-Keprta K.L. et al. (2022) Indigenous carbon-rich matter associated with unusual aqueous alteration features in Nakhla: Complex formation and preservation history. Geochim. et Cosmochim. Acta320, 41-78 (online).

Mars meteorite↗

NASA NTRS2006

The Need for High Spatial Resolution Multispectral Thermal Remote Sensing Data In Urban Heat Island Research

Although the study of the Urban Heat Island (UHI) effect dates back to the early 1800's when Luke Howard discovered London s heat island, it has only been with the advent of thermal remote sensing systems that the extent, characteristics, and impacts of the UHI have become to be understood. Analysis of the UHI effect is important because above all, this phenomenon can directly influence the health and welfare of urban residents. For example, in 1995, over 700 people died in Chicago due to heat-related causes. UHI s are characterized by increased temperature in comparison to rural areas and mortality rates during a heat wave increase exponentially with the maximum temperature, an effect that is exacerbated by the UHI. Aside from the direct impacts of the UHI on temperature, UHI s can produce secondary effects on local meteorology, including altering local wind patterns, increased development of clouds and fog, and increasing rates of precipitation either over, or downwind, of cities. Because of the extreme heterogeneity of the urban surface, in combination with the sprawl associated with urban growth, thermal infrared (TIR) remote sensing data have become of significant importance in understanding how land cover and land use characteristics affect the development and intensification of the UHI. TIR satellite data have been used extensively to analyze the surface temperature regimes of cities to help observe and measure the impacts of surface temperatures across the urban landscape. However, the spatial scales at which satellite TIR data are collected are for the most part, coarse, with the finest readily available TIR data collected by the Landsat ETM+ sensor at 60m spatial resolution. For many years, we have collected high spatial resolution (10m) data using an airborne multispectral TIR sensor over a number of cities across the United States. These high resolution data have been used to develop an understanding of how discrete surfaces across the urban environment (e.g., rooftops, pavements) interact from a surface-lower atmosphere energy flux perspective, to force the development of the UHI. Moreover, the airborne TIR sensor we used in our UHI studies was a multispectral sensor that had six channels in the 8-12pm range. The advantages of collecting multispectral TIR data became readily evident as a valuable tool for better calculation of unique surface thermal energy responses for urban materials over the 8-12 micrometer region, and also for getting a better handle on surface emissivity characteristics for these discrete surfaces. In this presentation, we will provide evidence on the virtues of how high spatial resolution multispectral TIR data can provide for better analysis of the UHI that cannot now be attained via TIR data obtained from satellites. Furthermore, we wish to provide compelling evidence on why future TIR satellite sensors should collect data at fine spatial resolutions (e.g. less than or equal to 30m) to better allow for measurement of surface thermal energy fluxes from discrete urban surfaces, and to better understand how surface fluxes from different urban materials in cities around the world in different climatic regimes, affect development of the UHI characteristics.

Quattrochi, Dale A.↗

NASA NTRS1988

Modeling of supersonic reacting flow fields

A detailed understanding of the scramjet combustor flow field is critical to the achievement of a successful design. Even though the combustor flow field is quite complex, it can be realistically viewed as a collection of spatially developing and reacting supersonic mixing layers that are initially discrete, but that ultimately merge into larger more complex zones. These mixing layers begin downstream of a set of fuel injectors that introduce gaseous hydrogen in both a parallel and transverse direction into a supersonic air stream entering from the engine inlet. The behavior of the initial portion of the combustor flow, in the mixing layers near the fuel injectors, appears to be most critical, since this is where the mechanism for efficient high speed mixing must be established to achieve the required degree of combustion downstream. Because of the structure of the flow field in this initial portion of the combustor, a single supersonic, spatially developing and reacting mixing layer serves as an excellent physical model for the overall flow field. Even though this reacting mixing layer flow is geometrically simple, it can still be made to retain all of the fluid mechanical and chemical complexities present in the actual combustor flow field.

Drummond, J. P.↗

NASA NTRS1990

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗

DOE OSTI2023

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Advances in hexagon mesh-based flow direction modeling

Watershed delineation and flow direction representation are the foundations of streamflow routing in spatially distributed hydrologic modeling. A recent study showed that hexagon-based watershed discretization has several advantages compared to the traditional Cartesian (latitude–longitude) discretization, such as uniform connectivity and compatibility with other Earth system model components based on unstructured mesh systems (e.g., oceanic models). Despite these advantages, hexagon-based discretization has not been widely adopted by the current generation of hydrologic models. One major reason is that there is no existing model that can delineate hexagon-based watersheds while maintaining accurate representations of flow direction across various spatial resolutions. In this study, we explored approaches such as spatial resampling and hybrid breaching-filling stream burning techniques to improve watershed delineation and flow direction representation using a newly developed hexagonal mesh watershed delineation model (HexWatershed). We applied these improvements to the Columbia River basin and performed 16 simulations with different configurations. The results show that (1) spatial resampling modulates flow direction around headwaters and provides an opportunity to extract subgrid information; and (2) stream burning corrects the flow directions in mountainous areas with complex terrain features.

58 GEOSCIENCES↗