PREDICTIONS OF THE STRUCTURE OF RADIATION-RESISTED SHOCK WAVES
Application of the continuum, inviscid-flow equations of gas dynamics to predict the effect of thermal radiation on the internal structure of a shock wave
SEARCH · Engineering Papers
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.
Application of the continuum, inviscid-flow equations of gas dynamics to predict the effect of thermal radiation on the internal structure of a shock wave
We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.
Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.
The yielding of soft materials is critical to many natural and industrial processes, yet experimental insights into microscopic aspects of yielding are limited. This study combines angle X-ray scattering, X-ray photon correlation spectroscopy, and in situ rheology (Rheo-SAXS-XPCS) with fast lubrication dynamics simulations to examine how interparticle interactions influence yielding in charged colloidal suspensions. By tuning attraction through salt addition, we compare repulsive and attractive systems under deformation. Repulsive suspensions yield uniformly with Andrade-like creep and minimal structural change. In contrast, attractive suspensions show complex behaviors, including shear banding, delayed yielding, and resolidification, governed by transient dynamics at shear band interfaces. These results directly link microscopic particle dynamics to macroscopic flow and demonstrate how interaction potentials control rheological behavior. This work offers a framework for designing soft materials with tailored properties for applications in coatings, food processing, drug delivery, and other technologies requiring precise mechanical control.
The major objectives of this study are to test the validity of quasi-geostrophic (QG) dynamics, compared to primitive equation (PE) dynamics for modeling the effect of cyclone waves on the larger scale flow, and to study the formation of frontal cyclones and the dynamics of occluded frontogenesis. In order to extend the realism of the initial conditions over those used in earlier work of this nature and in order to include strong surface frontal zones in the initial conditions, a horizontal resolution on the order of 100 km between gridpoints is needed. In order to allow for integrations of up to two weeks simulated time and to keep the logistics and cost of the runs feasible, the horizontal domain of the models must be reduced. Based on work concerning polar lows, the tropospheric static stability was reduced and was also allowed to vary in different tropospheric layers. Previous work used a 'standard atmosphere' temperature lapse rate. Stratospheric stability remains unchanged. The use of one of these 'reduced stability' cases, with no surface frontal zone present resulted in the discovery of normal mode perturbations.
Von Karman swirling flow is generated by the viscous pumping action of a solid disk spinning in a quiescent fluid media. When this spinning disk is ignited in an oxidizing environment, a flat diffusion flame is established adjacent to the disk, embedded in the boundary layer (see the preceding illustration). For this geometry, the conservation equations reduce to a system of ordinary differential equations, enabling researchers to carry out detailed theoretical models to study the effects of varying strain on the dynamics of diffusion flames. Experimentally, the spinning disk burner provides an ideal configuration to precisely control the strain rates over a wide range. Our original motivation at the NASA Glenn Research Center to study these flames arose from a need to understand the flammability characteristics of solid fuels in microgravity where slow, subbuoyant flows can exist, producing very small strain rates. In a recent work (ref. 1), we showed that the flammability boundaries are wider and the minimum oxygen index (below which flames cannot be sustained) is lower for the von Karman flow configuration in comparison to a stagnation-point flow. Adding a small forced convection to the swirling flow pushes the flame into regions of higher strain and, thereby, decreases the range of flammable strain rates. Experiments using downward facing, polymethylmethacrylate (PMMA) disks spinning in air revealed that, close to the extinction boundaries, the flat diffusion flame breaks up into rotating spiral flames (refs. 2 and 3). Remarkably, the dynamics of these spiral flame edges exhibit a number of similarities to spirals observed in biological systems, such as the electric pulses in cardiac muscles and the aggregation of slime-mold amoeba. The tail of the spiral rotates rigidly while the tip executes a compound, meandering motion sometimes observed in Belousov-Zhabotinskii reactions.
Computer flow simulation aided by dynamical systems analysis is used to investigate the kinematics of time-periodic vortex shedding past a two-dimensional circular cylinder in the context of the following general questions: (1) Is a dynamical systems viewpoint useful in the understanding of this and similar problems involving time-periodic shedding behind bluff bodies; and (2) Is it indeed possible, by adopting such a point of view, to complement previous analyses or to understand kinematical aspects of the vortex shedding process that somehow remained hidden in previous approaches. We argue that the answers to these questions are positive. Results are described.
Research in both artificial hearts and fluid dynamics benefits from computational studies. Algorithm that implements Navier-Stokes equations of flow extended to simulate flow of viscous, incompressible blood through articifial heart. Ability to compute details of such flow important for two reasons: internal flows with moving boundaries of academic interest in their own right, and many of deficiencies of artificial hearts attributable to dynamics of flow.
We report the first congruent integration of HPC, AI, and multiscale modeling (MSM) for solving a mainstream biomechanical problem of thrombogenesis involving 6 million particles at record molecular-scale resolutions in space and at simulation rates of milliseconds per day. The two supercomputers, the IBM Summit-like AiMOS and our University’s SeaWulf, are used for scalability analysis of, and production runs with, the LAMMPS with our customization and AI augmentation and they attained optimal simulation speeds of 3,077 µs/day and 266 µs/day respectively. The long-time and large scales simulations enable the first study of the integrated platelet flowing, flipping, aggregating dynamics in one dynamically-coupled production run. The platelets’ angular and translational speeds, membrane particles’ speeds, and the membrane stress distributions are presented for the analysis of platelets’ aggregations.
Community-scale building retrofits are not merely scaled-up versions of single-building retrofits. They involve complex challenges, such as reconciling individual interests with collective goals and managing the dynamic interplay between buildings through mechanisms like power grids and social connections. Internet of Things (IoT) connectivity holds the potential to leverage these interplays to balance individual and collective interests effectively in smart communities. One critical aspect of this interplay is information diffusion, which shapes how retrofit decisions spread among neighbors, influencing individual choices and ultimately impacting community-level retrofit outcomes. In other words, IoT-based smart devices automatically push tailored retrofit notifications to homeowners, which completely changes the format of information diffusion in the future. To investigate this influence by such information diffusion, the study used CityBES to simulate energy performance for different retrofits and applied an information diffusion model to analyze how decisions spread in a networked community of 192 buildings. The diffusion process was modeled on a weighted, directed network, capturing the dynamics of information flow and decision-making across 16 scenarios. Individual retrofit benefits were evaluated through payback years, while community-level retrofit outcomes were assessed using greenhouse gas (GHG) emission reductions. The results demonstrate that easier information diffusion among neighbors encourages households to prioritize retrofit measures that align with the majority’s optimal choices, even at the expense of individual financial benefits. In this case, such collective prioritization enhanced community-level retrofit performance, increasing GHG emission reductions by up to 29.4 %. However, this improvement came with trade-offs, as the average payback period for households extended by approximately 1.74 years. These findings highlight the potential of IoT-based information diffusion in future smart communities to coordinate individual interests with collective goals, ultimately accelerating community-level building retrofits.
The United States generates 42 Mt of plastic waste each year and is one of the biggest contributors to ocean plastic waste. Consequently, plastic has become synonymous with the linear economy, and many scholars are studying and proposing circular economy solutions to mitigate plastic pollution. Recycling has received much attention from both social sciences and engineering as a circular economy strategy, but no study has yet quantified how behavioral interventions could asymmetrically affect different populations. Here, this study combines agent-based modeling, material flow analysis, system dynamics, and life cycle assessment to assess the effect of four behavioral interventions on the collection rates of polyethylene terephthalate bottle waste, displaced virgin plastic manufacturing, and avoided greenhouse gas (GHG) emissions. Results show that, while behavioral interventions would require about 300–900 GJ of additional energy at end-of-life due to improved collection rates, they would avoid about 500–700 thousand metric tons of GHG emissions. Results also illustrate the importance of habits in disposal behaviors and show that different forms of interventions can be better adapted to particular social contexts than others. While the circular economy and its application to plastic waste should certainly not be restricted to recycling, this study demonstrates that improved collection rates and recycling technologies can contribute to reducing the amount of plastic waste polluting our oceans.
The failure of bulk metallic glasses (BMGs) during plastic deformation at room temperature is abrupt and instantaneous, while the analysis of precursor information based on avalanche events helps predict catastrophic failure. An acoustic emission (AE) signal can provide accurate precursor information for material failure, due to its sensitive and high fast calculation ability. In the current study, AE monitoring tests are carried out during uniaxial compression tests of BMGs at different strain rates. The AE experimental failure threshold, E max , is proposed on the basis of AE cumulative energy, which reflects the intensity of damage evolution at different loading conditions. Compared with the critical shear band velocity (CSBV) associated with stick-slip dynamics of serrated flow, E max is a more sensitive failure parameter since it is connected with the local microscopic changes that occur during the material response process. Here, the E max is obtained prior to reaching the CSBV since the calculation of these two avalanches analysis focuses on the different stages of shear band growth. In particular, AE events are related to the “dry” friction process in the first stage, however, the CSBV is responsible for the “viscous” glide in the second stage. Therefore, Emax is not affected by the complex interactions between the shear bands during the stick-slip process. The maximum avalanche of serrated flow, S max , is proposed as the experimental failure threshold, which depends on the applied strain rate as S max ~ $\dot{ε}$ –λ . According to the relationship of E max and S max , the theoretical failure threshold, E max , follows a criterion E max = 2545$\dot{ε}$ –λ - 4468, where λ is equivalent to 0.15 for this work. Finally, combining the different calculations and AE measurements, this model gives new insights to predict the deformation failure behavior of Zr-based BMGs.
Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.
Athermal semiflexible sheets dispersed in a fluid are simulated under simple shear flow, and the dynamical behavior of the sheets is found to depend strongly on initial orientation and the ratio of bending stress to viscous stress.
Collective properties of complex systems composed of many interacting components such as neurons in our brain can be modeled by artificial networks based on disordered systems. We show that a disordered neural network of superconducting loops with Josephson junctions can exhibit computational properties like categorization and associative memory in the time evolution of its state in response to information from external excitations. Superconducting loops can trap multiples of fluxons in many discrete memory configurations defined by the local free energy minima in the configuration space of all possible states. A memory state can be updated by exciting the Josephson junctions to fire or allow the movement of fluxons through the network as the current through them surpasses their critical current thresholds. Simulations performed with a lumped element circuit model of a 4-loop network show that information written through excitations is translated into stable states of trapped flux and their time evolution. Experimental implementation on a high-Tc superconductor YBCO-based 4-loop network shows dynamically stable flux flow in each pathway characterized by the correlations between junction firing statistics. Neural network behavior is observed as energy barriers separating state categories in simulations in response to multiple excitations, and experimentally as junction responses characterizing different flux flow patterns in the network. The state categories that produce these patterns have different temporal stabilities relative to each other and the excitations. This provides strong evidence for time-dependent (short-to-long-term) memories, that are dependent on the geometrical and junction parameters of the loops, as described with a network model.
Attosecond pulses in the soft x-ray (SXR) to water-window energy region offer the tools for creating and studying target-specific localized inner-shell electrons or holes in materials, enabling monitoring or controlling charge and energy flows in a dynamic system on attosecond timescales. Recently, a number of laboratories have reported generation of continuum harmonics in the hundred-electron-volt to kilovolt region with few-cycle long-wavelength mid-infrared lasers. These harmonics have the bandwidth to support pulses with durations of a few to a few-ten attoseconds. But harmonics generated in a gas medium have attochirps that cannot be fully compensated by materials over a broad spectral range; thus, realistically what are the typical shortest attosecond pulses that one can generate? To answer this question, it is essential that the temporal attosecond pulses be accurately characterized. By reanalyzing the soft x-ray attosecond metrology reported in three recent experiments using a newly developed broadband phase retrieval algorithm, we demonstrate here that the generated attosecond pulses in the first two papers have durations of about 60 as, longer than what they have reported. Similarly, the duration from the third experiment is retrieved to be about half of the 322-as upper limit cited in that work. We also introduce the autocorrelation (AC) of the streaking spectrogram. By comparing the ACs from the experiments and from the retrieved SXR pulses, the accuracy of the retrieved results can be directly visualized. Our retrieval method is fast and accurate, and it shall provide a powerful tool for the metrology of the emerging few-ten-attosecond pulses.
In this work, we propose a new formula that extracts the quantum Hall conductance from a single (2+1) D gapped wave function. The formula applies to general many-body systems that conserve particle number, and is based on the concept of modular flow, i.e., unitary dynamics generated from the entanglement structure of the wave function. The formula is shown to satisfy all formal properties of the Hall conductance: it is odd under time reversal and reflection, even under charge conjugation, universal and topologically rigid in the thermodynamic limit. Further evidence for relating the formula to the Hall conductance is obtained from conformal field theory arguments. Finally, we numerically check the formula by applying it to a noninteracting Chern band where excellent agreement is obtained.
This paper presents a stability analysis of parallel-connected bi-directional series-series resonant network wireless power transfer (WPT), optimized for Electric Vehicle (EV) charging and vehicle-to-grid (V2G) applications. The study addresses critical stability challenges in systems integrated with diverse distributed energy resources (DERs), including photovoltaics, fuel cells, wind turbines, energy storage systems, and the AC grid. The stability of such integrated DC grid systems is paramount for ensuring reliable operation, particularly under varying power flow conditions and dynamic interactions between parallel WPT systems. The analysis included system impedance characterization, state-space modeling, and open and closed-loop stability evaluations. The results demonstrated that the integration of a robust control architecture effectively mitigates instability risks and supports scalable, efficient operation. This work underscores the converter's adaptability and its potential for large-scale deployment in wireless EV charging infrastructures and integrated DC grid systems.