SEARCH · Engineering Papers

Results for “velocity model building”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1979

A global geochemical model for the evolution of the mantle

It is proposed that the upper mantle transition region, 220 to 670 km, is composed of eclogite which has been derived from primitive mantle by about 20 percent partial melting and that this is the source and sink of oceanic crust. The remainder of the upper mantle is garnet peridotite which is the source of continental basalts and hotspot magmas. This region is enriched in incompatible elements by hydrous and CO2 rich metasomatic fluids which have depleted the underlying layers in the L.I.L. elements and L.R.E.E. The volatiles make this a low-velocity, high attenuation, low viscosity region. The eclogite layer is internally heated and its controls the convection pattern in the upper mantle. Plate tectonics is intermittent. The continental thermal anomaly at a depth of 150-220 km triggers kimberlite and carbonatite activity, alkali and flood basalt volcanism, vertical tectonics and continental breakup. Hot spots remain active after the continents leave and build the oceanic islands. Mantle plumes rise from a depth of about 220 km. Midocean ridge basalts rise from the depleted layer below this depth. Material from this layer can also be displaced upwards by subducted oceanic lithosphere to form back-arc basins.

Anderson, D. L.↗

DOE OSTI2024

Qualification of Mixing Criteria by Computational Fluid Dynamics Modeling for the 325 Building Stack Revision

Additional ventilation capacity has been designed for the 325 Building filtered exhaust stack system. The four (4) existing main facility exhaust fans are past the end of their useful life. The fans are being replaced to provide additional exhaust capacity for future growth and to provide a more robust system. Stack operations will involve running various fan combinations at any given time. The air monitoring system of the existing stack previously was found to comply with the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard. Full-scale, three-dimensional computational fluid dynamics (CFD) modeling was used to evaluate the modified four-fan system for compliance with the ANSI/HPS N13.1-2021 standard, which essentially is equivalent in mixing requirements to the ANSI/HPS N13.1-1999 standard (and ANSI/HPS N13.1-2011). The four mixing criteria evaluated are 1) flow angle, 2) velocity, 3) gas tracer, and 4) particle tracer. In addition to the evaluating the modified four-fan system a temporary single fan stack configuration was also evaluated with CFD modeling. The temporary stack is planned to be used while the four-fan system is being modified. Modeling of the modified four-fan design and temporary ductwork showed that flow angles, velocity uniformity, gas tracer, and particle tracer were acceptable.

42 ENGINEERING↗

NASA NTRS2004

Transition from Pool to Flow Boiling: The Effect of Reduced Gravity

Applications of boiling heat transfer in space can be found in the areas of thermal management, fluid handling and control, power systems, on-orbit storage and supply systems for cryogenic propellants and life support fluids, and for cooling of electronic packages for power systems associated with various instrumentation and control systems. Recent interest in exploration of Mars and other planets, and the concepts of in-situ resource utiliLation on Mars highlights the need to understand the effect of gravity on boiling heat transfer at gravity levels varying from 1>= g/g(sub e) >=10(exp -6). The objective of the proposed work was to develop a mechanistic understanding of nucleate boiling and critical heat flux under low and micro-gravity conditions when the velocity of the imposed flow is small. For pool boiling, the effect of reduced gravity is to stretch both the length scale as well as the time scale for the boiling process. At high flow velocities, the inertia of the liquid determines the time and the length scales and as such the gravitational acceleration plays little role. However, at low velocities and at low gravity levels both liquid inertia and buoyancy are of equal importance. At present, we have little understanding of the interacting roles of gravity and liquid inertia on the nucleate boiling process. Little data that has been reported in the literature does not have much practical value in that it can not serve as a basis for design of heat exchange components to be used in space. Both experimental and complete numerical simulations of the low velocity, low-gravity nucleate boiling process were carried out. A building block type of approach was used in that first the growth and detachment process of a single bubble and flow and heat transfer associated with the sliding motion of the bubble over the heater surface after detachment was studied. Liquid subcooling and flow velocity were varied parametrically. The experiments were conducted at 1 g(sub e), while varying the orientation of surface with respect to the gravity vector. In the laboratory experiments, holographic interferometry was used to obtain data on velocity and temperature fields associated with a bubble prior to, and after detachment and during sliding motion. A test rig for conducting experiments in the KC-135 was developed, but experiments could not be conducted due to the unavailability of the aircraft prior to completion of the project. Numerical simulations modeling the micro and macro regions of the bubble were carried out in three dimensions. The results of the experiments were used to validate analytical/numerical models.

Dhir, Vijay K.↗

NASA NTRS1973

A boundary layer approach to the analysis of atmospheric motion over a surface obstruction

A boundary layer approach for the solution of the flow field induced over a two-dimensional surface obstruction, such as a building or other man-made structure, is proposed. Adopting a specific geometry in the form of a semi-elliptical cylinder, the characteristics of atmospheric shear flow over a rough terrain are coupled with the well-known boundary layer equations. Two approaches are presented to incorporate the pressure field and boundary conditions which exist within the large viscous region over the obstruction. The first considers a region in the immediate vicinity of the body in which the pressure distribution and outer boundary condition on the velocity are computed from potential theory for flow over the elliptical cylinder. The second approach considers a much larger region of influence, extending from the surface to the undisturbed flow at large heights above the obstruction. Methods which appear to provide an improved theoretical model of the flow over the ellipse, such as a technique for simulating the effect of the seperation regions upstream and downstream of the body, are also presented.

Frost, W.↗

DOE OSTI2021

Exploring Sensitivity of ICF Outputs to Design Parameters in Experiments Using Machine Learning

We report building a sustainable burn platform in inertial confinement fusion (ICF) requires an understanding of the complex coupling of physical processes and the effects that key experimental design changes have on implosion performance. While simulation codes are used to model ICF implosions, incomplete physics and the need for approximations deteriorate their predictive capability. Identification of relationships between controllable design inputs and measurable outcomes can help guide the future design of experiments and development of simulation codes, which can potentially improve the accuracy of the computational models used to simulate ICF implosions. In this article, we leverage developments in machine learning (ML) and methods for ML feature importance/sensitivity analysis to identify complex relationships in ways that are difficult to process using expert judgment alone. We present work using random forest (RF) regression for prediction of yield, velocity, and other experimental outcomes given a suite of design parameters, along with an assessment of important relationships and uncertainties in the prediction model. We show that RF models are capable of learning and predicting on ICF experimental data with high accuracy, and we extract feature importance metrics that provide insight into the physical significance of different controllable design inputs for various ICF design configurations. These results can be used to augment expert intuition and simulation results for optimal design of future ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A Lagrangian Study of the Transition from Shallow to Deep Convection Using ASR Observations and Les Simulations (Final Report)

With the support from this grant, we have made the following accomplishments. We have found that a formula for fractional entrainment rate, a key aspect for parameterization moist convection, that works for both shallow and deep convection, which provide a promising direction for unified treatment of entrainment for both shallow and deep convection. We also examined the formulation of vertical velocity and showed that the effect of non-hydrostatic pressure and entrainment on vertical velocity are not constants, as assumed in past studies, but can be computed using a simple model for the cloudy updrafts. These results also apply to both deep and shallow convection, which opens the door for a more physically based parameterization of the updraft vertical velocity. We have also conducted extensive analysis of the RACORO aircraft data and have published the results in Wei et al. (2020). We have made extensive effort scrutinizing the data in detail and examining the various hypotheses with the aid of the LES simulations. This was the first application of the modified Paluch diagram on modern aircraft data and the first comprehensive analysis of mixing in clouds using real data (instead of only LES output) in a very long time. These analyses enhanced the science return of the ASR aircraft measurements. This grant also supported a study of the collisions among cold pools and an effort to build parameterizations of moist convection using tools developed in control theory.

42 ENGINEERING↗

NASA NTRS2005

Dual-Band Deramp Radar Design for Ocean Current Measurements

A mission has been proposed to remotely measure ocean surface currents and surface wind velocities. It will provide the highest resolution and repeat time of these measurements to date for ocean current models with scientific and societal applications. A ground-based experimental radar unit is needed for proof of concept. The proposed experiment set up is to mount the radar on an oil rig to imitate satellite data acquisition. This summer, I completed the radar design. The design employs chirp/deramp topology with simultaneous transmit/receive channels. These two properties allow large system bandwidth, extended sample time, close range imaging, and low sampling rate. The radar operates in the Ku and Ka microwave bands, at 13.5 and 35.5 GHz, respectively, with a system bandwidth of 300 MHz. I completed the radar frequency analysis and research on potential components and antenna configurations. Subsequent work is needed to procure components, as well as to build, test, and deploy the radar.

Haynes, Mark S.↗

NASA NTRS1976

Flare build-up study; Proceedings of the Workshop, Falmouth, Mass., September 8-11, 1975

The papers deal with problems which might be common to solar flares and earth's magnetosphere, problems associated with the initial phases of the flare phenomenon, the acceleration processes that seem to occur in flares and the magnetosphere, as well as the buildup and storage of flare energy in magnetic-field structures. Topics include the active role of magnetic fields in providing flare energy, current-sheet models of solar flares, the role of plasma turbulence in flare development, similarities and differences between magnetospheric substorms and solar flares, observations of magnetic merging in earth's magnetotail during magnetospheric substorms, evidence for magnetic-energy storage in coronal active regions, the possible role of transition-zone instabilities in preflare energy buildup, and flare energy storage and deposition. Other papers discuss energy release through the interaction of coronal magnetic fields, photospheric electric currents as a source of flare energy, magnetic-energy buildup in the solar atmosphere, magnetic and velocity fields in an active region, flare onset at meter wavelengths, laboratory experiments on field-line reconnection, key problems in auroral flare processes, and the solar-physics Shuttle/Spacelab program. Individual items are announced in this issue.

Svestka, Z.↗

DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗

DOE OSTI2024

LIGHT-bgcArgo-1.0: using synthetic float capabilities in E3SMv2 to assess spatiotemporal variability in ocean physics and biogeochemistry

Since their advent over 2 decades ago, autonomous Argo floats have revolutionized the field of oceanography, and, more recently, the addition of biogeochemical and biological sensors to these floats has greatly improved our understanding of carbon, nutrient, and oxygen cycling in the ocean. While Argo floats offer unprecedented horizontal, vertical, and temporal coverage of the global ocean, uncertainties remain about whether Argo sampling frequency and density capture the true spatiotemporal variability in physical, biogeochemical, and biological properties. As the true distributions of, e.g., temperature or oxygen are unknown, these uncertainties remain difficult to address with Argo floats alone. Numerical models with synthetic observing systems offer one potential avenue to address these uncertainties. Here, we implement synthetic biogeochemical Argo floats into the Energy Exascale Earth System Model version 2 (E3SMv2), which build on the Lagrangian In Situ Global High-Performance Particle Tracking (LIGHT) module in E3SMv2 (E3SMv2-LIGHT-bgcArgo-1.0). Since the synthetic floats sample the model fields at model run time, the end user defines the sampling protocol ahead of any model simulation, including the number and distribution of synthetic floats to be deployed, their sampling frequency, and the prognostic or diagnostic model fields to be sampled. Using a 6-year proof-of-concept simulation, we illustrate the utility of the synthetic floats in different case studies. In particular, we quantify the impact of (i) sampling density on the float-derived detection of deep-ocean change in temperature or oxygen and on float-derived estimates of phytoplankton phenology, (ii) sampling frequency and sea-ice cover on float trajectory lengths and hence float-derived estimates of current velocities, and (iii) short-term variability in ecosystem stressors on estimates of their seasonal variability.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1995

The mosaic structure of plasma bulk flows in the Earth's magnetotail

Moments of plasma distributions observed in the magnetotail vary with different time scales. In this paper we attempt to explain the observed variability on intermediate timescales of approximately 10-20 min that result from the simultaneous energization and spatial structuring of solar wind plasma in the distant magnetotail. These processes stimulate the formation of a system of spatially disjointed. highly accelerated filaments (beamlets) in the tail. We use the results from large-scale kinetic modeling of magnetotail formation from a plasma mantle source to calculate moments of ion distribution functions throughout the tail. Statistical restrictions related to the limited number of particles in our system naturally reduce the spatial resolution of our results, but we show that our model is valid on intermediate spatial scales Delta(x) x Delta(z) equal to approximately 1 R(sub E) x 1000 km. For these spatial scales the resulting pattern, which resembles a mosaic, appears to be quite variable. The complexity of the pattern is related to the spatial interference between beamlets accelerated at various locations within the distant tail which mirror in the strong near-Earth magnetic field. Global motion of the magnetotail results in the displacement of spacecraft with respect to this mosaic pattern and can produce variations in all of the moments (especially the x-component of the bulk velocity) on intermediate timescales. The results obtained enable us to view the magnetotail plasma as consisting of two different populations: a tailward-Earthward system of highly accelerated beamlets interfering with each other, and an energized quasithermal population which gradually builds as the Earth is approached. In the near-Earth tail, these populations merge into a hot quasi-isotropic ion population typical of the near-Earth plasma sheet. The transformation of plasma sheet boundary layer (PSBL) beam energy into central plasma sheet (CPS) quasi-thermal energy occurs in the absence of collisions or noise. This paper also clarifies the relationship between the global scale where an MHD description might be appropriate and the lower intermediate scales where MHD fails and large-scale kinetic theory should be used.

Ashour-Abdalla, M.↗

DOE OSTI2023

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↗

NASA NTRS1976

Atmospheric flow over two-dimensional bluff surface obstructions

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

Bitte, J.↗

NASA NTRS1954

Design and Test of Mixed-flow Impellers IV : Experimental Results for Impeller Models MFI-1 and MFI-2 with Changes in Blade Height

Modifications A and B of impeller model MFI-1 and A, B, and C of impeller model MFI-2 were investigated experimentally in an attempt to determine what allowance in blade height should be made for boundary layer and viscous losses in an impeller designed for isentropic compressible flow. A gradual increase in blade height was arbitrarily made from inlet to outlet in anticipation of a gradual build-up of boundary layer. Apparently there was a rapid build-up of boundary layer near the inlet in the experimental case rather than a gradual one. Therefore, the proper allowance for boundary layer cannot be described from the data obtained. Decreasing the pressure gradient along the shroud by reducing the blade height allowance apparently did little to increase the overall efficiency. At the design speed of 1400 feet per second, the overall adiabatic efficiency was increased from 0.83 for the MFI-1A to 0.85 for the MFI-1B with reduction in height; however, it is indicated from the theoretical velocity distribution and outlet surveys that the increase was due to a change from decelerating to accelerating flow along the hub rather than from any change along the shroud. It is further indicated that the consequences of a thickened or separated boundary layer depend not only on the design velocity gradients but also on the shape of the passage.

Hamrick, Joseph T↗

DOE OSTI2022

Host-region parameters for an adjustable model for crustal earthquakes to facilitate the implementation of the backbone approach to building ground-motion logic trees in probabilistic seismic hazard analysis

The backbone approach to constructing a ground-motion logic tree for probabilistic seismic hazard analysis (PSHA) can address shortcomings in the traditional approach of populating the branches with multiple existing, or potentially modified, ground-motion models (GMMs) by rendering more transparent the relationship between branch weights and the resulting distribution of predicted accelerations. To capture epistemic uncertainty in a tractable manner, there are benefits in building the logic tree through the application of successive adjustments for differences in source, path, and site characteristics between the host region of the selected backbone GMM and the target region for which the PSHA is being conducted. The implementation of this approach is facilitated by selecting a backbone GMM that is amenable to such host-to-target adjustments for individual source, path, and site characteristics. The NGA-West2 GMM of Chiou and Youngs (CY14) has been identified as a highly adaptable model for crustal seismicity that is well suited to such adjustments. Rather than using generic source, path, and site characteristics assumed appropriate for the host region, the final suite of adjusted GMMs for the target region will be better constrained if the host-region parameters are defined specifically on the basis of their compatibility with the CY14 backbone GMM. To this end, making use of a recently developed crustal shear-wave velocity profile consistent with CY14, we present an inversion of the model to estimate the key source and path parameters, namely the stress parameter and the anelastic attenuation. With these outputs, the effort in constructing a ground-motion logic tree for any PSHA dealing with crustal seismicity can be focused primarily on the estimation of the target-region characteristics and their associated uncertainties. The inversion procedure can also be adapted for any application in which different constraints might be relevant.

Stafford, Peter J.↗

NASA NTRS1997

Dynamics Questions Associated with the AERCam Sprint Free-Flyer

The International Space Station will require the development of small robotic vehicles for such tasks as external inspection, monitoring of extravehicular activities (EVA's) and station build-up, and providing additional lighting at EVA worksites. The Autonomous EVA Robotic Camera (AERCam) family of free-flyers is currently being developed at NASA Johnson Space Center to perform these functions; the first member of the family is the AERCam Sprint vehicle. Many interesting dynamical questions are associated with the Sprint free-flyer. For instance, the reaction of a vehicle which is nearly spherically symmetric (such as Sprint) to a stuck-on thruster is significantly more complicated than that obtained for an idealized, perfectly spherical, spacecraft model. In particular, the real spacecraft will experience a form of forced nutation, with convergence towards either its major or minor principal axis, depending on both the applied torque and the mass properties of the vehicle. Furthermore, the body-fixed jet force vector may have a significant component along this principal axis, so giving rise to a considerable net linear acceleration of the spacecraft. The large velocity that can result is very important, as it may lead to collision with the nearby Orbiter, and is completely overlooked in the idealized analysis. This report will firstly briefly describe the stuck-on thruster dynamics of the real vehicle, and outline how the small products of inertia of the spacecraft determine the time constants of the motion. Secondly, the dynamical effects of a failed-off jet on the Sprint free-flyer will be described in more detail, and compared with the stuck-on thruster case. This will help to show whether the two malfunctions should be dealt with differently in flight. Finally, the stuck-on thruster detection software (known as the uncommanded motion algorithm) that is proposed to be flown on the Sprint vehicle will be analyzed, and all possible perturbation sources that may tend to give rise to false stuck-on thruster alarms quantified. It will be shown that false alarms can be triggered by discrete events such as a failed-off thruster, a glancing collision with Orbiter structure, or thruster saturation. They can also be triggered by a combination of errors introduced by constantly present sources such as thrust level errors, thruster misalignment, inertial cross-coupling, angular rate sensor noise, and structural vibrations induced by thruster firings. Based on this fact, the original plan to have the spacecraft enter safe mode whenever the uncommanded motion algorithm is triggered does not appear to be advisable. A better approach seems to be to provide the Sprint pilot with a warning whenever uncommanded motion is detected, and then allow the crew to determine whether safing the vehicle is appropriate.

Williams, Trevor↗

DOE OSTI2024

Momentum power spectrum of SDSS galaxies by massE cosmic ruler: 2.1×improvement in measure of growth rate

ABSTRACT Peculiar motion of galaxies probes the structure growth in the universe. In this study, we employ the galaxy stellar mass-binding energy (massE) relation with only two nuisance parameters to build the largest peculiar-velocity (PV) catalogue to date, consisting of 229 890 ellipticals from the main galaxy sample (MGS) of the Sloan Digital Sky Survey (SDSS). We quantify the distribution of the massE-based distances in individual narrow redshift bins (dz = 0.005), and then estimate the PV of each galaxy based on its offset from the Gaussian mean of the distribution. As demonstrated with the Uchuu-SDSS mock data, the derived PV and momentum power spectra are insensitive to accurate calibration of the massE relation itself, enabling measurements out to a redshift of 0.2, well beyond the current limit of z = 0.1 using other galaxy scaling laws. We then measure the momentum power spectrum and demonstrate that it remains almost unchanged if varying significantly the redshift bin size within which the distance is measured, as well as the intercept and slope of the massE relation, respectively. By fitting the spectra using the perturbation theory model with four free parameters, fσ8 is constrained to fσ8 = 0.459$^{+0.068}_{-0.069}$ over Δz = 0.02–0.2, 0.416$^{+0.074}_{-0.076}$ over Δz = 0.02–0.1, and 0.526$^{+0.133}_{-0.148}$ over Δz = 0.1–0.2. The error of fσ8 is 2.1 times smaller than that by the redshift space distortion (RSD) of the same sample. A Fisher matrix forecast illustrates that the constraint on fσ8 from the massE-based PV can potentially exceed that from the stage-IV RSD in late universe (z<0.5).

Shi, Yong (ORCID:0000000286146275)↗

DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence 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 the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{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 sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved 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 Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗