SEARCH · Engineering Papers

Results for “embedded boundary”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Embedded Ferroelectric Nanoclusters can drive Polarization Reversal in a Non-Ferroelectric Polar Film via the Proximity Effect

Heterogeneous nucleation from defects dominates the electric field required for polarization switching of ferroelectrics. Here, we consider the switching of a nominally non-switchable polar thin film of AlN due to the proximity effect arising from embedded ferroelectric nanoclusters of Al 1-x Sc x N. Using a Landau-Ginzburg-Devonshire thermodynamic approach and finite element modeling, we study the influence of nanocluster shape on polarization switching and domain nucleation emerging in AlN. The ferroelectric nanocluster boundary is modeled as a thin layer transitioning from Al 1-x Sc x N to AlN. We analyze the conditions under which polarization switching in the AlN film occurs at coercive fields significantly lower than its dielectric breakdown field. In the presence of spike-like Al 1-x Sc x N nanoclusters, the proximity effect enables switching of the spontaneous polarization in AlN and significantly reduces the corresponding coercive field. The internal field, which is depolarizing inside the AlN (due to its larger spontaneous polarization) and polarizing within the ferroelectric Al 1-x Sc x N nanoclusters (due to its smaller spontaneous polarization), lowers the potential barrier in the clusters and nucleates nanodomains at the Al 1-x Sc x N-AlN interface, forming localized regions of reversed polarization. Proximity effect can thus provide a pathway towards "thawing" previously "frozen" ferroelectrics through engineered nucleation for memory, actuation and optical technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗

DOE OSTI2022

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗

DOE OSTI2025

Temperature effects of ion irradiation on the nanostructural features in ductile-phase-toughened tungsten composites

Ductile-phase toughened tungsten (DPT W) composites have emerged as promising candidates for load-bearing components behind the plasma-facing tungsten armor in fusion reactors due to their enhanced thermomechanical properties. This study focuses on a composite consisting of W particles embedded in a ductile NiFeW solution matrix, hot-rolled to a thickness reduction of 87% (87R DPT W). Sequential irradiations with Ni2+ and He+ ions were performed to identical doses and helium concentrations at room temperature (RT) and 1273 K. Irradiation at RT produced no discernible nanostructural features due to the immobility of mono-vacancies, whereas cavity formation was observed at 973 K. At 1273 K, the W phase exhibited larger cavities, reduced cavity number density, and lower volumetric swelling compared to 973 K. Notably, nanosized NiFeW precipitates formed within the W phase at 1273 K, a phenomenon absent at 973 K. A new phase of cubic (NiFe)6W6C was also observed at the interphase boundary. In contrast, the NiFeW matrix showed no nanostructural changes at 1273 K, likely due to cavity dissociation. Separate irradiations at 1273 K indicated that Ni2+ ions induced precipitate formation in the W phase, while He+ ions exclusively caused cavity formation. The microstructure of 87R DPT W irradiated at RT and subsequently annealed at 1273 K closely resembled that of material irradiated directly at 1273 K. Like oxide-dispersion-strengthened steels, the observed nanoparticle-embedded W can inhibit dislocation propagation, potentially delaying the ductile-to-brittle transition temperature. These findings highlight the potential of NiFeW nanoparticle-reinforced W composites as irradiation-resistant materials for fusion reactors.

Jiang, Weilin (ORCID:0000000183028313)↗

DOE OSTI2022

Incorporating two-dimensional conduction modeling techniques into an energy simulation program: The EnergyPlus radiant system example

Due to the inherent complexity of multi-dimensional solutions for transient conduction in a multilayered building element, most heat balance-based algorithms assume that conduction occurs only in one dimension—through the building element between the inside and outside environments. Usually, since the length and width of the elements are much greater than the thickness of the element, this is a valid assumption. However, there are some cases such as hydronic radiant slabs where this assumption may not lead to an accurate estimate of the energy consumption of the system. In a hydronic radiant slab, water tubes are embedded in a slab, and hot or chilled water is circulated to meet any heating or cooling loads. The energy consumption of such a system will depend highly on the temperature of the water being sent to the slab (i.e., produced by the boiler or chiller). The temperature of the water that must be sent to the slab will depend on the spacing of the tubes, the multilayer construction of the building element, and the boundary conditions to which the construction is exposed. Thus, to accurately predict the energy consumption of such a system, a transient two-dimensional heat conduction model must be integrated with a heat balance-based simulation. This paper demonstrates the potential loss of accuracy associated with the one-dimensional radiant system models by providing a comparison of one- and two-dimensional solutions using EnergyPlus.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Hybrid Additive Manufacturing and Electric Field-Assisted Sintering of High-temperature Heat Exchangers via Sacrificial Channel Molds

A hybrid manufacturing approach, integrating additive manufacturing (AM) with powder methodology via electric field-assisted sintering (EFAS), was developed for the fabrication of high-temperature compact heat exchangers (CHX) from refractory metals. The methodology employed additively manufactured sacrificial channel molds (SCMs) as shapeholders for CHX channels, which were embedded in metal powders using EFAS. Following embedding, the SCMs were chemically dissolved to form the internal channel network. SCMs were fabricated using both digital light processing (DLP) and direct ink writing (DIW) from chemically reactive, calcium-based ceramic feedstocks with varying ratios of Al2O3 reinforcement. The microstructure, phase composition, and dissolution behavior of both as-printed and embedded SCMs were investigated. The shrinkage behavior of the SCMs embedded in refractory metals, as well as the interfacial characteristics between the SCMs and metal matrix, were studied. The results showed that the SCMs containing sufficient chemical reactive ceramics dissolved effectively before and after embedding. The as-printed SCMs retained the phase composition of their feedstocks, but the embedded SCMs containing calcium-based ceramics and Al2O3 exhibited the formation of calcium aluminates due to high temperature exposure during embedding. Most SCMs exhibited a cellular Al2O3 network filled with Ca-rich ceramics. Shrinkage after embedding was strongly dependent on SCM density, with lower density SCMs exhibiting greater shrinkage. A thin SCM-affected zone was observed at the metal matrix surface, characterized by increased porosity compared to the bulk matrix. This effect was attributed to infiltration of the SCM materials into powder particle boundaries under pressure, followed by their removal during dissolution. This study demonstrates the feasibility of manufacturing CHXs from hard-to-process refractory metals for use in harsh environments.

36 - MATERIALS SCIENCE↗

DOE OSTI2023

Extensive 3D mapping of dislocation structures in bulk aluminum

Abstract Thermomechanical processing such as annealing is one of the main methods to tailor the mechanical properties of materials, however, much is unknown about the reorganization of dislocation structures deep inside macroscopic crystals that give rise to those changes. Here, we demonstrate the self-organization of dislocation structures upon high-temperature annealing in a mm-sized single crystal of aluminum. We map a large embedded 3D volume ( $$100\times 300\times 300$$ 100 × 300 × 300 $$\upmu$$ μ m $$^3$$ 3 ) of dislocation structures using dark field X-ray microscopy (DFXM), a diffraction-based imaging technique. Over the wide field of view, DFXM’s high angular resolution allows us to identify subgrains, separated by dislocation boundaries, which we identify and characterize down to the single-dislocation level using computer-vision methods. We demonstrate how even after long annealing times at high temperatures, the remaining low density of dislocations still pack into well-defined, straight dislocation boundaries (DBs) that lie on specific crystallographic planes. In contrast to conventional grain growth models, our results show that the dihedral angles at the triple junctions are not the predicted 120 $$^{\circ }$$ ∘ , suggesting additional complexities in the boundary stabilization mechanisms. Mapping the local misorientation and lattice strain around these boundaries shows that the observed strain is shear, imparting an average misorientation around the DB of $$\approx$$ ≈ 0.003 to 0.006 $$^{\circ }$$ ∘ .

36 MATERIALS SCIENCE↗

DOE OSTI2025

Three-dimensional nucleation and growth of deformation twins in magnesium

At two-thirds the weight of aluminum, magnesium alloys have the potential to reduce the fuel consumption of transportation vehicles. These advancements depend on our ability to optimize the desirable versus undesirable effects of deformation twins, which are three-dimensional (3D) microstructural domains that form under mechanical stresses. Previously only characterized through surface or thin-film measurements, we here present 3D in situ characterization of deformation twinning inside an embedded grain over mesoscopic fields of view using dark-field x-ray microscopy supported by crystal plasticity finite element analysis. The results revealed the role of triple junctions on twin nucleation and the sequence and irregularity of twin growth and showed that twin-grain junctions, twin-twin junctions, and twin boundaries were the sites of localized dislocation accumulation.

Lee, Sangwon [Univ. of Michigan, Ann Arbor, MI (Un↗

DOE OSTI2025

Simulating Atmospheric Processes in Earth System Models and Quantifying Uncertainties With Deep Learning Multi‐Member and Stochastic Parameterizations

Abstract Deep learning is a powerful tool to represent subgrid processes in climate models, but many application cases have so far used idealized settings and deterministic approaches. Here, we develop stochastic parameterizations with calibrated uncertainty quantification to learn subgrid convective and turbulent processes and surface radiative fluxes of a superparameterization embedded in an Earth System Model (ESM). We explore three methods to construct stochastic parameterizations: (a) a single Deep Neural Network (DNN) with Monte Carlo Dropout; (b) a multi‐member parameterization; and (c) a Variational Encoder Decoder with latent space perturbation. We show that the multi‐member parameterization improves the representation of convective processes, especially in the planetary boundary layer, compared to individual DNNs. The respective uncertainty quantification illustrates that methods (b) and (c) are advantageous compared to a dropout‐based DNN parameterization regarding the spread of convective processes. Hybrid simulations with our best‐performing multi‐member parameterizations remained challenging and crash within the first days. Therefore, we develop a pragmatic partial coupling strategy relying on the superparameterization for condensate emulation. Partial coupling reduces the computational efficiency of hybrid Earth‐like simulations but enables model stability over 5 months with our multi‐member parameterizations. However, our hybrid simulations exhibit biases in thermodynamic fields and differences in precipitation patterns. Despite this, the multi‐member parameterizations enable improvements in reproducing tropical extreme precipitation compared to a traditional convection parameterization. Despite these challenges, our results indicate the potential of a new generation of multi‐member machine learning parameterizations leveraging uncertainty quantification to improve the representation of stochasticity of subgrid effects.

Behrens, Gunnar [Deutsches Zentrum für Luft‐ und R↗

DOE OSTI2026

Integrating science for water security governance

Hydrological extremes are intensifying globally, increasing the complexity of decisions required to ensure water security. Advances in hydrological science, modeling, and data systems have expanded the technical frontier of water research, yet uptake of scientific insights in policy and management decisions remains limited. This persistent science–policy gap is not primarily a failure of knowledge generation or robustness, but an institutional challenge shaped by how scientific and governance systems are organized, coordinated, and connected to support the effective use of scientific knowledge. These challenges are particularly pronounced in multi-level and transboundary water governance, where decisions span jurisdictions and require coordination across institutional and political boundaries. We synthesize research at the science–policy interface and evidence from water security initiatives to show how institutional arrangements, scientific tool development, and research practices enable or constrain the sustained use of scientific knowledge in water-security governance processes. Building on these insights, we develop ‘shared decision infrastructure’ as a framing to describe how scientific knowledge is embedded within the institutional, relational, and procedural arrangements that connect science to decision-making processes over time. We translate this framing into a practical intervention roadmap centered on institutional design, tool translation, sustained co-production, and outcome-oriented evaluation to support the integration of science into ongoing governance processes. By positioning science as shared decision infrastructure, the roadmap clarifies how researchers can design scientific efforts that support more coordinated, accountable, and adaptive water security decisions amid deepening uncertainty.

M whitney, Kristen [NASA Goddard Space Flight Cent↗

DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗

DOE OSTI2022

Summertime Marine Boundary Layer Cloud, Thermodynamic, and Drizzle Morphology over the Eastern North Atlantic: A Four-Year Study

Abstract Summertime remote sensor and in situ data from 2016 to 2019 collected at the ARM Eastern North Atlantic (ENA) Observatory are combined with aircraft measurements from the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign to quantify marine boundary layer (MBL) cloud, thermodynamic, and drizzle morphology in the region. A radar reflectivity–rainfall rate relationship ( Z – R ) is developed from aircraft data and 6-h cloud morphological regimes are identified from ENA data using a k -means algorithm driven by three independent inputs quantifying cloud thickness, drizzle intensity, and cloud field geometric complexity. Four separate MBL structural regimes representing non- or weakly drizzling single-layer stratocumulus, drizzling stratocumulus and cumulus-coupled stratocumulus, deep convection, and broken clouds embedded in northerly flow are identified. Single-layer stratocumulus is indicated when weak subtropical anticyclones are significantly west of the ENA site, and the MBL is cooler and drier than when drizzling and cumulus-coupled stratocumulus and broken clouds are observed. Drizzling and cumulus-coupled stratocumulus clouds are observed on the eastern flank of strong subtropical anticyclones in deep warm moist air masses with wind speeds exceeding 7 m s −1 and strong near-surface wind shear. Broken clouds exhibit strong wind shear near the inversion, while single-layer stratocumulus clouds have lower wind speeds and minimal shear. Net latent heat fluxes in the subcloud layer resulting from a combination of the ocean surface heat flux and evaporating drizzle average near zero over long periods in drizzling and cumulus-coupled stratocumulus. The ECMWF reanalysis version 5 (ERA5) is found to accurately represent single-layer stratocumulus properties, while producing significant discrepancies when drizzling stratocumulus and cumulus-coupled stratocumulus are observed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A high-energy and long-cycling lithium–sulfur pouch cell via a macroporous catalytic cathode with double-end binding sites

Lithium-sulfur batteries are attractive alternatives to lithium-ion batteries because of their high theoretical specific energy and natural abundance of sulfur. However, the practical specific energy and cycle life of Li-S pouch cells are significantly limited by the use of thin sulfur electrodes/flooded electrolytes and polysulfides shuttle/Li metal degradation, respectively. In this work we propose a cathode design concept to achieve good Li-S pouch cell performances. The cathode is composed of uniformly embedded ZnS nanoparticles and Co-N-C single-atom catalyst to form double-end biding sites inside a highly oriented macroporous host, which can effectively immobilize and catalytically convert polysulfides intermediates during cycling, thus eliminating the shuttle effect and lithium metal corrosion. The ordered macropores enhance ionic transport under high sulfur loading by forming sufficient catalyst/conductive support/electrolyte triple-phase boundaries. This prevents the formation of inactive sulfur (dead sulfur). Our cathode structure shows improved performances in a pouch cell configuration under high sulfur loading and lean electrolyte operation. An Ah-level pouch cell with only 100% lithium excess can deliver a cell specific energy of > 300 Wh kg –1 with a Coulombic efficiency >95% for 80 cycles.

25 ENERGY STORAGE↗

DOE OSTI2020

Does the rotational direction of a wind turbine impact the wake in a stably stratified atmospheric boundary layer?

Stably stratified atmospheric boundary layers are often characterized by a veering wind profile, in which the wind direction changes clockwise with height in the Northern Hemisphere. Wind-turbine wakes respond to this veer in the incoming wind by stretching from a circular shape into an ellipsoid. We investigate the relationship between this stretching and the direction of the turbine rotation by means of large-eddy simulations. Clockwise rotating, counterclockwise rotating, and non-rotating actuator disc turbines are embedded in wind fields of a precursor simulation with no wind veer and in wind fields with a Northern Hemispheric Ekman spiral, resulting in six combinations of rotor rotation and inflow wind condition. The wake strength, extension, width, and deflection depend on the interaction of the meridional component of Ekman spiral with the rotational direction of the actuator disc, whereas the direction of the disc rotation only marginally modifies the wake if no veer is present. The differences result from the amplification or weakening/reversion of the spanwise and the vertical wind components due to the effect of the superposed disc rotation. They are also present in the streamwise wind component of the wake and in the total turbulence intensity. In the case of an counterclockwise rotating actuator disc, the spanwise and vertical wind components increase directly behind the rotor, resulting in the same rotational direction in the whole wake while its strength decreases downwind. In the case of a clockwise rotating actuator disc, however, the spanwise and vertical wind components of the near wake are weakened or even reversed in comparison to the inflow. This weakening/reversion results in a downwind increase in the strength of the flow rotation in the wake or even a different rotational direction in the near wake in comparison to the far wake. The physical mechanism responsible for this difference can be explained by a simple linear superposition of a veering inflow with a Rankine vortex.

17 WIND ENERGY↗

DOE OSTI2024

Thermal partition function of $$ {J}_3{\overline{J}}_3 $$ deformed AdS3

Abstract We derive a compact formula for the one-loop, bosonic string partition function of Euclideanized$$ {J}_3{\overline{J}}_3 $$ J 3 J ¯ 3 deformedAdS 3 with periodic Euclidean time as an integral transform of the partition function of the undeformed EuclideanizedAdS 3 . Such a deformation is interpretable as an irrelevant “single-trace$$ T\overline{T} $$ T T ¯ deformation” of the boundary. We will do this by first establishing a formal procedure to compute a worldsheet torus zero point function for an exactly marginal$$ J\overline{J} $$ J J ¯ deformation of a sigma model with U(1) L × U(1) R global symmetry. We then describe how this procedure is implemented on SL(2,R) sigma model and its Euclidean continuation. Finally, we describe the embedding of the deformed SL(2,R) torus amplitude into critical string theory and interpret the result as the leading perturbative contribution to the thermal partition function of the deformed theory.

Physics↗

DOE OSTI2022

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗