Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scaling laws of complex systems”

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.

50 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effective Permeability of a Nuclear Fuel Assembly

This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Photocurrent Nanoscopy of Quantum Hall Bulk

Understanding nanoscale electronic and thermal transport of two-dimensional (2D) electron systems in the quantum Hall regime, particularly in the bulk insulating state, poses considerable challenges. One of the primary difficulties arises from the presence of chiral edge channels, whose transport behavior obscures the investigation of the insulating bulk. Using near-field optical and photocurrent nanoscopy, we probe real-space variations of the optical and thermal dynamics of graphene in the quantum Hall regime without relying on complex sample or electrode geometries. Near the charge-neutrality point, we detect strong optical and photothermal signals from resonant inter-Landau-level (LL) magnetoexciton excitations between the zeroth and ± 1 st LLs, which gradually weaken with increasing doping due to Pauli blocking. Interestingly, at higher doping levels and full-integer LL fillings, photothermal signals reappear across the entire sample over an approximately 10 − μ m scale, indicating unexpectedly long cooling lengths and nonlocal photothermal heating through the insulating bulk. This observation suggests that thermal conductivity persists for the localized states even as electronic transport is suppressed—a clear violation of the Wiedemann-Franz law. Our experiments provide novel insights into nanoscale thermal and electronic transport in incompressible 2D gases, highlighting the roles of magnetoexcitons and chiral edge states in the thermo-optoelectric dynamics of the Dirac quantum Hall state. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗
DOE OSTI2023

Breaking Rayleigh’s Law with Spatially Correlated Disorder to Control Phonon Transport

Controlling thermal transport in insulators and semiconductors is crucial for many technological fields such as thermoelectrics and thermal insulation, for which a low thermal conductivity (κ) is desirable. A major obstacle for realizing low κ materials is Rayleigh’s law, which implies that acoustic phonons, which carry most of the heat, are insensitive to scattering by point defects at low energy. We demonstrate, with large scale simulations on tens of millions of atoms, that isotropic long-range spatial correlations in the defect distribution can dramatically reduce phonon lifetimes of important low-frequency heat-carrying modes, leading to a large reduction of κ—potentially an order of magnitude at room temperature. We propose a general and quantitative framework for controlling thermal transport in complex functional materials through structural spatial correlations, and we establish the optimal functional form of spatial correlations that minimize κ. Finally, we end by briefly discussing experimental realizations of various correlated structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simulating dirty bosons on a quantum computer

Abstract Quantum computers hold the potential to unlock new discoveries in complex quantum systems by enabling the simulation of physical systems that have heretofore been impossible to implement on classical computers due to intractability. A system of particular interest is that of dirty bosons, whose physics highlights the intriguing interplay of disorder and interactions in quantum systems, playing a central role in describing, for instance, ultracold gases in a random potential, doped quantum magnets, and amorphous superconductors. Here, we demonstrate how quantum computers can be used to elucidate the physics of dirty bosons in one and two dimensions. Specifically, we explore the disorder-induced delocalized-to-localized transition using adiabatic state preparation. In one dimension, the quantum circuits can be compressed to small enough depths for execution on currently available quantum computers. In two dimensions, the compression scheme is no longer applicable, thereby requiring the use of large-scale classical state vector simulations to emulate quantum computer performance. In addition, simulating interacting bosons via emulation of a noisy quantum computer allowed us to study the effect of quantum hardware noise on the physical properties of the simulated system. Our results suggest that scaling laws control how noise modifies observables versus its strength, the circuit depth, and the number of qubits. Moreover, we observe that noise impacts the delocalized and localized phases differently. A better understanding of how noise alters the observed properties of the simulated system is essential for leveraging near-term quantum devices for simulation of dirty bosons, and indeed for condensed matter systems in general.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Discovery of Powerful Multivelocity Ultrafast Outflows in the Starburst Merger Galaxy IRAS 05189–2524 with XRISM

We observed the X-ray-bright ultraluminous infrared galaxy IRAS 05189−2524 with XRISM during its performance verification phase. The unprecedented energy resolution of the onboard X-ray microcalorimeter revealed complex spectral features at ∼7–9 keV, which can be interpreted as blueshifted Fe XXV/XXVI absorption lines with various velocity dispersions, originating from ultrafast outflow (UFO) components with multiple bulk velocities of ∼0.076c, ∼0.101c, and ∼0.143c. In addition, a broad Fe–K emission line was detected around ∼7 keV, forming a P Cygni profile together with the absorption lines. The onboard X-ray CCD camera revealed a 0.4–12 keV broadband spectrum characterized by a neutrally absorbed power-law continuum with a photon index of ∼2.3 and intrinsic flare-like variability on timescales of ∼10 ks, both of which are likely associated with near-Eddington accretion. We also found potential variability of the UFO parameters on a timescale of ∼140 ks. Using these properties, we propose new constraints on the outflow structure and suggest the presence of multiple outflowing regions on scales of about tens to 100 Schwarzschild radii, located within roughly 2000 Schwarzschild radii. Since both the estimated momentum and energy outflow rates of the UFOs exceed those of galactic molecular outflows, our results indicate that powerful, multivelocity UFOs are already well developed during a short-lived evolutionary phase following a major galaxy merger, characterized by intense starburst activity and likely preceding the quasar phase. This system is expected to evolve into a quasar, sustaining strong UFO activity and suppressing star formation in the host galaxy.

Astronomy and AstroPhysics↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗
DOE OSTI2022

Scaling of the strange-metal scattering in unconventional superconductors

Marked evolution of properties with minute changes in the doping level is a hallmark of the complex chemistry that governs copper oxide superconductivity as manifested in the celebrated superconducting domes and quantum criticality taking place at precise compositions. The strange-metal state, in which the resistivity varies linearly with temperature, has emerged as a central feature in the normal state of copper oxide superconductors. The ubiquity of this behaviour signals an intimate link between the scattering mechanism and superconductivity. However, a clear quantitative picture of the correlation has been lacking. In this study we report the observation of precise quantitative scaling laws among the superconducting transition temperature (T c ), the linear-in-T scattering coefficient (A 1 ) and the doping level (x) in electron-doped copper oxide La 2-x Ce x CuO 4 (LCCO). High-resolution characterization of epitaxial composition-spread films, which encompass the entire overdoped range of LCCO, has enabled us to systematically map its structural and transport properties with unprecedented accuracy and with increments of Δx = 0.0015. Additionally, we have uncovered the relations T c ~ (x c - x) 0.5 ~ (A 1 $\square$ ) 0.5 , where x c is the critical doping in which superconductivity disappears and A 1 $\square$ is the coefficient of the linear resistivity per CuO 2 plane. The striking similarity of the T c versus A 1 $\square$ relation among copper oxides, iron-based and organic superconductors may be an indication of a common mechanism of the strange-metal behaviour and unconventional superconductivity in these systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗
DOE OSTI2020

Adapting Existing Energy Planning, Simulation, and Operational Models for Resilience Analysis

Resilience is a topic receiving much attention in relation to energy systems, with particular attention being paid to the supply of electricity. As a result, research communities have proposed a plethora of candidate indicators and metrics for resilience, most of which remain immature at different scales and segments within the energy system. Given the complexity of resilience analyses and mitigation strategies, there is limited value in attempting to identify a single resilience metric, as no one metric can quantify resilience or its associated value for all stakeholders. Instead, a necessary focus of the research community should lie in implementing, testing, and validating resilience metrics and analysis approaches in energy sector models, which will be invaluable for informing resilience planning and investment decisions. Recognizing that implementing, testing, and validating resilience metrics are challenges that need to be addressed, the National Renewable Energy Laboratory (NREL) dedicated staff and time to researching how to effectively integrate resilience considerations into energy sector models and tools, as part of the Laboratory Directed Research and Development (LDRD) program. The overarching goal of the effort was to evaluate the data needs, methodologies, and outcomes - including consequences and/or changes in investment or operational decisions due to avoided consequences - based on resilience analysis in a range of existing tools. In particular, we selected five models originally built at NREL to explore non-resilience energy research questions to implement and exercise resilience metrics and analysis approaches.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Employing Long-Range Inductive Effects to Modulate Metal-to-Ligand Charge Transfer Photoluminescence in Homoleptic Cu(I) Complexes

Four Cu(I) bis(phenanthroline) photosensitizers formulated from a new ligand structural motif (Cu1-Cu4) coded according to their 2,9-substituents, were synthesized, structurally characterized, and fully evaluated using steady-state and time-resolved absorption and photoluminescence (PL) measurements as well as electrochemistry. The 2,9-disubstituted-3,4,7,8-tetramethyl-1,10-phenanthroline ligands feature the following six-membered ring systems prepared through photochemical synthesis: 4,4-dimethylcyclohexyl (1), tetrahydro-2H-pyran-4-yl (2), tetrahydro-2H-thiopyran-4-yl (3), and 4,4-difluorocyclohexyl (4). Universally, these Cu(I) metal-to-ligand charge transfer (MLCT) chromophores display excited state lifetimes on the microseconds time scale at room temperature, including the three longest-lived homoleptic cuprous phenanthroline excited states measured to date in deaerated CH 2 Cl 2 , τ = 2.5 to 4.3 μs. This series of molecules also feature high PL quantum efficiencies (ΦPL = 5.3 to 12% in CH 2 Cl 2 ). Temperature-dependent PL lifetime experiments confirmed that all these molecules exhibit reverse intersystem crossing and display thermally activated delayed PL from a 1 MLCT excited state lying slightly above the 3 MLCT state, 1050 – 1490 cm -1 . Ultrafast and conventional transient absorption measurements confirmed that the PL originates from the MLCT excited state, which remains sterically arrested, preventing an excessive flattening distortion even when dissolved in Lewis basic CH 3 CN. Combined PL and electrochemical data provided evidence that Cu1-Cu4 are highly potent photoreductants (E ox * = -1.73 to -1.62 V vs. Fc +/o in CH 3 CN), whose potentials are altered solely based on which heteroatoms or substituents are resident on the 2,9-appended ring derivatives. It is proposed that long-range electronic inductive effects are responsible for the systematic modulation observed in the PL spectra, excited state lifetimes, as well as the ground state absorption spectra and redox potentials. Cu1-Cu4 quantitatively follow the energy gap law, correlating well with structurally related cuprous phenanthrolines and are also shown to triplet photosensitize the excited states of 9,10-diphenylanthracene with bimolecular rate constants ranging from 1.61 to 2.82 × 10 8 M –1 s –1 . Furthermore, the ability to tailor both photophysical and electrochemical properties using long-range inductive effects imposed by the 2,9-ring platforms advocates new directions for future MLCT chromophore discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Entropy production in model colloidal suspensions under shear via the fluctuation theorem

Dissipative systems often exhibit novel and unexpected properties. This is, for instance, the case of simple liquids, which, when subjected to shear and after reaching a steady state, can exhibit a negative entropy production over finite length scales and timescales. This result, among others, is captured by nonequilibrium relations known as fluctuation theorems. Using nonequilibrium molecular dynamics simulations, we examine how, by fine-tuning the properties of the components of a complex fluid, we can steer the nonequilibrium response of the fluid. More specifically, we show how we control the nonequilibrium probability distribution for the shear stress and, in turn, how often states with a negative entropy production can occur. To achieve this, we start by characterizing how the size for the liquid matrix impacts the probability of observing negative entropy states, as well as the timescale over which these can be observed. We then measure how the addition of larger particles to this liquid matrix, i.e., simulating a model colloidal suspension, results in an increase in the occurrence of such states. As a result, this suggests how modifications in the composition of the mixture and in the properties of its components lead to an increase in the probability of observing states of negative entropy production and, thus, for the system to run in reverse.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗