Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Tropospheric aerosols over the western North Atlantic Ocean during the winter and summer deployments of ACTIVATE 2020: life cycle, transport, and distribution

The Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) is a NASA mission to characterize aerosol–cloud interactions over the western North Atlantic Ocean (WNAO). Such characterization requires understanding of life cycle, composition, transport pathways, and distribution of aerosols over the WNAO. This study uses the GEOS-Chem model to simulate aerosol distributions and properties that are evaluated against aircraft, ground-based, and satellite observations during the winter and summer field deployments in 2020 of ACTIVATE. Transport in the boundary layer (BL) behind cold fronts was a major mechanism for the North American continental outflow of pollution to the WNAO in winter. Turbulent mixing was the main driver for the upward transport of sea salt within and ventilation out of BL in winter. The BL aerosol composition was dominated by sea salt, which increased in the summer, followed by organics and sulfate. Aircraft in situ aerosol measurements provided useful constraints on wet scavenging in GEOS-Chem. The model generally captured observed features such as continental outflow, land–ocean gradient, and mixing of anthropogenic aerosols with sea salt. Model sensitivity experiments with elevated smoke injection heights to the mid-troposphere (versus within BL) better reproduced observations of smoke aerosols from the western US wildfires over the WNAO in the summer. Model analysis suggests strong hygroscopic growth of sea salt particles and their seeding of marine BL clouds over the WNAO (< 35° N). Future modeling efforts should focus on improving parameterizations for aerosol wet scavenging, implementing realistic smoke injection heights, and applying high-resolution models that better resolve vertical transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The suitability of differentiable, physics-informed machine learning hydrologic models for ungauged regions and climate change impact assessment

As a genre of physics-informed machine learning, differentiable process-based hydrologic models (abbreviated as δ or delta models) with regionalized deep-network-based parameterization pipelines were recently shown to provide daily streamflow prediction performance closely approaching that of state-of-the-art long short-term memory (LSTM) deep networks. Meanwhile, δ models provide a full suite of diagnostic physical variables and guaranteed mass conservation. Here, we ran experiments to test (1) their ability to extrapolate to regions far from streamflow gauges and (2) their ability to make credible predictions of long-term (decadal-scale) change trends. We evaluated the models based on daily hydrograph metrics (Nash–Sutcliffe model efficiency coefficient, etc.) and predicted decadal streamflow trends. For prediction in ungauged basins (PUB; randomly sampled ungauged basins representing spatial interpolation), δ models either approached or surpassed the performance of LSTM in daily hydrograph metrics, depending on the meteorological forcing data used. They presented a comparable trend performance to LSTM for annual mean flow and high flow but worse trends for low flow. For prediction in ungauged regions (PUR; regional holdout test representing spatial extrapolation in a highly data-sparse scenario), δ models surpassed LSTM in daily hydrograph metrics, and their advantages in mean and high flow trends became prominent. In addition, an untrained variable, evapotranspiration, retained good seasonality even for extrapolated cases. The δ models' deep-network-based parameterization pipeline produced parameter fields that maintain remarkably stable spatial patterns even in highly data-scarce scenarios, which explains their robustness. Combined with their interpretability and ability to assimilate multi-source observations, the δ models are strong candidates for regional and global-scale hydrologic simulations and climate change impact assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Global and regional modeling of clouds and aerosols in the marine boundary layer during VOCALS: the VOCA intercomparison

A diverse collection of models are used to simulate the marine boundary layer in the southeast Pacific region during the period of the October–November 2008 VOCALS REx (VAMOS Ocean Cloud Atmosphere Land Study Regional Experiment) field campaign. Regional models simulate the period continuously in boundary-forced free-running mode, while global forecast models and GCMs (general circulation models) are run in forecast mode. The models are compared to extensive observations along a line at 20° S extending westward from the South American coast. Most of the models simulate cloud and aerosol characteristics and gradients across the region that are recognizably similar to observations, despite the complex interaction of processes involved in the problem, many of which are parameterized or poorly resolved. Some models simulate the regional low cloud cover well, though many models underestimate MBL (marine boundary layer) depth near the coast. Most models qualitatively simulate the observed offshore gradients of SO 2 , sulfate aerosol, CCN (cloud condensation nuclei) concentration in the MBL as well as differences in concentration between the MBL and the free troposphere. Most models also qualitatively capture the decrease in cloud droplet number away from the coast. However, there are large quantitative intermodel differences in both means and gradients of these quantities. Many models are able to represent episodic offshore increases in cloud droplet number and aerosol concentrations associated with periods of offshore flow. Most models underestimate CCN (at 0.1% supersaturation) in the MBL and free troposphere. The GCMs also have difficulty simulating coastal gradients in CCN and cloud droplet number concentration near the coast. The overall performance of the models demonstrates their potential utility in simulating aerosol–cloud interactions in the MBL, though quantitative estimation of aerosol–cloud interactions and aerosol indirect effects of MBL clouds with these models remains uncertain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A box model for representing estuarine physical processes in Earth system models

Appropriately treating riverine freshwater discharge into the oceans in Earth system models is a challenging problem. Commonly, the river runoff is discharged into the ocean models with zero salinity and arbitrarily distributed either horizontally or vertically over several grid cells. Those approaches entirely neglect estuarine physical processes that modify river inputs before they reach the open ocean. Here in order to realistically represent riverine freshwater inputs in Earth system models, a physically based Estuary Box Model (EBM) is developed to parameterize the mixing processes in estuaries. The EBM represents the estuary exchange circulation with a two-layer box structure. It takes as input the river volume flux from the land surface model and the subsurface salinity at the estuary mouth from the ocean model. It delivers the estuarine outflow salinity and net volume flux into and out of the estuary to the ocean model. An offline test of the EBM forced with observed conditions for the Columbia River system shows good agreement with observations of outflow salinity and high-resolution simulations of the exchange flow volume flux. To illustrate the practicality of use of the EBM in an Earth system model, the EBM is implemented for all coastal grid cells with river runoff in the Community Earth System Model (CESM). Compared to the standard version of CESM, which treats runoff as an augmentation to precipitation, the EBM increases sea surface salinity and reduces stratification near river mouths. The EBM also leads to significant regional and remote changes in CESM ocean surface salinities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Measurements of pore-scale flow through apertures

Pore-scale aperture effects on flow in pore networks was studied in the laboratory to provide a parameterization for use in transport models. Four cases were considered: regular and irregular pillar/pore alignment with and without an aperture. The velocity field of each case was measured and simulated, providing quantitatively comparable results. Two aperture effect parameterizations were considered: permeability and transmission. Permeability values varied by an order of magnitude between the cases with and without apertures. However, transmission did not correlate with permeability. Despite having much greater permeability the regular aperture case permitted less transmission than the regular case. Moreover, both irregular cases had greater transmission than the regular cases, a difference not supported by the permeabilities. Overall, these findings suggest that pore-scale aperture effects on flow though a pore-network may not be adequately captured by properties such as permeability for applications that are interested in determining particle transport volume and timing.

58 GEOSCIENCES↗
DOE OSTI2021

Integrating crop growth models with remote sensing for predicting biomass yield of sorghum

Abstract Plant phenotypes are often descriptive, rather than predictive of crop performance. As a result, extensive testing is required in plant breeding programmes to develop varieties aimed at performance in the target environments. Crop models can improve this testing regime by providing a predictive framework to (i) augment field phenotyping data and derive hard-to-measure phenotypes and (ii) estimate performance across geographical regions using historical weather data. The goal of this study was to parameterize the Agricultural Production Systems sIMulator (APSIM) crop growth models with remote-sensing and ground-reference data to predict variation in phenology and yield-related traits in 18 commercial grain and biomass sorghum hybrids. Genotype parameters for each hybrid were estimated using remote-sensing measurements combined with manual phenotyping in West Lafayette, IN, in 2018. The models were validated in hybrid performance trials in two additional seasons at that site and against yield trials conducted in Bushland, TX, between 2001 and 2018. These trials demonstrated that (i) maximum plant height, final dry biomass and radiation use efficiency (RUE) of photoperiod-sensitive and -insensitive forage sorghum hybrids tended to be higher than observed in grain sorghum, (ii) photoperiod-sensitive sorghum hybrids exhibited greater biomass production in longer growing environments and (iii) the parameterized and validated models perform well in above-ground biomass simulations across years and locations. Crop growth models that integrate remote-sensing data offer an efficient approach to parameterize larger plant breeding populations.

Yang, Kai-Wei↗
DOE OSTI2013

Uncertainties in nuclear transition matrix elements of neutrinoless ββ decay

To estimate the uncertainties associated with the nuclear transition matrix elements M{sup (K)} (K=0ν/0N) for the 0{sup +} → 0{sup +} transitions of electron and positron emitting modes of the neutrinoless ββ decay, a statistical analysis has been performed by calculating sets of eight (twelve) different nuclear transition matrix elements M{sup (K)} in the PHFB model by employing four different parameterizations of a Hamiltonian with pairing plus multipolar effective two-body interaction and two (three) different parameterizations of Jastrow short range correlations. The averages in conjunction with their standard deviations provide an estimate of the uncertainties associated the nuclear transition matrix elements M{sup (K)} calculated within the PHFB model, the maximum of which turn out to be 13% and 19% owing to the exchange of light and heavy Majorana neutrinos, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Validation Tests of Resuspension Models for a Finite and Infinite Site

Dose assessment for deposited radionuclides often requires estimates of air concentrations that are derived from measured soil concentrations. For this, dose assessors typically use literature resuspension values that while empirically based, can vary by orders of magnitude making it difficult to provide accurate dose estimates. Despite the complexities of the physical processes involved in resuspension, the models generally used for dose assessment are relatively simplistic and rarely are the models validated for a specific site, thus making prediction of air concentrations or airborne emissions highly uncertain. Additionally, the size of the contaminated area can have an impact on downwind concentrations, yet literature values do not account for smaller-sized contaminated sites adding additional uncertainty. To test resuspension models for soil-bound radionuclides at finite and infinite spatial scales, measurements of soil and air concentrations are made at 1) a location downwind of a former outfall where plutonium-239 was released into the environment ( a finite site), and 2) uncontaminated locations where regional air sampling provides background measurements of naturally occurring uranium-238 in sampled dust (an infinite site). Measured air concentrations were compared to those predicted using the Resuspension Factor Model and the Mass Loading Model. An Area Factor was applied to the smaller contaminated site to account for dilution of dust from the contaminated site with dust originating from offsite locations. Here, results show that when properly parameterized to site conditions, resuspension models can predict air concentrations to within a factor of ten.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Building spatial symmetries into parameterized quantum circuits for faster training

Abstract Practical success of quantum learning models hinges on having a suitable structure for the parameterized quantum circuit. Such structure is defined both by the types of gates employed and by the correlations of their parameters. While much research has been devoted to devising adequate gate-sets, typically respecting some symmetries of the problem, very little is known about how their parameters should be structured. In this work, we show that an ideal parameter structure naturally emerges when carefully considering spatial symmetries (i.e. the symmetries that are permutations of parts of the system under study). Namely, we consider the automorphism group of the problem Hamiltonian, leading us to develop a circuit construction that is equivariant under this symmetry group. The benefits of our novel circuitstructure, called ORB, are numerically probed in several ground-state problems. We find a consistent improvement (in terms of circuit depth, number of parameters required, and gradient magnitudes) compared to literature circuit constructions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Improving Precipitation Estimation Using Convolutional Neural Network

Precipitation process is generally considered to be poorly represented in numerical weather/climate models. Statistical downscaling (SD) methods, which relate precipitation with model resolved dynamics, often provide more accurate precipitation estimates compared to model's raw precipitation products. Herein we introduce the convolutional neural network model to foster this aspect of SD for daily precipitation prediction. Specifically, we restrict the predictors to the variables that are directly resolved by discretizing the atmospheric dynamics equations. In this sense, our model works as an alternative to the existing precipitation-related parameterization schemes for numerical precipitation estimation. We train the model to learn precipitation-related dynamical features from the surrounding dynamical fields by optimizing a hierarchical set of spatial convolution kernels. We test the model at 14 geogrid points across the contiguous United States. Results show that provided with enough data, precipitation estimates from the convolutional neural network model outperform the reanalysis precipitation products, as well as SD products using linear regression, nearest neighbor, random forest, or fully connected deep neural network. Evaluation for the test set suggests that the improvements can be seamlessly transferred to numerical weather modeling for improving precipitation prediction. Based on the default network, we examine the impact of the network architectures on model performance. Also, we offer simple visualization and analyzing approaches to interpret the models and their results. Our study contributes to the following two aspects: First, we offer a novel approach to enhance numerical precipitation estimation; second, the proposed model provides important implications for improving precipitation-related parameterization schemes using a data-driven approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model

Abstract Characteristics of, and fundamental differences between, the radiative‐convective equilibrium (RCE) climate states following the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) protocols in the Community Atmosphere Model version 5 (CAM5) and version 6 (CAM6) are presented. This paper explores the characteristics of clouds, moisture, precipitation and circulation in the RCE state, as well as the tropical response to surface warming, in CAM5 and CAM6 with different parameterizations. Overall, CAM5 simulates higher precipitation rates that result in larger global average precipitation, despite lower outgoing longwave radiation compared to CAM6. Differences in the structure of clouds, particularly the amount and vertical location of cloud liquid, exist between the CAM versions and can, in part, be related to distinct representations of shallow convection and boundary layer processes. Both CAM5 and CAM6 simulate similar peaks in cloud fraction, relative humidity, and cloud ice, linked to the usage of a similar deep convection parameterization. These anvil clouds rise and decrease in extent in response to surface warming. More generally, extreme precipitation, aggregation of convection, and climate sensitivity increase with warming in both CAM5 and CAM6. This analysis provides a benchmark for future studies that explore clouds, convection, and climate in CAM with the RCEMIP protocols now available in the Community Earth System Model. These results are discussed within the context of realistic climate simulations using CAM5 and CAM6, highlighting the usefulness of a hierarchical modeling approach to understanding model and parameterization sensitivities to inform model development efforts.

Reed, Kevin A.↗
DOE OSTI2025

Reduce Revenue Versus Increase Expenditure: Fires and Plant Invasion Drive Soil Carbon Loss With Different Mechanisms in a Mediterranean Shrubland

Fires and plant invasions threaten Mediterranean ecosystems substantially, particularly in the context of changing climate. Our study utilized a data-model integration approach to assess the response of soil organic carbon (SOC) to fires and plant invasion under three Shared Socio-Economic Pathway (SSP) scenarios (SSP1-26, SSP2-45, and SSP5-85). We parameterized the CLM-Microbe model and then investigated the individual and interactive impacts of fires and plant invasion on soil C by comparing factorial simulations of initialization (fire/no wildfire in 2021), fire module on/off, and with and without plant invasion during 2023–2100 in a Mediterranean ecosystem. The simulations indicated a marked C loss due to the 2021 wildfire, projected fires, and plant invasion across all future climate scenarios. Specifically, the 2021 wildfire, projected fires, and plant invasion reduced the SOC (0–30 cm) by 0.12, 0.26, and 0.15 kg C m −2 under SSP1-26, 0.12, 0.30, and 0.12 kg C m −2 under SSP2-45, and 0.12, 0.24, and 0.13 kg C m −2 under SSP5-85, respectively. However, fires and plant invasion decreased SOC through distinct mechanisms. The effects of the 2021 wildfire occurred due to its negative legacy on soil microbial community and, thus, litter accumulation, suppressing the formation of soil carbon via decomposition. Influences of projected fires happen via consuming fuel and suppressing carbon input to soils. In contrast, the impacts of plant invasions were due to enhanced microbial respiration, leading to C loss. In conclusion, these findings emphasize the need for tailored C sequestration strategies considering the disparate effects of fires and plant invasions in the Mediterranean climate.

microbe↗
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↗
DOE OSTI2018

Shallow-to-Deep Transition of Continental Moist Convection: Cold Pools, Surface Fluxes, and Mesoscale Organization

Large-eddy simulation is used to investigate the effects of cold pools driven by rain evaporation on the shallow-to-deep convection transition over land. The physically consistent methodologies are developed to obtain a time-dependent reference ensemble without cold pools and to apply interactive surface heat fluxes without modeling of surface energy and water budgets. Three different simulation ensembles are contrasted. The reference ensemble, in the spirit of one-dimensional single-column models, eliminates cold pools by horizontally homogenizing negative buoyancy production due to rain evaporation. The additional ensembles complement the reference cold-pool-free ensemble by including cold pools and by applying either interactive or prescribed surface fluxes. Contrasting these ensembles suggests possible improvements of convection parameterization in large-scale models of weather and climate. Without cold pools, the reference ensemble preserves key features of buoyancy-driven cellular convection associated with a field of convective plumes, as assumed in a typical convection parameterization. With cold pools, a significant enhancement of surface heat and moisture fluxes and about an hour delay of their daily maximum is simulated. Cold pools enhance near-surface temperature and moisture standard deviations as well as maxima of the near-surface updraft velocity. They also lead to the reduction of cloud lateral entrainment, deeper vertical development of the cloud layer, and a few-times-larger accumulated surface precipitation. Interactive surface fluxes provide a damping mechanism that noticeably suppresses all these effects. Perhaps surprisingly, cold pools do not significantly change the cloud-base convective mass flux that approximately follows the evolution of surface heat fluxes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (MEME-CCS)

This effort employs and adapts an existing, rigorous multiscale electricity modeling platform at the National Renewable Energy Laboratory (NREL) to evaluate carbon capture and sequestration (CCS) and negative emissions technologies (NET) from the ARPA-E FLECCS program. NREL's modeling platform includes the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model, which projects future electricity generation mixes at sub-state-level resolution that are downscaled to the unit-level to enable hourly, zonal or nodal electricity production cost modeling in the PLEXOS model. The resulting hourly price data from PLEXOS is provided to technology developers under the ARPA-E FLECCS program to enable technology-specific economic analysis. The ReEDS-PLEXOS modeling suite is well-established for examining electric sector futures with high renewable energy penetrations, energy storage, electrification, and distributed generation. The key advancement proposed herein utilizes collaboration with CCS experts at the University of Wyoming and the FLECCS teams to create innovative methods for representing CCS and NET in the ReEDS and PLEXOS models. Expanded technology options, new operational parameterizations, and detailed data defining CO2 capture, transportation, and storage systems are being integrated into these models to allow an unprecedented combination of scope and resolution for exploring the future of CCS and NET. The resulting capabilities will take advantage of high-performance computing resources to permit wide-ranging scenario analysis to assess CCS and NET deployment under alternative CO2 prices, fossil fuel prices, electricity demand growth, and other electric sector characteristics. Final outcomes will include publicly available hourly grid operation and price data for any U.S. region of interest along with open-access capacity expansion tools for evaluating CCS/NET systems. These products will help an emerging CCS/NET industry in the United States by providing economics-driven guidance to technology developers while informing policy and investment decisions in the public and private sectors.

air capture↗
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↗