Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipation rate”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗
DOE OSTI2025

Cosmic recombination in the presence of primordial magnetic fields

Primordial magnetic fields (PMFs) may explain observations of magnetic fields on extragalactic scales. They are most cleanly constrained by measurements of cosmic microwave background radiation (CMB) anisotropies. Their effects on cosmic recombination may even be at the heart of the resolution of the Hubble tension. We present the most detailed analysis of the effects of PMFs on cosmic recombination to date. To this end we extend the public magneto-hydrodynamic code ENZO with a new cosmic recombination routine, Monte-Carlo simulations of Lyman-α photon transport, and a Compton drag term in the baryon momentum equation. The resulting code allows us, for the first time, to realistically predict the impact of PMFs on the cosmic ionization history and the clumping of baryons during cosmic recombination. Our results identify the importance of mixing of Lyman-α photons between overdense- and underdense- regions for small PMF strength. This mixing speeds up recombination beyond the speed-up due to clumping. We also investigate the effects of pecuilar flows on the recombination rate and find it to be small for small PMF strengths. For non-helical PMFs with a Batchelor spectrum we find a surprising dependency of results on ultra-violet magnetic modes. We further show that the increase in the ionization fraction at low redshift by hydrodynamic baryon heating due to PMF dissipation is completely compensated by the faster recombination from baryon clumping. In conclusion, the present study shall serve as a theoretical foundation for a future precise comparison of recombination with PMFs to CMB data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗
DOE OSTI2023

Dependence of the kinetic energy absorption capacity of bistable mechanical metamaterials on impactor mass and velocity

Using an alternative mechanism to dissipation or scattering, bistable structures and mechanical metamaterials have shown promise for mitigating the detrimental effects of impact by reversibly locking energy into strained material. Herein, we extend prior works on impact absorption via bistable metamaterials to computationally explore the dependence of kinetic energy transmission on the velocity and mass of the impactor, with strain rates exceeding 10 2 s –1 . We observe a large dependence on both impactor parameters, ranging from significantly better to worse performance than a comparative linear material. We then correlate the variability in performance to solitary wave formation in the system and give analytical estimates of idealized energy absorption capacity under dynamic loading. In addition, we find a significant dependence on damping accompanied by a qualitative difference in solitary wave propagation within the system. The complex dynamics revealed in this study offer potential future guidance for the application of bistable metamaterials to applications including human and engineered system shock and impact protection devices.

42 ENGINEERING↗
DOE OSTI2023

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗
DOE OSTI2022

Mitigation of plasma–wall interactions with low-Z powders in DIII-D high confinement plasmas

Experiments with low-Z powder injection in DIII-D high confinement discharges demonstrated increased divertor dissipation and detachment while maintaining good core energy confinement. Lithium (Li), boron (B), and boron nitride (BN) powders were injected in H-mode plasmas (I p = 1 MA, B t = 2 T, P NB = 6 MW, < n e > = 3.6–5.0 · 10 19 m –3 ) into the upper small-angle slot divertor for 2 s intervals at constant rates of 3–204 mg s –1 . The multi-species BN powders at a rate of 54 mg s –1 showed the most substantial increase in divertor neutral compression by more than an order of magnitude and lasting detachment with minor degradation of the stored magnetic energy W mhd by 5%. Rates of 204 mg s –1 of boron nitride powder further reduce edge localized mode-fluxes on the divertor but also cause a drop in confinement performance by 24% due to the onset of an n = 2 tearing mode. The application of powders also showed a substantial improvement of wall conditions manifesting in reduced wall fueling source and intrinsic carbon and oxygen content in response to the cumulative injection of non-recycling materials. Furthermore, the results suggest that low-Z powder injection, including mixed element compounds, is a promising new core-edge compatible technique that simultaneously enables divertor detachment and improves wall conditions during high confinement operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

High Power Test Results For A Cylindrical-Shell SiC Higher-Order-Mode DAMPER

The next high-current Electron Ion Collider (EIC) is a new accelerator to be built at Brookhaven National Laboratory (BNL) in collaboration with Thomas Jefferson National Accelerator Facility (TJNAF). In the EIC Electron Storage Ring (ESR), there will be beam currents up to 2.5 A, which will excite massive Higher-Order-Modes (HOM) power in the 17 single-cell 591 MHz Superconducting Radio Frequency (SRF) cavities. To damp the HOM power in the ESR SRF cavities is a challenge. A room temperature cylindrical shell shape Silicon Carbide (SiC) Beamline HOM Absorber (BLA) was chosen as the baseline design, due to its broadband and high-power capability, and previous demonstrations at other accelerator facilities, albeit at much lower power. Because the EIC BLA HOM power dissipation is significantly greater than the previous applications, it is imperative to carry out high power testing to determine the maximum device performance levels achievable for thermal transport, RF breakdown, and mechanical stress, prior to finalizing the design. A SiC HOM absorber with state-of-the-art geometry size was prototyped to verify shrink-fit technique, test outgassing rate, and high power handling capability. This paper presents the HOM damper’s prototyping and test results.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Interaction between cloud–radiation, atmospheric dynamics and thermodynamics based on observational data from GoAmazon 2014/15 and a cloud-resolving model

Observational meteorological data from the field experiment GoAmazon 2014/15 and data from numerical simulations with the cloud-resolving model (CRM) called the System for Atmospheric Modeling (SAM) are used to study the interaction between the cloudiness–radiation as well as the atmospheric dynamics and thermodynamics variables for a site located in the central Amazon region (–3.2°S, –60.6°W) during the wet and dry periods. The main aims are to (a) analyze the temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux as well as (b) to determine the relationship between the integrated cloud fraction, radiative fluxes and large-scale variable anomalies as a function of the previous day's average. The temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux from SAM simulations showed physical consistency with the observations from GoAmazon 2014/15. Shallow and deep convection clouds show to have a meaningful impact on radiation fluxes in the Amazon region during wet and dry periods. Anomalies of large-scale variables (relative to the previous day's average) are physically associated with cloud formation, evolution and dissipation. SAM consistently simulated these results, where the cloud fraction vertical profile shows a pattern very close to the observed data (cloud type). Additionally, the integrated cloud fraction and large-scale variable anomalies, as a function of the previous day's average, have a good correlation. These results suggest that the memory of the large-scale dynamics from the previous day can be used to estimate the cloud fraction as well as the water content, which is a variable of the cloud itself. In general, the SAM satisfactorily simulated the interaction between cloud–radiation as well as dynamic and thermodynamic variables of the atmosphere during the periods of this study, being able to obtain atmospheric variables that are impossible to obtain in an observational way.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

No Escaping Helium from 55 Cnc e

We search for escaping helium from the hot super-Earth 55 Cnc e by taking high-resolution spectra of the 1083 nm line during two transits using Keck/NIRSPEC. We detect no helium absorption down to a 90% upper limit of 250 ppm in excess absorption or 0.27 mÅ in equivalent width. This corresponds to a mass-loss rate of less than ∼10{sup 9} g s{sup −1} assuming a Parker wind model with a plausible exosphere temperature of 5000–6000 K, although the precise constraint is heavily dependent on model assumptions. We consider both hydrogen- and helium-dominated atmospheric compositions and find similar bounds on the mass-loss rate in both scenarios. Our hydrodynamical models indicate that if a lightweight atmosphere exists on 55 Cnc e, our observations would have easily detected it. Together with the nondetection of Lyα absorption by Ehrenreich et al., our helium nondetection indicates that 55 Cnc e either never accreted a primordial atmosphere in the first place or lost its primordial atmosphere shortly after the dissipation of the gas disk.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Strong-Drive Limits in Josephson Circuits: From Chaos to an Unbound-Resonance Threshold

Strong microwave drives enable fast measurement and parametric control in superconducting circuits but can induce transitions out of the intended low-energy manifold. We develop a unified description of strong-drive limits in flux- and charge-driven Josephson circuits across drive frequency and dc flux bias. Using classical phase-space analysis and Floquet--Markov simulations, we identify distinct low- and high-frequency mechanisms. At low frequency, we characterize bound-state resonances and separatrix chaos and find that the flux-drive chaos threshold depends strongly on dc flux bias. At high frequency, these mechanisms are suppressed, and the dissipative steady state transfers from the central bound-state sector to outer resonances formed by above-barrier running trajectories. The resulting unbound-resonance threshold is nearly independent of drive frequency and circuit parameters over the regime studied and is controlled primarily by dc flux bias. Coherent simulations show that parametric operation persists beyond this threshold, but at a reduced rate, setting an effective upper bound on the achievable operation speed. We derive analytical criteria for both thresholds, validate them numerically, and experimentally confirm the predicted dc-bias dependence of the low-frequency threshold in a flux-driven SQUID. We also determine the timescales for transfer into the unbound-resonance regime and relaxation back to the bound-state manifold after the drive is removed. Finally, we relate the stability limits to a complementary picture based on the junction critical current and extend the framework to multitone drives and inductively shunted circuits. Together, these results identify the mechanisms limiting strong driving and suggest routes to extend the stable operating range of Josephson circuits.

You, Xinyuan [Fermilab] (ORCID:0000000291789419)↗
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Work of Mechanical Degradation in Elongating Polymer Melts

Molecular dynamics simulations are used to study the mechanical degradation of well-entangled polymer melts during uniaxial extensional flow. Simulations measure the transient rise in extensional stresses and relate them to the molecular alignment and scission of chain backbones. Intermolecular entanglements couple chain scission in space and time, making degradation sensitive to deformation history and strain rate in ways not displayed by dilute polymer solutions. The rate of chain scission is nonmonotonic and peaks at strains corresponding to the maximum extensibility of entanglement segments but prior to the full extension of chain backbones. We measure a specific work per scission event w* and decompose it into separate contributions associated with chain alignment, chemical bond breaking, and scission-induced plasticity. We find chain scission in melts requires activating plastic dissipation that is multiple orders of magnitude larger than the chemical work required to break a covalent backbone bond. Our findings underscore the critical need to consider bulk polymer mechanics and rheology in designing efficient mechanical degradation and mechanochemical processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-power test results for a cylindrical-shell silicon carbide higher-order-mode damper

The next high-current Electron-Ion Collider (EIC) is a new accelerator to be built at Brookhaven National Laboratory in collaboration with Thomas Jefferson National Accelerator Facility. In the EIC Electron Storage Ring (ESR), there will be beam currents of up to 2.5 A, which will excite massive higher-order-mode (HOM) power in the 17 single-cell 591 MHz superconducting radio-frequency (SRF) cavities. Damping the HOM power in the ESR SRF cavities is a challenge. A room temperature cylindrical shell shape silicon carbide (SiC) beamline HOM absorber (BLA) was chosen as the baseline design, due to its broadband and high-power capability, and previous demonstrations at other accelerator facilities, albeit at much lower power. Because the EIC BLA HOM power dissipation is significantly greater than the previous applications, it is imperative to carry out high-power testing to determine the maximum device performance levels achievable for thermal transport, rf breakdown, and mechanical stress, prior to finalizing the design. A SiC HOM absorber with a state-of-the-art geometry size was prototyped to verify the shrink-fit technique, test outgassing rate, and high-power handling capability. This paper presents the HOM damper’s prototyping and test results. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulating Atomic Dark Matter in Milky Way Analogs

Dark sector theories naturally lead to multicomponent scenarios for dark matter where a subcomponent can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogs where the majority of dark matter is collisionless cold dark matter (CDM) but a subcomponent (6%) is strongly dissipative minimal atomic dark matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly standard model–like masses and couplings can cool to form a rotating “dark disk” with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark “clumps” (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ∼kiloparsec of the galaxy, affecting the galaxy’s star formation history and present-day baryonic and CDM distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Microwave assisted heating of plastic waste: Effect of plastic/susceptor (SiC) contacting patterns

We discuss plastic accumulation as a growing sustainability challenge. Most plastics are thermally inert at room temperature, requiring high temperatures for depolymerization. In the quest for developing sustainable plastic upcycling, we investigate silicon carbide (SiC) as microwave absorbing susceptors for reaching depolymerization temperatures. We utilize three SiC topologies – particles, hot-pressed plastic films containing SiC particles, and monoliths. SiC particles of >0.5 mm in diameter can provide heating rates approaching ~200 °C/min at modest powers. Still, physically mixed plastic-susceptor beds suffer from density differences, causing non-uniform contact during plastic softening and melting. Hot-pressed plastic films containing SiC particles allow uniform plastic-susceptor contact but suffer from slow heating as only a single layer of susceptors exists between polymer layers. SiC monoliths mitigate these challenges by preferentially dissipating the microwaves and rapidly heating the plastics independently of polymeric properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microbubble dynamics and heat transfer in boiling droplets

Dissipating large heat fluxes from a surface is critically important in numerous industrial and natural applications. Boiling based spray cooling and surface texturing are two of the most promising methods being investigated to address this problem. Although our understanding on these topics has significantly improved over past decades, critical gaps remain in the knowledgebase stymieing the realization of their full potential. As an example, while bubble growth in pool boiling have been investigated in detail, comparatively little is known about how the bubbles evolve inside boiling drops. In the present work, we have investigated for the first time, the microbubble dynamics inside water droplets boiling on superhydrophilic textured substrates using high-speed X-ray phase contrast imaging (XRPCI). Our observations show that the transient bubble density variation follows similar characteristics irrespective of the texture spacing at a given surface temperature. For an example microstructure, we found that the number of discrete bubbles on the surface decreases as temperature is increased although their growth rate increases. Here, we observe that bubble growth is highly non-uniform during the lifetime of a drop on the surface. Initially, bubbles grow under diffusion-limited regime, but at later times they grow as ~t 1.45 due to combined effects of coalescence and evaporation. In some conditions, we found that bubbles shrink dramatically af- ter the initial growth spurt presumably due to severe quenching of the surface, and migration of bubbles on the surface. Using the bubble sizes, for the first time we analyzed the heat flux removed by a single bubble and also by all the bubbles at a given time. We find that the highest dissipation through latent- heat component (~600 W/cm 2 ) occurs just in the beginning and thereafter it decreases. We expect that our findings and the analysis would guide further work on the topic and will aid in the overarching goal of engineering surfaces that are more efficient in boiling heat transfer.

42 ENGINEERING↗