Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Particle Acceleration, Magnetic Field Generation and Emission from Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), supernova remnants, and Galactic microquasar systems usually have power-law emission spectra. Fermi acceleration is the mechanism usually assumed for the acceleration of particles in astrophysical environments. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that particle acceleration occurs within the downstream jet, rather than by the scattering of particles back and forth across the shock as in Fermi acceleration. Shock acceleration' is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different spectral properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. We will review recent PIC simulations of relativistic jets and try to make a connection with observations.

Nishikawa, K.-I.↗
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↗
NASA NTRS2003

Crystallization Physics in Biomacromolecular Systems

The crystals are built of molecules of protein, nucleic acid and their complexes, like viruses, approx. 5x10(exp 3)+ 3x10(exp 6) Da in weight and 2 + 20 nm in effective diameter. This size strongly exceeds action range of molecular forces and makes a big difference with inorganic crystals. Intermolecular contacts form patches on the biomacromolecular surface. Each patch may occupy only a small percent of the whole surface and vary from polymorph to polymorph of the same protein. Thus, under different conditions (pH, solution chemistry, temperature, any area on the macromolecular surface may form a contact. The crystal Young moduli, E approx. equals 0.1 + 0.5 GPa are more than 10 times lower than that of inorganics and the biomolecules themselves. Water within biocrystals (30-70%) is unable to flow unless typical deformation time is longer than approx. 10(exp -5)s. This explains the discrepancy between light scattering and static measurements of E. Nucleation and Growth requires typically concentrations exceeding the equilibrium ones up to 100 times - because of the new size scale results in 10 - 10(exp 3) times lower kinetic coefficients than that needed for inorganic solution growth. All phenomena observed in the latter occur with protein crystallization and are even better studied by AFM. Crystals are typically facetted. Among unexpected findings of general significance are - net molecular exchange flux at kinks is much lower than that expected from supersaturation, steps with low (< approx. 10(exp -2)) kink density at steps follow Gibbs-Thomson law only at very low supersaturations, step segment growth rate may be independent of step energy. Crystal perfection is a must of biocrystallization to achieve the major goal to find 3-D atomic structure of biomacromolecules by x-ray diffraction. Poor diffraction resolution (> 3Angstrom) makes crystallization a bottleneck for structural biology. All defects typical of small molecule crystals are found in biocrystals, but the defects responsible for poor resolution are not identified. Conformational changes are one of them. Biocrystallization in microgravity reportedly results in 20% cases of better crystals. The mechanism of how lack of convection can do this is still not clear. Lower supersaturation, self-purification &om preferentially trapped homologous impurities and step bunching are viable hypotheses.

Chernov, A. A.↗
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↗
NASA NTRS1990

Rapid near-optimal trajectory generation and guidance law development for single-stage-to-orbit airbreathing vehicles

General problems associated with on-board trajectory optimization, propulsion system cycle selection, and with the synthesis of guidance laws were addressed for an ascent to low-earth-orbit of an air-breathing single-stage-to-orbit vehicle. The NASA Generic Hypersonic Aerodynamic Model Example and the Langley Accelerator aerodynamic sets were acquired and implemented. Work related to the development of purely analytic aerodynamic models was also performed at a low level. A generic model of a multi-mode propulsion system was developed that includes turbojet, ramjet, scramjet, and rocket engine cycles. Provisions were made in the dynamic model for a component of thrust normal to the flight path. Computational results, which characterize the nonlinear sensitivity of scramjet performance to changes in vehicle angle of attack, were obtained and incorporated into the engine model. Additional trajectory constraints were introduced: maximum dynamic pressure; maximum aerodynamic heating rate per unit area; angle of attack and lift limits; and limits on acceleration both along and normal to the flight path. The remainder of the effort focused on required modifications to a previously derived algorithm when the model complexity cited above was added. In particular, analytic switching conditions were derived which, under appropriate assumptions, govern optimal transition from one propulsion mode to another for two cases: the case in which engine cycle operations can overlap, and the case in which engine cycle operations are mutually exclusive. The resulting guidance algorithm was implemented in software and exercised extensively. It was found that the approximations associated with the assumed time scale separation employed in this work are reasonable except over the Mach range from roughly 5 to 8. This phenomenon is due to the very large thrust capability of scramjets in this Mach regime when sized to meet the requirement for ascent to orbit. By accounting for flight path angle and flight path angle rate in construction of the flight path over this Mach range, the resulting algorithm provides the means for rapid near-optimal trajectory generation and propulsion cycle selection over the entire Mach range from take-off to orbit.

Calise, A. J.↗
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↗
NASA NTRS1999

The Isothermal Dendritic Growth Experiment (IDGE): USMP-4 One-Year-Report

Dendrites describe the tree-like crystal morphology commonly assumed in many material systems--particularly in metals and alloys that freeze from supercooled or supersaturated melts. There remains a high level of engineering interest in dendritic solidification because of the role of dendrites in the determination of cast alloy microstructures. Microstructure plays a key role in determining the physical properties of cast or welded products. In addition, dendritic solidification provides an example of non-equilibrium physics and is one of the simplest non-trivial examples of dynamic pattern formation, where an amorphous melt, under simple starting conditions, evolves into a complex ramified microstructure. Although it is well-known that dendritic growth is controlled by the transport of latent heat from the moving solid-melt interface as the dendrite advances into a supercooled melt, an accurate, and predictive model has not been developed. Current theories consider: 1) the transfer of heat or solute from the solid-liquid interface into the melt, and 2) the interfacial crystal growth and growth selection physics for the interface. However, the effects of gravity-induced convection on the transfer of heat from the interface prevent either element from being adequately tested solely under terrestrial conditions. The Isothermal Dendritic Growth Experiment (IDGE) constituted a series of three NASA-supported microgravity experiments, all of which flew aboard the space shuttle, Columbia. This experimental space flight series was designed and operated to grow and record dendrite solidification in the absence of gravity-induced convective heat transfer, and thereby produce a wealth of benchmark-quality data for testing solidification scaling laws. The data collection from the on-orbit phase of the IDGE flight series is now complete. We are currently completing analyses and moving towards final data archiving.

Glicksman, M. E.↗
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↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗
NASA NTRS2012

Development of a Robust and Efficient Parallel Solver for Unsteady Turbomachinery Flows

The traditional design and analysis practice for advanced propulsion systems relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. In the design of advanced propulsion systems, CFD plays a major role in defining the required performance over the entire flight regime, as well as in testing the sensitivity of the design to the different modes of operation. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment. The turbomachinery simulation capability presented here is being developed in a computational tool called Loci-STREAM [1]. It integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci [2] which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective is to be able to routinely simulate problems involving complex geometries requiring large unstructured grids and complex multidisciplinary physics. An immediate application of interest is simulation of unsteady flows in rocket turbopumps, particularly in cryogenic liquid rocket engines. The key components of the overall methodology presented in this paper are the following: (a) high fidelity unsteady simulation capability based on Detached Eddy Simulation (DES) in conjunction with second-order temporal discretization, (b) compliance with Geometric Conservation Law (GCL) in order to maintain conservative property on moving meshes for second-order time-stepping scheme, (c) a novel cloud-of-points interpolation method (based on a fast parallel kd-tree search algorithm) for interfaces between turbomachinery components in relative motion which is demonstrated to be highly scalable, and (d) demonstrated accuracy and parallel scalability on large grids (approx 250 million cells) in full turbomachinery geometries.

West, Jeff↗
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↗