Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multilevel models”

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 199 records · Page 11

What Multilevel Parallel Programs do when you are not Watching: A Performance Analysis Case Study Comparing MPI/OpenMP, MLP, and Nested OpenMP

With the current trend in parallel computer architectures towards clusters of shared memory symmetric multi-processors, parallel programming techniques have evolved that support parallelism beyond a single level. When comparing the performance of applications based on different programming paradigms, it is important to differentiate between the influence of the programming model itself and other factors, such as implementation specific behavior of the operating system (OS) or architectural issues. Rewriting-a large scientific application in order to employ a new programming paradigms is usually a time consuming and error prone task. Before embarking on such an endeavor it is important to determine that there is really a gain that would not be possible with the current implementation. A detailed performance analysis is crucial to clarify these issues. The multilevel programming paradigms considered in this study are hybrid MPI/OpenMP, MLP, and nested OpenMP. The hybrid MPI/OpenMP approach is based on using MPI [7] for the coarse grained parallelization and OpenMP [9] for fine grained loop level parallelism. The MPI programming paradigm assumes a private address space for each process. Data is transferred by explicitly exchanging messages via calls to the MPI library. This model was originally designed for distributed memory architectures but is also suitable for shared memory systems. The second paradigm under consideration is MLP which was developed by Taft. The approach is similar to MPi/OpenMP, using a mix of coarse grain process level parallelization and loop level OpenMP parallelization. As it is the case with MPI, a private address space is assumed for each process. The MLP approach was developed for ccNUMA architectures and explicitly takes advantage of the availability of shared memory. A shared memory arena which is accessible by all processes is required. Communication is done by reading from and writing to the shared memory.

Jost, Gabriele↗

MFNets: Multifidelity data-driven networks for Bayesian learning and prediction

This paper presents a multifidelity uncertainty quantification framework called MFNets. We seek to address three existing challenges that arise when experimental and simulation data from different sources are used to enhance statistical estimation and prediction with quantified uncertainty. Specifically, we demonstrate that MFNets can (1) fuse heterogeneous data sources arising from simulations with different parameterizations, e.g simulation models with different uncertain parameters or data sets collected under different environmental conditions; (2) encode known relationships among data sources to reduce data requirements; and (3) improve the robustness of existing multi-fidelity approaches to corrupted data. MFNets construct a network of latent variables (LVs) to facilitate the fusion of data from an ensemble of sources of varying credibility and cost. These LVs are posited as explanatory variables that provide the source of correlation in the observed data. Furthermore, MFNets provide a way to encode prior physical knowledge to enable efficient estimation of statistics and/or construction of surrogates via conditional independence relations on the LVs. We highlight the utility of our framework with a number of theoretical results which assess the quality of the posterior mean as a frequentist estimator and compare it to standard sampling approaches that use single fidelity, multilevel, and control variate Monte Carlo estimators. We also use the proposed framework to derive the Monte Carlo-based control variate estimator entirely from the use of Bayes rule and linear-Gaussian models -- to our knowledge the first such derivation. Finally, we demonstrate the ability to work with different uncertain parameters across different models.

97 MATHEMATICS AND COMPUTING↗

High-Efficiency Flexible Multilevel Photon Sieves by Single-Step Laser-Based Fabrication and Optical Analysis

Over the past several decades, the need for high-resolution, high-efficiency, lightweight, high contrast focusing optics has continued to increase due to their applications in fields such as astronomy, spectroscopy, free-space optical communications, defense, and remote sensing. In recent years, photon sieve planar diffractive optics have been developed on flexible, lightweight polyimide substrates. However, transmission efficiencies have continuously been very low (~1-4%), thus impeding the widespread use of photon sieves in practical applications. Here, we present a flexible, lightweight, N = 4 level phase photon sieve with over 25% transmission efficiency; nearly triple that of any other photon sieve reported thus far. In this study, the photon sieve was fabricated via a novel pulsed laser ablation method using an ultraviolet source, and the total time to fabricate a ~3 sq cm sample was tens of seconds. Theoretical analysis of the photon sieve was carried out via the FDTD method, and was in very good agreement with experimental results. We also calculated via FDTD modeling the behavior of higher step (N = 8, 16, 32) photon sieves for further enhanced efficiencies, and show a fundamental limit on photon sieve efficiency of 70% in the limit of increasing step number. This multilevel photon sieve represents a new step in high-resolution diffractive optics, moving towards efficiencies suitable for widespread applications.

Julian, Matthew N.↗

Finite-element analysis of nonlinear conduction problems subject to moving fields

Through the use of a space-time warp, specialized moving finite elements are developed that can be employed to generate a nonlinear heat conduction model for situations involving traveling boundary and heat generation fields superposed on an initial state. To facilitate the solution of the resulting nonlinear finite-element formulation, a multilevel heuristic iterative solution strategy is developed. In order to demonstrate the versatility and accuracy of the moving elements and their associated nonlinear solution strategy, the results of several numerical experiments are presented.

Padovan, J.↗

Interrelated structures of the transport shock and collisional relaxation layer in a multitemperature, multilevel ionized gas

The gas dynamic structures of the transport shock and the downstream collisional relaxation layer are evaluated for partially ionized monatomic gases. Elastic and inelastic collisional nonequilibrium effects are taken into consideration. Three electronic levels are accounted for in the microscopic model of the atom. Nonequilibrium processes with respect to population of levels and species plus temperature are considered. By using an asymptotic technique the shock morphology is found on a continuum flow basis. The asymptotic procedure gives two distinct layers in which the nonequilibrium effects to be considered are different. A transport shock appears as the inner solution to an outer collisional relaxation layer in which the gas reaches local equilibrium. A family of numerical examples is displayed for different flow regimes. Argon and helium models are used in these examples.

Vinolo, A. R.↗

Indefinite and bidirectional near-infrared nanocrystal photoswitching

Materials whose luminescence can be switched by optical stimulation drive technologies ranging from superresolution imaging, nanophotonics, and optical data storage, to targeted pharmacology, optogenetics, and chemical reactivity. Furthermore, these photoswitchable probes, including organic fluorophores and proteins, can be prone to photodegradation and often operate in the ultraviolet or visible spectral regions. Colloidal inorganic nanoparticles can offer improved stability, but the ability to switch emission bidirectionally, particularly with near-infrared (NIR) light, has not, to our knowledge, been reported in such systems. Here, we present two-way, NIR photoswitching of avalanching nanoparticles (ANPs), showing full optical control of upconverted emission using phototriggers in the NIR-I and NIR-II spectral regions useful for subsurface imaging. Employing single-step photodarkening and photobrightening, we demonstrate indefinite photoswitching of individual nanoparticles (more than 1,000 cycles over 7 h) in ambient or aqueous conditions without measurable photodegradation. Critical steps of the photoswitching mechanism are elucidated by modelling and by measuring the photon avalanche properties of single ANPs in both bright and dark states. Unlimited, reversible photoswitching of ANPs enables indefinitely rewritable two-dimensional and three-dimensional multilevel optical patterning of ANPs, as well as optical nanoscopy with sub-Å localization superresolution that allows us to distinguish individual ANPs within tightly packed clusters.

36 MATERIALS SCIENCE↗

Parallelized modelling and solution scheme for hierarchically scaled simulations

This two-part paper presents the results of a benchmarked analytical-numerical investigation into the operational characteristics of a unified parallel processing strategy for implicit fluid mechanics formulations. This hierarchical poly tree (HPT) strategy is based on multilevel substructural decomposition. The Tree morphology is chosen to minimize memory, communications and computational effort. The methodology is general enough to apply to existing finite difference (FD), finite element (FEM), finite volume (FV) or spectral element (SE) based computer programs without an extensive rewrite of code. In addition to finding large reductions in memory, communications, and computational effort associated with a parallel computing environment, substantial reductions are generated in the sequential mode of application. Such improvements grow with increasing problem size. Along with a theoretical development of general 2-D and 3-D HPT, several techniques for expanding the problem size that the current generation of computers are capable of solving, are presented and discussed. Among these techniques are several interpolative reduction methods. It was found that by combining several of these techniques that a relatively small interpolative reduction resulted in substantial performance gains. Several other unique features/benefits are discussed in this paper. Along with Part 1's theoretical development, Part 2 presents a numerical approach to the HPT along with four prototype CFD applications. These demonstrate the potential of the HPT strategy.

Padovan, Joe↗

Operation and Control of a Back to Back Modular Multilevel Converter System for Grid Forming Application with Advanced Grid Support Functionalities

The main focus of this paper is to investigate the possibility of using modular multilevel converter based back-to-back system for grid forming applications. A decentralized hierarchical control architecture with a modified local controller based on nonlinear techniques for grid forming application have been investigated. Grid functionalities based on IEEE 1547-2018 for grid forming converters have been utilized. The local controllers ensure distortion free balanced sinusoidal output voltage on the grid forming side with unbalanced or nonlinear loading condition. The grid following side ensures unity power factor currents under unbalanced grid voltage condition. The dc bus voltage control ensures oscillation free voltage under balanced conditions and control of the average voltage under unbalanced condition. Efficacy of the overall system is verified by modeling the system in MATLAB/Simulink and PLECS domain and the most important case studies are presented.

advanced grid support functionalities↗

Multigrid-Reduction-in-Time for the Rotating Shallow Water Equations

We consider multilevel time-parallel methods for the numerical solution of the rotating shallow water equations. In particular, the multigrid-reduction-in-time (MGRIT) algorithm is used for the parallel time integration. An asymptotic model is used at the coarse levels while the full model is employed at the finer levels. The asymptotic model is well-suited for highly oscillatory partial differential equations like the rotating shallow water equations because it can accurately and stably take the required large time-steps on coarse levels. Our work exploits the flexibility of the MGRIT algorithm in terms of the number of levels and relaxation schemes to show some computational benefits, especially with respect to FCF-relaxation and data reuse.

97 MATHEMATICS AND COMPUTING↗

Parallelized solvers for heat conduction formulations

Based on multilevel partitioning, this paper develops a structural parallelizable solution methodology that enables a significant reduction in computational effort and memory requirements for very large scale linear and nonlinear steady and transient thermal (heat conduction) models. Due to the generality of the formulation of the scheme, both finite element and finite difference simulations can be treated. Diverse model topologies can thus be handled, including both simply and multiply connected (branched/perforated) geometries. To verify the methodology, analytical and numerical benchmark trends are verified in both sequential and parallel computer environments.

Padovan, Joe↗

Scalable multilevel Monte Carlo methods exploiting parallel redistribution on coarse levels

Here, we study an element agglomeration coarsening strategy that requires data redistribution at coarse levels when the number of coarse elements becomes smaller than the number of MPI processes used on the finest level. The overall procedure generates coarse elements (general unstructured unions of fine grid elements) within the framework of element-based algebraic multigrid methods (or AMGe) studied previously. The AMGe-generated coarse spaces have the ability to exhibit approximation properties of the same order as the fine-level spaces since by construction they contain the piecewise polynomials of the same order as on the fine level. These approximation properties are key for the successful use of AMGe in multilevel solvers for nonlinear partial differential equations as well as for multilevel Monte Carlo (MLMC) simulations. The ability to coarsen without being constrained by the number of MPI processes, as described in the present paper, allows to improve the scalability of these solvers as well as the overall MLMC method. The paper illustrates this latter fact with detailed scalability study of MLMC simulations applied to model Darcy equations with a stochastic log-normal permeability field.

AMGe↗

Analyzing Tropical Waves Using the Parallel Ensemble Empirical Model Decomposition Method: Preliminary Results from Hurricane Sandy

In this study, we discuss the performance of the parallel ensemble empirical mode decomposition (EMD) in the analysis of tropical waves that are associated with tropical cyclone (TC) formation. To efficiently analyze high-resolution, global, multiple-dimensional data sets, we first implement multilevel parallelism into the ensemble EMD (EEMD) and obtain a parallel speedup of 720 using 200 eight-core processors. We then apply the parallel EEMD (PEEMD) to extract the intrinsic mode functions (IMFs) from preselected data sets that represent (1) idealized tropical waves and (2) large-scale environmental flows associated with Hurricane Sandy (2012). Results indicate that the PEEMD is efficient and effective in revealing the major wave characteristics of the data, such as wavelengths and periods, by sifting out the dominant (wave) components. This approach has a potential for hurricane climate study by examining the statistical relationship between tropical waves and TC formation.

PEEMD↗

Large space structures controls research and development at Marshall Space Flight Center: Status and future plans

Work performed in Large Space Structures Controls research and development program at Marshall Space Flight Center is described. Studies to develop a multilevel control approach which supports a modular or building block approach to the buildup of space platforms are discussed. A concept has been developed and tested in three-axis computer simulation utilizing a five-body model of a basic space platform module. Analytical efforts have continued to focus on extension of the basic theory and subsequent application. Consideration is also given to specifications to evaluate several algorithms for controlling the shape of Large Space Structures.

Buchanan, H. J.↗

Estimating posterior quantity of interest expectations in a multilevel scalable framework

Scalable approaches for uncertainty quantification are necessary for characterizing prediction confidence in large-scale subsurface flow simulations with uncertain permeability. To this end we explore a multilevel Monte Carlo approach for estimating posterior moments of a particular quantity of interest, where we employ an element-agglomerated algebraic multigrid (AMG) technique to generate the hierarchy of coarse spaces with guaranteed approximation properties for both the generation of spatially correlated random fields and the forward simulation of Darcy's law to model subsurface flow. In both these components (sampling and forward solves), we exploit solvers that rely on state-of-the-art scalable AMG. To illustrate the applicability of this approach, numerical tests are performed on two 3D examples-a unit cube and an egg-shaped domain with an irregular boundary-where the scalability of each simulation as well as the scalability of the overall algorithm are demonstrated.

97 MATHEMATICS AND COMPUTING↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct Power Control of Back to Back Modular Multilevel Converter with Advanced Grid Support Functions for Grid Forming Application

Possibility of using back-to-back modular multilevel converter system with a direct power control architecture has been investigated in this paper. The rectifier side is controlled to be in grid following mode whereas the inverter side is controlled in grid forming mode. Both the sides’ local controllers use Lyapunov energy function based architecture to accomplish the objective of active and reactive power control as well as maintaining a fixed user defined voltage over the dc bus on the rectifier side. Therefore, the proposed direct power control architecture has embedded dc bus control architecture to maintain fixed user defined value of the capacitors. Similar architecture is implemented on the inverter side to accomplish either voltage control during grid forming or power control during grid following. Advanced grid support functionalities based on IEEE-1547-2018 is implemented and utilized to generate the reference values for either the rectifier or the inverter sides. The overall system is modeled based on MATLAB/Simulink and PLECS domain and various important case studies to verify the efficacy of the overall system has been presented.

direct power control (DPC)↗

Bosonic Pair Production and Squeezing for Optical Phase Measurements in Long-Lived Dipoles Coupled to a Cavity

Here, we propose to simulate bosonic pair creation using large arrays of long-lived dipoles with multilevel internal structure coupled to an undriven optical cavity. Entanglement between the atoms, generated by the exchange of virtual photons through a common cavity mode, grows exponentially fast and is described by two-mode squeezing of effective bosonic quadratures. The mapping between an effective bosonic model and the natural spin description of the dipoles allows us to realize the analog of optical homodyne measurements via straightforward global rotations and population measurements of the electronic states, and we propose to exploit this for quantum-enhanced sensing of an optical phase (common and differential between two ensembles). We discuss a specific implementation based on Sr atoms and show that our sensing protocol is robust to sources of decoherence intrinsic to cavity platforms. Our proposal can open unique opportunities for next-generation optical atomic clocks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗