Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical hydrodynamics”

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 181 records · Page 10

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatio-temporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Furthermore, our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of baryonic feedback on HSC-Y1 weak lensing non-Gaussian statistics

Baryonic feedback is a major systematic in weak lensing cosmology. Its most studied effect is the suppression of the lensing power spectrum, a second-order statistic, on small scales. Motivated by the growing interest in statistics beyond the second order, we investigate the effect of baryons on lensing non-Gaussian statistics and the resulting biases in the matter clustering amplitude S 8 = σ 8 Ω m / 0.3 . We focus on the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) data which, with its high source number density, closely resembles those expected from the upcoming Euclid and Rubin Legacy Survey of Space and Time. We study four non-Gaussian statistics of convergence maps—peak counts, minimum counts, the probability distribution function, and the scattering transform coefficients—in addition to the usual power spectrum. We first estimate the biases in S 8 using mock observations built from the IllustrisTNG and BAHAMAS hydrodynamical simulations and theoretical models built from dark-matter-only simulations. We find up to 1 σ bias in S 8 when the smallest scales (2 arcmin) and the highest feedback level are considered. We then analyze the HSC-Y1 data and compare the S 8 obtained for each statistic with different smoothing scales or scale cuts. As we expect that baryons mostly affect the small scales, comparing the results obtained from including and excluding small scales can indicate the level of impact from baryons. With HSC data, we find only minor ( ≤ 0.5 σ ) differences in S 8 for all statistics, even when considering very small scales (2 arcmin). Our results suggest that the effect of baryons is insignificant at the level of HSC-Y1 down to 2 arcmin for all statistics examined here, or it is canceled by other scale-dependent systematics.

79 ASTRONOMY AND ASTROPHYSICS↗

Viscous Dissipation in One-Dimensional Quantum Liquids

We develop a theory of viscous dissipation in one-dimensional single-component quantum liquids at low temperatures. Such liquids are characterized by a single viscosity coefficient, the bulk viscosity. We show that for a generic interaction between the constituent particles this viscosity diverges in the zerotemperature limit. In the special case of integrable models, the viscosity is infinite at any temperature, which can be interpreted as a breakdown of the hydrodynamic description. In conclusion, our consideration is applicable to all single-component Galilean- invariant one-dimensional quantum liquids, regardless of the statistics of the constituent particles and the interaction strength.

36 MATERIALS SCIENCE↗

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗

Experimenting with the GMAO 4D Data Assimilation

The Global Modeling and Assimilation Office (GMAO) has been working to promote its prototype four-dimensional variational (4DVAR) system to a version that can be exercised at operationally desirable configurations. Beyond a general circulation model (GeM) and an analysis system, traditional 4DV AR requires availability of tangent linear (TL) and adjoint (AD) models of the corresponding GeM. The GMAO prototype 4DVAR uses the finite-volume-based GEOS GeM and the Grid-point Statistical Interpolation (GSI) system for the first two, and TL and AD models derived ITom an early version of the finite-volume hydrodynamics that is scientifically equivalent to the present GEOS nonlinear GeM but computationally rather outdated. Specifically, the TL and AD models hydrodynamics uses a simple (I-dimensional) latitudinal MPI domain decomposition, which has consequent low scalability and prevents the prototype 4DV AR ITom being used in realistic applications. In the near future, GMAO will be upgrading its operational GEOS GCM (and assimilation system) to use a cubed-sphere-based hydrodynamics. This versions of the dynamics scales to thousands of processes and has led to a decision to re-derive the TL and AD models for this more modern dynamics, thus taking advantage of a two-dimensional MPI decomposition and improved scalability properties. With the aid of the Transformation of Algorithms in FORTRAN (l'AF) automatic adjoint generation tool and some hand-coding, a version of the cubed-sphere-based TL and AD models, with a simplified vertical diffusion scheme, is now available, enabling multiple configurations of standard implementations of 4DV AR in GEOS. Concurrent to this development, collaboration with the National Centers for Environmental Prediction (NCEP) and the Earth System Research Laboratory (ESRL) has allowed GMAO to implement a hybrid-ensemble capability within the GEOS data assimilation system. Both 3Dand 4D-ensemble capabilities are presently available thus allowing GMAO to now evaluate the performance and benefit of various ensemble and variational assimilation strategies. This presentation will cover the most recent developments taking place at GMAO and show results from various comparisons from traditional techniques to more recent ensemble-based ones.

Todling, R.↗

Crystallization and collapse in relativistically degenerate matter

In this paper, it is shown that a mass density limit exists beyond which the relativistically degenerate matter would crystallize. The mass density limit, found here, is quite analogous to the mass limit predicted by Chandrasekhar for a type of compact stars called white dwarfs (M{sub Ch} Asymptotically-Equal-To 1.43 Solar Mass). In this study, the old problem of white dwarf core collapse, which has been previously investigated by Chandrasekhar using hydrostatic stability criteria, is revisited in the framework of the quantum hydrodynamics model by inspection of the charge screening at atomic scales in the relativistic degeneracy plasma regime taking into account the relativistic Fermi-Dirac statistics and electron interaction features such as the quantum statistical pressure, Coulomb attraction, electron exchange-correlation, and quantum recoil effects. It is revealed that the existence of ion correlation and crystallization of matter in the relativistically degenerate plasma puts a critical mass density limit on white dwarf core region. It is shown that a white dwarf star with a core mass density beyond this critical limit can undergo the spontaneous core collapse (SCC). The SCC phenomenon, which is dominantly caused by the electron quantum recoil effect (interference and localization of the electron wave function), leads to a new exotic state of matter. In such exotic state, the relativistic electron degeneracy can lead the white dwarf crystallized core to undergo the nuclear fusion and an ultimate supernova by means of the volume reduction (due to the enhanced compressibility) and huge energy release (due to the increase in cohesive energy), under the stars huge inward gravitational pressure. Moreover, it is found that the SCC phenomenon is significantly affected by the core composition (it is more probable for heavier plasmas). The critical mass density found here is consistent with the values calculated for core density of typical white dwarf stars.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reactive Hydrodynamics in Rotating Spherical and Cylindrical Geometry

In turbulent spray combustion among many complex interactions between local flow structures called turbulent eddies and droplets are those associated with rotation of droplets. In general, for a complete statistical description of turbulent sprays, consideration of at least four degrees of freedom respectively associated with translational, rotational, vibrational (pulsational), and internal motions of the droplet are needed. Clearly the interactions between all degrees of freedom of the droplets and those for the gaseous background field will be exceedingly complex. For example, one type of interaction between the translational and the rotational velocity of droplets results in droplet helicity, H(d) = w(d).v(d), the significance of which in turbulent spray combustion is yet to be recognized. The role of droplet rotation in turbulent spray combustion modeling and its impact on the evaporation of liquid fuel droplets was recently investigated. Also, the impact of rotation on combustion of solid particles such as is encountered in pulverized coal combustion has been emphasized. The problem of viscous flow around a rotating sphere discussed above also occurs in other areas of physical sciences such as astrophysics and geophysics. Consequently, the subject has been addressed in many classical as well as more recent investigations. According to these investigations, the rotation of a rigid sphere in an otherwise quiescent, unconfined environment results in the motion of the fluid towards the poles. The polar flows from the northern and southern hemispheres move along helical trajectories towards the equatorial plane. Eventually, the polar flows collide at the equatorial plane, thus producing a sheet of rotating fluid that is radially ejected outward on this plane. Therefore, a droplet induces a strained flow field as a result of its rotation. Since the spatial extent of equatorial jets could easily exceed many droplet diameters, interactions between neighboring droplets are enhanced as a result of their rotation. Also, the equatorial jet substantially alters the spherical geometry of the diffusion flame surface that surrounds a rotating droplet. The objective of the research is to gain more knowledge about the hydrodynamics within and around rotating spherical and cylindrical body of fluid, and the behavior of diffusion or premixed flame surfaces that could surround such symmetric body of rotating fluid.

Sohrab, Siavash H.↗

TURBULENT AMPLIFICATION AND STRUCTURE OF THE INTRACLUSTER MAGNETIC FIELD

We compare DNS calculations of homogeneous isotropic turbulence with the statistical properties of intracluster turbulence from the Matryoshka Run and find remarkable similarities between their inertial ranges. This allowed us to use the time-dependent statistical properties of intracluster turbulence to evaluate dynamo action in the intracluster medium, based on earlier results from a numerically resolved nonlinear magneto-hydrodynamic turbulent dynamo. We argue that this approach is necessary (a) to properly normalize dynamo action to the available intracluster turbulent energy and (b) to overcome the limitations of low Re affecting current numerical models of the intracluster medium. We find that while the properties of intracluster magnetic field are largely insensitive to the value and origin of the seed field, the resulting values for the Alfvén speed and the outer scale of the magnetic field are consistent with current observational estimates, basically confirming the idea that the magnetic field in today’s galaxy clusters is a record of its past turbulent activity.

79 ASTRONOMY AND ASTROPHYSICS↗

Shape and flow fluctuations in ultracentral Pb + Pb collisions at the energies available at the CERN Large Hadron Collider

In ultracentral heavy-ion collisions, anisotropic hydrodynamic flow is generated by density fluctuations in the initial state rather than by geometric overlap effects. For a given centrality class, the initial fluctuation spectrum is sensitive to the method chosen for binning the events into centrality classes. We show that sorting events by total initial entropy or by total final multiplicity yields event classes with equivalent statistical fluctuation properties, in spite of viscous entropy production during the fireball evolution. With this initial entropy-based centrality definition we generate several classes of ultracentral Pb + Pb collisions at CERN Large Hadron Collider energies and evolve the events using viscous hydrodynamics with nonzero shear but vanishing bulk viscosity. Comparing the predicted anisotropic flow coefficients for charged hadrons with CMS data we find that both the Monte Carlo Glauber (MC-Glb) and Monte Carlo Kharzeev-Levin-Nardi (MC-KLN) models produce initial fluctuation spectra that are incompatible with the measured final anisotropic flow power spectrum, for any choice of the specific shear viscosity. In spite of this failure, we show that the hydrodynamic model can qualitatively explain, in terms of event-by-event fluctuations of the anisotropic flow coefficients and flow angles, the breaking of flow factorization for elliptic, triangular, and quadrangular flow measured by the CMS experiment. For elliptic flow, this factorization breaking is large in ultracentral collisions. Finally, we conclude that the bulk of the experimentally observed flow factorization breaking effects are qualitatively explained by hydrodynamic evolution of initial-state fluctuations, but that their quantitative description requires a better understanding of the initial fluctuation spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled Langmuir oscillations in 2-dimensional quantum plasmas

In this work, we present a hydrodynamic model to study the coupled quantum electron plasma oscillations (QEPO) for two dimensional (2D) degenerate plasmas, which incorporates all the essential quantum ingredients such as the statistical degeneracy pressure, electron-exchange, and electron quantum diffraction effect. Effects of diverse physical aspects like the electronic band-dispersion effect, the electron exchange-correlations and the quantum Bohm-potential as well as other important plasma parameters such as the coupling parameter (plasma separation) and the plasma electron number-densities on the linear response of the coupled system are investigated. By studying three different 2D plasma coupling types, namely, graphene-graphene, graphene-metalfilm, and metalfilm-metalfilm coupling configurations, it is remarked that the collective quantum effects can influence the coupled modes quite differently, depending on the type of the plasma configuration. It is also found that the slow and fast QEPO frequency modes respond very differently to the change in plasma parameters. Current findings can help in understanding of the coupled density oscillations in multilayer graphene, graphene-based heterojunctions, or nanofabricated integrated circuits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Probing cosmic velocities with the pairwise kinematic Sunyaev-Zel’dovich signal in DESI Bright Galaxy Sample DR1 and ACT DR6

We present a measurement of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) signal using the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Sample (BGS) Data Release 1 (DR1) galaxy sample overlapping with the Atacama Cosmology Telescope (ACT) CMB temperature map. Our analysis makes use of 1.6 million galaxies with stellar masses log⁡ 𝑀 ⋆ /𝑀 ⊙ >10, and we explore measurements across a range of aperture sizes (2.1′ <𝜃 ap <3.5′) and stellar mass selections. This statistic directly probes the velocity field of the large-scale structure, a unique observable of cosmic dynamics and modified gravity. In particular, at low redshifts, this quantity is especially interesting, as deviations from General Relativity are expected to be largest. Notably, our result represents the highest-significance low-redshift (𝑧 ∼ 0.3) detection of the kSZ pairwise effect yet. In our most optimal configuration (𝜃 ap =3.3′, log⁡ 𝑀 ⋆ >11), we achieve a 5⁢𝜎 detection. Assuming that an estimate of the optical depth and galaxy bias of the sample exists via e.g., external observables, this measurement constrains the fundamental cosmological combination 𝐻 0 ⁡𝑓⁡𝜎$^2_8$. A key challenge is the degeneracy with the galaxy optical depth. We address this by combining CMB lensing, which allows us to infer the halo mass and galaxy population properties, with hydrodynamical simulation estimates of the mean optical depth, $\bar{𝜏}$ . We stress that this is a proof-of-concept analysis; with BGS DR2 data we expect to improve the statistical precision by roughly a factor of two, paving the way toward robust tests of modified gravity with kSZ-informed velocity-field measurements at low redshift.

Hadzhiyska, Boryana [Institute of Astronomy; Kavli↗

Modelling the Lyman-α forest with Eulerian and SPH hydrodynamical methods

ABSTRACT We compare two state-of-the-art numerical codes to study the overall accuracy in modelling the intergalactic medium and reproducing Lyman-α forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modelling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite-volume method. The second code, CRK-HACC , uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modelled with particles, to treat the intergalactic gas. We compare the convergence behaviour of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable ≲1 per cent agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in CRK-HACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, z ∼ 5, as the Lyman-α forest mostly comes from lower-density regions that are intrinsically poorly sampled by SPH particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Chiral Symmetry Breaking in Crystal Growth: Is Hydrodynamic Convection Relevant?

The effects of mechanical stirring on nucleation and chiral symmetry breaking have been investigated for a simple inorganic molecule, sodium chlorate (NaClO3). In contrast to earlier findings, our experiment suggests that the symmetry breaking may have little to do with hydrodynamic convection. Rather the effect can be reasonably accounted for by mechanical damage to incipient crystals. The catastrophic events, creating numerous small 'secondary' crystals, produce statistical domination of one chiral species over the other. Our conclusion is supported by a number of observations using different mixing mechanisms.

Martin, B.↗

Instability and dynamics of two nonlinearly coupled intense laser beams in a quantum plasma

We consider nonlinear interactions between two relativistically strong laser beams and a quantum plasma composed of degenerate electron fluids and immobile ions. The collective behavior of degenerate electrons is modeled by quantum hydrodynamic equations composed of the electron continuity, quantum electron momentum (QEM) equation, as well as the Poisson and Maxwell equations. The QEM equation accounts the quantum statistical electron pressure, the quantum electron recoil due to electron tunneling through the quantum Bohm potential, electron-exchange, and electron-correlation effects caused by electron spin, and relativistic ponderomotive forces (RPFs) of two circularly polarized electromagnetic (CPEM) beams. The dynamics of the latter are governed by nonlinear wave equations that include nonlinear currents arising from the relativistic electron mass increase in the CPEM wave fields, as well as from the beating of the electron quiver velocity and electron density variations reinforced by the RPFs of the two CPEM waves. Furthermore, nonlinear electron density variations associated with the driven (by the RPFs) quantum electron plasma oscillations obey a coupled nonlinear Schroedinger and Poisson equations. The nonlinearly coupled equations for our purposes are then used to obtain a general dispersion relation (GDR) for studying the parametric instabilities and the localization of CPEM wave packets in a quantum plasma. Numerical analyses of the GDR reveal that the growth rate of a fastest growing parametrically unstable mode is in agreement with the result that has been deduced from numerical simulations of the governing nonlinear equations. Explicit numerical results for two-dimensional (2D) localized CPEM wave packets at nanoscales are also presented. Possible applications of our investigation to intense laser-solid density compressed plasma experiments are highlighted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gaussian distributions, Jacobi group, and Siegel-Jacobi space

Let N be the space of Gaussian distribution functions over ℝ, regarded as a 2-dimensional statistical manifold parameterized by the mean μ and the deviation σ. In this paper, we show that the tangent bundle of N, endowed with its natural Kähler structure, is the Siegel-Jacobi space appearing in the context of Number Theory and Jacobi forms. Geometrical aspects of the Siegel-Jacobi space are discussed in detail (completeness, curvature, group of holomorphic isometries, space of Kähler functions, and relationship to the Jacobi group), and are related to the quantum formalism in its geometrical form, i.e., based on the Kähler structure of the complex projective space. This paper is a continuation of our previous work [M. Molitor, “Remarks on the statistical origin of the geometrical formulation of quantum mechanics,” Int. J. Geom. Methods Mod. Phys. 9(3), 1220001, 9 (2012); M. Molitor, “Information geometry and the hydrodynamical formulation of quantum mechanics,” e-print arXiv (2012); M. Molitor, “Exponential families, Kähler geometry and quantum mechanics,” J. Geom. Phys. 70, 54–80 (2013)], where we studied the quantum formalism from a geometric and information-theoretical point of view.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reconstructing 3D asymmetries in laser-direct-drive implosions on OMEGA

Three-dimensional reconstruction algorithms have been developed which determine the hot-spot velocity, hot-spot apparent ion temperature distribution, and fuel areal density distribution present in laser-direct-drive inertial con finement fusion implosions on the OMEGA laser. These reconstructions rely on multiple independent measurements of the neutron energy spectrum emitted from the fusing plasma. Measurements of the neutron energy spectrum on OMEGA are made using a suite of quasi-orthogonal neutron time-of-flight detectors and a magnetic recoil spectrometer. These spectrometers are positioned strategically around the OMEGA target chamber to provide unique 3-D measurements of the conditions of the fusing hot spot and compressed fuel near peak compression. Here, the uncertainties involved in these 3-D reconstructions are discussed and are used to identify a new nTOF diagnostic line of sight which when built will reduce the uncertainty in the hot-spot apparent ion temperature distribution from 0.7 to < 0.4 keV.

47 OTHER INSTRUMENTATION↗

Equation of state at finite densities for QCD matter in nuclear collisions

We construct the QCD equation of state at finite chemical potentials including net baryon, electric charge, and strangeness, based on the conserved charge susceptibilities determined from lattice QCD simulations and the equation of state of the hadron resonance gas model. For the application to relativistic heavy-ion collisions we consider the situation of strangeness neutrality and matter with a fixed electric charge-to-baryon ratio, resembling that of heavy nuclei. The importance of finite electric charge and strangeness chemical potentials for particle production in heavy-ion collisions is demonstrated using hydrodynamic simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗