Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 793 records · Page 44

Nanoscale Imaging of Magnetotransport around a Circular 𝑝−𝑛 Junction in Graphene

Magnetoresistance studies of 2D systems are often shaped by the motion of electrons that occupy spatially confined wave functions, such as topological edge modes and disorder-induced bound states. Directly probing how such states form and behave in situ , under applied currents, provides a clear way of connecting microscopic physics to the macroscopic transport response. Here, in this Letter, scanning tunneling potentiometry is used to probe the local, current-induced electrochemical potential of carriers in graphene near circular 𝑝−𝑛 junctions in an out-of-plane magnetic field ranging from 0 to 1.4 T. These measurements provide detailed information about the motion of carriers at the nanometer scale, revealing how it evolves with increasing field. The electrochemical potential displays distinct patterns, such as dipoles, spirals, and concentric disks in weak, moderate, and high fields, respectively. The size and orientation of these patterns can be used to understand how local carrier dynamics change in different transport regimes as well as to directly extract physical parameters such as the electron mean free path and cyclotron diameter.

36 MATERIALS SCIENCE↗

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↗

Low-Energy Backgrounds in Solid-State Phonon and Charge Detectors

Solid-state phonon and charge detectors probe the scattering of weakly interacting particles, such as dark matter and neutrinos, through their low recoil thresholds. Recent advancements have pushed sensitivity to eV-scale energy depositions, uncovering previously unseen low-energy excess backgrounds. While some arise from known processes such as thermal radiation, luminescence, and stress, others remain unexplained. This review examines these backgrounds, their possible origins, and parallels to low-energy effects in solids, an understanding of which is essential for interpreting particle interactions at and below the eV scale.

coherent neutrino–nucleus scattering↗

Low-Energy Backgrounds in Solid-State Phonon and Charge Detectors

Solid-state phonon and charge detectors probe the scattering of weakly interacting particles, such as dark matter and neutrinos, through their low recoil thresholds. Recent advancements have pushed sensitivity to eV-scale energy depositions, uncovering previously-unseen low-energy excess backgrounds. While some arise from known processes such as thermal radiation, luminescence, and stress, others remain unexplained. This review examines these backgrounds, their possible origins, and parallels to low-energy effects in solids. Their understanding is essential for interpreting particle interactions at and below the eV-scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Scattering of acoustic waves into Tollmien-Schlichting waves by small streamwise variations in surface geometry

By using the triple-deck scaling of Stewartson (1969) and Messiter (1970) it is shown that small but relatively sudden surface geometry variations that produce only very weak static pressure variations can nevertheless produce strong, i.e. O(1), coupling between an externally imposed acoustic disturbance and a spatially growing Tollmien-Schlichting wave. The analysis provides a qualitative explanation of the Leehey and Shapiro (1979) boundary-layer receptivity measurements and is in good quantitative agreement with the Aizin and Poliakov (1979) experiment. It may also explain why small 'trip wires' can promote early transition.

Goldstein, M. E.↗

The comet/solar wind transition region at Giacobini-Zinner

An account is given of the electron density, temperature and flow speed measurements made during the encounter with Comet Giacobini-Zinner by the Los Alamos plasma electron experiment on the ICE spacecraft. Between about 70,000 and 120,000 km from the nucleus of the comet, ICE found a region in which the solar wind flow speed decreased and the temperature increased by factors of about two. This transition region was characterized by large fluctuations in the plasma parameters and by highly variable electron velocity distributions. Electron temperature and density variations through the transition region reveal that ICE never crossed a short-scalelength bow shock during the encounter, although the scalelength for the gross transition in plasma properties is not incompatible with a shock dominated by the length scales of cometary ions. However, many of the electron distributions in the transition region and sheath are similar to those seen behind weak collisionless shocks elsewhere in the heliosphere. A model of the comet/solar wind interaction is suggested in which a standing shock exists sunward of the ICE trajectory, but, due to large variations in the upstream conditions, it is highly variable and perhaps only intermittent along the flanks.

Thomsen, M. F.↗

Ion flow at comet Halley

The three-dimensional positive ion analyzer aboard the Giotto spacecraft has been used to study the interaction between protons and alpha-particles in the solar wind and positive ions from comet Halley. Although the first impression of the overall structure is that the plasma flow evolves smoothly as the nucleus is approached, three sharp transitions of relatively small amplitude can be identified on both the inbound and outbound legs of the trajectory. The outermost one, at about one million km from the nucleus, appears to be a multiple crossing of a weak bow shock. The innermost one, at 80,000 km, is the boundary where the flowing plasma becomes depleted. On a microscopic scale, the turbulence created by the interaction between the two ion populations extends to a distance of several million km from the nucleus. At Giotto's closest approach to the nucleus, the plasma produced around the spacecraft by dust and gas impacts was much more energetic than had been expected.

Johnstone, A.↗

A multidimensional model for aerosols - Description of computational analogs

The numerical algorithms used to simulate the advection, diffusion, sedimentation, coagulation, and condensational growth of atmospheric aerosols are described. The model can be used in one, two, or three spatial dimensions. The continuity equation in a generalized horizontal and vertical coordinate system is developed, which allows the model to be quickly adapted to a wide variety of dynamical models of global or regional scale. Algorithms are developed to treat the various physical processes, and the results of simulations are presented, which show the strengths and weaknesses of these algorithms. Although the emphasis is on the modeling of aerosols, the work is also applicable to the simulations of the transport of gases.

Toon, O. B.↗

Timing analysis of binary X ray pulsars observed by HEAO 1

Timing analysis of the pointed observations of 5 selected binary X ray pulsars, in two categories of disk-fed and wind-fed sources, by HEAO 1 A-2 is reported. The power spectral analysis was performed on the data in the frequency range from approx. 1 MHz to 6.25 Hz. The coherent signal of the pulsation, the continuum of the power spectrum, varies in time and differs among sources. Quasi-Periodic Oscillation (QPO) is probably related to a fast spinning but weakly magnetized neutron star in the low mass X ray binaries. QPO was searched for in this frequency range to see if scaling laws exist among these two systems which may have possessed different order of magnitude of magnetic field strengths and the inner disk radii. One possible QPO is found centering at 0.062 Hz in 4U0115+63 during a flare.

Soong, Y.↗

Linear mixing rule in screened binary ionic mixtures

The validity of the linear mixing rule is examined for the following two cases (1) when the response of the electron gas is taken into account in the effective ionic interaction and (2) when finite-temperature effects are included in the dielectric response of the electrons, i.e., when the ions interact with both temperature- and density-dependent screened Coulomb potentials. It is found that the linear mixing rule remains valid when the electron response is taken into account in the interionic potential at any density, even though the departure from linearity can reach a few percent for the asymmetric mixtures in the region of weak degeneracy for the electron gas. A physical explanation of this behavior is proposed which is based on a simple additional length scale.

Chabrier, G.↗

Extreme ultraviolet spectroscopy of Capella

The X-ray active binary system Capella was observed with a moderate-resolution extreme ultraviolet spectrograph from 200 to 330 A. Two low-level emission features were detected. One most likely is geocoronal He II 304 A emission, while the other probably originates from the corona of Capella. The weak stellar emission at 304 A is in direct conflict with predictions of the intrinsic stellar He II flux based on standard scaling arguments but is consistent with the only previous observation of Capella in the EUV. The most plausible explanation for the lack of stellar 304 A emission is a warm wind from the active G0 III star.

Green, James C.↗

Localization of Ad-Hoc Lunar Constellations in Communication Failure Modes for Distributed Spacecraft Autonomy

As lunar missions increase in complexity inspired by NASA’s Artemis Program, they will require reliable and sufficient capability of the Position, Navigation, and Timing (PNT) system to support their scientific objectives. In addition, NASA's Commercial Lunar Payload Services (CLPS) program initiates the proliferation of public and private exploration partnerships using small satellites from commercial and private organizations, expanding traditionally confined low Earth orbit to be used for missions beyond geosynchronous orbit (Zucherman et al., 2022). Therefore, the Lunar PNT system is also required to provide navigation services compatible with the smaller platforms being sent by the public and private sectors, like CubeSats. However, traditional approaches to deep space missions’ navigation based on ground radio facilities have difficulties in providing sufficient support for the increasing number of users and communication at a distance from the Earth (Kaplev et al., 2022). In particular, the existing Lunar navigation technologies such as weak signal global positioning system (GPS) and deep space network (DSN) are not able to ensure operations of the upcoming small-scale Lunar missions due to their limitations in localization performance as well as capacity aspects. Another way to provide Lunar PNT service is to create a dedicated Lunar global navigation satellite system (GNSS) constellation, like GNSS systems on Earth. Space agencies like NASA, ESA, and JAXA are now developing the lunar communications relay and navigation systems (LCRNS) and Lunar navigation satellite systems (LNSS). In their systems, satellites will be deployed in moon orbits to provide the communication, positioning, navigation, and timing (CPNT) service at the lunar south pole region where the Artemis base camp will be expected (Murata et al., 2022). Meanwhile, common challenges considered in lunar PNT research arise from poor geometry of the terrestrial GNSS satellites when seen from the lunar user, highly perturbed lunar orbits, and limitations in power, size, and cost of the equipment on lunar satellites (Iiyama et al., 2023). It is also not clear if there will be enough Lunar users to support the cost and resources this would require as the Low-cost surface missions may not be able to support the large power, mass, and weight requirements that these navigation solutions entail (Niemoeller et al., 2022). As an alternative, existing Lunar science and exploration assets could be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions (Hagenau et al., 2021). Introducing the non-dedicated and ad-hoc Lunar navigation constellation gives a way to provide PNT services on-demand. The non-dedicated swarm assets of Lunar constellations are designed to localize themselves with minimal interaction with Earth by adding cooperative autonomous localization to lunar missions, freeing up valuable bandwidth and ground segment resources. An autonomous localization of Lunar constellations is based on the concept of the decentralized PNT system with a distributed extended Kalman filter (DEKF) approach to state estimation for minimal onboard operating costs. In the distributed data processing algorithm, computation is broken down and assigned to each satellite, resulting in a considerably decreased computational amount while maintaining the accuracy of the orbit ephemeris and clock offsets as the result of centralized data processing (Wen et al., 2019). The DEKF requires spacecraft to perform two-way ranging operations with each other to communicate simultaneously, leveraging neighbor two-way intersatellite link (ISL) measurements such as pseudoranges to, and relative velocities between, visible satellites as sensor values (Frank et al., 2021). The Lunar autonomous PNT simulation (LAPS) demonstrated the feasibility of orbital asset localization among ad-hoc Lunar small-sat constellations based on the DEKF in Hagenau et al. (2021) and evaluated the matching algorithm proposed by Frank et al. (2021) in scheduling position estimation updates. In previous papers, all assets and measurements are assumed to be always available without consideration of the impact of intermittent and permanent communication failure. This study presents localization performance with increasing levels of network degradation for swarm assets and users to demonstrate the robustness of the decentralized Lunar PNT service in more realistic scenarios. Main issues arising from communication failure include spacecraft permanent or transient loss, antenna failures, message delays, etc. We tested four possible reasons for network degradation for 7 days in 21 satellites frozen with an altitude of 5500 km, evenly spaced around 3 circular, 40 inclination orbital planes where each spacecraft has two directional antennas. As anchor nodes with an independent estimate of their position are required in the DEKF approach, two ground nodes in each pole and one node in the gateway were implemented in the simulation. First, the most probable failure scenario involves the loss of a single spacecraft due to solar interference and technical malfunctions of the assets. Losing the availability of a single spacecraft means losing the two-way ISL measurement of the asset in the DEKF update. In order to provide the best possible quality of PNT service with limited time and resources, the distributed Lunar constellations must schedule the communication activities. The scheduler leverages mixed-integer linear programming (MILP) for the coordination and scheduling of the desired “as-needed” localization service (Niemoeller et al., 2022). We assume the scheduler has completely excluded the spacecraft information before the DEKF update in the failure scenario. When a random spacecraft has been turned off at a specific time, the robustness of the autonomous Lunar PNT system is evaluated. The simulation results give an 11.5% degradation in median position accuracy compared to the idealized performance excluding the asset loss. Second, a large number of assets may vanish due to major hardware problems or meteor strikes around the moon. A multiple spacecraft loss can degrade the localization performance very fast by losing the communication ability to do cross-plane measurements and in-plane measurements in a 3-plane constellation. When the matching-based scheduler is aware of ISL availability, we investigate a large number of in-plane and cross-plane asset vanishments both in close proximity and equally spaced throughout the orbital plane. According to the simulations, the loss of in-plane measurements gives 40.2% degradation while cross-plane measurements degrade 50.5% of asset localization performance among available assets. Therefore, it is concluded that cross-plane measurements are more important in improving the position estimation accuracy. Third, spacecraft failure information can be lost due to the internal message delay, resulting in the DEKF update scheduler to solve the matching problem with unavailable assets. The DEKF update cycle is comprised of network setup, communication, and computations where a global broadcast network and a 2-way ISL network setup take 6 minutes in total (Frank et al., 2021). Once the broadcast network successfully transmits and receives information, a random spacecraft may lose its availability right before solving the matching problem. This means the matching solution is no longer optimal, resulting in degradation in the localization performance. A numerical assessment shows the matching-based scheduler with knowing failure holds 11.5% of position accuracy degradation, whereas the scheduler without knowing failure gives 34% degraded localization performance without asset loss. Fourth, a transient loss of a single or multiple spacecraft may occur due to their antenna outages. After losing the two-way ISL availability for a few DEKF update cycles, the availability of spacecraft can easily be recovered as their states have been independently updated using measurements from anchor nodes. It is likely that the longer failure will result in worse localization performance. We have tested the transient failure of a random single asset for 30 min in the simulation, which is losing 3 update cycles in the DEKF system. From the simulation results, the position accuracy has been degraded to 4.84% which is better than the degraded localization performance of 11.5% from the permanent loss scenario among available assets. In conclusion, the autonomous Lunar PNT system based on the DEKF approach shows the ability to maintain resilience and robustness in the possible communication failure scenarios, ensuring that localization accuracy is preserved across various network degradation and outages. Future studies on investigating user localization performance near the South Pole and the broadcast network system will be continued in the following months.

Yeji Kim↗

COMMUNITY AND INFRASTRUCTURE ADAPTION TO CLIMATE CHANGE

(SF-23-046) The CIACC tool labels a large number of documents using a weak-supervision-based NLP technique that makes use of semantic similarity between categories and documents. This method allows the labeling process to scale across multiple GPU's in a computing cluster at the Argonne Leadership Computing Facility.

Verner, Duane↗

First Detection of the GI-type of Intrinsic Alignments of Galaxies Using the Self-calibration Method in a Photometric Galaxy Survey

Weak gravitational lensing is one of the most promising cosmological probes to constrain dark matter, dark energy, and the nature of gravity at cosmic scales. Intrinsic alignments (IAs) of galaxies have been recognized as one of the most serious systematic effects facing gravitational lensing. Such alignments must be isolated and removed to obtain a pure lensing signal. Furthermore, the alignments are related to the processes of galaxy formation, so their extracted signal can help in understanding such formation processes and improving their theoretical modeling. We report in this Letter the first detection of the gravitational shear–intrinsic shape (GI) correlation and the intrinsic shape–galaxy density (Ig) correlation using the self-calibration method in a photometric redshift survey. These direct measurements are made from the KiDS-450 photometric galaxy survey with a significance of 3.65σ in the third bin for the Ig correlation, and 3.51σ for the GI cross-correlation between the third and fourth bins. The self-calibration method uses the information available from photometric surveys without needing to specify an IA model and will play an important role in validating IA models and IA mitigation in future surveys such as the Rubin Observatory Legacy Survey of Space and Time, Euclid, and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗

Interpreting Sunyaev–Zel’dovich observations with MillenniumTNG: mass and environment scaling relations

ABSTRACT Sunyaev–Zel’dovich (SZ) measurements can dramatically improve our understanding of the intergalactic medium and the role of feedback processes in galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, τ, and the Compton-y parameter, Y, receive large line-of-sight contributions that can be removed effectively by applying a compensated aperture photometry filter. In contrast with other τ probes (e.g. X-rays and fast radio bursts), the kinematic SZ-inferred τ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line of sight. Additionally, we perform fits to the Y–M and τ–M scaling relations and report best-fitting parameters adopting the smoothly broken power law formalism. We note that subgrid physics modelling can broaden the error bar on these by 30 per cent for intermediate-mass haloes (${\sim }10^{13} \, {\rm M}_{\odot }$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than haloes in observations, which can bias the inferred profiles by ∼20 per cent for L* galaxies.

Hadzhiyska, Boryana (ORCID:0000000223123121)↗

Non-Gaussianity in the weak lensing correlation function likelihood – implications for cosmological parameter biases

ABSTRACT We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in 1D marginal distributions of shear two-point correlation functions in simulated weak lensing data. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in 1D posteriors of Ωm and σ8 due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data (100 deg2). We also present a systematic approach to constructing approximate multivariate likelihoods with 1D parametric functions by assuming independence or more flexible non-parametric multivariate methods after decorrelating the data points using principal component analysis (PCA). While the use of PCA does not modify the non-Gaussianity of the multivariate likelihood, we find empirically that the 1D marginal sampling distributions of the PCA components exhibit less skewness and kurtosis than the original shear correlation functions. Modelling the likelihood with marginal parametric functions based on the assumption of independence between PCA components thus gives a lower limit for the biases. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of ∼5.

79 ASTRONOMY AND ASTROPHYSICS↗

Joining and Assembly of Bulk Metallic Glass Composites Through Capacitive Discharge

Bulk metallic glasses (BMGs), a class of amorphous metals defined as having a thickness greater than 1 mm, are being broadly investigated by NASA for use in spacecraft hardware. Their unique properties, attained from their non-crystalline structure, motivate several game-changing aerospace applications. BMGs have low melting temperatures so they can be cheaply and repeatedly cast into complex net shapes, such as mirrors or electronic casings. They are extremely strong and wear-resistant, which motivates their use in gears and bearings. Amorphous metal coatings are hard, corrosion-resistant, and have high reflectivity. BMG composites, reinforced with soft second phases, can be fabricated into energy-absorbing cellular panels for orbital debris shielding. One limitation of BMG materials is their inability to be welded, bonded, brazed, or fastened in a convenient method to form larger structures. Cellular structures (which can be classified as trusses, foams, honeycombs, egg boxes, etc.) are useful for many NASA, commercial, and military aerospace applications, including low-density paneling and shields. Although conventional cellular structures exhibit high specific strength, their porous structures make them challenging to fabricate. In particular, metal cellular structures are extremely difficult to fabricate due to their high processing temperatures. Aluminum honeycomb sandwich panels, for example, are used widely as spacecraft shields due to their low density and ease of fabrication, but suffer from low strength. A desirable metal cellular structure is one with high strength, combined with low density and simple fabrication. The thermoplastic joining process described here allows for the fabrication of monolithic BMG truss-like structures that are 90% porous and have no heat-affected zone, weld, bond, or braze. This is accomplished by welding the nodes of stacked BMG composite panels using a localized capacitor discharge, forming a single monolithic structure. This removes many complicated and costly fabrication steps. Moreover, the cellular structures detailed in this work are among the highest- strength and most energy-absorbent materials known. This implies that a fabricated structure made from these materials would have unequaled mechanical properties compared to other metal foams or trusses. The process works by taking advantage of the electrical properties of the matrix material in the metal-matrix composite, which in this case is a metallic glass. Due to the random nanoscale arrangement of atoms (without any grain boundaries), the matrix glass exhibits a near-constant electrical resistivity as a function of temperature. By placing the composite panels between two copper electrode plates and discharging a capacitor, the entire matrix of the panel can be heated to approximately 700 C in 10 milliseconds, which is above the alloy s solidus but below the liquidus. By designing the geometry of the panels into the shape of an egg box, the electrical discharge localizes only in the tips of each pyramidal cell. By applying a forging load during discharge, the nodes of the panels can be fused together into a single piece, which then dissipates heat through radiation back into a glassy state. This means that two panels can be metallurgically fused into one panel with no heat-affected zone, creating a seamless connection between panels. During the process, the soft metal particles (dendrites) that are uniformly distributed in the glassy matrix to increase the toughness are completely unaffected by the thermoplastic joining. The novelty is that a truss (or foam-like) structure can be formed with excellent energy- absorbing capabilities without the need for machining. The technique allows for large-scale fabrication of panels, well-suited for spacecraft shields or military vehicle door panels. Crystalline metal cellular structures cannot be fabricated using the thermoplastic joining technique described here. If metal panels were te assembled into a cellular structure, they would either have to be welded, brazed, bonded, or fastened together, creating a weak spot in the structure at each connection. Welded parts require a welding material to be added to the joint and exhibit a soft and weak heat-affected zone. Brazing and bonding do not form a metallurgical joint and thus exhibit low strengths, especially when the panels are pulled apart and fasteners require high-stress-concentration holes to be drilled. No equivalent rapid heating method exists for assembling metal panels together into cellular structures, and thus, those parts must be foamed, machined, or investment cast if they are to form a monolithic structure. If the crystalline panels were to be joined using capacitive discharge, as with a spot welder, their bond would be very weak, and the panels would have to be extremely thin. In contrast, the strength of joined BMG parts has been demonstrated to have strength comparable to the parent material. This technique opens up the possibility of using large-scale BMG hardware in spacecraft, military, or commercial applications.

Hofmann, Douglas C.↗

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗