Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data predictive control”

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,207 records · Page 67

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗

Permeability and Induced Polarization of Mudstones

Electrical measurements can be used to estimate hydraulic properties such as permeability ( k ) in sedimentary rocks. Previous work has focused on sandstones, siltstones, and carbonates, while investigations on mudstones have rarely been reported. In this study, we report on electrical geophysical measurements for 23 mudstone samples using an experimental approach designed to reliably saturate these low permeability mudstones. The modified Hagen-Poiseuille model linking permeability to the formation factor ( F ) and an effective pore radius ( r ) provides an excellent fit to the data set with a near-constant pore radius, indicating that the effective porosity (1/ F ) is the controlling factor on k . In these samples, the surface area normalized to pore volume ( S por ), frequently used in permeability estimation models, varies by 1–2 orders of magnitude and is thus not a reliable proxy of the inverse effective hydraulic radius. The formation factor also exerts the primary control on induced polarization (IP) parameters, whereas Spor shows no relation to the IP parameters. A strong linear relationship is found between IP parameters (imaginary conductivity and normalized chargeability) and surface conductivity, although the proportionality factor is significantly lower than those observed in more permeable rocks and sediments. Apparent relationships between the polarization strength-derived and time constant-derived geophysical length scales and the effective hydraulic radius appear to be driven by variations in the electrochemical parameters (i.e., specific polarizability and diffusion coefficient). Overall, these findings emphasize that predicting hydraulic properties from electrical measurements in fine-grained rocks remains challenging and requires further investigation into the electrochemical properties involved.

58 GEOSCIENCES↗

Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas

This research is a concentrated and coordinated experimental and physics analysis program by a small U.S. team, focused on the topical area of “Long Pulse Control” as stated in the proposal solicitation. The overall goal of the research is to produce world-class analysis regarding forecasting and avoidance of disruptions in a long-pulse, superconducting tokamak at high stability parameters (e.g. normalized beta) with direct, quantifiable figures of merit (i.e. plasma disruption rate - “disruptivity”) and its reduction based on improved stability physics and control. The research targets the backbone of the problem – the basic, main issues – to produce this quantitative analysis by: (1) Creating and applying a new physics-based paradigm to understand how to eliminate disruptions: disruption event characterization, and event chain analysis – a direct manifestation of the approach suggested in the DOE Transients Report, (2) Ensuring that the components required to produce high-quality measurements and analysis needed for disruption event characterization for the tokamak in question (here, KSTAR) exist (e.g. pressure profile peaking, high-accuracy pitch angle data, compensated low frequency toroidal mode number n = 1 mode measurement) to create a primarily physics-based disruption forecasting capability, and (3) Addressing disruption avoidance for the first time in KSTAR by attacking primary disruption ‘event chains’, utilizing the most important associated plasma measurements (dynamic n = 1 field). The focus of this work will be avoidance of macroscopic instabilities, especially global modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-informed machine-learning model of temperature evolution under solid phase processes

We model temperature dynamics during Shear Assisted Proccess Extrusion (ShAPE), a solid phase process that plasticizes feedstock with a rotating tool and subsequently extrudes it into a consolidated tube, rod, or wire. Control of temperature is critical during ShAPE processing to avoid liquefaction, ensure smooth extrusion, and develop desired material properties in the extruded products. Accurate modeling of the complicated thermo-mechanical feedbacks between process inputs, material temperature, and heat generation presents a significant barrier to predictive modeling and process design. In particular, connecting micro-structural scale mechanisms of heat generation to macro-scale predictions of temperature can become computationally intractable. In this work we use a neural network (NN) model of heat generation to bridge this gap, by combining it with a simplified model of the temperature dynamics due to conduction and convection to capture the macro scale evolution of temperature. We inform the construction of the NN heat generation model using crystal plasticity simulations at the micro-structural scale to model the effects of process inputs on generation of heat. We achieved close fits of the temperature dynamics model to a diverse experimental data-set. Further, the relationships learned by the NN model between process inputs and heat generation showed qualitative agreement with those predicted by crystal plasticity simulations.

36 MATERIALS SCIENCE↗

Startup Physics Testing of Advanced Reactors

This report documents the startup physics testing from initial fuel loading through ascension to full power for past advanced reactor startup physics testing programs. The review includes an assessment of what nuclear physics data was measured, why this data was measured, how was the measurement made, and the agreement with predictive reactor performance calculations of that time. The purpose of this review is to establish historical precedence for test inclusion for future advanced reactors planned for demonstration at the National Reactor Innovation Center (NRIC). The historical review includes reactor designs considered to be significantly different from current light-water reactor designs, or use simplified, inherent, passive, or other innovative means to accomplish their safety functions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

C3PU:Real-Time Predictive Charge Control Software for Battery Management Systems (Code for Charge Control and Predictive Unit) [SWR-21-48]

C3PU sets up a model predictive control formulation for the optimal charging control of battery packs and can be deployed for real-time operation on an embedded microprocessor. C3PU uses electrochemical and thermal models of battery packs to predict the charging trajectories of the battery pack over a time horizon under different operating conditions. It then selects the optimal charging trajectory such that a pre-defined objective is minimized, which in the present state of the code is to minimize the charging duration of the battery pack. The optimal charging trajectory is obtained by solving an underlying mathematical problem. However, C3PU is flexible to incorporate other charging objectives. The underlying battery models can be swapped as well, as long as they follow certain mathematical properties. Mathematically, optimal control problems (in this case, optimally controlling the charging current of the battery pack) are computationally expensive. C3PU implements advanced numerical techniques, namely pseudo-spectral optimization, to reduce the computational burden of the underlying problem to solve. This allows for the problem to be solved in real-time in an embedded system. Once the optimal charging trajectory is computed over a time horizon, C3PU can package such data and send it out via appropriate communication protocols.

Mishra, Partha↗

Integrated parameter and process learning for hydrologic and biogeochemical modules in Earth System Models

Focus area: Primary focal area #2; secondary focal area #3: Learning about parameters and processes of land surface hydrologic and biogeochemical models in Earth System models by integrating machine learning, physics, and big data. Science challenges: How do we maximally leverage big-data observations to improve hydrobiogeochemical process description and parameterization so that such modules more realistically capture hydrologic and vegetation responses and feedbacks under the future climate? For example, how can we leverage physics, limited observations of vegetation and streamflow to better estimate evapotranspiration, and, relatedly, net primary productivity, especially for drought areas? Vegetation plays a critical role in regional and global water cycles; however, existing vegetation models have failed to predict vegetation response to droughts (McDowell & Xu, 2017) , arctic greening (Keenan & Riley, 2018) , and critical transitions between forest and savanna (Hirota et al., 2011) . These studies suggest that when we build process-based models (PBM) parameterized from regional and global plant traits, we tend to poorly describe plant adaptation and local-scale competition processes. The models and their associated parameters assigned for different regions in the world are not capturing essential heterogeneity in vegetation responses at finer spatial scales. Many parameters of the land surface models control hydrology and vegetation dynamics at the same time. The heterogeneity in vegetation response is a function of (i) plant type, (ii) plant size, (iii) competition and succession, (iv) environmental controls, and (v) local variations due to the unique ecological community that are very difficult to describe (e.g., the size of gaps resulting from fire that facilitated the coexistence of pioneering species). In the demographic models, only factors (i) and (iv) were captured, and plant types were generally described only by leaf phenology and climate zones. With current demographic models, we generally consider more traits to define plant types (i) and calibrate these traits to consider factors (ii), (iii) and (iv); however, it is substantially challenging to scale to regional and global simulations due to trait variations across space (Ali et al., 2016). Moreover, it has been noted that hillslope processes, including ridge-to-valley flow and sunny vs. shady slopes are primary organizers of water, energy, and vegetation (Clark et al., 2015; Fan et al., 2019) . Although gradual improvements in the hydrologic model component in earth system models may reduce this error (at a remarkably slow pace), the long-term, gradual impact of hydrology on plant traits are not well captured. Recent work showed that the hydrologic controls exerted by groundwater and lateral flow are primary regulators of rooting depth (Fan et al., 2017) . Such hydrologic controls have seldom been reflected in vegetation model parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Characterization of the In Vivo Deuteration of Native Phospholipids by Mass Spectrometry Yields Guidelines for Their Regiospecific Customization

Customization of deuterated biomolecules is vital for many advanced biological experiments including neutron scattering. However, because it is challenging to control the proportion and regiospecificity of deuterium incorporation in live systems, often only two or three synthetic lipids are mixed together to form simplistic model membranes. This limits the applicability and biological accuracy of the results generated with these synthetic membranes. Despite some limited prior examination of deuterating Escherichia coli lipids in vivo, this approach has not been widely implemented. In this report an extensive mass spectrometry-based profiling of E. coli phospholipid deuteration states with several different growth media was performed, and a computational method to describe deuterium distributions with a one-number summary is introduced. The deuteration states of 36 lipid species were quantitatively profiled in 15 different growth conditions, and tandem mass spectrometry was used to reveal deuterium localization. Regressions were employed to enable the prediction of lipid deuteration for untested conditions. Small-angle neutron scattering was performed on select deuterated lipid samples, which validated the deuteration states calculated from the mass spectral data. Based on these experiments, guidelines for the design of specifically deuterated phospholipids are described. This unlocks even greater capabilities from neutron-based techniques, enabling experiments that were formerly impossible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonparametric extensions of nuclear equations of state: Probing the breakdown scale of relativistic mean-field theory

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. Here, we further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

Equations of state of nuclear matter↗

Artificial Intelligence for Accelerating Nuclear Applications, Science, and Technology

Artificial intelligence (AI) and machine learning (ML) methods have had significant impacts in science and technology in recent years. These methods for generating models from datasets or logic-based algorithms that emulate aspects of human performance can similarly accelerate the fields of nuclear applications, science, and technology toward the IAEA goals of contributing to peace, health, and prosperity. In order to accomplish advances with AI in general and ML in particular across these fields, IAEA can play a significant role by establishing, hosting and curating centralised resources, including databases, adhering to FAIR (findable, accessible, interoperable and reusable) principles and Open Science best practices, providing stewardship of data sharing, supporting training efforts and development of relevant workforces, as well as enabling connections among the scientific, technology, mathematics, AI and ethics communities. Many areas can benefit from the use of AI in the realm of nuclear applications. In human health, these areas include clinical research, epidemiology, nutrition, medical imaging, radiotherapy and education of health professionals. AI-based tools are also being used to facilitate different clinical tasks in imaging, computer-assisted diagnosis in mammography and lung cancer screening programmes, and dose prediction in nuclear medicine procedures. ML methods in particular may also increase the efficiency and accuracy of the analysis of computerised tomography and dual-energy absorptiometry scans for body composition and bone analysis. The application of AI methods to nuclear and related technologies in food and agriculture can lead to significant advances and improved efficiency in the optimisation of agricultural production, food product development, management of supply chains, food safety and food authenticity control. In the water and environmental sector, AI can help inform policies to mitigate the world’s water problems. The application of AI techniques to hydrology and environmental sciences is expected to improve patterns identification and enable model predictions under a changing climate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S 0 ( T ) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S 0 ( T ) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S 0 ( T ) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Science & Technology - Other Topics↗

A Method for Projecting Cloud Shadows Onto a Central Receiver Field to Predict Receiver Damage

This work demonstrates methods of mapping high-spatial-resolution direct normal irradiance (DNI) data from satellites, Total Sky Imagers (TSIs), and analogous data sources onto a heliostat field for characterizing the spatial and temporal variation of the incident flux on a central receiver tower during cloud transient events. The mapping methods are incorporated into an optical software module that interfaces with CoPylot–SolarPILOT’s python API– to provide computationally efficient optical simulation of the heliostat field and the solar power tower. Eventually, this optical model will be incorporated into optimization models whereby a plant operator can understand the effects of cloud transient events on overall power production and receiver lifetime due to creep-fatigue damage and therefore make better informed decisions about receiver shutdown events. By more accurately modelling the effects of cloud events on receiver flux maps, this work may determine the magnitude and frequency of thermal cycling on receiver tubes and panels using actual or realistic cloud shapes instead of averaged DNI values–which may undercount the total cycle number. This work may also prevent unnecessary plant shutdowns due to overly precautionary control strategies and characterize the relative impact of various cloud types on receiver life. We plan to eventually integrate this methodology into the System Advisor Model (SAM) to improve performance model accuracy during periods of cloudiness. In this paper, we demonstrate generating DNI maps and mapping them to a solar field in CoPylot using 10 m resolution data from publicly available Sentinel-2 satellite data over the Crescent Dunes plant.

Mullin, Matthew↗

Characterization of MoS 2 films via simultaneous grazing incidence X-ray diffraction and grazing incidence X-ray fluorescence (GIXRD/GIXRF)

Physical vapor deposited (PVD) molybdenum disulfide (nominal composition MoS 2 ) is employed as a thin film solid lubricant for extreme environments where liquid lubricants are not viable. The tribological properties of MoS 2 are highly dependent on morphological attributes such as film thickness, orientation, crystallinity, film density, and stoichiometry. These structural characteristics are controlled by tuning the PVD process parameters, yet undesirable alterations in the structure often occur due to process variations between deposition runs. Nondestructive film diagnostics can enable improved yield and serve as a means of tuning a deposition process, thus enabling quality control and materials exploration. Grazing incidence X-ray diffraction (GIXRD) for MoS 2 film characterization provides valuable information about film density and grain orientation (texture). However, the determination of film stoichiometry can only be indirectly inferred via GIXRD. The combination of density and microstructure via GIXRD with chemical composition via grazing incidence X-ray fluorescence (GIXRF) enables the isolation and decoupling of film density, composition, and microstructure and their ultimate impact on film layer thickness, thereby improving coating thickness predictions via X-ray fluorescence. We have augmented an existing GIXRD instrument with an additional X-ray detector for the simultaneous measurement of energy-dispersive X-ray fluorescence spectra during the GIXRD analysis. This combined GIXRD/GIXRF analysis has proven synergetic for correlating chemical composition to the structural aspects of MoS 2 films provided by GIXRD. We present the usefulness of the combined diagnostic technique via exemplar MoS 2 film samples and provide a discussion regarding data extraction techniques of grazing angle series measurements.

MoS2↗

Wind turbine gearbox fault prognosis using high-frequency SCADA data

Condition-based maintenance using routinely collected Supervisory Control and Data Acquisition (SCADA) data is a promising strategy to reduce downtime and costs associated with wind farm operations and maintenance. New approaches are continuously being developed to improve the condition monitoring for wind turbines. Development of normal behaviour models is a popular approach in studies using SCADA data. This paper first presents a data-driven framework to apply normal behaviour models using an artificial neural network approach for wind turbine gearbox prognostics. A one-class support vector machine classifier, combining different error parameters, is used to analyse the normal behaviour model error to develop a robust threshold to distinguish anomalous wind turbine operation. A detailed sensitivity study is then conducted to evaluate the potential of using high-frequency SCADA data for wind turbine gearbox prognostics. The results based on operational data from one wind turbine show that, compared to the conventionally used 10-min averaged SCADA data, the use of high-frequency data is valuable as it leads to improved prognostic predictions. High-frequency data provides more insights into the dynamics of the condition of the wind turbine components and can aid in earlier detection of faults.

17 WIND ENERGY↗

Galaxy pairs in the Sloan Digital Sky Survey – XV. Properties of ionized outflows

ABSTRACT Powerful outflows are thought to play a critical role in galaxy evolution and black hole growth. We present the first large-scale systematic study of ionized outflows in paired galaxies and post-mergers compared to a robust control sample of isolated galaxies. We isolate the impact of the merger environment to determine if outflow properties depend on merger stage. Our sample contains ∼4000 paired galaxies and ∼250 post-mergers in the local universe (0.02 ≤ z ≤ 0.2) from the Sloan Digital Sky Survey Data Release 7 (SDSS DR 7) matched in stellar mass, redshift, local density of galaxies, and [O iii] λ5007 luminosity to a control sample of isolated galaxies. By fitting the [O iii] λ5007 line, we find ionized outflows in ∼15 per cent of our entire sample. Outflows are much rarer in star-forming galaxies compared to active galactic nuclei (AGNs), and outflow incidence and velocity increase with [O iii] λ5007 luminosity. Outflow incidence is significantly elevated in the optical + mid-infrared selected AGN compared to purely optical AGN; over 60 per cent show outflows at the highest luminosities ($L_{\mathrm{[OIII]~\lambda 5007}}\, \gtrsim$ 1042 erg s−1), suggesting mid-infrared AGN selection favours galaxies with powerful outflows, at least for higher [O iii] λ5007 luminosities. However, we find no statistically significant difference in outflow incidence, velocity, and luminosity in mergers compared to isolated galaxies, and there is no dependence on merger stage. Therefore, while interactions are predicted to drive gas inflows and subsequently trigger nuclear star formation and accretion activity, when the power source of the outflow is controlled for, the merging environment has no further impact on the large-scale ionized outflows as traced by [O iii] λ5007.

79 ASTRONOMY AND ASTROPHYSICS↗

C{sub 2}N{sub 2} Vertical Profile in Titan’s Stratosphere

In this paper, we present the first measurements of the vertical distribution of cyanogen (C{sub 2}N{sub 2}) in Titan’s lower atmosphere at different latitudes and seasons, using Cassini’s Composite Infrared Spectrometer far-infrared data. We also study the vertical distribution of three other minor species detected in our data: methylacetylene (C{sub 3}H{sub 4}), diacetylene (C{sub 4}H{sub 2}), and H{sub 2}O, in order to compare them to C{sub 2}N{sub 2}, but also to get an overview of their seasonal and meridional variations in Titan’s lower stratosphere from 85 km to 225 km. We measured an average volume mixing ratio of C{sub 2}N{sub 2} of 6.2±0.8×10{sup −11} at 125 km at the equator, but poles exhibit a strong enrichment in C{sub 2}N{sub 2} (up to a factor 100 compared to the equator), greater than what was measured for C{sub 3}H{sub 4} or C{sub 4}H{sub 2}. Measuring C{sub 2}N{sub 2} profiles provides constraints on the processes controlling its distribution, such as bombardment by Galactic Cosmic Rays which seem to have a smaller influence on C{sub 2}N{sub 2} than predicted by photochemical models.

79 ASTRONOMY AND ASTROPHYSICS↗

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗