Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Discrete sources”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Source Discretization Near Curvilinear Mesh Refinements in One Dimeions (Summer Internship Final Report)

Lawrence Livermore National Lab developed the seismic wave code SW4 for earthquake simulations, and has recently added the ability to specify mesh-refinement interfaces in the curvilinear mesh. The primary purpose of this is to allow for high resolution near the surface, with relatively fewer total points in the mesh leading to shorter computation time. Presently, the code does not allow for point-source forcing terms (commonly used in seismic simulations) to be located near these curvilinear mesh refinement interfaces. The goal of this internship was to investigate the discretization of point-source forcing terms in the seismic wave equation when the sources are located near the interfaces between two curvilinear meshes. Ultimately this will be incorporated into SW4 to allow for more realistic simulations to be carried out in this faster implementation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Theoretical and numerical studies of inverse source problem for the linear parabolic equation with sparse boundary measurements

We consider the inverse source problem in the parabolic equation, where the unknown source possesses the semi-discrete formulation. Theoretically, we prove that the flux data from any nonempty open subset of the boundary can uniquely determine the semi-discrete source. This means the observed area can be extremely small, and that is the reason we call it sparse boundary data. For the numerical reconstruction, we formulate the problem from the Bayesian sequential prediction perspective and conduct the numerical examples which estimate the space-time-dependent source state by state. To better demonstrate the method’s performance, we solve two common multiscale problems from two models with a long source sequence. The numerical results illustrate that the inversion is accurate and efficient.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Constraints on the origin of the radio synchrotron background via angular correlations

ABSTRACT The origin of the radio synchrotron background (RSB) is currently unknown. Its understanding might have profound implications in fundamental physics or might reveal a new class of radio emitters. In this work, we consider the scenario in which the RSB is due to extragalactic radio sources and measure the angular cross-correlation of Low-Frequency Array (LOFAR) images of the diffuse radio sky with matter tracers at different redshifts, provided by galaxy catalogues and cosmic microwave background lensing. We compare these measured cross-correlations to those expected for models of RSB sources. We find that low-redshift populations of discrete sources are excluded by the data, while higher redshift explanations are compatible with available observations. We also conclude that at least 20 per cent of the RSB surface brightness level must originate from populations tracing the large-scale distribution of matter in the Universe, indicating that at least this fraction of the RSB is of extragalactic origin. Future measurements of the correlation between the RSB and tracers of high-redshift sources will be crucial to constraining the source population of the RSB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2022

The channel-source hypothesis: Empirical evidence for in-channel sourcing of dissolved organic carbon to explain hysteresis in a headwater mountain stream

Catchment hydrologists have long puzzled over the question: How can catchments rapidly generate storm flows and pulses of solutes in response to storm events? Conceptual models viewing catchments as composed of discrete source areas generating flow at unique time scales and with unique chemical characteristics have been used to explain the observed changes in flow and water chemistry. Surprisingly, those conceptual models usually do not treat the stream channel as one of the potential source areas. We propose the channel-source hypothesis in which the stream itself should be considered as a potential source with the same rigour as other contributing areas. We pose this in the spirit of the scientific use of the word: a hypothesis is not a proven idea but rather a provisional supposition serving as the basis for further study. We suggest that the channel should be considered as a potential source for dissolved organic carbon (DOC). Channels store substantial amounts of organic matter, and stream ecologists have long studied stream carbon cycling. From those studies, we know that leaching and decomposition can generate DOC from particulate organic carbon (POC). Further, POC is stored in channel ‘dead-zones’—regions of low flow velocity - that can be activated as flow velocity increases, thus releasing accumulated DOC during storms. All catchments are different; there is no reason to assume that channel sources are always important, in every catchment, in every storm. Thus, the channel-source hypothesis does not replace existing conceptual models. Instead, it adds another potential mechanism that may explain DOC dynamics observed in streams. The channel-source hypothesis has substantial implications for catchment studies examining sources of DOC in stream water or using DOC as a tracer to determine the locations of, and proportional contributions of, different source areas for streamflow generation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Granular Flow

LIGGGHTS-INL is a capability-extended version of LIGGGHTS. LIGGGHTS is an Open Source Discrete Element Method Particle Simulation Software (https://www.cfdem.com/media/DEM/docu/Manual.html). LIGGGHTS stands for LAMMPS improved for general granular and granular heat transfer simulations. LAMMPS is a classical molecular dynamics simulator. It is widely used in the field of Molecular Dynamics (https://lammps.sandia.gov/doc/Manual.html). Thanks to physical and algorithmic analogies, LAMMPS is a very good platform for DEM simulations. LAMMPS offers a GRANULAR package to perform these kinds of simulations. LIGGGHTS aims to improve those capability with the goal to apply it to industrial applications.

Xia, Yidong↗
DOE OSTI2023

Biomass Feedstock Flow Modeling at Exascale

Variable biomass feedstock handling in hoppers and screw conveyors can cause upstream failures and disrupt downstream conversion operations. Understanding how particle size, shape, and moisture content affect feedstock flowability is important for efficient biorefinery operation. We present the use of our open-source discrete element method (DEM) solver for predicting flows in large scale hoppers and assess its performance on state-of-the-art high performance computing platforms that use heterogenous computing architectures.

bioenergy↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comparative Investigation of Current-Source Inverters using SiC Discrete Devices and Power Modules

The purpose of this paper is to investigate the impact of the SiC device packages on the commutation performance characteristics of current-source inverters (CSIs). The parasitic components in the CSI current commutation loop between the two-phase legs and output capacitors have a significant impact on the high-frequency performance of the SiC devices. To meet the elevated current requirements of high-power CSIs, it is often necessary to connect multiple discrete devices in parallel which increases the current commutation loop length. The selection of compact high-power SiC MOSFET and Schottky diode modules instead of discrete devices can be highly desirable to reduce the loop inductance and improve the system performance and power density. Two CSI benchtop prototype units, one with SiC discrete devices and the other with power modules, have been designed and tested, and the performances of the two CSIs are compared. The CSI with SiC power modules significantly reduces the inverter volume and lowers the parasitic inductance by 60% and the voltage ripple amplitude by 20% compared to the CSI unit with discrete devices.

42 ENGINEERING↗
DOE OSTI2025

SPADES (Scalable Parallel Discrete Events Simulation) [SWR-24-99]

SPADES (Solver for PArallel Discrete Event Simulation) is an open-source parallel discrete event simulation (PDES) package built on the AMReX library. Targeted at solving discrete event systems in parallel, this software package aims to be performance portable and scalable on heterogeneous computing architectures, e.g., graphic processing units (GPU). SPADES implements optimistic synchronization with rollback through an implementation of the Time Warp algorithm. An alternative conservative synchronization approach is also implemented using the Lower Bound on Incoming Time Stamp. In our implementation, logical processes are represented as cells in a grid and event messages are represented as particles. SPADES supports various parallel decomposition strategies, including the use of the Message Passing Interface (MPI) and OpenMP threading. All major GPU architectures (e.g., Intel, AMD, NVIDIA) are supported through the use of performance portability functionalities implemented in AMReX. The SPADES software is released in NREL Software Record SWR-24-99 “SPADES (Scalable Parallel Discrete Events Simulation)”.

Henry de Frahan, Marc [National Renewable Energy L↗
DOE OSTI2020

High Order Implicit Residual-Based Spatial Discretization Error Estimation for S N Neutron Transport

This work demonstrates our novel residual source spatial discretization error estimator (LeR/TEAD) for a DGFEM-1 discretization and assesses it along with two contemporary estimators, Ragusa and Wang's h -refinement estimator (RW) and Duo, Azmy, and Zikatanov's explicit residual-based estimator (DAZ), on a suite of Method of Manufactured Solutions (MMS) 2D problems and three realistic problem geometries. LeR/TE-AD is attractive because it directly estimates the local error in the angular flux, as opposed to a mere indicator of the error's behavior, on the same mesh and method order as the original numerical solution, thus typically being less computationally intensive than a refinement-based method. On the MMS suite, LeR/TE-AD consistently displayed a reduced performance versus its DGFEM-0 results in terms of accuracy and precision metrics, though it was not typically grossly inaccurate. This is attributed to the irregularities in the true solution across singular characteristics limiting the local accuracy of the numerical flux solution, leading to poor derivative approximations used in the residual approximations. The error transport problem then spreads the error in the residual to nearby cells, causing a greater degree of imprecision that did not afflict DAZ or RW. In testing the estimators on realistic problem geometries, however, LeR/TE-AD fared better. In practice, the true error is much larger in non-idealized geometries like in MMS, and a superlinear true solution means that RW and DAZ are not beneficially biased for DGFEM-1 error estimation. LeR/TE-AD was typically first or second in accuracy, primarily competing with RW, but the latter usually consumed 2-4 times the computational time as LeR/TE-AD, and requires a solution with four times as many unknowns. Furthermore, RW and LeR/TE-AD can be used to compute direct estimates of the error in any quantity of interest that is based on the angular ux solution, such as the fission rate density in a fuel pin, whereas DAZ requires a heuristic extension due to its norm-based nature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗
DOE OSTI2023

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Thermal Management of Discretized Heaters Using CuW Microchannel Heat Sinks and FC3283

Single-phase cooling using microchannel heat sinks (MCHS) has become a popular approach for overcoming the thermal challenges associated with high-powered microelectronic devices. Thermal management is one of the largest barriers to higher power densities in electronics and frequently limits overall device performance. The implementation of forced convective cooling via single-phase liquid cooling in MCHS reduces the thermal resistance resulting in lower device temperatures at high-power conditions, which can decrease the package size and extend the lifespan of devices. The goal of the work described in this article was to investigate practical cooling solutions for laser diode bars. This study examined the effectiveness of a copper tungsten (CuW) microchannel heat sink paired with a dielectric coolant (FC3283) for dissipating both discrete and uniform heat fluxes up to 600 W/cm2 across a 0.25-cm2 surface area through a numerical and experimental study. CuW was chosen as the MCHS material because it is thermal expansion matched to GaAs, which is a common laser diode substrate. FC3283 serves as a dielectric coolant that is compatible with power electronics cooling. The microchannels utilized in this work were approximately $365$ $μ$ m deep and $100$ $μ$ m wide resulting in a hydraulic diameter of $160$ $μ$ m. The study investigated the cooling performance across a variety of applied power loads and flow rates to determine optimal effectiveness. The resulting thermal resistance ranged from 0.15 cm2 K/W at the highest flow rate to 0.26 cm2 K/W at the lowest flow rate. Further, the resulting thermal performance from this study demonstrates the importance of considering discrete heat sources separately from uniform heat sources and proved that the unique combination of CuW microchannels with FC3283 can be a promising cooling option toward future advancements of laser diode bars and other high-power microelectronics.

42 ENGINEERING↗
DOE OSTI2023

A Hybrid Physics-Based, Data-Driven Approach to Model Damage Accumulation in Corrosion of Polymeric Adhesives

The degradation of polymeric adhesives in dissimilar joints is a menace to electric and structural components of the automotive industry. In polymeric adhesives, failure often occurs due to damage accumulated from discrete environmental sources—thermal degradation, oxidation, hydrolysis, radiation, and moisture—that are combined with progressive damage mechanisms, such as aging and fatigue. Because most damage mechanisms that lead to degradation act in parallel, the estimation of the lifetime of a polymeric material becomes a significant reliability issue, no study currently exists that can simultaneously consider even two of these phenomena. Most current efforts are directed toward a single-phenomenon, single-model approach in which a model is developed to cover one specific damage mechanism. While this approach is advantageous in exploring single mechanisms, it is strongly inconclusive for degradation that occurs by multiple damage mechanisms. Our team plans to take a fresh perspective by designing a hybrid constitutive/data-driven approach that allows models of different damage mechanisms to be integrated into one model and then used to explore the degradation-induced failure of polymeric adhesives. The theoretical-computational part of this work is divided into four research components as described in the Approach section below.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First Direct Observations of Propagation of Discrete Chorus Elements From the Equatorial Source to Higher Latitudes, Using the Van Allen Probes and Arase Satellites

Whistler mode chorus waves have recently been established as the most likely candidate for scattering relativistic electrons to produce the electron microbursts observed by low altitude satellites and balloons. These waves would have to propagate from the equatorial source region to significantly higher magnetic latitude in order to scatter electrons of these relativistic energies. This theoretically proposed propagation has never been directly observed. We present the first direct observations of the same discrete rising tone chorus elements propagating from a near equatorial (Van Allen Probes) to an off-equatorial (Arase) satellite. The chorus is observed first on the more equatorial satellite and is found to be more oblique and significantly attenuated at the off-equatorial satellite. Further, this is consistent with the prevailing theory of chorus propagation and with the idea that chorus must propagate from the equatorial source region to higher latitudes. Ray tracing of chorus at the observed frequencies confirms that these elements could be generated parallel to the field at the equator, and propagate through the medium unducted to Van Allen Probes A and then to Arase with the observed time delay, and have the observed obliquity and intensity at each satellite.

79 ASTRONOMY AND ASTROPHYSICS↗