SEARCH · Engineering Papers

Results for “Coefficient of performance”

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,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Removing flat directions in standard model EFT fits: How polarized electron-ion collider data can complement the LHC

We study the potential of future electron-ion collider (EIC) data to probe four-fermion operators in the Standard Model effective field theory (SMEFT). The ability to perform measurements with both polarized electron and proton beams at the EIC provides a powerful tool that can disentangle the effects from different SMEFT operators. We compare the potential constraints from an EIC with those obtained from Drell-Yan data at the Large Hadron Collider. We show that EIC data play an important complementary role since they probe combinations of Wilson coefficients not accessible through available Drell-Yan measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Design of a Single Layer Metamaterial for Pressure Offloading of Transtibial Amputees

While using a prosthesis, transtibial amputees may experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients by using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells which can be 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh 3rd order model. The 3rd order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from: 1) a common homogenous silicone liner, 2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and 3) a metamaterial solution with optimized Yeoh 3rd order coefficients. When compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Verification of RESRAD-OFFSITE Code (V.4)

This report documents the verification of RESRAD-OFFSITE Version 4.0 and describes, where necessary, the verification of the following: • The data comprising the standard dose and risk coefficient libraries in the RESRAD database files Master_dcf_ICRP07.mdb and Master_dcf_2k.mdb. • The extraction and transfer of the data from the selected database file to the computational code by the RESRAD-OFFSITE 4.0 interface, ResOWin.exe. • The different processes that are modeled by the main computational code in RESRAD OFFSITE 4.0, ResOMain.exe. • The data displayed in the graphical and text reports. Many verifications were performed as part of the quality assurance quality control program associated with the development and release of RESRAD-OFFSITE 4.0, namely: • developer testing, • internal independent testing, and • release testing. Some were also performed in response to questions from users regarding the performance of the code. The main text of the report focuses on summarizing a subset of those tests, both independent and developer tests that verified the computations performed by the code. The verifications included in this report served as the basis for the development of the release tests of the computational executables and provided the quantitative results to be compared with the code output. The input and output interfaces and the data transfers between the various executables of the code were tested while performing the verification testing. They were tested intentionally during release testing. This report also provides some basic information to help in understanding the activities that were verified. The report: • outlines the components of RESRAD-OFFSITE 4.0 and the interconnections between these components, • outlines the processes modeled by the computational code, • provides summary figures and tables to offer confirmation of the verification of the computational components of the code, • reproduces the verifiers’ reports, if available, in individual appendices, • refers to the previous verification report (Yu et al. 2011) for more details about some of the verifications, and • reproduces the test cases and the testers’ reports from the release testing in individual appendices, when possible.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

CFD Modeling of Concentration Polarization in High-Performance Flat-Sheet Modules

Membrane-based gas separation is a promising, regeneration-free route for industrial CO₂ capture, but its performance is often limited by concentration polarization (CP). Classical 1-D models neglect CP and thus overpredict fluxes, especially for today’s high-permeance, high-selectivity membranes. Therefore, we built a 3-D CFD model that resolves local species transport, boundary-layer resistance, permeate-side effects, and pressure drop in flat-sheet modules. Using the Concentration Polarization Coefficient (CPC) and the overall Effective Membrane CP Index (EAC), the model quantifies how geometry and operating conditions shape CP intensity. Simulations show CP becomes pronounced for advanced membranes, with feed-channel height emerging as the dominant geometric lever: shrinking the gap mitigates CP but raises pressure drop, forcing a trade-off between mass-transfer enhancement and energy cost. These insights emphasize the need for CFD-level analysis in next-generation module design, where 1-D treatments can severely underestimate required membrane area and misguide scale-up.

carbon capture↗

DOE OSTI2024

Benchmark Calculation for BEAVRS Cycles 1 and 2 Using the SCALE 6.3/Polaris–PARCS v3.4.2 Code Package

This study performed the benchmark calculation for BEAVRS to validate the SCALE 6.3/Polaris–PARCS v3.4.2 code package with the ENDF/B-VII.1 56-group library, comparing the simulated results with the measured data. The benchmark results will be used in evaluating uncertainties of the SCALE/Polaris–PARCS code package for pressurized water reactor physics analysis for key nuclear parameters such as reactivity, control bank worth, temperature coefficients, and pin and assembly power peaking factors. This report details plant and fuel design specifications and input data for SCALE/Polaris as well as GenPMAXS and PARCS. Additionally, detailed information is provided for all the input and output files produced for the benchmark calculations. The benchmark results are summarized herein so that they can be used in evaluating uncertainties for key nuclear parameters with other benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Enhancing isotope mixing in U-10Mo downblend castings with electromagnetic stirring

Uranium alloyed with 10 wt percent molybdenum (U-10Mo) is the proposed fuel for use in the United States’ high-performance research reactors. The U-10Mo fuel is fabricated by a two-step casting process that downblends highly enriched uranium with depleted uranium (DU) and/or natural uranium (NU) and subsequently alloys the resulting high-assay low-enriched uranium (HALEU) with Mo. Currently, a two-step casting process is needed to meet the ingots’ U-235 enrichment homogeneity specifications. This work demonstrates the ability to provide more homogeneous U-235 distributions during downblend casting by using zirconia for crucibles rather than graphite. When graphite is used as a crucible, the electromagnetic field produced by the induction heater couples directly with the graphite. Zirconia is nonconducting in this study, and therefore does not couple with the field. The induction field couples directly with the metal, and causes electromagnetic stirring (EMS) in the molten metal pool. Eight downblend casting experiments were carried out in this work. Four were performed with zirconia crucibles and four with graphite crucibles. The U-235 enrichment was measured at nine discrete points in each casting using laser ablation multi-collector inductively coupled mass spectrometry (LA-MC-ICP-MS). The enrichment homogeneity was approximately eight times better in zirconia crucible castings than in graphite crucible castings as measured by the plate enrichment range and enrichment coefficient of variation. The results show a single casting step could be used to meet USHPRR enrichment specifications if electromagnetic stirring is present during casting. A related investigation was carried out to determine whether it is appropriate to measure enrichment on as-cast ingots or whether heat treated specimens should be used to accurately characterize enrichment. The results show that the as-cast microstructure plays a significant role in enrichment homogeneity. It is recommended that all enrichment measurements be done after the homogenization heat treatment typical of U-10Mo processing.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effect of inlet throttling on thermohydraulic instability in a large scale water-based RCCS: An experimental study

The objective of the present experimental study is to investigate the effect of inlet throttling on the thermohydraulic stability of a large scale water-based Reactor Cavity Cooling System (RCCS). The test was performed using the water-based Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne, which represented a ½ axial scale and 12.5° sector slice of the full scale Framatome 625 MW t SC-HTGR RCCS concept. A two-phase steady state was first established through direct condensate refill, followed by increased inlet throttling over 10 stages, corresponding to a loss coefficient K over the range of 0.05–653. With the inlet throttling gradually increased, the system experienced a unique transition process between stabilization and destabilization. Through a stability analysis, three instability mechanisms were identified in the present test, including a compound mechanism due to both natural circulation oscillations (NCOs) and density wave oscillations (DWOs), Type-II DWOs, and geysering.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Lithium inventory tracking as a non-destructive battery evaluation and monitoring method

Tracking the active lithium (Li) inventory in an electrode shows the true state of a Li battery, akin to a fuel gauge for an engine. However, non-destructive Li inventory tracking is currently unavailable. Here, in this work, we used the theoretical capacity of a transition metal oxide to convert capacity into a Li inventory analysis. The Li inventory in electrodes was tracked reliably to show how battery formulations and test methods affect performance. Contrary to capacity, Li inventory tracking reveals stoichiometric variations near the electrode–electrolyte interface. Verifiable results rationalized differences in measurements, clarifying and reducing interferences from cell formulations and experimental manipulations. By tracing four variables from formation to end-of-life, we characterize electrode and cell performance with a thermodynamic framework. Accurate rationalization of subtle differences in Li inventory utilization promises precise battery engineering, evaluation, failure analysis and risk mitigation. The method could be applicable from cell design optimization and fabrication to battery management, improving battery performance and reliability.

25 ENERGY STORAGE↗

DOE OSTI2022

Non-Equilibrium Doping of Poly-Si:Ga/SiOx Passivating Contacts Solar Cells

In this contribution, we show that Ga-hyperdoped poly-Si/SiOx contacts are realized by pulsed laser annealing and post-hydrogenation. The motivation behind this work is to improve the inferior performance of boron-doped poly- Si/SiOx contacts when compared to the phosphorus-doped contacts. Here we used Ga as a novel p-type dopant to replace B, as B is known to accumulate at the interface between the oxide and the crystalline silicon and cause recombination. Due to a high diffusivity of Ga in oxide and a much higher segregation coefficient (Si:SiOx), Ga is shown to be immune to the oxide pileup phenomenon. However, Ga has a low solid solubility limit in Si (5x10^19 cm-3), which is not enough to form a low resistivity contact. To overcome this problem, we explored pulsed laser annealing via excimer laser by varying the laser energy densities and pulses. The passivation and electrical properties are studied systematically using photoluminescence (PL), secondary ion mass spectrometry (SIMS), and Van der Pauw-Hall measurements. We show an active Ga doping concentration above 1x10^20 cm^-3, which is above the solid solubility limit of Ga in Si. Good passivation quality is demonstrated with the highest implied Voc (iVoc) of 721 mV on single side polished samples. Furthermore, a low diode contact resistivity of ~22 m..omega.. cm^2 is obtained using Ti/Ag contacts, with minimal metal-to-semiconductor contact resistivity of 0.9 m..omega.. cm^2.

non-equilibrium doping↗

DOE OSTI2025

Measurements of attenuation and scattering properties of water-based liquid scintillator

Water-based liquid scintillator (WbLS) is a hybrid detector medium which has been proposed as a fill material for future large-volume rare event searches, including ν − detectors. This family of scintillating suspensions promises waterlike attenuation and scattering behavior, while offering scintillation light yield for sub-Cherenkov-threshold events. These features may allow for improved vertex and energy resolution, and concomitant improvements in background rejection and particle identification. While subscale measurements of light yield, timing and pulse-shape have been performed in several WbLS formulations, detailed measurements of the attenuation and scattering properties of WbLS remain a critical, unresolved step along the pathway to deployment. In pursuit of a better understanding of these parameters, a “long-arm” attenuation and scattering instrument has been developed at Lawrence Livermore National Laboratory, dubbed LASE (Livermore Attenuation and Scattering Experiment). Optical property measurements have been performed using LASE in DI water, Gd-water, WbLS and Gd-WbLS. Measurements of the optical properties of WbLS (1% LAB-PPO scintillator) provided by Brookhaven National Laboratory demonstrate an attenuation minimum of 4.87 × 10 − 4 ± 1.42 × 10 − 4 cm − 1 at 450 nm, while Gd-loaded WbLS (1% LAB-PPO scintillator, 0.1% Gd) demonstrated a minimum attenuation coefficient of 4.77 × 10 − 4 ± 1.27 × 10 − 4 cm − 1 at the same wavelength. Measurements of scattering in WbLS show a scattering coefficient of 3.39 × 10 − 4 ± 1.98 × 10 − 5 cm − 1 at the 450 nm attenuation minimum, while Gd-WbLS has a scattering coefficient of 2.8 × 10 − 4 ± 1.63 × 10 − 5 cm − 1 at that wavelength. These scattering and attenuation coefficients are significantly larger than measured for DI water at similar wavelengths ( 8.8 × 10 − 5 cm − 1 ± 3.8 × 10 − 5 attenuation at 430 nm, 2.2 × 10 − 5 ± 1.28 × 10 − 6 cm − 1 scattering). Though the WbLS attenuation and scattering coefficients measured are lower than the corresponding values published for ultrapure oils and LAB-PPO liquid scintillator, the material tested attenuates significantly more than water across the blue and green portions of the visible spectrum. This attenuation is markedly stronger at wavelengths under 430 nm. Published by the American Physical Society 2025

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Cluster dynamics simulation of xenon diffusion during irradiation in UO 2

Diffusion of fission gas in UO 2 nuclear fuel impacts several important performance metrics, such as fission gas release, swelling, and thermal conductivity. Current empirical models of fission gas release have significant uncertainty, some of which derives from the bulk diffusion rate and its dependence on, for example, fuel chemistry and irradiation. In this work, we have applied the previously-developed Free Energy Cluster Dynamics (FECD) methodology in the code Centipede to calculate xenon cluster concentrations in UO 2 under intrinsic (high temperature) and irradiation-enhanced (intermediate temperature) conditions in order to develop a model of the xenon diffusion coefficient based on the atomic scale mechanisms responsible for transport. While the diffusion mechanism for xenon in UO 2 is adequately described by the Xe + U 2 O vacancy cluster for intrinsic conditions, a similar process is not capable of capturing measured in-pile fission gas diffusivity at intermediate temperatures. Therefore, a different diffusion mechanism must dominate under this regime. Using calculated atomistic data, we have shown that irradiation-enhanced diffusion at intermediate temperatures occurs via the larger Xe + U 4 O y vacancy clusters, which have lower migration barriers and increase in concentration by several orders of magnitude compared to intrinsic conditions. This mechanism is enabled by the increased uranium vacancy concentration under irradiation due to Frenkel pair production. In addition, the fast migration of uranium interstitials with two attached oxygen interstitials lowers the total uranium interstitial concentration through reactions with sinks. This allows the extended defects, such as Xe + U 4 O y vacancy clusters, to maintain high concentrations by limiting annihilation with attached vacancies. Furthermore, predictions using the Xe + U 4 O y diffusion mechanism are in good agreement with experiment, albeit with some differences in the Arrhenius slope, which we believe may be related to either experimental or model parameter uncertainty. Lastly, an analytical expression suitable for application in fuel performance simulations was derived to capture the predictions of the Centipede simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Measurement of Convective Heat Transfer Coefficients With Supercritical CO 2 Using the Wilson-Plot Technique

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in a smooth tube and compares the results with published correlations for validation. The paper also describes the Heat Exchange and Experimental Testing (HEET) rig recently designed and built at the U.S. Department of Energy’s (DoE’s) National Energy Technology Laboratory (NETL) in Morgantown, WV. Here, the Wilson-plot technique used for measuring the heat transfer coefficients is described along with the data reduction process. The Wilson-plot technique was chosen as the basis for the design of NETL’s HEET rig. Advantages of the Wilson-plot technique include the (1) ability to measure high convective heat transfer coefficients accurately, (2) ability to measure average heat transfer coefficient for complicated heat exchange geometries like those produced using additive manufacturing, (3) ability to measure heat transfer coefficients on both sides of a heat exchanger independently, and (4) simplicity of experimental setup. Capabilities of the HEET rig include pressure to 24 MPa (3500 psig), temperature to 538 °C (1000 °F), mass flow rate to 1.5 kg/s (3 lb/s), and Re to 500,000. The rig is designed to operate with pure CO 2 or a mixture of CO 2 and up to 10% N2 by volume to study the impact of gas mixtures typical of direct-fired sCO 2 power cycles on the convective heat transfer and pressure drop. Preliminary tests in the HEET rig were performed with smooth stainless-steel tube and pure CO 2 , and the results were compared with published correlations for Nusselt number (Nu) and friction factor. Over a Reynolds number (Re) range from 58,000 to 228,000, measured Nu was compared to predictions using the Dittus and Boelter equation (Kreith and Bohn, 1993, “Principles of Heat Transfer, West Publishing Company”) within 5% and measured friction factors were compared to predictions using the McAdams correlation (“McAdams, 1954, “Heat Transmission,” 3rd ed., McGraw Hill, New York)” for smooth tube to be within 5%.

42 ENGINEERING↗

DOE OSTI2020

Quantum Spectral Methods for Differential Equations

Recently developed quantum algorithms address computational challenges in numerical analysis by performing linear algebra in Hilbert space. Such algorithms can produce a quantum state proportional to the solution of a d-dimensional system of linear equations or linear differential equations with complexity poly(logd). While several of these algorithms approximate the solution to within ϵ with complexity poly(log(1/ϵ)), no such algorithm was previously known for differential equations with time-dependent coefficients. In this work, we develop a quantum algorithm for linear ordinary differential equations based on so-called spectral methods, an alternative to finite difference methods that approximates the solution globally. Using this approach, we give a quantum algorithm for time-dependent initial and boundary value problems with complexity poly(logd, log(1/ϵ)).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2021

The mobility of small, over-pressurized helium bubbles in tungsten at 2000 K

Fusion reactor environments inevitably lead to the formation of high-pressure helium bubbles whose nucleation, growth, and diffusion strongly impact the performance of plasma-facing components. Here this research describes a diffusion mechanism of over-pressurized bubbles via a sequence of Frenkel pair nucleation, self-interstitial migration, and Frenkel pair annihilation. Molecular dynamics was used to simulate the diffusion of small bubbles in tungsten at 2000 K with helium-per-vacancy ratios in the range of 4.5 to 7. The diffusion coefficients are calculated and their dependence on helium content, number of vacancies, and number of attached self-interstitials is characterized. It is found that bubbles are most mobile when the nucleation/annihilation rates of Frenkel pairs are nearly equal and when the bubbles nucleate and annihilate a single self-interstitial. All bubbles experience a peak diffusivity, which can be as high as 10 -11 m 2 /s decreasing with bubble size. The calculated diffusion coefficients provide valuable insight into the mobility of small, high-pressure bubbles, and can be used as input parameters in mesoscale models to improve predictions of plasma-surface interactions. (LA-UR-21-21881)

36 MATERIALS SCIENCE↗

DOE OSTI2025

Solvent Transport in Disordered and Dynamic Membrane Pores: Implications for Reverse Osmosis and Nanofiltration Membranes

Pressure-driven separations with nanoporous membranes, such as reverse osmosis and nanofiltration, play a vital role in addressing water scarcity and enabling resource recovery. Understanding water or solvent transport in membrane pores is essential for advancing membrane separation technologies. A key question in transport modeling is to establish a relationship between solvent permeability and membrane porous structure properties, such as porosity or pore size. The nano- and subnanometer pores in polymeric membranes such as reverse osmosis and nanofiltration membranes are highly tortuous and dynamically connected, which challenges the conventional methods of transport modeling. This study addresses this challenge by developing a theoretical framework to describe solvent transport through membranes with dynamic and disordered porous structure. Specifically, we propose a lattice model to describe the pore network, while preserving the viscous nature of solvent permeation. We further establish a relationship between solvent permeability and membrane porosity or pore size, which is validated by molecular dynamics simulations and experimental data. By integrating this relationship into the solution-friction model, we define pore connectivity and local friction coefficient to quantify the impact of pore structure on solvent permeability. Our analysis highlights the dominant influence of pore connectivity on the permeability of reverse osmosis and nanofiltration membranes, particularly when the pore size approaches the dimensions of solvent molecules. Overall, this study provides critical insights into water and solvent transport mechanisms in nanoporous membranes, opening the door for strategies to substantially enhance membrane performance.

membranes↗

DOE OSTI2024

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗