Search NASA⌕ Search

SEARCH · Search NASA

Results for “avoidance”

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 721 records · Page 40

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗

Peridynamic modeling of cementitious materials for nuclear waste management

Radioactive and hazardous waste generated from fuel processing plants, nuclear reactors, and hospitals, requires effective management strategies. Cementitious materials are widely applied for these needs, serving as structural materials, reactive barriers, or waste forms. In these applications cracking poses a significant risk to performance, driven by inconsistent shrinkage of the components and the varying strength and permeability of their interfaces. Traditional modeling approaches face challenges in representing the complex fracture behavior of cementitious materials due to the reliance on spatial derivatives and difficulties with mesh generation. Here, to overcome these limitations, we employ peridynamics, a novel continuum mechanics formulation that uses integrals to describe mechanical equilibrium, avoiding discontinuities associated with traditional methods. Through incorporation of a bi-linear softening model and quasistatics, an experimentally validated model for Portland cement concrete samples was created. Mechanical parameters, including compressive strength and elastic modulus were validated and variation due to aggregate packing was evaluated. Additionally, sensitivity analysis of the peridynamic parameters to the bulk material properties was established. The results lay the groundwork for evaluating the impact of unique conditions of cementitious waste forms that can be assessed to improve the reliability of waste management strategies.

Aggregates↗

Thickness threshold study of polyaniline-based lightning strike protection coating for carbon/glass fiber reinforced polymer composites

Different polyaniline (PANI) based lightning strike protection (LSP) coating thicknesses were tested against simulated lightning strike tests to evaluate the threshold thickness of the coating. A series of lightning strike tests were performed on carbon/glass fiber reinforced polymer (CFRP/GFRP) composites subjected to simulated lightning currents ranging from 40 to 100 kA. Five LSP-coated specimens were considered herein: two CFRP and three GFRP composites. A minimum of ~ 345 μm thick PANI-based LSP coating on CFRP and GFRP composites under a 40 kA peak current retained 83 % and 100 % residual strength, respectively. Additionally, GFRP specimens tested against 60 and 100 kA peak currents exhibited 100 % residual strength. Finally, this study proposes that a PANI-based LSP coating is suitable for protecting GFRP composites, while the careful design is necessary for CFRP composites to avoid any catastrophic failure.

36 MATERIALS SCIENCE↗

A critical verification of beam and shell models of wind turbine blades

Ever-increasing wind turbine size has challenged predictive capabilities on several fronts. Here, to address part of the blade structural modeling uncertainty, a systematic model fidelity comparison study was conducted on commonly used finite elements. pyNuMAD was utilized to create beam, shell, and solid models of a 100 m long blade undergoing large static deflections. The solid model avoided the use of layered-solid elements by resolving core and facesheet layers. An unprecedented model with 73.7 million elements revealed insights that have never been possible from prior experimental and numerical studies. As compared to the solid element model, the tip deflection from the shell and beam model was found to be about 2% and 4.3% too low, respectively. The twist from the beam model was found to be about 5.6% too high, while the twist from shell model was 24% too low, though improvement was demonstrated with mesh refinement. The beam model adhesive stresses were more accurate than the shell model. Out-of-plane stresses were of great significance near geometric and material discontinuities, and neither the shell nor beam model captured these effects well. Failure predictions from beam, shell, or layered-solid models are unlikely to be reliable at trailing edges, adhesives, ply-drops, spar-cap boundaries.

17 WIND ENERGY↗

Historic and modern nuclear graphite impurities: Pathways to improved waste strategies

Graphite has been used in large volumes as a structural material and neutron moderator since the earliest days of nuclear fission. However, no international consensus exists on the disposal of irradiated graphite, leaving much of the historic radioactive graphite inventory in interim vault or silo storage. With several new graphite-moderated reactors planned or under construction, the issue of graphite waste management is becoming increasingly urgent. This paper reviews and quantifies impurities in both historic and modern nuclear graphite, with emphasis on nitrogen—responsible for much of the 14 C inventory—and chlorine, which plays a critical role in repository performance and design. Modern graphites, benefitting from stringent quality-control measures developed for non-nuclear industries, meet or exceed the ASTM Ultra-High Purity nuclear standards, even without halide purification. Both chlorine and nitrogen concentrations have declined over time. For chlorine, identified as a key impurity influencing U.S. waste repository design, we propose a target of 0.1 appm in as-fabricated billets as a reasonable benchmark. Nitrogen sources are traced throughout the graphite production process, with surface and bulk concentrations characterized for all materials studied. Modern graphites commonly exhibit nitrogen levels below 5 appm, with values approaching 1 appm achievable. Using such reduced-nitrogen grades is critical to keeping graphite-induced radioactivity below the greater-than-Class-C waste threshold, thereby avoiding disposal cost penalties of nearly an order of magnitude.

Chlorine↗

Efficient Fourier transforms for transverse momentum dependent distributions

Hadron production at low transverse momenta in semi-inclusive deep inelastic scattering can be described by transverse momentum dependent (TMD) factorization. This formalism has also been widely used to study the Drell-Yan process and back-to-back hadron pair production in $e^+e^-$ collisions. These processes are the main ones for extractions of TMD parton distribution functions and TMD fragmentation functions, which encode important information about nucleon structure and hadronization. One of the most widely used TMD factorization formalism in phenomenology formulates TMD observables in coordinate $b_\perp$-space, the conjugate space of the transverse momentum. The Fourier transform from $b_\perp$-space back into transverse momentum space is sufficiently complicated due to oscillatory integrands that it requires a careful and computationally intensive numerical treatment in order to avoid potentially large numerical errors. Within the TMD formalism, the azimuthal angular dependence is analytically integrated and the two-dimensional $b_\perp$ integration reduces to a one-dimensional integration over the magnitude $b_\perp$. In this paper we develop a fast numerical Hankel transform algorithm for such a $b_\perp$-integration that improves the numerical accuracy of TMD calculations in all standard processes. Libraries for this algorithm are implemented in Python 2.7 and 3, C++, as well as FORTRAN77. All packages are made available open source.

97 MATHEMATICS AND COMPUTING↗

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

De-noising drift chambers in CLAS12 using convolutional auto encoders

Modern Nuclear Physics experimental setups run experiments with higher beam intensity resulting in increased noise in detector components used for particle track reconstruction. Increased uncorrelated signals (noise) result in decreased particle reconstruction efficiency. In this paper, we investigate the usage of Machine Learning, specifically Convolutional Neural Network Auto-Encoders (CAE), for de-noising raw hits from drift chambers in the CLAS12 detector. To the best of our knowledge, this is the first time CAE is employed to perform such an operation in this field. During the de-noising phase, it is important to remove as much noise as possible while retaining the valid hits to avoid losing crucial information about the experiment. Here, we show that using CAE, it is possible to remove noise hits while retaining up to 94% of valid tracks for a beam current of 110nA while for lower beam currents (45-55nA), we get up to 98% efficiency. Studies on experimental conditions with increasing noise show that CAE performs better than conventional tracking algorithms in isolating hits belonging to tracks. Specifically, the de-noising algorithm results in tracking efficiency improvements greater than 15%, in real data production procedures with nominal conditions, and up to two times better efficiency in synthetically generated data with high luminosity conditions (90-110nA), indicating that machine learning can lead to significantly shorter times for conducting physics experiments.

97 MATHEMATICS AND COMPUTING↗

HiPACE++: A portable, 3D quasi-static particle-in-cell code

Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropriate memory management, HiPACE++ demonstrates orders-of-magnitude speedups on modern scientific GPUs over CPU-only implementations. The 2D transverse operations are performed on a single GPU, avoiding time-consuming communications. The longitudinal parallelization is done through temporal domain decomposition, enabling near-optimal strong scaling from 1 to 512 GPUs. HiPACE++ is a modular, open-source code enabling efficient modeling of plasma accelerators from laptops to state-of-the-art supercomputers.

97 MATHEMATICS AND COMPUTING↗

GenASiS Mathematics: Object-oriented manifolds, operations, and solvers for large-scale physics simulations (version 2)

We report GenASiS Mathematics provides modern Fortran classes furnishing extensible object-oriented functionality for the solution of fields governed by selected partial differential equations. The initial release included extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon. This revision - Version 2 of Mathematics - includes significant reorganization and streamlining of these classes, higher-order reconstruction by a different method, a Poisson solver, coarsening to avoid Courant time step limitations near coordinate singularities, and the offloading of computational kernels to GPUs.

97 MATHEMATICS AND COMPUTING↗

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

General solution to the Kohn–Luttinger nonconvergence problem

A simple, but general solution is proposed for the Kohn–Luttinger problem, i.e., the nonconvergence of the finite-temperature many-body perturbation theory with its zero-temperature counterpart as temperature is lowered to zero under some circumstances. How this nonconvergence can be avoided by altering the reference wave function is illustrated numerically by using up to the fifth order of the perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗

An editorial to the Special Issue on “Severe climate Risks”

The history of this Special Issue (https://www.sciencedirect.com/special-issue/10JD7LNJNQ0) indirectly dates back to the early 1990s, when the signature of the United Nations Framework Convention on Climate Change kicked-off an international political process based on one overarching and foundational principle: to avoid “dangerous anthropogenic interference with the climate system” at the global level. More than three decades later, such a principle remains central, though complementary aims made their way through the climate negotiation process, such as the importance of ensuring equity and justice, to give just one example here. Scientific knowledge also considerably progressed and we know more about the range of risks that climate change imposes and will continue to impose to the biosphere and humankind, worldwide and at all territorial levels. It is also clear that societal responses to these risks —“climate adaptation” as we know it— are increasingly happening, but definitely not at the pace of climate risk trends (Berrang-Ford et al., 2021, Erisken et al., 2021, Olazabal and Ruiz De Gopegui, 2021, Magnan et al., 2023a, Reckien et al., 2023, UNEP, 2023). As a result, concerns have emerged over the recent years in both the scientific and policy arenas around the idea that societies may not be able to address all climate risks, and that limits to adaptation and induced residual risks need to be considered more seriously. Such concerns further highlight the continuing importance of the imperative to minimise dangerous anthropogenic interference with the climate system, at any scale. But what does “dangerous interference” mean? How can we decide that we are entering the “dangerous” space, compared to a broader range of climate risks that would qualify as problematic but not necessarily “dangerous”? Who should make such a decision? Which conditions drive risk severity over time, including in the future? And what would be the environmental, economic, social and cultural implications of prioritising some climate risks over others? The Intergovernmental Panel on Climate Change (IPCC) was a pioneer in addressing such questions through the development of the “Key Risks” framing that describes those climate risks having the potential to become dangerous or “severe” over the course of this century (Schellnhuber et al., 2006, Schneider et al., 2007, Oppenheimer et al., 2014, O’Neill et al., 2022). The Fifth and Sixth assessment cycles (AR5 and AR6) went a step further by identifying about 120 Key Risks across regions and sectors, and clustering them into 8 “Representative Key Risks” covering a range of geographical systems (low-lying coasts, and to terrestrial and ocean ecosystems), sectors (critical infrastructure, living standards, human health, food security, and water security) and human dimensions (peace and mobility) (Oppenheimer et al., 2014, O’Neill et al., 2022). This Special Issue was born of the efforts of a range of authors, during the development of the IPCC AR6 main Assessment Report between 2019 and 2022, to characterise Key Risks and Representative Key Risks, and advance knowledge on what shapes “severe climate risks” conceptually as well as in the real-world. The series of papers forming this Special Issue is not intended to cover the topic exhaustively, but rather to give readers an overview through the following narrative: defining “severe climate risks” is highly challenging (Magnan et al., 2023b), but knowledge is expanding on the driving climate hazards (Tebaldi et al., 2023) and their implications on geographical systems, sectors and human dimensions, using here food security (Mirzabaev et al., 2023), human mobility (Gilmore et al., 2024) and peace (Buhaug et al., 2023) as illustrative examples. The overall intention is to support especially decision-makers, whatever the scale or sector considered, in asking severity-driven questions to identify sector- and context-specific “priority” risks from climate change.

54 ENVIRONMENTAL SCIENCES↗

Comparative techno-economic analysis of synthetic renewable natural gas production via reactive CO 2 capture and conversion

Reactive CO 2 capture and conversion (RCC) is an emerging carbon management strategy that integrates CO 2 capture and conversion and avoids intermediate CO 2 purification. In this study, we design an RCC process to capture atmospheric CO 2 and react it with renewable hydrogen to produce synthetic renewable natural gas (SRNG), which serves as a carbon-neutral energy source and a chemical form of long-duration renewable energy storage. We assess the technological potential of RCC through process modeling, techno-economic, carbon footprint, and sensitivity analyses. Our findings demonstrate that RCC offers energy savings and comparable cost to separated capture and conversion processes. The cost is dominated by renewable hydrogen and material replacement cost. SRNG produced via RCC is competitive with existing low-carbon natural gas technologies and presents a promising low-cost option for long-duration energy storage. This work highlights the potential for deploying RCC technologies within a circular carbon economy and the scientific and technical challenges that must be overcome for material and technology developers.

03 NATURAL GAS↗

Quantifying the reductions in mortality from air-pollution by cancelling new coal power plants

Deep decarbonization paths to the 1.5°C or 2°C temperature stabilization futures require a rapid reduction in conventional coal-fired power plants, but countries are currently building 223 GW of new coal power capacity and plan to build 377 GW more in the next decade. Coal-fired plants are also a major contributor to air pollution related health impacts. Here, we couple an integrated human-earth system model (GCAM) with an air quality model (TM5-FASST) to examine regional health co-benefits from cancelling new coal-fired plants worldwide. Additionally, we find that cancelling all new proposed projects would decrease air pollution related premature mortality between 101,388 - 213,205 deaths (2-5%) in 2030, and 213,414 - 373,054 (5-8%) in 2050, globally, but heavily concentrated in developing Asia. Furthermore, we estimate that strengthening the climate target from 2°C to 1.5°C would avoid 326,351 additional mortalities in 2030, of which 251,011 (75%) are attributable to the incremental coal plant shutdown.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗