Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative dynamics”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Study of solid molecular deuterium D 2 growth under gas pressure

The injection of high-speed cryogenic pellets made of frozen hydrogen-isotopes, represents to date the most effective method to fuel magnetically confined thermonuclear fusion plasmas. Additionally, the injection of very large pellets composed of cryogenic solid of some suitable impurity (typically a noble-gas such as H 2 , Ne, or H 2 /Ne, D 2 /Ne mixtures), shattered in relatively small fragments just before entering the plasma, seems to be the most promising method to reduce the damage risks for the plasma-facing components in case of a plasma disruption. This technology, known as "Shattered Pellet Injection" (SPI), allows to spread out the plasma energy and mitigate possible damage to the in-vessel components, as well as to densify the plasma to suppress the formation of runaway electrons, and/or dissipate their energy. Several techniques to produce and launch cryogenic pellets have been investigated in the past decades. "Pipe gun" injectors are reliable and relatively simple devices are still commonly used today. They make use of single- or two-stage pneumatic light-gas guns to accelerate the pellet at high speeds. In these injectors, the cryogenic pellets are formed “ in situ ” (i.e., inside the launching barrel), by de-sublimating them directly from the gas phase, i.e., at temperatures and pressures below those of the triple point. The simplest case is pure deuterium pellets ($\mathcal{T}$ < 18.7 K, $\mathcal{P}$ < 171.3 hPa). The production of good quality solid deuterium, capable of withstanding the mechanical stress during the acceleration of the pellets, is a key issue. To this end the phase transition of deuterium from gas to solid (and vice versa) is modeled with extensive molecular-dynamics (MD) simulations. Furthermore, the solid growth from the gas phase is simulated in an ample range of temperatures and pressures, to find the best compromise between growth velocity and mechanical properties of the resulting solid system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of Two Models for the Angular Structure of the Outflows Producing the Swift/XRT “Larger-angle Emission” of Gamma-Ray Bursts

The quasi-instantaneous emission from a relativistic surface endowed with a Lorentz factor that decreases away from the outflow symmetry axis can naturally explain the three phases observed by Swift X-Ray Telescope (XRT) in gamma-ray bursts (GRBs) and their afterglows (GRB tail, afterglow plateau, and postplateau) based only on the angular change of the relativistic Doppler boost across the outflow surface. We develop further the analytical formalism of the “larger-angle emission” model for the case of “n-exponential” outflows (where the Lorentz factor Γ dependence of the angular location θ is Γ ~ exp{-(θ/θ c ) n }), and compare its ability to account for the X-ray emission of XRT afterglows relative to that of “power-law” outflows (Γ ∼ θ −g ). Power-law outflows yield longer afterglow plateaus, followed by slower postplateau flux decays than n-exponential outflows, features which may be used in identifying which type of angular structure is at work in a given afterglow. Identifying the Γ(θ) angular structure that accommodates XRT light curves is slightly complicated by the fact that the afterglow X-ray light curve is also determined by how two characteristics of the comoving-frame emission spectrum (peak energy $E'_p$ and peak intensity $i'_p$) change with the angular location or, equivalently, with the Lorentz factor. Here, we assume power-law Γ dependences of those spectral characteristics and find that, unlike power-law outflows, n-exponential outflows cannot account for plateaus with a temporal dynamical range larger than 100 (2 dex in logarithmic space). To capture all the information contained in XRT afterglow measurements (0.3–10 keV unabsorbed flux and effective spectral slope), we calculate 0.3 and 10 keV light curves using a broken-power-law emission spectrum of peak energy and low- and high-energy slopes that are derived from the effective slope measured by XRT. This economical peak energy determination is found to be consistent with the results of more expensive spectral fits. The angular distributions of the Lorentz factor, comoving frame peak energy, and peak intensity (Γ(θ), $E'_p$(θ), $i'_p$(θ)) constrain the (yet-to-be determined) convolution of various features of the production of relativistic jets by solar-mass black holes and of their propagation through the progenitor/circumburst medium, while the $E'_p$(Γ) and $i'_p$(Γ) dependences may constrain the GRB dissipation mechanism and the GRB emission process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

An effective description of charge diffusion and energy transport in a charged plasma from holography

We discuss the physics of sound propagation and charge diffusion in a plasma with non-vanishing charge density. Our analysis culminates the program initiated in to construct an open effective field theory of low-lying modes of the stress tensor and charge current in such plasmas. We model the plasma holographically as a Reissner-Nordström-AdS d+1 black hole, and study linearized fluctuations of longitudinally polarized scalar gravitons and photons in this background. We demonstrate that the perturbations can be decoupled and repackaged into the dynamics of two designer scalars, whose gravitational coupling is modulated by a non-trivial dilatonic factor. The holographic analysis allows us to isolate the phonon mode from the charge diffusion mode, and identify the combination of currents that corresponds to each of them. We use these results to obtain the real-time Gaussian effective action, which includes both the retarded response and the associated stochastic (Hawking) fluctuations, accurate to quartic order in gradients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Unraveling excitation energy transfer assisted by collective behaviors of vibrations

Abstract We investigate how collective behaviors of vibrations such as cooperativity and interference can enhance energy transfer in a nontrivial way, focusing on an example of a donor–bridge–acceptor trimeric chromophore system coupled to two vibrational degrees of freedom. Employing parameters selected to provide an overall uphill energy transfer from donor to acceptor, we use numerical calculations of dynamics in a coupled exciton–vibration basis, together with perturbation-based analytics and calculation of vibronic spectra, to identify clear spectral features of single- and multi-phonon vibrationally-assisted energy transfer (VAET) dynamics, where the latter include up to six-phonon contributions. We identify signatures of vibrational cooperation and interference that provide enhancement of energy transfer relative to that obtained from VAET with a single vibrational mode. We observe a phononic analogue of two-photon absorption, as well as a novel heteroexcitation mechanism in which a single phonon gives rise to simultaneous excitation of both the trimeric system and the vibrational degrees of freedom. The impacts of vibrations and of the one- and two-phonon VAET processes on the energy transfer are seen to be quite different in the weak and strong site–vibration coupling regimes. In the weak coupling regime, two-phonon processes dominate, whereas in the strong coupling regime up to six-phonon VAET processes can be induced. The VAET features are seen to be enhanced with increasing temperature and site–vibration coupling strength, and are reduced in the presence of dissipation. We analyze the dependence of these phenomena on the explicit form of the chromophore–vibration couplings, with comparison of VAET spectra for local and non-local couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dissecting the contributions to non-photochemical quenching in a land plant under fluctuating light

Photosynthetic organisms have evolved multiple non-photochemical quenching (NPQ) processes, providing photoprotection by safely dissipating excess excitation energy. These processes involve various molecular players functioning on overlapping timescales from seconds to days, making it challenging to isolate and quantify their individual kinetics. In this study, we perform whole-leaf chlorophyll fluorescence lifetime and xanthophyll concentration measurements on wild-type and various newly characterized NPQ mutants of Nicotiana benthamiana, a vascular land plant. Based on these measurements, we construct a fluorescence lifetime-based quantitative kinetic model that disentangles individual photoprotection components and, when integrated additively, accurately predicts wild-type and mutant NPQ behaviors under various light-dark regimes. Additionally, the model quantifies the per-molecule quenching effectiveness of various xanthophylls and the contributions of six quenching components (qE V , qE A , qE Z , qE L , qZ, and qI) across different genotypes. It also suggests improved overall quenching efficiency at specific VDE:ZEP:PsbS overexpression stoichiometries, aligning with previous studies and supporting translational efforts to optimize photoprotection and enhance crop yields under dynamic light environments.

Lam, Lam [University of California, Berkeley, CA (↗
DOE OSTI2025

Enhancing Thermal Transport in Polymeric Composites Via Engineered Noncovalent Filler–Polymer Interactions

Understanding thermal transport mechanisms in polymeric composites allows us to expand the boundaries of thermal conductivity in them, either increasing it for more efficient heat dissipation or decreasing it for better thermal insulation. But, these mechanisms are not fully understood. Systematic experimental investigations remain limited. Practical strategies to tune the interfacial thermal resistance (ITR) between fillers and polymers and the thermal conductivity of composites remain elusive. Here, we studied the thermal transport in representative polymer composites, using polyethylene (PE) or polyaniline (PANI) as matrices and graphite as fillers. PANI, with aromatic rings in its backbone, interacts with graphite through strong noncovalent π–π stacking interactions, whereas PE lacks such interactions. We can then quantify how π–π stacking interactions between graphite and polymers enhance thermal transport in composites. PE/graphite and PANI/graphite composites with the same 1.5% filler volume fractions show a ∼22.82% and ∼34.85% enhancement in thermal conductivity compared to pure polymers, respectively. Calculated ITRs in PE/graphite and PANI/graphite are ~6 x 10 -8 m 2 KW -1 and ~1 x 10 -8 m 2 KW -1 , respectively, highlighting how π–π stacking interactions reduce ITR. Molecular dynamics (MD) simulations suggest that π–π stacking interactions between PANI chains and graphite surfaces enhance alignment of PANI's aromatic rings with graphite surfaces. This allows more carbon atoms from PANI chains to interact with graphite surfaces at a shorter distance compared to PE chains. Finally, our work indicates that tuning the π–π stacking interactions between polymers and fillers is an effective approach to reduce the ITR and enhance the thermal conductivity of composites.

atoms↗
DOE OSTI2026

How earthquakes organize stress

Stress is not uniform in the Earth. Therefore, we must use natural experiments to measure the distribution of stresses and related quantities, rather than single values. For instance, dynamic triggering shows that faults are uniformly distributed over their loading cycles in Southern California. The probability that a fault ruptures across a barrier measures the in situ energy distribution. Fault roughness reflects the distribution of strength. These natural experiments produce observable distributions that are surprisingly consistent and suggest some degree of self-organization in the Earth’s crust. Once established, the functional form of the distributions can be used to track changes in response to earthquakes as well as to distinguish fundamentally different fault systems. Transient fault locking before stress release in laboratory experiments can be interpreted as a consequence of self-organization of fault stress. The robust self-organization of multiple variables in earthquake systems suggests that the most consequential mechanical outcome of earthquakes may be the redistribution of stress and the strain energy associated with it. The low friction on a fault during seismic slip as inferred by temperature measurements of the Tohoku earthquake is consistent with dissipation playing a secondary role to this redistribution process. Through stress redistribution and interaction, subduction zone faults tend to synchronize, perhaps due to their geometric simplicity, while the continental system of Southern California cannot synchronize, perhaps due to the complexity of the fault network. Earthquakes organize stress in the crust and produce a suite of well-defined, consistent distributions.

earthquakes↗
DOE OSTI2024

Bilayer Ion Trap Design for 2D Arrays

Junctions are fundamental elements that support qubit locomotion in two-dimensional ion trap arrays and enhance connectivity in emerging trapped-ion quantum computers. In surface ion traps they have typically been implemented by shaping radio frequency (RF) electrodes in a single plane to minimize the disturbance to the pseudopotential. However, this method introduces issues related to RF lead routing that can increase power dissipation and the likelihood of voltage breakdown. Here, we propose and simulate a novel two-layer junction design incorporating two perpendicularly rotoreflected (rotated, then reflected) linear ion traps. The traps are vertically separated, and create a trapping potential between their respective planes. The orthogonal orientation of the RF electrodes of each trap relative to the other provides perpendicular axes of confinement that can be used to realize transport in two dimensions. While this design introduces manufacturing and operating challenges, as now two separate structures have to be precisely positioned relative to each other in the vertical direction and optical access from the top is obscured, it obviates the need to route RF leads below the top surface of the trap and eliminates the pseudopotential bumps that occur in typical junctions. Here in this paper the stability of idealized ion transfer in the new configuration is demonstrated, both by solving the Mathieu equation analytically to identify the stable regions and by numerically modeling ion dynamics. Our novel junction layout has the potential to enhance the flexibility of microfabricated ion trap control to enable large-scale trapped-ion quantum computing.

42 ENGINEERING↗
DOE OSTI2024

Dual Mechanism for Transient Capacitance Anomaly in Improper Ferroelectrics

Negative capacitance (NC) effects in ferroelectrics can potentially break fundamental limits of power dissipation known as “Boltzmann tyranny.” However, the origin of transient NC of ferroelectrics, which is attributed to two different mechanisms involving free-energy landscape and nucleation, is under intense debate. Here, we report the coexistence of transient NC and an S-shaped anomaly during the switching of ferroelectric hexagonal ferrites capacitor in an RC circuit. The early-stage NC arises from the nucleation process, while the late-stage S-shaped anomaly corresponds to a nascent NC associated with the free-energy landscape. The entire waveform can be reproduced using a hybrid model that simultaneously incorporates these two mechanisms. These results highlight the multivariable free-energy landscape of hexagonal ferrites that enables an abrupt change of the internal field and demonstrate that the two mechanisms are not mutually exclusive, resolving the long-standing debate. In conclusion, the behavior of the S-shaped anomaly also provides a pathway to extract parameters of free-energy landscape and switching dynamics.

Capacitance↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
DOE OSTI2023

Hybrid Simulation and Quasi-linear Theory of Bi-Kappa Proton Instabilities

The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropy T ⊥ > T ∥ (where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κ power-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Thermalization of nuclear matter in heavy-ion collisions at Fermi energies

We analyze the time evolution of the kinetic properties of nuclear matter produced in heavy-ion collisions at Fermi energies. The collision system is simulated using Constrained Molecular Dynamics (CoMD) transport calculations whose output is the isospin, position, and momentum of the nucleons. Focusing on central 35 A·MeV 40 Ca+ 40 Ca collisions we utilize this information to extract localized momentum distributions in volume elements of 8 fm 3 and time steps of 5 fm/$c$. We then parameterize the single-particle momentum distributions with thermally motivated fit functions in the local rest frame of each cell. While the transverse-momentum distributions are well reproduced by thermal ones, the longitudinal ones carry a marked imprint of the initial nuclear motion which we capture by introducing a centroid motion into our fit functions. In particular, we find that Fermi distributions yield significantly better fits than Boltzmann ones, a consequence of the Pauli blocking implemented in CoMD. From the fits we extract the time dependence of the thermodynamic and collective properties of the excited nuclear medium. We find that the transverse temperature gradually rises to about 6 MeV, which is accompanied by a dissipation of the initial centroid motion of the incoming nuclei which vanishes at about 100 fm/c after initial impact. In conclusion, we are therefore able to track the transition of beam energy into random kinetic energy for nucleons, suggesting a three-dimensional equilibration of energy in the late stages of the collision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantum dynamics in strongly driven random dipolar magnets

The random dipolar magnet LiHo x Y 1– x F 4 enters a strongly frustrated regime for small Ho 3+ concentrations with x < 0.05. In this regime, the magnetic moments of the Ho 3+ ions experience small quantum corrections to the common Ising approximation of LiHo x Y 1– x F 4 , which lead to a $\mathbb{Z}_2$-symmetry breaking and small, degeneracy breaking energy shifts between different eigenstates. Here we show that destructive interference between two almost degenerate excitation pathways burns spectral holes in the magnetic susceptibility of strongly driven magnetic moments in LiHo x Y 1– x F 4 . Such spectral holes in the susceptibility, microscopically described in terms of Fano resonances, can already occur in setups of only two or three frustrated moments, for which the driven level scheme has the paradigmatic Λ shape. For larger clusters of magnetic moments, the corresponding level schemes separate into almost isolated many-body Λ schemes, in the sense that either the transition matrix elements between them are negligibly small or the energy difference of the transitions is strongly off-resonant to the drive. This enables the observation of Fano resonances, caused by many-body quantum corrections to the common Ising approximation also in the thermodynamic limit. Finally, we discuss its dependence on the driving strength and frequency as well as the crucial role that is played by lattice dissipation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Simulation of 15% and 50% Thermal Power Dispatch to an Industrial Facility Using a Flexible Generic Full-Scope Pressurized Water Reactor Plant Simulator

Nuclear power plants in the United States are increasingly challenged to compete in wholesale electricity markets due to the low electricity costs and increasingly dynamic grid conditions from competing generation sources. An alternative market for nuclear power is industrial facilities that can use the thermal and/or electrical power generated by a nuclear power plant to offset the economic losses incurred by electricity market challenges. A generic pressurized water reactor (PWR) simulator was used to show the results of a basic design for a generic thermal power extraction system and tests were run using a set of procedures to show what happens when a nuclear power plant transitions from full electrical power dispatch to 15% and 50% thermal power dispatch. This type of operation leads to losses in turbine performance efficiency due to the deviation from the design operating point, but because the thermal power is also used by the industry load without conversion losses, the combined thermal efficiency of the PWR increases. For the 15% case, the thermal efficiency increased from 32% to 41.9%, while for the 50% case, the efficiency increased up to 60.1%. In addition, for 50% thermal power dispatch, the power dissipated by the condenser decreased from approximately 2000 to approximately 1300 MW (thermal), indicating a substantially diminished impact on the environment in terms of releasing heat into the cooling water reservoir.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Distinctive features of structural evolution and thermodynamic response in wide-bandgap semiconductors driven by intense electronic excitation

Radiation-tolerant material selection requires balancing lattice rigidity, defect dynamics, and electronic stability, as shown by covalent SiC outperforming ionic Ga 2 O 3 and GaN under extreme environments. Responding to intense electronic excitation, irradiation-driven phase segregation (β → δ/κ in Ga 2 O 3 ) and core–shell track (disordered structure in GaN), accompanied by elemental redistribution, contrastingly, exceptional radiation tolerance manifested by comparatively minimal lattice distortion (0.17 % strain variation) was demonstrated in SiC. These differential responses are primarily attributed to two fundamental mechanisms: (i) thermodynamic driving forces governing defect migration and phase separation, and (ii) the synergistic effects of robust covalent bonding composition coupled with efficient defect recombination processes. Here, the stronger electron–phonon (e-ph) coupling in Ga 2 O 3 (4.34 × 1018 W m −3 K −1 ) and GaN (3.55 × 10 18 W m −3 K −1 ) enhances lattice energy deposition, triggering thermal spikes (ΔT ≫ T m ) and structural transition behaviors, whereas weaker e-ph coupling in SiC (3.69 × 10 18 W m −3 K −1 ), relatively high thermodynamic parameters and efficient energy dissipation suppress thermal spikes to maintaining lattice integrity. The photoresponse degradation driven by enhanced radiative recombination is dominant in N-doped SiC, while V-doped systems achieve defect-mediated photoconduction optimization characterized by abrupt current transitions, matching fluorescence yield evolutions, and directly connecting defect engineering to optoelectronic performance.

Intense electronic excitation↗
DOE OSTI2024

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗
DOE OSTI2023

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗