Modeling of rarefied ionized flows with the QUIPS and velocity-space hybrid methods.
Abstract not provided.
Oblapenko, Georgii↗
SEARCH · Search NASA
Results for “Velocity Model”
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.
A single-point modeling approach for the intercomparison and evaluation of ozone dry deposition across chemical transport models (Activity 2 of AQMEII4)
A primary sink of air pollutants and their precursors is dry deposition. Dry deposition estimates differ across chemical transport models, yet an understanding of the model spread is incomplete. Here, we introduce Activity 2 of the Air Quality Model Evaluation International Initiative Phase 4 (AQMEII4). We examine 18 dry deposition schemes from regional and global chemical transport models as well as standalone models used for impact assessments or process understanding. We configure the schemes as single-point models at eight Northern Hemisphere locations with observed ozone fluxes. Single-point models are driven by a common set of site-specific meteorological and environmental conditions. Five of eight sites have at least 3 years and up to 12 years of ozone fluxes. The interquartile range across models in multiyear mean ozone deposition velocities ranges from a factor of 1.2 to 1.9 annually across sites and tends to be highest during winter compared with summer. No model is within 50 % of observed multiyear averages across all sites and seasons, but some models perform well for some sites and seasons. For the first time, we demonstrate how contributions from depositional pathways vary across models. Models can disagree with respect to relative contributions from the pathways, even when they predict similar deposition velocities, or agree with respect to the relative contributions but predict different deposition velocities. Both stomatal and nonstomatal uptake contribute to the large model spread across sites. Our findings are the beginning of results from AQMEII4 Activity 2, which brings scientists who model air quality and dry deposition together with scientists who measure ozone fluxes to evaluate and improve dry deposition schemes in the chemical transport models used for research, planning, and regulatory purposes.
54 ENVIRONMENTAL SCIENCES↗
Multi-fidelity information fusion with concatenated neural networks
Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.
42 ENGINEERING↗
Uniting Surface Properties With Hydrodynamic Roughness in Shallow Overland Flow Models
Describing flow resistance from the properties of an underlying surface is a challenge in surface hydrology. Runoff models must specify a resistance formulation or “roughness scheme”—describing the functional relationship between flow resistance and flow depth/velocity—and its parameters. Uncertainty in runoff predictions derives from both the selected roughness scheme (e.g., Darcy Weisbach, Manning's, or laminar flow equations), and its parameterization with a roughness coefficient (e.g., Manning's n). Both choices are informed by model calibration to data, usually discharge, and, if available, velocity. In this study, a Saint Venant Equation-based runoff model is calibrated to discharge and velocity data from 112 rainfall simulator experiments. The results are used to identify the optimal roughness scheme among four widely-used options for each experiment, and to explore whether surface properties can be used to select the optimal roughness scheme and its coefficient. Among the tested roughness schemes, a transitional flow equation provided the best fit to the plurality of experiments. The most suitable roughness scheme for a given experiment was not related to measured surface properties. Regression models predicted the calibrated roughness coefficients with adjusted values between 0.48 and 0.54, depending on the roughness scheme used. Litter cover was the best predictor of the roughness coefficient, followed by soil cover and average canopy gap size. The results suggest that selection of an optimal roughness scheme based on surface properties alone remains difficult, but that once a scheme is selected, roughness coefficients can be estimated from surface properties.
54 ENVIRONMENTAL SCIENCES↗
Investigation of turbulent inflow specification in Euler–Lagrange simulations of mid-field spray
The process of atomization of a liquid jet by a parallel high-speed gas stream results in a spray, whose downstream development is of considerable interest to several applications. The round jet spray can be spatially divided into (i) a near-field (near-nozzle) region of liquid atomization and (ii) a downstream mid-field region of fully-dispersed droplets. In order to accurately model mid-field droplet dispersion, this work aims at developing a rigorous and robust injection model for Euler–Lagrange spray simulations. Results from experiments are used to obtain the relevant droplet number density, size distribution, and mean and standard deviation velocity distributions of the injection model, systematically in a step-by-step process. Two-phase large eddy simulations are performed by stochastically generating the Lagrangian droplets at the inlet of the mid-field region. Number flux, diameter distribution, mean velocity, and other time-averaged statistics at several downstream locations are shown to agree well with the corresponding experimental data.
42 ENGINEERING↗
Non-Virialized Axion Search Sensitive to Doppler Effects in the Milky Way Halo
The Axion Dark Matter eXperiment (ADMX) has previously excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions between 680-790 MHz under the assumption that the dark matter is described by the isothermal halo model. However, the precise nature of the velocity distribution of dark matter is still unknown, and alternative models have been proposed. We report the results of a non-virialized axion search over the mass range 2.81-3.31 μeV, corresponding to the frequency range 680-800 MHz. This analysis marks the most sensitive search for non-virialized axions sensitive to Doppler effects in the Milky Way Halo to date. Accounting for frequency shifts due to the detector's motion through the Galaxy, we exclude cold flow relic axions with a velocity dispersion of order 10^-7 c with 95% confidence.
79 ASTRONOMY AND ASTROPHYSICS↗
Nonvirialized axion search sensitive to Doppler effects in the Milky Way halo
The Axion Dark Matter eXperiment (ADMX) has previously excluded Dine-Fischler-Srednicki- Zhitnisky (DFSZ) axions between 680-790 MHz under the assumption that the dark matter is described by the isothermal halo model. However, the precise nature of the velocity distribution of dark matter is still unknown, and alternative models have been proposed. Here, we report the results of a non-virialized axion search over the mass range 2.81–3.31 µeV, corresponding to the frequency range 680–800 MHz. This analysis marks the most sensitive search for non-virialized axions sensitive to Doppler effects in the Milky Way Halo to date. Accounting for frequency shifts due to the detector’s motion through the Galaxy, we exclude cold flow relic axions with a velocity dispersion of $\mathscr{O}$(10 -7 )c with 95% confidence.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Tidal Stream Energy Resource Characterization in the Salish Sea
Tidal stream energy holds great potential as a renewable energy source in regions of the world where tidal currents are strong and predictable. The tidal stream energy resource of a region is strongly controlled by its tidal wave characteristics, the local bathymetry, and coastal geometric features. The Salish Sea, a large estuary on the Pacific Northwest coast, represents a great tidal stream resource because of its strong tidal currents in many tidal channels. However, the tidal resource in the Salish Sea has not been systematically characterized, primarily because of its large area and complex bathymetry and coastlines. This paper presents a modeling study conducted to characterize the tidal energy resource of the Salish Sea based on a high-resolution three-dimensional tidal hydrodynamic model of the Salish Sea, which was extensively validated using data derived from 10 tidal gauge stations and 132 Acoustic Doppler Current Profiler stations. Model validation results indicated the Salish Sea hydrodynamic model is skillful in simulating tidal wave propagation and velocity distributions in the Salish Sea. Based on model results, a total of 16 tidal channels with strong tidal currents, 9 in the San Juan Islands and Rosario Strait regions and 7 in the Puget Sound, were identified as hotspots for potential tidal energy development. Velocity probability distributions and exceedance curves of cross-channel average velocity magnitudes were calculated at all 16 channels based on the recommendations of the International Electrotechnical Commission Technical Specifications for tidal energy resource characterization. The tidal energy resource at the 16 hotspots was also characterized using power density distributions and kinetic energy fluxes. The ranking of the kinetic energy fluxes suggested that Admiralty Inlet, Rosario Strait, and Middle Channel are the top three hotspot sites in the Salish Sea for tidal stream energy development. The study demonstrated the need for a high-resolution 3-D modeling framework for accurate simulation of tidal currents in large complex estuarine systems, in the context of tidal energy resource characterization.
tidal stream energy, resource characterization, Sa↗
Extensions to models of the galaxy–halo connection
ABSTRACT We explore two widely used empirical models for the galaxy–halo connection, subhalo abundance matching (SHAM) and the halo occupation distribution (HOD), and compare them with the hydrodynamical simulation IllustrisTNG (TNG) for multiple statistics quantifying the galaxy distribution at $n_{\rm gal}\approx 1.3\times 10^{-3}\, ({\rm Mpc}\,h^{-1})^{-3}$. We observe that in their most straightforward implementations, both models fail to reproduce the two-point clustering measured in TNG. We find that SHAM models that use the relaxation velocity, Vrelax, and the peak velocity, Vpeak, perform best, and match the clustering reasonably well, although neither captures adequately the one-halo clustering. Splitting the total sample into sub-populations, we discover that SHAM overpredicts the clustering of high-mass, blue, star-forming, and late-forming galaxies and underpredicts that of low-mass, red, quiescent, and early-forming galaxies. We also study various baryonic effects, finding that subhaloes in the dark-matter-only simulation have consistently higher values of their SHAM-proxy properties than their full-physics counterparts. We then consider a 2D implementation of the HOD model augmented with a secondary parameter (environment, velocity anisotropy, σ2Rhalf-mass, and total potential) tuned so as to match the two-point clustering of the IllustrisTNG galaxies on large scales. We analyse these galaxy populations adopting alternative statistical tools such as galaxy–galaxy lensing, void–galaxy cross-correlations, and cumulants of the density field, finding that the hydrodynamical galaxy distribution disfavours σ2Rhalf-mass and the total potential as secondary parameters, while the environment and velocity anisotropy samples are consistent with full physics across all statistical probes examined. Our results demonstrate the power of examining multiple statistics for determining the secondary parameters that are vital for understanding the galaxy–halo connection.
79 ASTRONOMY AND ASTROPHYSICS↗
Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF
Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.
Smolentsev, Sergey↗
Origin of spin–orbit misalignments: the microblazar V4641 Sgr
ABSTRACT Of the known microquasars, V4641 Sgr boasts the most severe lower limit (>52°) on the misalignment angle between the relativistic jet axis and the binary orbital angular momentum. Assuming the jet and black hole spin axes coincide, we attempt to explain the origin of this extreme spin–orbit misalignment with a natal kick model, whereby an aligned binary system becomes misaligned by a supernova kick imparted to the newborn black hole. The model inputs are the kick velocity distribution, which we measure customized to V4641 Sgr, and the immediate pre/post-supernova binary system parameters. Using a grid of binary stellar evolution models, we determine post-supernova configurations that evolve to become consistent with V4641 Sgr today and obtain the corresponding pre-supernova configurations by using standard prescriptions for common envelope evolution. Using each of these potential progenitor system parameter sets as inputs, we find that a natal kick struggles to explain the origin of the V4641 Sgr spin–orbit misalignment. Consequently, we conclude that evolutionary pathways involving a standard common envelope phase followed by a supernova kick are highly unlikely for V4641 Sgr. An alternative interpretation is that the jet axis does not reliably trace the black hole spin axis. Our results raise concerns about compact object merger statistics gleaned from binary population synthesis models, which rely on unverified prescriptions for common envelope evolution and natal kicks. We also challenge the spin–orbit alignment assumption routinely invoked to measure black hole spin magnitudes.
79 ASTRONOMY AND ASTROPHYSICS↗
The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation
We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.
42 ENGINEERING↗
The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation
We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.
42 ENGINEERING↗
Velocity Independent Constraints on Spin-Dependent DM-Nucleon Interactions from IceCube and PICO
Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 $\hbox {C}_3\hbox {F}_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (JCAP 1509: 052, 2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma _{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm {cm}^2$ ($6 \times 10^{-38} \mathrm {cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1 TeV annihilating into $\tau ^+ \tau ^-$ ($b\bar{b}$) with a local density of $\rho _{\mathrm{DM}} = 0.3~\mathrm {GeV/cm}^3$. The constraints scale inversely with $\rho _{\mathrm{DM}}$ and are independent of the DM velocity distribution.
79 ASTRONOMY AND ASTROPHYSICS↗
Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation
Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.
54 ENVIRONMENTAL SCIENCES↗
Radiation and heat transport in divergent shock–bubble interactions
Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinematic Evidence for Bipolar Ejecta Flows in the Galactic Supernova Remnant W49B
W49B is a unique Galactic supernova remnant with centrally peaked, “bar”-like ejecta distribution, which was once considered evidence for a hypernova origin that resulted in a bipolar ejection of the stellar core. However, chemical abundance measurements contradict this interpretation. Closely connected to the morphology of the ejecta is its velocity distribution, which provides critical details for understanding the explosion mechanism. We report the first ever observational constraint on the kinematics of the ejecta in W49B using the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM). Using XRISM/Resolve, we measured the line-of-sight velocity traced by the Fe Heα emission, which is the brightest feature in the Resolve spectrum, to vary by ±300 km s −1 with a smooth east-to-west gradient of a few tens of kilometers per second per parsec along the major axis. Similar trends in the line-of-sight velocity structure were found for other Fe-group elements Cr and Mn, traced by the Heα emission, and also for intermediate-mass elements Si, S, Ar, and Ca, traced by the Lyα emission. The discovery of the east–west gradient in the line-of-sight velocity, together with the absence of a twin-peaked line profile or enhanced broadening in the central region, clearly rejects the equatorially expanding disk model. In contrast, the observed velocity structure suggests bipolar flows reminiscent of a bipolar explosion scenario. An alternative scenario would be a collimation of the ejecta by an elongated cavity sculpted by bipolar stellar winds.
High resolution spectroscopy↗
Optimal white-noise stochastic forcing for linear models of turbulent channel flow
In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ( $Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components.
Mechanics↗