SEARCH · Engineering Papers

Results for “Dynamic model reduction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A nudged hybrid analysis and modeling approach for realtime wake-vortex transport and decay prediction

Here, we put forth a long short-term memory (LSTM) nudging framework for the enhancement of reduced order models (ROMs) of fluid flows utilizing noisy measurements for air traffic improvements. Toward emerging applications of digital twins in aviation, the proposed approach allows for constructing a realtime predictive tool for wake-vortex transport and decay systems. We build on the fact that in realistic application, there are uncertainties in initial and boundary conditions, model parameters, as well as measurements. Moreover, conventional nonlinear ROMs based on Galerkin projection (GROMs) suffer from imperfection and solution instabilities, especially for advection-dominated flows with slow decay in the Kolmogorov n-width. In the presented LSTM nudging (LSTM-N) approach, we fuse forecasts from a combination of imperfect GROM and uncertain state estimates, with sparse Eulerian sensor measurements to provide more reliable predictions in a dynamical data assimilation framework. We illustrate our concept by solving the two-dimensional vorticity transport equation. We investigate the effects of measurements noise and state estimate uncertainty on the performance of the LSTM-N behavior. We also demonstrate that it can sufficiently handle different levels of temporal and spatial measurement sparsity, and offer a huge potential in developing next-generation digital twin technologies for aerospace applications.

42 ENGINEERING↗

DOE OSTI2019

Machine Learning-Guided Approach for Studying Solvation Environments

Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, in this work, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1–, and 2–, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Multi-Electron Transfer at H-Terminated p-Si Electrolyte Interfaces: Large Photovoltages under Inversion Conditions

Here, photovoltages for hydrogen terminated p-Si(111) in an acetonitrile electrolyte were quantified with methyl viologen [1,1'-(CH 3 ) 2 -4,4'-bipyridinium](PF 6 ) 2 , abbreviated MV 2+ , and [Ru(bpy) 3 ](PF 6 ) 2 , where bpy is 2,2'-bipyridine, that respectively undergo two and three one-electron transfer reductions. The reduction potentials, E°, of the two MV 2+ reductions occurred at energies within the forbidden bandgap, while the three [Ru(bpy) 3 ] 2+ reductions occurred within the continuum of conduction band states. Bandgap illumination resulted in reduction that were more positive than that measured with a degenerately doped n + -Si demonstrative of a photovoltage, Vph, that increased in the order: MV 2+/+ (260 mV) < MV +/0 (400 mV) < Ru 2+/+ (530 mV) ~ Ru +/0 (540 mV) ~ Ru 0/– (550 mV). Pulsed 532 nm excitation generated electron-hole pairs whose dynamics were nearly constant under depletion conditions and increased markedly as the potential was raised or lowered. A long wavelength absorption feature assigned to conduction band electrons provided additional evidence for the presence of an inversion layer. Collectively, the data reveal that the most optimal photovoltage, as well at the longest electron-hole pair lifetime and the highest surface electron concentration, occur when E° lies energetically within the unfilled conduction band states where an inversion layer is present. The bell-shaped dependence for electron-hole pair recombination with the surface potential was predicted by the time-honored SRH model providing a clear indication that this interface provides access to all four bias conditions, i.e., accumulation, flat band, depletion, and inversion. The implications of these findings for photocatalysis applications and solar energy conversion are discussed.

14 SOLAR ENERGY↗

DOE OSTI2023

Data-driven modeling to enhance municipal water demand estimates in response to dynamic climate conditions

Altered precipitation and temperature patterns from a changing climate will affect supply, demand, and overall municipal water system operations throughout the arid western U.S. While supply forecasts leverage hydrological models to connect climate influences with surface water availability, demand forecasts typically estimate water use independent of climate and other externalities. Stemming from an increased focus on seasonal water demand management, we use the Salt Lake City, Utah municipal water system as a test bed to assess model accuracy versus complexity trade-offs between simple climate-independent econometric-based models and complex climate-sensitive data-driven models to average to extreme wet and dry climate conditions—representative of a new climate normal. Here, the climate-independent model displayed low performance during extreme dry conditions with predictions exceeding 90% and 40% of the observed monthly and seasonal volumetric demands, respectively, which we attribute to insufficient model complexity. The climate-sensitive models displayed greater accuracy in all conditions, with an ordinary least squares model demonstrating a measurable reduction in prediction bias (3.4% vs. -27.3%) and RMSE (74.0 lpcd vs. 294 lpcd) compared to the climate-independent model. The climate-sensitive workflow increased model accuracy and characterized climate-demand interactions, demonstrating a novel tool to enhance water system management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Optimal management of electric hotel loads in mild hybrid heavy duty truck

Here, the problem of engine idling for heavy-duty trucks has been under study for decades with Auxiliary Power Units (APUs) and Truck Stop Electrification (TSE) as the most compelling solutions. With the electrification of trucks approaching feasibility in terms of cost-effective technology, hybridization offers another “degree of freedom” to tackle the problem. Here, this work aims at exploiting a battery pack of a 48 V mild-hybrid heavy-duty truck to store sufficient onboard energy for powering the auxiliary loads during the hoteling. This problem is not trivial, as the battery packs typically cannot recover the entire energy required through regeneration alone; hence an optimal energy management strategy needs to be employed to charge the battery through the engine during drive operation. This strategy optimizes powertrain performances among the four modes: (i) Engine of Coasting (EOC), (ii) Regeneration by braking, (iii) Regeneration by engine, and (iv) engine idling. This paper presents the development of a Dynamic Programming (DP) framework that employs a multi-objective cost function to minimize the fuel consumption and maximize the regeneration using the above-mentioned four modes. A typical heavy-duty truck drive cycle is used to represent the drive phase, with mandatory hoteling stops as per regulations. A comprehensive powertrain model is developed using validated components’ model. The DP employs two state variables: battery State-of-Charge (SOC) and engine mode, and three control inputs: (i) the engine ON–OFF state, (ii) clutch engagement state, and (iii) power request at the Electric Machine (EM) for calculating optimal SOC trajectory. The framework also tackles rapid engine ON–OFF scenarios to avoid the challenges associated with DP and the compromises in fuel cost with those approaches. Finally, the effectiveness of the proposed framework is tested for potential fuel savings on two different battery packs by performing the full cycle simulations. The results show 6.47% of fuel consumption reduction as compared to traditional APU-based heavy-duty truck.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2023

Circulating Current Reduction of Back-to-Back MMC with Advanced Grid-Support Functionalities for Medium-Voltage Applications

A Lyapunov energy function based circulating current reduction scheme for a back to back modular multilevel converter based system for a microgrid application has been proposed in this digest. In order to ensure no circulating or zero sequence currents within the modules, the dynamics from one leg of the converter to the other are evaluated. A Lyapunov energy function is then developed and utilizing the circulating current dynamics of the converter the overall control law is implemented. In this digest a decentralized hierarchical control architecture has been implemented for the main current control; where the outer loop or the central controller implements the higher level grid interconnection functions as per IEEE 1547-2018 standards and the inner loop based on synchronous reference frame implements the current control architecture. The inner loops' outputs are the commanded converter voltages either on either side which are augmented by the Lyapunov energy function based circulating current suppression architecture. To evaluate the efficacy of the overall system, the overall system is modeled around a 300kVA back to back SiC multilevel modular converter used as a direct-connected medium voltage (MV) grid intertie based on MATLAB/Simulink domain. Several important case studies are presented in this paper to prove the effectiveness of the proposed architecture. case study results are presented.

circulating current↗

DOE OSTI2024

Modeling Diurnal and Annual Ethylene Generation from Solar-Driven Electrochemical CO 2 Reduction Devices

Integrated solar fuels devices for CO 2 reduction (CO 2 R) are a promising technology class towards achieving net-negative carbon emissions. Designing integrated CO 2 R solar fuels devices requires careful co-design of electrochemical and photovoltaic components as well as consideration of the diurnal and seasonal effects of solar irradiance, temperature, and other meteorological factors expected for ‘on-sun’ deployment. Here, using a photovoltaic-electrochemical (PV-EC) platform, we developed a temperature and potential-dependent diurnal and annual model using experimental CO 2 R performance of Cu-based electrocatalysts, local meteorological data from the National Solar Radiation Database (NSRD), and modeled performance of commercial c-Si PVs. We simulated diurnal product outputs with and without the effects of ambient temperature to determine gaseous product temperature sensitivity. From these outputs, we observed seasonal variation in gaseous product generation, with up to two-fold increases in ethylene productivity between the Winter and Summer, analyzed the consequences of dynamic cloud coverage, and identified periods where device cooling/heating mechanisms could be implemented to maximize ethylene generation. Finally, we modeled the annual ethylene generation for a scaled 1 MW solar farm at three different locations (Beijing, CN; Sydney, AUS; Barstow, CA) to determine the consequences of local meteorological climates on PV-EC CO 2 R product output, recording a maximum ethylene output of 18.5 tonne/yr at Barstow. Overall, this model presents a critical tool for streamlining the translation of experimental solar-driven electrochemical research to real-world implementation.

Yap, Kyra M. K.↗

DOE OSTI2024

Personalized whole-brain neural mass models reveal combined Aβ and tau hyperexcitable influences in Alzheimer's disease.

Neuronal dysfunction and cognitive deterioration in Alzheimer's disease (AD) are likely caused by multiple pathophysiological factors. However, mechanistic evidence in humans remains scarce, requiring improved non-invasive techniques and integrative models. We introduce personalized AD computational models built on whole-brain Wilson-Cowan oscillators and incorporating resting-state functional MRI, amyloid-β (Aβ) and tau-PET from 132 individuals in the AD spectrum to evaluate the direct impact of toxic protein deposition on neuronal activity. This subject-specific approach uncovers key patho-mechanistic interactions, including synergistic Aβ and tau effects on cognitive impairment and neuronal excitability increases with disease progression. The data-derived neuronal excitability values strongly predict clinically relevant AD plasma biomarker concentrations (p-tau217, p-tau231, p-tau181, GFAP) and grey matter atrophy obtained through voxel-based morphometry. Furthermore, reconstructed EEG proxy quantities show the hallmark AD electrophysiological alterations (theta band activity enhancement and alpha reductions) which occur with Aβ-positivity and after limbic tau involvement. Microglial activation influences on neuronal activity are less definitive, potentially due to neuroimaging limitations in mapping neuroprotective vs detrimental activation phenotypes. Mechanistic brain activity models can further clarify intricate neurodegenerative processes and accelerate preventive/treatment interventions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

Optimisation of the Kaplan hydropower system via PID 2 and digital twin

Here, this paper proposes a proportional–integral-double–derivative (PID 2 ) optimisation method for the Kaplan hydropower system by building a digital twin. The study first uses one multilayer perceptron (MLP) to model the hydroturbine dynamic and then adopts three connected MLPs to model the generator dynamic, both in an open-loop fashion. Inspired by stochastic distribution control (SDC) theory, we regard the training of the turbine's neural network model as a process control problem, and we propose minimising entropy loss to update the network parameters. The next step is to build the digital twin by connecting the neural network models with a PID 2 controller and a lead-lag exciter and run the whole model in a closed-loop fashion. After that, a binary search approach is applied to optimise the PID 2 parameters based on the obtained digital twin model. The simulation results show that the proposed method can reduce the mean square tracking error by more than 90%. Furthermore, the method is extended to jointly optimise the PID 2 controller and excitation system gains through multiobjective optimisation, leveraging Pareto frontier analysis to balance active power and voltage tracking performance. Simulation results confirm the effectiveness of the proposed method, achieving a 83.46% reduction in relative mean square error of active power, a 47.13% reduction in terminal voltage tracking error, and an 82.78% improvement in the overall scalarized objective.

Hydropower system↗

DOE OSTI2020

Coordinated Multiscale Modeling and Synthesis of Novel Nanostructured Composite Membranes for Solar Fuels Generation (Final Report)

Generation of carbon-based liquid fuels from reduction of CO 2 using solar energy requires the development of improved membranes that provide good ionic conductivity and mechanical properties while minimizing the crossover of gases and reduction products. Using closely coordinated multiscale modeling, synthesis and characterization studies we propose to investigate and develop novel nanostructured composite membranes with applications to conversion of carbon dioxide to storable chemical fuels. The envisioned proton exchange membranes (PEMs) will be based on a support of self- assembling polymer-modified nanoparticles with controllable porosity and tortuosity. The use of a nanoparticle support will provide control of transport and mechanical properties of the membranes. Tuning the copolymer modification of nanoparticles (polymer architectures, brush grafting density) will allow to optimize the ion (proton) transport via the Grotthuss mechanism through narrow hydrated channels while minimizing crossover of gases and reaction products due to the high density of the grafted brush and use of a non-ionic conducting glassy polymers to fill spaces between polymer-grafted nanoparticles. In Phase I synthesis of the copolymers and their attachment to the nanoparticle support will be guided by multiscale simulations that include atomistic molecular dynamics and continuum-level transport simulations. The former will include reactive and non-reactive atomistic simulations that will provide key insight into nanoscale polymer morphology and ion/molecular transport mechanisms need to optimize the copolymer structure. The latter will provide insight into the role of nanoparticle self-assembly on the global transport of ionic and molecular species. The proton conduction and CO 2 /methanol permeability of the membranes will be experimentally characterized as well as simulated and will be correlated with morphology and polymer structure. Prototypes of optimal membranes will be fabricated and characterized at the end of the project. Successful development of improved nanostructured composite PEMs and demonstration of the simulation-guided materials-by-design paradigm will result in our ability to design and synthesize membranes for a wide variety of solar fuels generation and related electrochemical applications.

14 SOLAR ENERGY↗

NASA NTRS2017

Genetic Algorithm-Guided, Adaptive Model Order Reduction of Flexible Aircrafts

This paper presents a methodology for automated model order reduction (MOR) of flexible aircrafts to construct linear parameter-varying (LPV) reduced order models (ROM) for aeroservoelasticity (ASE) analysis and control synthesis in broad flight parameter space. The novelty includes utilization of genetic algorithms (GAs) to automatically determine the states for reduction while minimizing the trial-and-error process and heuristics requirement to perform MOR; balanced truncation for unstable systems to achieve locally optimal realization of the full model; congruence transformation for "weak" fulfillment of state consistency across the entire flight parameter space; and ROM interpolation based on adaptive grid refinement to generate a globally functional LPV ASE ROM. The methodology is applied to the X-56A MUTT model currently being tested at NASA/AFRC for flutter suppression and gust load alleviation. Our studies indicate that X-56A ROM with less than one-seventh the number of states relative to the original model is able to accurately predict system response among all input-output channels for pitch, roll, and ASE control at various flight conditions. The GA-guided approach exceeds manual and empirical state selection in terms of efficiency and accuracy. The adaptive refinement allows selective addition of the grid points in the parameter space where flight dynamics varies dramatically to enhance interpolation accuracy without over-burdening controller synthesis and onboard memory efforts downstream. The present MOR framework can be used by control engineers for robust ASE controller synthesis and novel vehicle design.

Numerical Analysi↗

NASA NTRS2016

Modeling and Simulation of Radiative Compressible Flows in Aerodynamic Heating Arc-Jet Facility

Numerical simulations of an arc heated flow inside NASA's 20 [MW] Aerodynamics heating facility (AHF) are performed in order to investigate the three-dimensional swirling flow and the current distribution inside the wind tunnel. The plasma is considered in Local Thermodynamics Equilibrium(LTE) and is composed of Air-Argon gas mixture. The governing equations are the Navier-Stokes equations that include source terms corresponding to Joule heating and radiative cooling. The former is obtained by solving an electric potential equation, while the latter is calculated using an innovative massively parallel ray-tracing algorithm. The fully coupled system is closed by the thermodynamics relations and transport properties which are obtained from Chapman-Enskog method. A novel strategy was developed in order to enable the flow solver and the radiation calculation to be preformed independently and simultaneously using a different number of processors. Drastic reduction in the computational cost was achieved using this strategy. Details on the numerical methods used for space discretization, time integration and ray-tracing algorithm will be presented. The effect of the radiative cooling on the dynamics of the flow will be investigated. The complete set of equations were implemented within the COOLFluiD Framework. Fig. 1 shows the geometry of the Anode and part of the constrictor of the Aerodynamics heating facility (AHF). Fig. 2 shows the velocity field distribution along (x-y) plane and the streamline in (z-y) plane.

Pasma flows↗

DOE OSTI2022

Efficient hybrid explicit-implicit learning for multiscale problems

Splitting method is a powerful method to handle application problems by splitting physics, scales, domain, and so on. Many splitting algorithms have been designed for efficient temporal discretization. Here, in this paper, our goal is to use temporal splitting concepts in designing machine learning algorithms and, at the same time, help splitting algorithms by incorporating data and speeding them up. We propose a machine learning assisted splitting scheme which improves the efficiency of the scheme meanwhile preserves the accuracy. We consider a recently introduced multiscale splitting algorithms, where the multiscale problem is solved on a coarse grid. To approximate the dynamics, only a few degrees of freedom are solved implicitly, while others explicitly. This splitting concept allows identifying degrees of freedom that need implicit treatment. In this paper, we use this splitting concept in machine learning and propose several strategies. First, the implicit part of the solution can be learned as it is more difficult to solve, while the explicit part can be computed. This provides a speed-up and data incorporation for splitting approaches. Secondly, one can design a hybrid neural network architecture because handling explicit parts requires much fewer communications among neurons and can be done efficiently. Thirdly, one can solve the coarse grid component via PDEs or other approximation methods and construct simpler neural networks for the explicit part of the solutions. We discuss these options and implement one of them by interpreting it as a machine translation task. This interpretation of the splitting scheme successfully enables us using the Transformer since it can perform model reduction for multiple time series and learn the connection between them. We also find that the splitting scheme is a great platform to predict the coarse solution with insufficient information of the target model: the target problem is partially given and we need to solve it through a known problem which approximates the target. Our machine learning model can incorporate and encode the given information from two different problems and then solve the target problems. We conduct four numerical examples and the results show that our method is stable and accurate.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A combined SECM and electrochemical AFM approach to probe interfacial processes affecting molecular reactivity at redox flow battery electrodes

Redox flow batteries are attractive technologies for grid energy storage since they use solutions of redox-active molecules that enable a superior scalability and the decoupling of power and energy density. However, the reaction mechanisms of the redox active components at RFB electrodes are complex, and there is currently a pressing need to understand how interfacial processes impact the kinetics and operational reversibility of RFB systems. Here, we developed a combined electrochemical imaging methodology rooted in scanning electrochemical microscopy (SECM) and atomic force microscopy (AFM) for exploring the impact of electrode structure and conditioning on the electron transfer properties of model redox-active dialkoxybenzene derivatives, 2,5-di-tert-butyl-1,4-bis(2-methoxyethoxy)benzene (C1) and 2,3-dimethyl-1,4-dialkoxybenzene (C7). Using AFM and secondary-ion mass spectrometry (SIMS), we observed the formation of interfacial films with distinct mechanical properties compared to those of cleaved graphitic surfaces, and exclusively during reduction of electrogenerated radical cations. These films had an impact on the median rate and distribution of the electron transfer rate constant at the basal plane of multilayer and single layer graphene electrodes, displaying kinetically-limited values that did not yield the activation expected per the Butler–Volmer model with a transfer coefficient ~0.5. These changes were dependent on redoxmer structure: furthermore, SECM showed strong attenuation of C7 kinetics by a surface layer on MLG and SLG, while C1 kinetics were only affected by SLG. SECM and AFM results together show that these limiting films operate exclusively on the basal plane of graphite, with the edge plane showing a relative insensitivity to cycling and operation potential. Ultimately, this integrated electrochemical imaging methodology creates new opportunities to understand the unique role of interfacial processes on the heterogeneous reactivity of redoxmers at electrodes for RFBs, with a future role in elucidating phenomena at high active concentrations and spatiotemporal variations in electrode dynamics.

25 ENERGY STORAGE↗

DOE OSTI2023

The Coastal Carbon Library and Atlas: Open source soil data and tools supporting blue carbon research and policy

Abstract Quantifying carbon fluxes into and out of coastal soils is critical to meeting greenhouse gas reduction and coastal resiliency goals. Numerous ‘blue carbon’ studies have generated, or benefitted from, synthetic datasets. However, the community those efforts inspired does not have a centralized, standardized database of disaggregated data used to estimate carbon stocks and fluxes. In this paper, we describe a data structure designed to standardize data reporting, maximize reuse, and maintain a chain of credit from synthesis to original source. We introduce version 1.0.0. of the Coastal Carbon Library, a global database of 6723 soil profiles representing blue carbon‐storing systems including marshes, mangroves, tidal freshwater forests, and seagrasses. We also present the Coastal Carbon Atlas, an R‐shiny application that can be used to visualize, query, and download portions of the Coastal Carbon Library. The majority (4815) of entries in the database can be used for carbon stock assessments without the need for interpolating missing soil variables, 533 are available for estimating carbon burial rate, and 326 are useful for fitting dynamic soil formation models. Organic matter density significantly varied by habitat with tidal freshwater forests having the highest density, and seagrasses having the lowest. Future work could involve expansion of the synthesis to include more deep stock assessments, increasing the representation of data outside of the U.S., and increasing the amount of data available for mangroves and seagrasses, especially carbon burial rate data. We present proposed best practices for blue carbon data including an emphasis on disaggregation, data publication, dataset documentation, and use of standardized vocabulary and templates whenever appropriate. To conclude, the Coastal Carbon Library and Atlas serve as a general example of a grassroots F.A.I.R. (Findable, Accessible, Interoperable, and Reusable) data effort demonstrating how data producers can coordinate to develop tools relevant to policy and decision‐making.

Holmquist, James R.↗

DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗

DOE OSTI2022

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗