Search NASA⌕ Search

SEARCH · Search NASA

Results for “Problem structuring methods”

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

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

At least 181 records · Page 10

Sparse Bayesian mass mapping with uncertainties: hypothesis testing of structure

ABSTRACT A crucial aspect of mass mapping, via weak lensing, is quantification of the uncertainty introduced during the reconstruction process. Properly accounting for these errors has been largely ignored to date. We present a new method to reconstruct maximum a posteriori (MAP) convergence maps by formulating an unconstrained Bayesian inference problem with Laplace-type l1-norm sparsity-promoting priors, which we solve via convex optimization. Approaching mass mapping in this manner allows us to exploit recent developments in probability concentration theory to infer theoretically conservative uncertainties for our MAP reconstructions, without relying on assumptions of Gaussianity. For the first time, these methods allow us to perform hypothesis testing of structure, from which it is possible to distinguish between physical objects and artefacts of the reconstruction. Here, we present this new formalism, and demonstrate the method on simulations, before applying the developed formalism to two observational data sets of the Abell 520 cluster. Initial reconstructions of the Abell 520 catalogues reported the detection of an anomalous ‘dark core’ – an overdense region with no optical counterpart – which was taken to be evidence for self-interacting dark matter. In our Bayesian framework, it is found that neither Abell 520 data set can conclusively determine the physicality of such dark cores at $99{{\ \rm per\ cent}}$ confidence. However, in both cases the recovered MAP estimators are consistent with both sets of data.

Price, M. A.↗

Environmental Durability Assessment and Development of Subcomponent Test Methods for Advanced Materials in Marine Energy Applications

Marine energy developers are beginning to rely more heavily on advanced composite materials for their load-bearing structures. Understanding how these materials respond and degrade in harsh marine environments is a complex problem. Much research has been performed at smaller laboratory scales to understand the fundamental physical and chemical degradation mechanisms, but little structural testing research has been performed at larger scales to understand how these physical and chemical interactions scale and how they should be accounted for in the design process. The Advanced Materials project is a multilaboratory research collaboration that addresses these issues with the intention of reducing risks and barriers to marine energy developers wanting to adopt advanced composite materials. The National Laboratory of the Rockies has been playing a leading role in developing larger-scale testing infrastructure and test methods to validate complex numerical degradation models and standardize testing procedures to do so. This paper reviews past structural testing activities, how they have grown into current research activities, and how they will continue to evolve to directly address industry needs. This paper covers key research areas - large-scale water absorption and subsequent structural validation testing, qualification of bolted and adhesive connections, and combined submerged fatigue testing at increasing scales. It provides insights into the complex requirements for testing infrastructure and test method development to effectively qualify materials and certify marine energy structures for harsh marine environments. This research program paves the way for better guidance and standardization within the industry when adopting advanced composite materials.

16 TIDAL AND WAVE POWER↗

Environmental Durability Assessment and Development of Subcomponent Test Methods for Advanced Materials in Marine Energy Applications: Preprint

Marine energy developers are beginning to rely more heavily on advanced composite materials for their load bearing structures. Understanding how these materials respond and degrade in harsh marine environments is a complex problem. Much research has been performed at smaller laboratory scales to understand the fundamental physical and chemical degradation mechanisms, but little structural testing research has been performed at larger scales to understand how these physical and chemical interactions scale and how they should be accounted for in the design process. The Advanced Materials project is a multilaboratory research collaboration that has been addressing these issues with the intention of reducing risks and barriers to marine energy developers wanting to adopt advanced composite materials. The National Laboratory of the Rockies has been playing a leading role in developing larger scale testing infrastructure and test methods to validate complex numerical degradation models and standardize testing procedures to do so. This paper provides a review of past structural testing activities and how they have grown into current research activities and will continue to evolve into the future to directly address industry needs. This paper covers key research areas large-scale water absorption and subsequent structural validation testing, qualification of bolted and adhesive connections, and conducting combined submerged fatigue testing at increasing scales. This paper provides insights into the complex requirements for testing infrastructure and test method development to effectively qualify materials and certify marine energy structures for harsh marine environments. This research program paves the way for better guidance and standardization within the industry when adopting advanced composite materials.

16 TIDAL AND WAVE POWER↗

A Novel Framework for Optimizing Ramping Capability of Hybrid Energy Storage Systems

Hybrid Energy Storage System has been widely applied in aerospace, electric vehicle, and microgrid applications. The advantages are that they include complimentary technologies with both high power and energy capabilities. HESS have the potential to be useful to the bulk power systems, for example to increase the value of energy produced by variable generation resource through enabling participation in ancillary service markets. Actualizing the benefits of HESSs requires optimizing the combined ramping capability of aggregated HESS resources. To provide quality ancillary service, source of HESSs locating at multiple sites of variable generation resource must be optimally sized and cohesively controlled. Optimizing aggregated resources can be formulated as an optimization problem. Successfully solving this problem not only requires using effective solution methods but also depends on accurately and rapidly setting the necessary parameters in the problem formula. This article proposes a novel framework with double-layer structure to solve the optimization problem of aggregating ramping capability. Program developed on the upper layer focuses on solving the optimization of aggregating ramping capability among multiple HESSs and is compatible with most existing optimization algorithms. Program on lower layer of the framework targets at optimizing the local control of a single HESS based on a thorough analytics of HESS operation strategy presented in this article. Result of local optimization is also provided to upper layer program for updating the parameters in formula of optimization problem. Real time hardware-in-the-loop test is conducted to verify the performance of the optimization framework developed. The work presented is expected to provide guidance for implementing this framework in practical operation.

Luo, Yushen↗

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗

Computationally Efficient Decompositions of Oblique Projection Matrices

Oblique projection matrices arise in problems in weighted least squares, signal processing, and optimization. While these matrices can be potentially very large, their low-rank structure can be exploited for efficient computation. Here, we propose fast and scalable algorithms for computing their eigendecomposition and singular value decomposition (SVD). Numerical experiments that compare our proposed approaches to existing methods, including randomized SVD, are presented. In addition, we test their accuracy on linear systems from equality constrained optimization problems.

97 MATHEMATICS AND COMPUTING↗

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗

Computational design of metamaterials with self contact

Inverse homogenization in combination with contact modeling, topology optimization and shape optimization is used to design metamaterials with optimized macroscopic response. The homogenization assumes length scale separation which allows the non-linear macroscopic behavior to be obtained by analyzing a single unit cell in a lattice structure. Self contact in the unit cell, which is modeled using a third medium contact method, is leveraged to obtain a complex homogenized response. The inverse homogenization problem is initially formulated as a topology optimization problem, where the macroscopic stress–strain behavior is tuned to our liking. However, it is well known that boundary phenomena are difficult to model in topology optimization and that interface modeling is crucial to accurately analyze contact. For that reason, the boundary representation of the topology optimized design is extracted and used as initial design in a subsequent shape optimization. The behaviors of our designs are verified by performing rigorous post-processing analyzes using conforming meshes and conventional contact formulations.

42 ENGINEERING↗

Deep Learning Segmentation of Complex Features in Atomic-Resolution Phase-Contrast Transmission Electron Microscopy Images

Phase-contrast transmission electron microscopy (TEM) is a powerful tool for imaging the local atomic structure of materials. TEM has been used heavily in studies of defect structures of two-dimensional materials such as monolayer graphene due to its high dose efficiency. However, phase-contrast imaging can produce complex nonlinear contrast, even for weakly scattering samples. It is, therefore, difficult to develop fully automated analysis routines for phase-contrast TEM studies using conventional image processing tools. For automated analysis of large sample regions of graphene, one of the key problems is segmentation between the structure of interest and unwanted structures such as surface contaminant layers. In this study, we compare the performance of a conventional Bragg filtering method with a deep learning routine based on the U-Net architecture. We show that the deep learning method is more general, simpler to apply in practice, and produces more accurate and robust results than the conventional algorithm. We provide easily adaptable source code for all results in this paper and discuss potential applications for deep learning in fully automated TEM image analysis.

47 OTHER INSTRUMENTATION↗

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗

Sparse thermal data for cellular automata modeling of grain structure in additive manufacturing

Grain growth in the wake of the melt pool formed during alloy-based additive manufacturing (AM) is complex and multifaceted, depending on parameters governing heat transport, fluid flow, and solidification itself. Cellular automata (CA) models have proven effective in providing computationally efficient and physically sound predictions of grain structure for several AM problems, but their efficiency is tied to the performance of heat transport models. CA models use only a small portion of the problem's temperature data (near the moving melt pool boundary), and much of the CA calculations do not affect the final result due to re-melting of material. Coupling of and communication between heat transport and solidification models, and eliminating operations irrelevant towards final grain structure prediction, will be necessary for using these methods for efficient simulation of large parts. In this work, we introduce a procedure of decoupling the CA from temperature field simulation, using files of relevant temperature data written by the heat transport model. This approach is validated against the standard coupling approach using data obtained through the computational fluid dynamics software OpenFOAM. Negligible differences are seen in grain size, volume, and texture distributions for multilayer simulation of test problems, while the quantity of temperature data for these test problems was reduced by four orders of magnitude (from 100s of GB to 10s of MB) and the code performance sped up by a factor of around 50. Furthermore, variability in microstructure as a function of cell size, substrate, time step, and nucleation parameters is studied, and it is found that cell sizes less than or equal to 1.67 μm and sufficiently small time steps yield Αstatistically equivalent microstructures. Finally, a potential use case for this CA approach—the layer-wise convergence in grain structure starting from extremes in initial grain size—is examined. This approach's ability to simulate expected trends in nucleation and epitaxial grain growth for large regions of microstructure, simulated independently of heat transport models themselves, should prove useful for investigation of various microstructure uncertainties and prediction of part-scale experimental results.

36 MATERIALS SCIENCE↗

Implementation of Sacrificial Support Structures for Hybrid Manufacturing of Thin Walls

Thin-walled features can be difficult to produce with traditional machining methods which often rely on excess stock material for stiffness. This challenge is increased in hybrid manufacturing where the feature is already near net shape before machining. Significant workpiece deflection can result in poor geometric and surface finish tolerances on the finished part. A potential solution to this problem is to implement sacrificial support structures to the as-printed geometry. The supports are then machined away during the finishing portion of the hybrid process. In the present work, several different design parameters for these sacrificial supports were evaluated to determine their impact on the quality of representative thin wall geometry samples. The angle, height, and spacing of triangular support structures were varied for each sample and then machined and examined. The addition of these supports relative to an unsupported configuration provided a deflection reduction of around 0.2 mm. Surface roughness was improved by approximately 1.5 µm. Increasing values of support height were found to correspond to reduced wall deflection. Similarly, decreasing values of support angle and support spacing improved geometric accuracy. Efficiency comparisons showed that increases in print time corresponded to rapidly diminishing gains in geometric accuracy but continued to improve surface roughness. Implications for hybrid finishing of additively manufactured thin-walled structures is briefly discussed.

36 MATERIALS SCIENCE↗

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Atomic-scale Imaging of PGM-free Catalyst Active Sites by 30 keV 4D-STEM

Platinum group metal (PGM)-free catalysts have attracted a large amount of attention due to their potential for enabling low-cost, commercially viable hydrogen fuel cells and electrolyzers [1]. Degradation of cell performance remains a significant challenge, however, currently limiting implementation of devices that utilize PGM-free materials [2]. Controlling degradation in these materials requires a better understanding of the atomic structure of active sites, thought to be FeN4 structures within a graphitic carbon lattice, which would enable more accurate prediction of the associated degradation pathways. The exact structural arrangement of active sites in these materials is still debated since a range of structures have been computationally predicted and methods to directly validate these models are still needed [1,3-5].While scanning transmission electron microscopy (STEM) has provided initial glimpses into the nature of the proposed active sites, detailed evaluation of the structure of these sites is challenging since they involve defects, edges, and nitrogen dopants in the graphitic lattice, which are susceptible to beam damage [6,7]. Conventional STEM imaging and spectroscopy methods exacerbate this problem with dose-inefficiency or nonideal contrast characteristics for imaging light and heavy elements simultaneously. A method that minimizes damage while generating dose-efficient, easily-interpretable contrast for both light and heavy elements is therefore needed to facilitate imaging of sensitive active site atomic structure.Here, we demonstrate direct atomic-scale imaging of PGM-free catalyst active sites by low-voltage four-dimensional (4D)-STEM. To accomplish this, we pair a 30 keV aberration-corrected probe, which minimizes knock-on damage, with a fast pixelated detector that has optimal performance at low beam energies [8]. We show how this setup enables relatively simple center-of-mass (CoM) techniques [9,10] to produce images with an increased signal-to-noise ratio (SNR) and light-element contrast over conventional imaging modes, allowing the entire active site structure to be imaged at the atomic scale. Moreover, we show how electron ptychography [11,12] enables images with further improved characteristics to be obtained, for example by minimizing residual aberrations, which provides a more accurate structural representation. The increased understanding that these low-voltage 4D-STEM techniques will provide about the atomic structure of PGM-free active sites and their associated degradation pathways will facilitate rational design of next-generation materials, promoting development of low-cost hydrogen fuel cells, electrolyzers, and other energy conversion devices [13]. References:[1] U Martinez et al., Adv. Mater. 31 (2019), p. 1806545.[2] Y Shao et al., Adv. Mater. 31 (2019), p. 1807615.[3] J Kneebone et al., J. Phys. Chem. C 121, (2017), p.16283.[4] T Mineva et al., ACS Catal. 9 (2019), p. 9359.[5] J Li et al., Nat. Catal. 4 (2021), p. 10.[6] T Susi et al., ACS Nano 6 (2012), p. 8837.[7] P Zelenay and DJ Myers, US Department of Energy Hydrogen and Fuel Cells Program 2017 Annual Merit Review and Peer Evaluation Meeting, Washington, DC (2017).[8] H Ryll et al., J. Inst. 11 (2016), p. P04006.[9] K Müller et al., Nat. Commun. 5 (2014), p. 5653.[10] I Lazic et al., Ultramicroscopy 160 (2016), p. 265.[11] H Yang et al., Nat. Commun. 7 (2016), p. 12532.[12] Y Jiang et al., Nature 559 (2018), p. 343. [13] Research sponsored by the Hydrogen and Fuel Cell Technologies Office, Office of Energy Efficiency and Renewable Energy, US Department of Energy (DOE). Research was conducted at the Center for Nanophase Materials Sciences, which is a DOE Office of Science User Facility.

Zachman, Michael↗

Transfer of Millimeter‐Scale Strained Multiferroic Epitaxial Thin Films on Rigid Substrates via an Epoxy Method Producing Magnetic Property Enhancement

The demonstration of epitaxial thin film transfer has enormous potential for thin film devices free from the traditional substrate epitaxy limitations. However, large-area continuous film transfer remains a challenge for the commonly reported polymer-based transfer methods due to bending and cracking during transfer, especially for highly strained epitaxial thin films. In this work, a new epoxy-based, rigid transfer method is used to transfer films from an SrTiO 3 (STO) growth substrate onto various new substrates, including those that will typically pose significant problems for epitaxy. An epitaxial multiferroic Bi 3 Fe 2 Mn 2 O x (BFMO) layered supercell (LSC) material is selected as the thin film for this demonstration. The results of surface and structure studies show an order of magnitude increase in the continuous area of transferred films when compared to previous transfer methods. The magnetic properties of the BFMO LSC films are shown to be enhanced by the release of strain in this method, and ferromagnetic resonance is found with an exceptionally low Gilbert damping coefficient. The large-area transfer of this highly strained complex oxide BFMO thin film presents enormous potential for the integration of many other multifunctional oxides onto new substrates for future magnetic sensors and memory devices.

bismuth oxides↗

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗