Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance map”

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

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

At least 1,657 records · Page 92

Weakly Supervised Minirhizotron Image Segmentation with MIL-CAM

We present a multiple instance learning class activation map (MIL-CAM) approach for pixel-level minirhizotron image segmentation given weak image-level labels. Minirhizotrons are used to image plant roots in situ. Minirhizotron imagery is often composed of soil containing a few long and thin root objects of small diameter. The roots prove to be challenging for existing semantic image segmentation methods to discriminate. In addition to learning from weak labels, our proposed MILCAM approach re-weights the root versus soil pixels during analysis for improved performance due to the heavy imbalance between soil and root pixels. Furthermore, the proposed approach outperforms other attention map and multiple instance learning methods for localization of root objects in minirhizotron imagery.

97 MATHEMATICS AND COMPUTING↗

Airframe Noise from a Hybrid Wing Body Aircraft Configuration

A high fidelity aeroacoustic test was conducted in the NASA Langley 14- by 22-Foot Subsonic Tunnel to establish a detailed database of component noise for a 5.8% scale HWB aircraft configuration. The model has a modular design, which includes a drooped and a stowed wing leading edge, deflectable elevons, twin verticals, and a landing gear system with geometrically scaled wheel-wells. The model is mounted inverted in the test section and noise measurements are acquired at different streamwise stations from an overhead microphone phased array and from overhead and sideline microphones. Noise source distribution maps and component noise spectra are presented for airframe configurations representing two different approach flight conditions. Array measurements performed along the aircraft flyover line show the main landing gear to be the dominant contributor to the total airframe noise, followed by the nose gear, the inboard side-edges of the LE droop, the wing tip/LE droop outboard side-edges, and the side-edges of deployed elevons. Velocity dependence and flyover directivity are presented for the main noise components. Decorrelation effects from turbulence scattering on spectral levels measured with the microphone phased array are discussed. Finally, noise directivity maps obtained from the overhead and sideline microphone measurements for the landing gear system are provided for a broad range of observer locations.

Hutcheson, Florence V.↗

Towards Generic Parallel Programming in Computer Science Education with Kokkos

Parallel patterns, views, and spaces are promising abstractions to capture the programmer's intent as well as the contextual information that can be used by an underlying runtime to efficiently map software to parallel hardware. These abstractions can be valuable in cases where an algorithm must accommodate requirements of code and performance portability across hardware architectures and vendor programming models. Kokkos is a parallel programming model for host- and accelerator architectures that relies on these abstractions and targets these requirements. It consists of a pure C++ interface, a specification, and a programming library. The programming library exposes patterns and types and maps them to an underlying abstract machine model. The abstract machine model offers a generic view of parallel hardware. While Kokkos is gaining popularity in large-scale HPC applications at some DOE laboratories, we believe that the implemented concepts are of interest to a broader audience including academia as they may contribute to a generic, vendor, and architecture-independent education of parallel programming. In this work, we give an insight into the design considerations of this programming model and list important abstractions. Further, we document best practices obtained from giving virtual classes on Kokkos and give pointers to resources that the reader may consider valuable for a lecture on generic parallel programming for students with preexisting knowledge on this matter.

Ciesko, Jan↗

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.↗

Integrating Spatial and Ethnographic Methods for Resilience Research: A Thick Mapping Approach for Hurricane Maria in Puerto Rico

Hurricane Maria left unprecedented impacts on Puerto Rican communities, leaving some without infrastructure services and unable to communicate with family for several months. Here, to understand the forms of community-level resilience that emerged while hard infrastructure systems took time recover, this article (1) abductively explores resilience as an emergent phenomenon of complex adaptive systems; (2) identifies subsequent forms of social capital, local adaptive capacities, and manifestations of quantifiable variables, such as infrastructure performance, in community experiences; and (3) demonstrates a framework to integrate disparate methodologies for resilience assessments via a multiplicity of mappings of space and place. We combine ethnographic and geospatial methods into an interactive GeoApp for analysis using participant-coded narratives and a series of geospatial indicators as a thick map. Thick mapping facilitates quantitative and qualitative data analysis at several scales, while enabling qualitative query of collected narratives. Results highlight local innovation, community bonding and bridging, and nuances in the role of public institutions as emergent elements of resilience. The thick map shows how top-down assessments can be augmented by thick data and how multiple framings can be anchored in the same system or place. These findings are important to inform and integrate community-oriented and technocentric solutions toward resilience-enhancing measures.

54 ENVIRONMENTAL SCIENCES↗

Domain connectivity among systems of overset grids

This progress report describes research performed to achieve algorithmic improvements in the use of overset grids for the computation of unsteady flowfields about geometrically complex and moving component configurations. A hole-map computational method is detailed that reduces computational demand and corresponding demand on human resources. Results indicate that the method is both efficient and flexible for use in complex simulation cases. Tests were performed on shuttle orbiter/external tank, wing/store, and V-22 tilt-rotor configurations. The number of simulated points and the corresponding time on a SGI 4D-210 workstation are presented.

Meakin, Robert L.↗

Parasitic Interference in Nulling Interferometry

Nulling interferometry aims to detect faint objects close to bright stars. Its principle is to produce a destructive interference along the line of sight so that the stellar flux is rejected, while the flux of the off-axis source can be transmitted. In practice, various instrumental perturbations can degrade the nulling performance. Any imperfection in phase, amplitude or polarization produces a spurious flux that leaks to the interferometer output and corrupts the transmitted off-axis flux. One of these instrumental perturbations is the crosstalk phenomenon, which occurs because of multiple parasitic reflections inside transmitting optics, and/or diffraction effects related to beam propagation along finite size optics. It can include a crosstalk of a beam with itself, and a mutual crosstalk between different beams. This can create a parasitic interference pattern, which degrades the intrinsic transmission map - or intensity response - of the interferometer. In this context, we describe how this instrumental effect impairs the performance of a Bracewell interferometer. A simple formalism is developed to derive the corresponding modified intensity response of the interferometer, as a function of the two parameters of interest: the crosstalk level (or contamination rate) and the phase shift between the primary and secondary - parasitic - beams. We then apply our mathematical approach to a few scientific cases, both analytically and using the GENIESIM simulation software, adapted to handle coherent crosstalk. Our results show that a coherent crosstalk level of about 1 per cent implies a 20 per cent drop of the signal-to-noise ratio at most. Careful attention should thus be paid to reduce the crosstalk level inside an interferometric instrument and ensure an instrumental stability that provides the necessary sensitivity through calibration procedures.

circumstellar matter↗

A Code Generator for High-Performance Tensor Contractions on GPUs

Tensor contractions are higher dimensional generalizations of matrix-matrix multiplication. They form the compute-intensive core for many applications in computational science and data science. Although high-performance library implementations are widely available for dense matrix-matrix multiplication on GPUs, the same is not true for contraction of dense tensors. In this paper, we describe a high- performance GPU code generator for arbitrary tensor contractions. It exploits domain-specific properties about data reuse in tensor contractions to devise an effective code generation schema, coupled with an effective model-driven search to determine parameters for mapping of computation to threads and staging of data through the GPU memory hierarchy. Experimental evaluation using a set of tensor contraction benchmarks demonstrates performance improvement and/or significantly reduced code generation time over other state-of-the-art tensor contraction libraries and code generators.

Kim, Jinsung↗

Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies

This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of a supercharged direct-injection stratified-charge rotary combustion engine

A zero-dimensional thermodynamic performance computer model for direct-injection stratified-charge rotary combustion engines was modified and run for a single rotor supercharged engine. Operating conditions for the computer runs were a single boost pressure and a matrix of speeds, loads and engine materials. A representative engine map is presented showing the predicted range of efficient operation. After discussion of the engine map, a number of engine features are analyzed individually. These features are: heat transfer and the influence insulating materials have on engine performance and exhaust energy; intake manifold pressure oscillations and interactions with the combustion chamber; and performance losses and seal friction. Finally, code running times and convergence data are presented.

Bartrand, Timothy A.↗

Field validation of data-driven BSDF and peak extraction models for light-scattering fabric shades

Shading and daylighting systems affect cooling, heating, and lighting energy use by modulating solar radiation through the building façade. Characterizing shading systems holistically and accurately helps designers and engineers evaluate shading systems to achieve energy and non-energy performance goals. These complex fenestration systems can be modeled using Bidirectional Scattering Distribution Functions (BSDF), which map incident radiation to hemispherical distributions of outgoing radiation. Data-driven, tabulated BSDFs are derived from interpolated goniophotometer measured data, then sampled during the raytracing calculation. A peak extraction (PE) algorithm was developed to circumvent limits in BSDF angular resolution, where the specular peak is extracted during simulation by evaluating the BSDF in the through direction and surrounding region. The objective of this study was to validate this measurement and modeling workflow using field monitored data from a full scale testbed with eleven installed fabrics of different weaves, openness factors, and colors and assess the accuracy of the workflow under different adaptation and contrast conditions. Test conditions were limited to clear sky conditions with the sun in the field of view. Results showed that, for tensor tree datasets, vertical illuminance, solar luminance (2.5° apex), and daylight glare probability (DGP) were predicted to within a mean bias error (MBE) error of -456 lx (-12.3%), -3.46e5 (-38.4%), and -0.042 (-7.8%) when full PE occurred. With a binary classification of glare/ no glare, DGP was predicted accurately with a true positive rate of 0.98 and true negative rate of 1.0 using tensor tree data and less accurately with Klems BSDF data, particularly for cases of no glare. The workflow may be of insufficient accuracy to distinguish borderline performance between fabrics using the four-point glare scale, particularly under low adaptation, high contrast daylit conditions. Errors were due to reductions in peak shape and intensity across the BSDF interpolation and data reduction workflow. Future work is needed to better preserve measurement fidelity during interpolation and sampling, which in turn will improve PE performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mapping microhabitats of lignocellulose decomposition by a microbial consortium

The leaf-cutter ant fungal garden ecosystem is a naturally evolved model system for efficient plant biomass degradation. Degradation processes mediated by the symbiotic fungus Leucoagaricus gongylophorus are difficult to characterize due to dynamic metabolisms and spatial complexity of the system. Herein, we performed microscale imaging across 12-µm-thick adjacent sections of Atta cephalotes fungal gardens and applied a metabolome-informed proteome imaging approach to map lignin degradation. This approach combines two spatial multiomics mass spectrometry modalities that enabled us to visualize colocalized metabolites and proteins across and through the fungal garden. Spatially profiled metabolites revealed an accumulation of lignin-related products, outlining morphologically unique lignin microhabitats. Metaproteomic analyses of these microhabitats revealed carbohydrate-degrading enzymes, indicating a prominent fungal role in lignocellulose decomposition. Integration of metabolome-informed proteome imaging data provides a comprehensive view of underlying biological pathways to inform our understanding of metabolic fungal pathways in plant matter degradation within the micrometer-scale environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

International Low-Earth-Orbit Spacecraft Materials Test Program Initiated for Better Prediction of Durability and Performance

Spacecraft in low Earth orbit (LEO) are subjected to many components of the environment, which can cause them to degrade much more rapidly than intended and greatly shorten their functional life. The atomic oxygen, ultraviolet radiation, and cross contamination present in LEO can affect sensitive surfaces such as thermal control paints, multilayer insulation, solar array surfaces, and optical surfaces. The LEO Spacecraft Materials Test (LEO-SMT) program is being conducted to assess the effects of simulated LEO exposure on current spacecraft materials to increase understanding of LEO degradation processes as well as to enable the prediction of in-space performance and durability. Using ground-based simulation facilities to test the durability of materials currently flying in LEO will allow researchers to compare the degradation evidenced in the ground-based facilities with that evidenced on orbit. This will allow refinement of ground laboratory test systems and the development of algorithms to predict the durability and performance of new materials in LEO from ground test results. Accurate predictions based on ground tests could reduce development costs and increase reliability. The wide variety of national and international materials being tested represent materials being functionally used on spacecraft in LEO. The more varied the types of materials tested, the greater the probability that researchers will develop and validate predictive models for spacecraft long-term performance and durability. Organizations that are currently participating in the program are ITT Research Institute (USA), Lockheed Martin (USA), MAP (France), SOREQ Nuclear Research Center (Israel), TNO Institute of Applied Physics (The Netherlands), and UBE Industries, Ltd. (Japan). These represent some of the major suppliers of thermal control and sensor materials currently flying in LEO. The participants provide materials that are exposed to selected levels of atomic oxygen, vacuum ultraviolet radiation, contamination, or synergistic combined environments at the NASA Lewis Research Center. Changes in characteristics that could affect mission performance or lifetime are then measured. These characteristics include changes in mass, solar absorptance, and thermal emittance. The durability of spacecraft materials from U.S. suppliers is then compared with those of materials from other participating countries. Lewis will develop and validate performance and durability prediction models using this ground data and available space data. NASA welcomes the opportunity to consider additional international participants in this program, which should greatly aid future spacecraft designers as they select materials for LEO missions.

Rutledge, Sharon K.↗

Comparison of level I land cover classification accuracy for MSS and AVHRR data

The capabilities of the Advanced Very-High-Resolution Radiometer (AVHRR) for land-cover mapping were investigated by comparing the accuracy of land-cover information for the Washington, DC area derived from NOAA-7 AVHRR data with that from Landsat Multispectral Scanner Subsystem (MSS) data. Unsupervised level I land-cover classifications were performed for MSS and AVHRR data sets collected on July 11, 1981. A detailed accuracy assessment was conducted based on ground data delineated on 12 U.S. Geological Survey 7-5 min series topographic maps. These results produced overall land-cover classification accuracies of 71.9 and 76.8 per cent for AVHRR and MSS, respectively. While the accuracies for predominant categories were similar for both sensors, land-cover discrimination for less commonly occurring and/or spatially heterogeneous categories was improved with the MSS data set. The AVHRR, however, performed as well as or better than the MSS in classifying large homogeneous areas. The application of AVHRR data with its lower processing cost and more frequent worldwide coverage appears promising for regional land-cover mapping.

Gervin, J. C.↗

Observing Planets and Small Bodies in Sputtered High Energy Atom (SHEA) Fluxes

The evolution of the surfaces of bodies unprotected by either strong magnetic fields or thick atmospheres in the Solar System is caused by various processes, induced by photons, energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, with implications for their evolution. Ion precipitation produces neutral atom releases into the exosphere through ion sputtering, with velocity distribution extending well above the particle escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA The evolution of the surfaces of bodies unprotected by either strong magnetic fields or thick atmospheres in the Solar System is caused by various processes, induced by photons, energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, with implications for their evolution. Ion precipitation produces neutral atom releases into the exosphere through ion sputtering, with velocity distribution extending well above the particle escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA and suggests the need for quantitative results from laboratory simulations and molecular physic modeling in order to understand SHEA data from planetary missions. In the Appendix, referenced computer simulations using existing sputtering data are reviewed.

Milillo, A.↗

SHAD: the Scalable High-performance Algorithms and Data-structures Library

The unprecedented amount of data that needs to be processed in emerging data analytics applications poses novel challenges to industry and academia. Scalability and high performance become more than a desirable feature because, due to the scale and the nature of the problems, they draw the line between what is achievable and what is unfeasible. In this paper, we propose SHAD, the Scalable High-performance Algorithms and Data-structures library. SHAD adopts a modular design that confines low level details and promotes reuse. SHAD's core is built on an Abstract Runtime Interface which enhances portability and identifies the minimal set of features of the underlying system required by the framework. The core library includes common data-structures such as: Array, Vector, Map and Set. These are designed to accommodate significant amount of data which can be accessed in massively parallel environments, and used as building blocks for SHAD extensions, i.e. higher level software libraries. We have validated and evaluated our design with a performance and scalability study of the core components of the library. We have validated the design flexibility by proposing a Graph Library as an example of SHAD extension, which implements two different graph data-structures; we evaluate their performance with a set of graph applications. Experimental results show that the approach is promising in terms of both performance and scalability. On a distributed system with 320 cores, SHAD Arrays are able to sustain a throughput of 65 billion operations per second, while SHAD Maps sustain 1 billion of operations per second. Algorithms implemented using the Graph Library exhibit performance and scalability comparable to a custom solution, but with smaller development effort.

Castellana, Vito G.↗

Simplified ISCCP Cloud Regimes for Evaluating Cloudiness in CMIP5 Models

We take advantage of the ISCCP simulator implemented in many CMIP5 models, in order to introduce a credible method to compare model cloud output with corresponding ISCCP observations based on the cloud regime (CR) concept. A novel simplified CR derivation method is introduced that exploits the co-variations of three variables, cloud optical thickness, cloud top pressure and cloud fraction (, pc, CF). Following evaluation criteria established in a previous paper of ours, we assess model cloud simulation performance based on how well the simplified CRs are simulated in terms of similarity of centroids, global values and map correlations of relative frequency-of-occurrence (RFO), and long term total clouds amounts (as opposed to average instantaneous cloud fractions). Mirroring prior results, modeled clouds tend to be too thick optically and not as extensive as in observations. Averaging across models does not help in this or other evaluation aspects, except for cloud geographical locations. A good model for one criterion is usually found to be also a good performer for other criteria despite plenty of exceptions for the criterion of centroid similarity. On the other hand, CRs that include convectively generated high-altitude clouds are not as well simulated here compared to the previous study, but other regimes containing low clouds with high CF improve. Models that have previously performed well distinguish themselves again here, but improvements from other models are also found. The simplified cloud evaluation method thus proves to be more forgiving while still good enough to reveal model weaknesses.

global climate models↗

A reduced basis approach for PDEs on parametrized geometries based on the shifted boundary finite element method and application to a Stokes flow

Here we propose a model order reduction technique integrating the Shifted Boundary Method (SBM) with a POD-Galerkin strategy. This approach allows to deal with complex parametrized domains in an efficient and straightforward way. The impact of the proposed approach is threefold. First, problems involving parametrizations of complex geometrical shapes and/or large domain deformations can be efficiently solved at full-order by means of the SBM. This unfitted boundary method permits to avoid remeshing and the tedious handling of cut cells by introducing an approximate surrogate boundary. Second, the computational effort is reduced by the development of a Reduced Order Model (ROM) technique based on a POD-Galerkin approach. Third, the SBM provides a smooth mapping from the true to the surrogate domain, and for this reason, the stability and performance of the reduced order basis are enhanced. This feature is the net result of the combination of the proposed ROM approach and the SBM. Similarly, the combination of the SBM with a projection-based ROM gives the great advantage of an easy and fast to implement algorithm considering geometrical parametrization with large deformations. The transformation of each geometry to a reference geometry (morphing) is in fact not required. These combined advantages will allow the solution of PDE problems more efficiently. We illustrate the performance of this approach on a number of two-dimensional Stokes flow problems.

42 ENGINEERING↗