Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,585 records · Page 88

Experimental validation of multiphysics model simulations of the thermal response of a cement clinker rotary kiln at laboratory scale

Abstract An increasing demand for buildings, transportation systems and civil infrastructure development has driven expansion of cement consumption world‐wide, producing a significant increase in related global energy demand. With approximately 7% of the world‐wide industrial energy consumption (10.7 exajoules [EJ]), the cement industry is the third most energy intensive industrial processes and a key component for concrete, the most consumed composite material in the global construction industry. In cement manufacturing, the cement kiln accounts for most of the energy consumption in the production process. As the heart of a cement plant, the cement kiln is where the kiln feed primarily containing calcium oxide (CaO), silica (SiO 2 ), alumina (Al 2 O 3 ), and iron (Fe 2 O 3 ) are thermally and chemically transformed into clinker minerals. The presented work developed a multiphysics model, designed and built a laboratory‐scale rotary cement clinker kiln, and produced cement clinker at laboratory‐scale. The model was developed to study the interaction between the various thermal, fluid dynamic and chemical interactions involved in the sintering process used to form Portland cement clinker in an effort to reduce energy use. The analytical model was validated through experimental testing using a unique laboratory‐scale rotary cement kiln developed during the investigation. Also demonstrated was the feasibility of producing clinker at laboratory scale. This modeling and lab scale tests were designed to better understand the clinker sintering process so that operational and quality decisions can be made to optimize energy consumption without compromising cement clinker quality. The computational fluid dynamics modeling was developed in COMSOL Multiphysics 6.0. The characteristics of the combustion fluid flow, concentration of species, temperature and heat transfer were studied for a turbulent flow of methane (CH 4 ) gas and oxygen (O 2 ). Theory suggests that heat transfer impacts the cement production process but the multiphysics model more accurately describes the convection, conduction, and radiant heat transfer in the kilning process and thus allows for a better understanding of the energy exchange driving the chemical reactions that produce Portland cement. Clinker minerals were formed because of appropriate burning conditions implemented during experimental model validation.

Tabares, Juan David↗

Predicting Power for a Scaled‐up Non‐Newtonian Biomass Slurry

Abstract High‐solids biomass slurries exhibit non‐Newtonian behavior with a yield stress and require high power input for mixing. The goals were to determine the effect of scale and geometry on power number P 0 , and estimate the power for mixing a pretreated biomass slurry in a 3.8 million L hydrolysis reactor of conventional design. A lab‐scale computational fluid dynamics model was validated against experimental data and then scaled up. A pitched‐blade turbine and A310 hydrofoil were tested for various geometric arrangements. Flow was transitional; laminar and turbulence models resulted in equivalent P 0 which increased with scale. The ratio of impeller diameter to tank diameter affected P 0 for both impellers, but impeller clearance to tank diameter affected P 0 only for the A310. At least 2 MW is required to operate at this scale.

Russ, David C.↗

Applying cumulative effects to strategically advance large-scale ecosystem restoration

Programs restoring large terrestrial-aquatic ecosystems are expected to increase in area and number over the next decade. Factors ensuring long-term effectiveness at the program scale differ from those at the individual restoration project scale, necessitating a new synthesis. Here we contextualize hydrological, geomorphological, and biological effects of large-scale ecosystem restoration activities within a cumulative-effects paradigm. Evidence of eight modes of cumulative effects benefiting species and ecosystems was elucidated by reviewing restoration efforts in the greater Florida Everglades, Gulf Coast, lower Columbia River and estuary, Puget Sound estuary, Missouri River, northeastern coastal states, San Francisco Bay and Sacramento Delta. Additive and nonlinear cumulative effects of interacting restoration projects propagated across landscape and regional geographic scales, and throughout trophic levels, affecting population structures. Applying a cumulative-effects paradigm to planning, implementing, and evaluating large-scale restoration will avoid countervailing effects and capitalize on positive feedback, advancing evidence-based programs to recover priority species and ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Implications from Monitoring Gopher Tortoises at Two Spatial Scales

A problem that conservation biologists face is how to monitor species, given both resource limitations and the inherent challenges of assessing long-term demographic processes. We assessed gopher tortoise (Gopherus polyphemus) abundance at a landscape scale and at the scale of 3 local populations within the Conecuh National Forest (CNF).USA between 1991 and 2017. Landscape-level data were collected 26 from line transect distance sampling arranged uniformly across the CNF and collected during a single season (2011); data for local populations were generated from long-term mark-recapture of individuals at three sites selected based on prior knowledge of high density at each. At a landscape scale, we estimated 5,242 (95% CI = 3,538-7,768) tortoises occurred across the approximately 34,000-ha forest, yielding a density of 0.14-0.32 tortoises/ ha. These low densities across the landscape suggest that, on average, management activities across the property have not allowed tortoise populations to retain social structure needed for long-term persistence. The three local populations, however, contained 25-60 individuals and densities of 1.9–6.9 tortoises/ha. Over the study period, populations at two sites were stable and the third experienced significant population growth. Mean annual survival of individuals was 0.89 and invariant across size classes. Altogether, line transect distance sampling is important for assessing landscape-scale abundance of tortoises but may fail to detect local clusters of high-density sites important for population persistence. Our mark-recapture efforts at the local scale revealed that small populations on these high-density sites can exhibit long-term stability or growth even though they do not meet current established criteria for viability. Improved models that incorporate immigration and emigration and better reflect dynamics of peripheral populations would assist in determining how such populations best contribute to species recovery and regional conservation targets.

59 BASIC BIOLOGICAL SCIENCES↗

A micromorphic filter for determining stress and deformation measures from direct numerical simulations of lower length scale behavior

A micromorphic filter is presented for the extraction of stress and deformation measures from underlying direct numerical simulation (DNS) of microstructural mechanical response. The filter is consistent with the micromorphic continuum theory of Eringen and Suhubi (1964), and allows the interrogation of a DNS to generate higher order constitutive models. The filter adopts aspects of the overlap coupling techniques for atomistic-continuum scale-bridging in order to update the nodal degrees of freedom of the finite-element-based filter. Here, two methods for the extraction of quantities of interest from the underlying DNS are presented. The first is a least squares fitting technique which is appropriate for finite-element-style macro-scale domains, and the second is a variationally based approach which generalizes the method to any macro-scale domain which has a variational form. No restrictions on the form of the micro-scale simulation are inherent to the approach, provided that the required quantities (stress, position, etc.) can be provided. The filter is applied to a simple homogeneous medium undergoing homogeneous deformation to demonstrate the recovery of a classical continuum response, what we call the null case. We then apply the framework to the same homogeneous medium with a central and edge penny-shaped crack to show micromorphic effects. We demonstrate that varying the size of the filter causes the micromorphic behaviors to arise and then dissipate as the filter becomes smaller than the crack's length scale. Additionally we show the effect of filter size on a chiral element network along with a plate perforated with evenly spaced and perturbed holes through the thickness. The micromorphic effects are shown to be present for the chiral elements under any loading but become most pronounced when a single element is located within each averaging domain. Micromorphic effects become pronounced for the plate upon the onset of localization.

97 MATHEMATICS AND COMPUTING↗

Light dark matter, naturalness, and the radiative origin of the electroweak scale

We study classically scale invariant models in which the Standard Model Higgs mass term is replaced in the Lagrangian by a Higgs portal coupling to a complex scalar field of a dark sector. We focus on models that are weakly coupled with the quartic scalar couplings nearly vanishing at the Planck scale. The dark sector contains fermions and scalars charged under dark SU(2) × U(1) gauge interactions. Radiative breaking of the dark gauge group triggers electroweak symmetry breaking through the Higgs portal coupling. Requiring both a Higgs boson mass of 125.5 GeV and stability of the Higgs potential up to the Planck scale implies that the radiative breaking of the dark gauge group occurs at the TeV scale. We present a particular model which features a long-range abelian dark force. The dominant dark matter component is neutral dark fermions, with the correct thermal relic abundance, and in reach of future direct detection experiments. The model also has lighter stable dark fermions charged under the dark force, with observable effects on galactic-scale structure. Collider signatures include a dark sector scalar boson with mass ≲ 250 GeV that decays through mixing with the Higgs boson, and can be detected at the LHC. As a result, the Higgs boson, as well as the new scalar, may have significant invisible decays into dark sector particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Higgs coupling measurements and the scale of new physics

A primary goal of present and future colliders is measuring the Higgs couplings to Standard Model (SM) particles. Any observed deviation from the SM predictions for these couplings is a sign of new physics whose energy scale can be bounded from above by requiring tree-level unitarity. In this paper, we extend previous work on unitarity bounds from the Higgs cubic coupling to Higgs couplings to vector bosons and top quarks. We find that HL-LHC measurements of these couplings compatible with current experimental bounds may point to a scale that can be explored at the HL-LHC or a next-generation collider. Our approach is completely model-independent: we assume only that there are no light degrees of freedom below the scale of new physics, and allow arbitrary values for the infinitely many couplings beyond the SM as long as they are in agreement with current measurements. We also extend and clarify the methodology of this analysis, and show that if the scale of new physics is above the TeV scale, then the deviations can be described by the leading higher-dimension gauge invariant operator, as in the SM effective field theory.

43 PARTICLE ACCELERATORS↗

A Two-length Scale Turbulence Model for Single-phase Multi-fluid Mixing

A two-length scale, second moment turbulence model (Reynolds averaged Navier-Stokes, RANS) is proposed to capture a wide variety of single-phase flows, spanning from incompressible flows with single fluids and mixtures of different density fluids (variable density flows) to flows over shock waves. The two-length scale model was developed to address an inconsistency present in the single-length scale models, e.g. the inability to match both variable density homogeneous Rayleigh-Taylor turbulence and Rayleigh-Taylor induced turbulence, as well as the inability to match both homogeneous shear and free shear flows. The two-length scale model focuses on separating the decay and transport length scales, as the two physical processes are generally different in inhomogeneous turbulence. This allows reasonable comparisons with statistics and spreading rates over such a wide range of turbulent flows using a common set of model coefficients. The specific canonical flows considered for calibrating the model include homogeneous shear, single-phase incompressible shear driven turbulence, variable density homogeneous Rayleigh-Taylor turbulence, Rayleigh-Taylor induced turbulence, and shocked isotropic turbulence. The second moment model shows to compare reasonably well with direct numerical simulations (DNS), experiments, and theory in most cases. The model was then applied to variable density shear layer and shock tube data and shows to be in reasonable agreement with DNS and experiments. Additionally, the importance of using DNS to calibrate and assess RANS type turbulence models is highlighted.

04 OIL SHALES AND TAR SANDS↗

Effect of Scale Spallation During Coiling on the Electrochemical and Pickling Behavior of a Hot-Rolled Dual-Phase Steel

Tertiary oxide scale consists of magnetite and hematite, and in dual-phase steel strips were observed after the exit of final rolling stand in a thin slab casting rolling (TSCR) mill. Microscopic observation confirmed that these scales cracked and spalled off from many places of surface during coiling causing a non-uniform oxide layer on the steel surface. The present work discusses the synergistic effect of such non-uniform oxide layer and run-out-table cooling water on the electrochemical behavior of the dual-phase steel. The electrochemical response was assessed with the help of dynamic and linear polarization tests considering cooling water as electrolyte. The presence of non-uniform scale produced high cathodic (scaled zone) to anodic (spalled-off zone) area ratio that eventually led to extensive pitting corrosion in the present steel. Furthermore, the pickling behavior of the steel was also studied in detail. It was observed that the steel surface, with spalled-off area of about one-tenth of total surface area, corroded preferentially 3-4 times faster at spalled-off area than that of the uniformly pickled surface. The characterization study was performed using scanning electron microscopy, confocal laser 3D scanning microscopy, x-ray diffraction and Raman spectroscopy in order to establish the reasons for pitting corrosion due to the scale spallation.

36 MATERIALS SCIENCE↗

The Effect of CO 2 Pressure on Chromia Scale Microstructure at 750°C

To understand and model performance in supercritical CO 2 (sCO 2 ) for high-efficiency, concentrating solar power (CSP) and fossil energy power cycles, reaction rates are compared at 750°C in 0.1 MPa CO 2 and 30 MPa sCO 2 as well as laboratory air as a baseline on structural materials such as Ni-based alloy 625. Due to the thin reaction products formed even after 5000 h, scanning transmission electron microscopy was used to study the Cr-rich surface oxide scale. The scales formed in CO 2 and sCO 2 had a much finer grain size with more voids observed in CO 2 . However, the observations on alloy 625 were complicated by Mo and Nb-rich precipitates in the adjacent substrate and Al internal oxidation. To simplify the system, a binary Ni-22Cr alloy was exposed for 1000 h in similar environments. After exposure in sCO 2 , there was an indication of carbon segregation detected on the Cr 2 O 3 grain boundaries. After exposure in air, metallic Ni precipitates were observed in the scale that were not observed in the scale formed on alloy 625. The scale formed in air on a second Ni-22Cr model alloy with Mn and Si additions did not contain Ni precipitates, suggesting caution when drawing conclusions from model alloys.

36 MATERIALS SCIENCE↗

On modeling subgrid-scale macro-structures in narrow twisted channels

Porosity-based subgrid topography models often fail to capture the effects of subgrid-scale topographic features in the interior of coarse grid cells. Existing approaches that modify bottom roughness or a drag coefficient are inadequate for macro-structures (large emergent or submerged obstacles) in subgrid-scale narrow twisted channels. Such structures partially block the cross-sectional area and provide enhanced topographic dissipation – effects that are not well represented by a drag coefficient that scales on a coarse-grid cell-averaged velocity and the cell volume. The relative alignment between mesh and flow further complicates this problem as it makes the subgrid model sensitive to mesh design. In the present study, three new approaches for simulating subgrid-scale macro-structures in narrow channels are proposed. The interior partial-blocking effect of structures is modeled as reduction of grid face-area. The sheltering of flow volumes around obstacles, which leads to topographic dissipation, is modeled by reducing the cell volume in the momentum equation (only). A mesh-shift procedure is designed to optimize mesh alignment for identifiable subgrid features. Combining the three subgrid methods improves the approximation of surface elevation and in-channel flow rate with a coarse-grid model. Tests are conducted for channelized flow using both synthetic domains and real marsh topography. Finally, the new methods reduce the overall mesh dependency of the subgrid model and provides stronger physical connection between effects of macro-structures and their geometry at coarse grid scales.

42 ENGINEERING↗

Spaced out: An economic framework to explore the impacts of PV panel spacing on large-scale farming in Colorado

CONTEXT Agrivoltaic systems co-locate solar technologies with agricultural operations on an integrated plot of land and potentially provide benefits to both energy and agricultural systems. To date, large-scale (>5-MW) agrivoltaic projects in the United States have been limited to grazing and ecovoltaic applications, raising questions about the impact and scalability of agrivoltaic crop systems. Many agrivoltaic designs raise the height of the solar panels to accommodate agricultural practices while keeping energy density high. However, raising the panels results in increased photovoltaic (PV) development costs, which often are higher than the economic returns of crop production underneath the panels. This leads to unfavorable project economics and the need for other agrivoltaic solutions than raising panels. OBJECTIVE To explore other solutions, we perform an initial feasibility analysis for an agrivoltaic solution that can integrate with large-scale farming practices by increasing the row spacing in between panels. Increased PV row spacing is a low-cost approach for scaling agrivoltaics to accommodate crop production and this spacing can be tailored to required crop equipment for different regions. Increasing row spacing will reduce the power density (PV installed per acre), but in areas that are not land limited, these agrivoltaic designs could be economically feasible. Our analysis establishes a framework for a feasibility analysis for where and with what crops spaced out panel agrivoltaic solutions might be economical. METHODS Using a case study for large-scale agriculture crops in Colorado, we establish a framework for wide-row agrivoltaic economic feasibility analysis. We utilized the System Advisor Model to calculate technoeconomic metrics to compare different row spacing solutions and capture tradeoffs of these system designs. RESULTS AND CONCLUSIONS We find that, in some circumstances, wider row agrivoltaic solutions that allow for continued mechanized crop production can provide economic benefits over a traditional utility-scale PV system. For most crops examined in this analysis, roughly $\$$200/acre in agricultural profit justified spacing out the panels to at least 31.7 ft. to accommodate agrivoltaic configurations versus PV only configurations. Additionally, opportunities for increased agricultural revenue with agrivoltaic systems allow PV project economics to tolerate a larger range of CAPEX variability while remaining economically viable relative to the PV only configurations. SIGNIFICANCE This framework can be adapted for a wide variety of crops and regions and allows for examination of economically favorable sites for future agrivoltaic systems that utilize different configuration and expand opportunities for agrivoltaics.

14 SOLAR ENERGY↗

Scale-dependent enhancement of productivity and stability in xenic Nannochloropsis cultures

Although bacterial amendments can enhance algal productivity and stability, the development of algal-bacteria consortia for commercial-scale utilization is limited. Here, for this work, we used an established high-throughput approach to generate algal-bacteria consortia, and tested consortia performance at spatial scales from microplates to 320 L raceway ponds. We used both lab and field-reared strains of Nannochloropsis oceanica to build consortia. In some experiments, we imposed environmental perturbations to test the ability of bacteria to enhance algal culture stability. In repeated assays at the scale of well plates, flasks, and bioreactors, strong effects of bacterial amendments on N. oceanica were observed. These effects were most dramatic when cultures experienced stressors such as temperature perturbations or removal of CO 2 augmentation. Isolates that were advanced for field testing included species in the genera Algoriphagus, Oceanicaulis, and Marinobacter. When consortia were generated in the field, positive effects of bacterial amendments were not observed. The amended bacteria were outcompeted, and bacterial community composition across treatments converged after the first grow out. These results highlight the complexity of using consortia in open systems, where interactions between the existing bacterial community, inoculated bacteria, and changing environmental conditions are layered upon other differences in scale and cultivation regimes. Moreover, functionally redundant bacteria are likely present in the field. Following this work, we hypothesize that tight interactions (e.g., obligate relationships between partners) will scale more predictably to outdoor systems. We suggest relying on true synthetic ecological approaches in which the relationships between bacteria and algae partners are well understood, or synthetic ecological approaches coupled with high throughput approaches to design and test consortia. We also recommend future work to examine the effect of algae-bacteria inoculation ratios on productivity and stability, track dynamics of partners through time, and manage ponds to retain beneficial symbioses.

59 BASIC BIOLOGICAL SCIENCES↗

A compact 341 model at TeV scale

We build a gauge model based on the SU(3){sub c}⊗SU(4){sub L}⊗U(1){sub X} symmetry where the scalar spectrum needed to generate gauge boson and fermion masses has a smaller scalar content than usually assumed in literature. We compute the running of its abelian gauge coupling and show that a Landau pole shows up at the TeV scale, a fact that we use to consistently implement those fermion masses that are not generated by Yukawa interactions, including neutrino masses. This is appropriately achieved by non renormalizable effective operators, suppressed by the Landau pole scale. Also, SU(3){sub c}⊗SU(3){sub L}⊗U(1){sub N} models embedded in this gauge structure are bound to be strongly coupled at this same energy scale, contrary to what is generally believed, and neutrino mass generation is rather explained through the same effective operators used in the larger gauge group. Besides, their nice features, as the existence of cold dark matter candidates and the ability to reproduce the observed standard model Higgs-like phenomenology, are automatically inherited by our model. Finally, our results imply that this model is constrained to be observed or discarded soon, since it must be realized at the currently probed energy scale in LHC. - Highlights: • We build a SU(4){sub L}⊗U(1){sub N} electroweak model with a reduced scalar spectrum. • The existence of a Landau pole at TeV scale is investigated. • Fermion masses in this model, including neutrinos, emerge from effective operators. • Imposition of an extra Z{sub 2} symmetry avoids fast proton decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparing and contrasting In-Vial and full-scale systems for sparging volatile analytes

Here, in-vial sparging was demonstrated as an effective, practical alternative to a full-scale sparging system for supporting the analysis of volatile constituents. Using elemental mercury (Hg 0 ) and toluene as representative purgeable analytes, the mass removal for various sparge configurations was measured and a reduced order model was developed and validated. In the primary experiments, Hg 0 in the sparge gas was trapped on activated carbon or gold, thermally desorbed, and quantified using atomic absorption or atomic fluorescence spectroscopy. Toluene experiments using the same in-vial sparge apparatus and sparge parameters were performed to demonstrate the applicability of the reduced order model to a broad range of compounds. Toluene removal was tracked by measuring the remaining toluene in sparged aliquots using Ultraviolet-visible (UV–Vis) spectroscopy. For the sparging, flow rates varied from 25 to 75 mL/min for periods from 0 to 30 min. Sparge performance, mass removal as a function of time, and sparge gas volume were measured for both in-vial and full-scale systems. A model based on dimensionless Henry's Law coefficient, normalized sparge gas volume, and fractional extent of equilibrium matched the experimental data for both compounds and provides a practical tool for future applications. For the conditions tested in this study, the calibrated model indicated that the sparge gas in the in-vial system reached approximately 33% of its equilibrium value before exiting the water surface, while a full-scale system reached approximately 100%. The tests validated the quality, reproducibility, and predictability of sparging performance for both full scale and in-vial sparge systems. Related factors such as waste generation, worker risk, and labor were also assessed. Full scale sparge systems provide the advantage of lower detection levels due to larger sample volume, while the in-vial sparge systems provide advantages for most other factors; including automatability, reducing secondary wastes, lessening the need to clean and check the sparge apparatus, and lowering labor and costs. The data and associated reduced order model support continued development and deployment of in-vial sparge platforms as a practical option for analysis of purgeable analytes such as volatile organic compounds and volatile metals/organometallics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A multiscale design method using interpretable machine learning for phononic materials with closely interacting scales

Manipulating the dispersive characteristics of vibrational waves is beneficial for many applications, e.g., high-precision instruments. architected hierarchical phononic materials have sparked promise tunability of elastodynamic waves and vibrations over multiple frequency ranges. In this article, hierarchical unit-cells are obtained, where features at each length scale result in a band gap within a targeted frequency range. Our novel approach, the ‘‘hierarchical unit-cell template method,’’ is an interpretable machine-learning approach that uncovers global unit-cell shape/topology patterns corresponding to predefined band-gap objectives. A scale-separation effect is observed where the coarse-scale band-gap objective is mostly unaffected by the fine-scale features despite the closeness of their length scales, thus enabling an efficient hierarchical algorithm. Moreover, the hierarchical patterns revealed are not predefined or self-similar hierarchies as common in current hierarchical phononic materials. Furthermore, our approach offers a flexible and efficient method for the exploration of new regions in the hierarchical design space, extracting minimal effective patterns for inverse design in applications targeting multiple frequency ranges.

Architected materials↗

Understanding collective human movement dynamics during large-scale events using big geosocial data analytics

Conventional approaches for modeling human mobility pattern often focus on human activity and movement dynamics in their regular daily lives and cannot capture changes in human movement dynamics in response to large-scale events. With the rapid advancement of information and communication technologies, many researchers have adopted alternative data sources (e.g., cell phone records, GPS trajectory data) from private data vendors to study human movement dynamics in response to large-scale natural or societal events. Big geosocial data such as georeferenced tweets are publicly available and dynamically evolving as real-world events are happening, making it more likely to capture the real-time sentiments and responses of populations. However, precisely-geolocated geosocial data is scarce and biased toward urban population centers. In this research, we developed a big geosocial data analytical framework for extracting human movement dynamics in response to large-scale events from publicly available georeferenced tweets. The framework includes a two-stage data collection module that collects data in a more targeted fashion in order to mitigate the data scarcity issue of georeferenced tweets; in addition, a variable bandwidth kernel density estimation(VB-KDE) approach was adopted to fuse georeference information at different spatial scales, further augmenting the signals of human movement dynamics contained in georeferenced tweets. To correct for the sampling bias of georeferenced tweets, we adjusted the number of tweets for different spatial units (e.g., county, state) by population. To demonstrate the performance of the proposed analytic framework, we chose an astronomical event that occurred nationwide across the United States, i.e., the 2017 Great American Eclipse, as an example event and studied the human movement dynamics in response to this event. Finally, this analytic framework can easily be applied to other types of large-scale events such as hurricanes or earthquakes.

54 ENVIRONMENTAL SCIENCES↗

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗