Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale features”

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 199 records · Page 11

Adaptive Mesh Refinement Simulations for Turbulent Reacting Flow

With the increased availability of exascale computing hardware, detailed simulations of realistic devices can be performed at practically relevant time and length scales. Insights into the multiscale driving mechanisms in compressible reacting flow systems with complex geometry, such as combustors, can be used for design optimization and technology improvements. However, to effectively perform these simulations, advanced numerical algorithms must be used to maintain solution accuracy without incurring undue computational costs. PeleC, part of the Pele suite of codes, leverages block-structured adaptive mesh refinement (AMR) through the AMReX library to capture fine-scale flow features in compressible reacting flows. In this talk, we discuss recent improvements to the numerical algorithms, particularly in regard to describing flows at complex boundary structures, and PeleC's performance on exascale computing hardware. We will demonstrate that PeleC is well-suited for modern, extreme-scale, heterogenous compute platforms.

combustion↗

Parallel Multiscale Algorithms for Astrophysical Fluid Dynamics Simulations

Our goal is to develop software libraries and applications for astrophysical fluid dynamics simulations in multidimensions that will enable us to resolve the large spatial and temporal variations that inevitably arise due to gravity, fronts and microphysical phenomena. The software must run efficiently on parallel computers and be general enough to allow the incorporation of a wide variety of physics. Cosmological structure formation with realistic gas physics is the primary application driver in this work. Accurate simulations of e.g. galaxy formation require a spatial dynamic range (i.e., ratio of system scale to smallest resolved feature) of 104 or more in three dimensions in arbitrary topologies. We take this as our technical requirement. We have achieved, and in fact, surpassed these goals.

Norman, Michael L.↗

Computed Tomography Scanning and Geophysical Measurements of Appalachian Basin Core from the Jones and Laughlin #1 Well, Beaver County, PA

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize Appalachian Basin core from Beaver County, Pennsylvania. The primary impetus of this work is a collaboration between the U.S. Department of Energy (DOE) and the Pennsylvania Geological Survey to characterize and make publicly available core information from the Onondaga-Huntersville formations of the Appalachian Basin. This stratigraphic well and the core data produced in this report will aid in understanding the structural complexities of the Onondaga-Huntersville formations. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/jonesandlaughlin1well. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with imaging performed with the NETL medical CT scanner over the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for further study. Targeted higher resolution CT scanning of select sections was performed with NETL’s micro-CT scanner. The combination of methods used provides a multiscale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

58 GEOSCIENCES↗

Development of Efficient Process for Manufacturing of Thermoplastic Composites with Tailored Properties (CRADA 511)

This is a collaborative effort between Battelle Memorial Institute as manager and operator of Pacific Northwest National Laboratory (PNNL) and ESI North America Inc. (“ESI” or “Participant”) to apply computation and data analytics to the challenge of light weighting with a focus on the battery enclosures of electric vehicles (EVs). EVs use heavy batteries to increase range and power. A complex-shaped battery enclosure is required to meet a host of challenging performance requirements. The ability to virtually develop composite parts such as battery enclosure with tailored properties to meet required performance will be highly valuable to the automotive industry. However, efficient simulation of composite-manufacturing processes remains a challenging issue since simulation involves multiscale models in space and time, highly non-linear and anisotropic behavior, strongly coupled multi-physics, and complex geometries. This work will advance the state of the art by reducing the computational burden of composite optimization by using simulation data from a limited number of configurations off-line and then developing a reduced order model (ROM) using data analytics and machine learning (ML). Develop a data driven approach to link features of the material and manufacturing processes to the mechanical properties of thermoplastic composite parts.

42 ENGINEERING↗

Paper or Plastic? Multiscale Material Handling Properties of Two Model Municipal Solid Waste Streams

Purpose: Municipal Solid Waste (MSW) is a potentially valuable sustainable feedstock for fuel and chemical production due to its carbon-rich content and low cost. This study aims to assess the material handling properties of paperand plastic-rich MSW feedstocks to mitigate equipment failure and processing downtime. Methods: The material handling properties of crumbled MSW feedstocks were measured using apowder rheometer with mass flow hopper calculations to assess handling performance. Inverse gas chromatography was use to measure the surface energy differences between feedstocks. Electron microscopy and Raman spectroscopy was used to evaluate microscale features that may contribute to material handling differences. Results: Plastic and paper rich feedstocks crumbled to a nominal 2 mm particle size were observed to have similar flow and handling characteristics with reasonable hopper outlets. 2 mm plastic rich crumbles, with their higher bulk density, exhibited superior flow performance. By contrast, 4 mm material required significantly larger hopper outlets, indicating poor flowability. Paper rich and 4 mm plastic rich samples displayed broad particle size distributions, which contributed to particle interlocking, jamming, and other flow issues. Electron microscopy revealed that plastic rich samples were significantly smoother, enhancing their flowability compared to the rougher, paper rich materials. Conclusions: This study establishes critical material handling baselines for processing MSW as a viable feedstock for fuel and chemical production. The findings highlight the importance of optimizing particle size and feedstock composition to improve flowability and handling performance.

09 BIOMASS FUELS↗

WUS256: An Adjoint Waveform Tomography Model of the Crust and Upper Mantle of the Western United States for Improved Waveform Simulations

Abstract We report a new model (WUS256) of radially anisotropic seismic wavespeeds of the crust and upper mantle of the western United States (WUS) obtained from adjoint waveform tomography for the purpose of improving synthetic waveform fits to observed data. WUS256 is based on inversion of over 94,000 waveforms from 72 earthquakes recorded by nearly 3,400 stations. We started with the SPiRaL global model (Simmons et al., 2021, https://doi.org/10.1093/gji/ggab277 ) and waveforms in the period band of 50–120 s. We followed a conservative multiscale inversion approach with eight stages and 256 total inversion iterations which enabled monotonic misfit reduction to 20‐s minimum‐period waves. WUS256 relied on time‐frequency (TF) phase misfits and a trust region limited memory Broyden–Fletcher–Goldfarb–Shanno (L‐BFGS) optimization. Hessian‐vector products were used to qualitatively assess model resolution. Results indicate that WUS256 has good coverage of the continental regions to depths of about 150 km and is able to resolve features on lateral scales of about 200 km. We quantify waveform fits by the reduction in TF and normalized amplitude difference misfits between WUS256 and the SPiRaL starting model. WUS256 significantly improves waveform fits with misfit reduction 64% for both inversion and validation data sets compared to the SPiRaL starting model and shows even better fits compared to other models. Waveform fits illustrate that WUS256 reproduces body‐waves, fundamental mode surface waves as well as late arriving dispersed and/or scattered short period surface waves. The improvement in waveform fit indicates that WUS256 can be used to reproduce path effects on regional complete waveforms and moment tensor inversions.

58 GEOSCIENCES↗

Multiscale electric-field imaging of polarization vortex structures in PbTiO 3 /SrTiO 3 superlattices

In ferroelectric heterostructures, the interaction between intrinsic polarization and the electric field generates a rich set of localized electrical properties. The local electric field is determined by several connected factors, including the charge distribution of individual unit cells, the interfacial electromechanical boundary conditions, and chemical composition of the interfaces. However, especially in ferroelectric perovskites, a complete description of the local electric field across micro-, nano-, and atomic-length scales is missing. Here, by applying four-dimensional scanning transmission electron microscopy (4D STEM) with multiple probe sizes matching the size of structural features, we directly image the electric field of polarization vortices in (PbTiO 3 ) 16 /(SrTiO 3 ) 16 superlattices and reveal different electric field configurations corresponding to the atomic scale electronic ordering and the nanoscale boundary conditions. The separability of two different fields probed by 4D STEM offers the possibility to reveal how each contributes to the electronic properties of the film.

36 MATERIALS SCIENCE↗

On the use of a multigrid-reduction-in-time algorithm for multiscale convergence of turbulence simulations

Simulations of turbulent flow present challenges in terms of accuracy and affordability on modern highly-parallel computer architectures. A multigrid-reduction-in-time algorithm is used to provide a framework for separately evolving different scales of turbulence and for parallelizing the temporal domain, thereby increasing the concurrency. It is hypothesized that the space–time locality of the small scales of turbulence can be used to circumvent difficulties in applying temporal multigrid to flows dominated by inertial physics. For algorithms that fall well short of spectral accuracy (fourth-order is used in this work) attention must be paid to the accuracy of features on scales transferred between multigrid levels. Numerical experiments were performed using implicit large-eddy simulation. Results from applying the approach to an infinite-Reynolds number Taylor–Green flow and a double-shear flow at a Reynolds number of 11650 provide strong evidence that the approach has merit. The multigrid-reduction-in-time framework can be used to parallelize the temporal domain of a high-Reynolds-number turbulent flow and permit independent convergence of different scales. Establishing this foundation allows for future research in reducing the wall-clock time to solve turbulent flows while retaining the same accuracy as sequential solvers. In conclusion, current performance results from parallelizing the temporal domain are not competitive with those from sequential-in-time methods.

97 MATHEMATICS AND COMPUTING↗

Hierarchical median narrow band for level set segmentation of cervical cell nuclei

This paper presents a novel hierarchical nuclei segmentation algorithm for isolated and overlapping cervical cells based on a narrow band level set implementation. Our method applies a new multiscale analysis algorithm to estimate the number of clusters in each image region containing cells, which turns into the input to a narrow band level set algorithm. We assess the nuclei segmentation results on three public cervical cell image databases. Overall, our segmentation method outperformed six state-of-the-art methods concerning the number of correctly segmented nuclei and the Dice coefficient reached values equal to or higher than 0.90. We also carried out classification experiments using features extracted from our segmentation results and the proposed pipeline achieved the highest average accuracy values equal to 0.89 and 0.77 for two-class and three-class problems, respectively. Furthermore, these results demonstrated the suitability of the proposed segmentation algorithm to integrate decision support systems for cervical cell screening.

47 OTHER INSTRUMENTATION↗

Lower Hybrid Drift Waves During Guide Field Reconnection

Generation and propagation of lower hybrid drift wave (LHDW) near the electron diffusion region (EDR) during guide field reconnection at the magnetopause is studied with data from the Magnetospheric Multiscale mission and a theoretical model. Inside the current sheet, the electron beta (β e ) determines which type of LHDW is excited. Inside the EDR, where the electron beta is high (β e ~5), the long-wavelength electromagnetic LHDW is observed propagating obliquely to the local magnetic field. In contrast, the short-wavelength electrostatic LHDW, propagating nearly perpendicular to the magnetic field, is observed slightly away from the EDR, where β e is small (~0.6). These observed LHDW features are explained by a local theoretical model, including effects from the electron temperature anisotropy, finite electron heat flux, electrostatics, and parallel current. The short-wavelength LHDW is capable of generating significant drag force between electrons and ions.

79 ASTRONOMY AND ASTROPHYSICS↗

A reliable workflow for improving nanoscale X-ray fluorescence tomographic analysis on nanoparticle-treated HeLa cells

Abstract Scanning X-ray fluorescence (XRF) tomography provides powerful characterization capabilities in evaluating elemental distribution and differentiating their inter- and intra-cellular interactions in a three-dimensional (3D) space. Scanning XRF tomography encounters practical challenges from the sample itself, where the range of rotation angles is limited by geometric constraints, involving sample substrates or nearby features either blocking or converging into the field of view. This study aims to develop a reliable and efficient workflow that can (1) expand the experimental window for nanoscale tomographic analysis of local areas of interest within a laterally extended specimen, and (2) bridge 3D analysis at micrometer and nanoscales on the same specimen. We demonstrate the workflow using a specimen of HeLa cells exposed to iron oxide core and titanium dioxide shell (Fe3O4/TiO2) nanocomposites. The workflow utilizes iterative and multiscale XRF data collection with intermediate sample processing by focused ion beam (FIB) sample preparation between measurements at different length scales. Initial assessment combined with precise sample manipulation via FIB allows direct removal of sample regions that are obstacles to both incident X-ray beam and outgoing XRF signals, which considerably improves the subsequent nanoscale tomography analysis. This multiscale analysis workflow has advanced bio-nanotechnology studies by providing deep insights into the interaction between nanocomposites and single cells at a subcellular level as well as statistical assessments from measuring a population of cells.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗

Explicit physics-informed neural networks for nonlinear closure: The case of transport in tissues

In upscaling methods, closures for nonlinear problems present a well-known challenge. While a number of theoretical methods have been proposed for handling such closures, nonlinearities still remain a significant obstacle for many problems. In this work, we use a combination of formal upscaling and data-driven machine learning for explicitly closing a nonlinear transport and reaction process in multiscale tissues. The classical effectiveness factor model is used to formulate the macroscale reaction kinetics. We train a multilayer perceptron network using training data generated by direct numerical simulations over microscale examples. Once trained, the network is used in an algorithm for numerically solving the upscaled (coarse-grained) differential equation describing mass transport and reaction in two example tissues. The network is described as being explicit in the sense that the network is trained using macroscale concentrations and gradients of concentration as components of the feature space rather than incorporating them as part of a constraint in the optimization process. Network training and solutions to the macroscale transport equations were computed for two different tissues. The two tissue types (brain and liver) exhibit markedly different geometrical complexity and spatial scale (cell size and sample size). The upscaled solutions for the average concentration are compared with numerical solutions derived from the microscale concentration fields by a posteriori averaging. There are three outcomes of this work of particular note. 1) Our overall approach results in an upscaled nonlinear PDE. The PDE is closed using a neural network, and our approach results in the definition of the classical effectiveness factor for effecting closure. 2) We identify particular source terms for the closure problem that are important for representing the structure of the closure. These source terms involve macroscale concentrations and their gradients. We adopt these source terms to use as explicit features in the learning algorithm. We find the trained networks that include the macroscale source terms generate models that are able to predict the correction factor with increased fidelity over those that do not. 3) We find that the trained network exhibits good generalizability, and it is able to predict the effectiveness factor with high fidelity for realistically-structured tissues despite the significantly different scale and geometrical complexity of the two example tissue types. This latter result emphasizes our purposeful connection between conventional averaging methods with the use of machine learning for closure; this contrasts with some machine learning methods for upscaling where the exact form of the macroscale equation remains unknown.

97 MATHEMATICS AND COMPUTING↗

River Delta Morphotypes Emerge From Multiscale Characterization of Shorelines

Abstract Delta shoreline structure has long been hypothesized to encode information on the relative influence of fluvial, wave, and tidal processes on delta formation and evolution. We introduce here a novel multiscale characterization of shorelines by defining three process‐informed morphological metrics. We show that this characterization yields self‐emerging classes of morphologically similar deltas, that is, delta morphotypes, and also predicts the dominant forcing of each morphotype. Then we show that the dominant forcings inferred from shoreline structure generally align with those estimated via relative sediment fluxes, while positing that misalignments arise from spatiotemporal heterogeneity in deltaic sediment fluxes not captured in their estimates. The proposed framework for shoreline characterization advances our quantitative understanding of how shoreline features reflect delta forcings, and may aid in deciphering paleoclimate from images of ancient deposits and projecting delta morphologic response to changes in sediment fluxes.

Vulis, L.↗

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗

Generation of Electron Whistler Waves at the Mirror Mode Magnetic Holes: MMS Observations and PIC Simulation

The Magnetospheric Multiscale mission has observed electron whistler waves at the center and at the edges of magnetic holes in the dayside magnetosheath. The magnetic holes are nonlinear mirror structures since their magnitude is anticorrelated with particle density. In this article, we examine the growth mechanisms of these whistler waves and their interaction with the host magnetic hole. In the observations, as magnetic holes develop and get deeper, an electron population gets trapped and develops a temperature anisotropy favorable for whistler waves to be generated. In addition, the decrease in magnetic field magnitude and the increase in density reduce the electron resonance energy, which promotes the electron cyclotron resonance. To investigate this process, we used expanding box particle-in-cell simulations to produce the mirror instability, which then evolve into magnetic holes. The simulation shows that whistler waves can be generated at the center and edges of magnetic holes, which reproduces the primary features of the MMS observations. The simulation shows that the electron temperature anisotropy develops in the center of the magnetic hole once the mirror instability reaches its nonlinear stage of evolution. The plasma is then unstable to whistler waves at the minimum of the magnetic field structures. In the saturation regime of mirror instability, when magnetic holes are developed, the electron temperature anisotropy appears at the edges of the holes and electron distributions become more isotropic at the magnetic field minimum. At the edges, the expansion of magnetic holes decelerates the electrons, which leads to temperature anisotropies.

Magnetospheric Multiscale↗

Impact of Material and Architecture Model Parameters on the Failure of Woven Ceramic Matrix Composites (CMCs) via the Multiscale Generalized Method of Cells

It is well known that failure of a material is a locally driven event. In the case of ceramic matrix composites (CMCs), significant variations in the microstructure of the composite exist and their significance on both deformation and life response need to be assessed. Examples of these variations include changes in the fiber tow shape, tow shifting/nesting and voids within and between tows. In the present work, the effects of many of these architectural parameters and material scatter of woven ceramic composite properties at the macroscale (woven RUC) will be studied to assess their sensitivity. The recently developed Multiscale Generalized Method of Cells methodology is used to determine the overall deformation response, proportional elastic limit (first matrix cracking), and failure under tensile loading conditions. The macroscale responses investigated illustrate the effect of architectural and material parameters on a single RUC representing a five harness satin weave fabric. Results shows that the most critical architectural parameter is weave void shape and content with other parameters being less in severity. Variation of the matrix material properties was also studied to illustrate the influence of the material variability on the overall features of the composite stress-strain response.

Liu, Kuang C.↗