SEARCH · Search NASA

Results for “Conformable gridding”

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

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Operating-Envelopes-Aware Decentralized Welfare Maximization for Energy Communities: Preprint

We propose an operating-envelope-aware, prosumer-centric, and efficient energy community that aggregates individual and shared community distributed energy resources downstream of a regulated distribution system operator's (DSO) net energy metering revenue meter. Due to the elevated risk of grid constraint violations and to ensure safe network operation, the DSO imposes dynamic export and import limits, known as dynamic operating envelopes, on end-users' revenue meters. Given the operating envelopes, the proposed community market mechanism maximizes the community's social welfare in a decentralized fashion while every community member abides by its own operating envelopes. We show that the proposed market mechanism conforms with the cost-causation principle and guarantees community members a surplus level no less than their maximum surplus when they autonomously face the DSO. Lastly, a numerical study is implemented to showcase and compare the community's welfare under the proposed operating-envelopes-aware mechanisms to others, including the welfare of customers under the DSO's regime.

distributed energy resources aggregation↗

NASA NTRS1990

A kinematical/numerical analysis of rotor-stator interaction noise

In this study, the unsteady, thin-layer Navier-Stokes equations are solved using a system of patched grids for a rotor-stator configuration of an axial turbine. The study examines the plurality of spinning modes that are present in such an interaction. The propagation of these modes is analyzed and appropriate grid spacing chosen in the far upstream and downstream regions to attenuate reflections from the computational boundaries. In addition, radiating boundary conditions are implemented based on the farfield acoustical behavior of the flow field. Results in the form of pressure amplitudes and the spectra of turbine tones are presented. Numerical results and experimental data are compared wherever possible. The numerical results are also shown to conform with the predictions of a kinematical analysis of the flowfield.

Rangwalla, Akil A.↗

DOE OSTI2020

Will Solar Panels Help When the Power Goes Out? Planning for PV Resilience

PV systems, like all energy systems, must chart a difficult course in order to improve resilience. They need to simultaneously withstand extreme conditions that cause grid power disruptions, be configured in a way that allows electricity to be delivered to meet the PV user's needs, and overcome disruptions arising from solar resource intermittency. Certain considerations reveal themselves as key to a PV system's resilience benefits. Choice of hardware is critical, as is site selection for the PV array. Voluminous evidence in a hurricane context demonstrates that durable hardware, properly maintained and sited, leads to PV systems surviving extreme weather with their functionality intact. Structurally weak or poorly maintained PV equipment, on the other hand, has demonstrated failures during extreme weather events. PV system stakeholders should ensure that best practices are followed for hardware procurement and siting. It is also important to understand that PV systems are not one-size-fits-all. While incorporating islanding capabilities and energy storage into the system's electrical topology will typically enhance resilience, the user's needs will determine the optimally resilient conformation of any given system. These factors - appropriate hardware, siting, maintenance, and system design - are the primary drivers of PV resilience. PV systems possess high potential for resilient energy generation due to their distributed, abundant, and free fuel supply, but the technology itself does not inherently lead to resilient outcomes. Realizing PV's potential for resilience requires consideration of these determinative factors.

14 SOLAR ENERGY↗

DOE OSTI2022

Variable resolution Poisson-disk sampling for meshing discrete fracture networks

Here, we present the near-Maximal Algorithm for Poisson-disk Sampling (nMAPS) to generate point distributions for variable resolution Delaunay triangular and tetrahedral meshes in two and three-dimensions, respectively. nMAPS consists of two principal stages. In the first stage, an initial point distribution is produced using a cell-based rejection algorithm. In the second stage, holes in the sample are detected using an efficient background grid and filled in to obtain a near-maximal covering. Extensive testing shows that nMAPS generates a variable resolution mesh in linear run time with the number of accepted points. We demonstrate nMAPS capabilities by meshing three-dimensional discrete fracture networks (DFN) and the surrounding volume. The discretized boundaries of the fractures, which are represented as planar polygons, are used as the seed of 2D-nMAPS to produce a conforming Delaunay triangulation. The combined mesh of the DFN is used as the seed for 3D-nMAPS, which produces conforming Delaunay tetrahedra surrounding the network. Under a set of conditions that naturally arise in maximal Poisson-disk samples and are satisfied by nMAPS, the two-dimensional Delaunay triangulations are guaranteed to only have well-behaved triangular faces. While nMAPS does not provide triangulation quality bounds in more than two dimensions, we found that low-quality tetrahedra in 3D are infrequent, can be readily detected and removed, and a high-quality balanced mesh is produced.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Lifting MGARD: Construction of (pre)wavelets on the interval using polynomial predictors of arbitrary order

MGARD (MultiGrid Adaptive Reduction of Data) is an algorithm for compressing and refactoring scientific data, based on the theory of multigrid methods. The core algorithm is built around stable multilevel decompositions of conforming piecewise linear $C^0$ finite element spaces, enabling accurate error control in various norms and derived quantities of interest. In this work, we extend this construction to arbitrary order Lagrange finite elements $\mathbb{Q}_p$, $p \geq 0$, and propose a reformulation of the algorithm as a lifting scheme with polynomial predictors of arbitrary order. Additionally, a new formulation using a compactly supported wavelet basis is discussed, and an explicit construction of the proposed wavelet transform for uniform dyadic grids is described.

Reshniak, Viktor [Oak Ridge National Laboratory (O↗

NASA NTRS1998

Simulation of Conformal Spiral Slot Antennas on Composite Platforms

During the course of the grant, we wrote and distributed about 12 reports and an equal number of journal papers supported fully or in part by this grant. The list of reports (title & abstract) and papers are given in Appendices A and B. This grant has indeed been instrumental in developing a robust hybrid finite element method for the analysis of complex broadband antennas on doubly curved platforms. Previous to the grant, our capability was limited to simple printed patch antennas on mostly planar platforms. More specifically: (1) mixed element formulations were developed and new edge-based prisms were introduced; (2) these elements were important in permitting flexibility in geometry gridding for most antennas of interest; (3) new perfectly matched absorbers were introduced for mesh truncations associated with highly curved surfaces; (4) fast integral algorithms were introduced for boundary integral truncations reducing CPU time from O(N-2) down to O(N-1.5) or less; (5) frequency extrapolation schemes were developed for efficient broadband performance evaluations. This activity has been successfully continued by NASA researchers; (6) computer codes were developed and extensively tested for several broadband configurations. These include FEMA-CYL, FEMA-PRISM and FEMA-TETRA written by L. Kempel, T. Ozdemir and J. Gong, respectively; (7) a new infinite balun feed was designed nearly constant impedance over the 800-3000 MHz operational band; (8) a complete slot spiral antenna was developed, fabricated and tested at NASA Langley. This new design is a culmination of the projects goals and integrates the computational and experimental efforts. this antenna design resulted in a U.S. patent and was revised three times to achieve the desired bandwidth and gain requirements from 800-3000 MHz.

Volakis, J. L.↗

DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal photon sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal electron sources

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1995

Algorithms for high aspect ratio oriented triangulations

Grid generation plays an integral part in the solution of computational fluid dynamics problems for aerodynamics applications. A major difficulty with standard structured grid generation, which produces quadrilateral (or hexahedral) elements with implicit connectivity, has been the requirement for a great deal of human intervention in developing grids around complex configurations. This has led to investigations into unstructured grids with explicit connectivities, which are primarily composed of triangular (or tetrahedral) elements, although other subdivisions of convex cells may be used. The existence of large gradients in the solution of aerodynamic problems may be exploited to reduce the computational effort by using high aspect ratio elements in high gradient regions. However, the heuristic approaches currently in use do not adequately address this need for high aspect ratio unstructured grids. High aspect ratio triangulations very often produce the large angles that are to be avoided. Point generation techniques based on contour or front generation are judged to be the most promising in terms of being able to handle complicated multiple body objects, with this technique lending itself well to adaptivity. The eventual goal encompasses several phases: first, a partitioning phase, in which the Voronoi diagram of a set of points and line segments (the input set) will be generated to partition the input domain; second, a contour generation phase in which body-conforming contours are used to subdivide the partition further as well as introduce the foundation for aspect ratio control, and; third, a Steiner triangulation phase in which points are added to the partition to enable triangulation while controlling angle bounds and aspect ratio. This provides a combination of the advancing front/contour techniques and refinement. By using a front, aspect ratio can be better controlled. By using refinement, bounds on angles can be maintained, while attempting to minimize the number of Steiner points.

Posenau, Mary-Anne K.↗

NASA NTRS1984

Tabulations of ambient ozone data obtained by GASP (Global Air Sampling Program) airliners, March 1975 to July 1979

Tabulations are given of GASP ambient ozone mean, standard deviation, median, 84th percentile, and 98th percentile values, by month, flight level, and geographical region. These data are tabulated to conform to the temporal and spatial resolution required by FAA Advisory Circular 120-38 (monthly by 2000 ft in altitude by 5 deg in latitude) for climatological data used to show compliance with cabin ozone regulations. In addition seasonal x 10 deg latitude tabulations are included which are directly comparable to and supersede the interim GASP ambient ozone tabulations given in appendix B of FAA-EE-80-43 (NASA TM-81528). Selected probability variations are highlighted to illustrate the spatial and temporal variability of ambient ozone and to compare results from the coarse and fine grid analyses.

Jasperson, W. H.↗

NASA NTRS1990

Obtaining an equivalent beam

In modeling a complex structure the researcher was faced with a component that would have logical appeal if it were modeled as a beam. The structure was a mast of a robot controlled gantry crane. The structure up to this point already had a large number of degrees of freedom, so the idea of conserving grid points by modeling the mast as a beam was attractive. The researcher decided to make a separate problem of of the mast and model it in three dimensions with plates, then extract the equivalent beam properties by setting up the loading to simulate beam-like deformation and constraints. The results could then be used to represent the mast as a beam in the full model. A comparison was made of properties derived from models of different constraints versus manual calculations. The researcher shows that the three-dimensional model is ineffective in trying to conform to the requirements of an equivalent beam representation. If a full 3-D plate model were used in the complete representation of the crane structure, good results would be obtained. Since the attempt is to economize on the size of the model, a better way to achieve the same results is to use substructuring and condense the mast to equivalent end boundary and intermediate mass points.

Butler, Thomas G.↗

DOE OSTI2020

The arbitrary‐order virtual element method for linear elastodynamics models: convergence, stability and dispersion‐dissipation analysis

Abstract We design the conforming virtual element method for the numerical approximation of the two‐dimensional elastodynamics problem. We prove stability and convergence of the semidiscrete approximation and derive optimal error estimates under h ‐ and p ‐refinement in both the energy and the L 2 norms. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including nonconvex cells up to order four in the h ‐refinement setting. Exponential convergence is also experimentally observed under p ‐refinement. Finally, we present a dispersion‐dissipation analysis for both the semidiscrete and fully discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion‐dissipation properties.

Antonietti, Paola F.↗

DOE OSTI2021

A 24-Channel Digitizer With a JESD204B-Compliant Serial Interface for High-Speed Detectors

In this report a 24-channel application-specific integrated circuit (ASIC) for the readout of high-speed CMOS active pixel sensors for charged particle detection is presented. The chip comprises 24 preamplifiers, 24 distinct 12-bit, 25 MSPS Pipelined analog-to-digital converters (ADCs) with self-calibration, an internal phase-locked loop (PLL), and a 3 Gbps serial interface that conforms to the JESD204B standard. To simplify interfacing with a variety of sensors, the ASIC also includes an automatic offset calibration loop. The high level of integration of the ASIC reduces overall system cost and area, and exploiting the signal characteristics of the image sensor allows the ADC to be optimized for reduced power dissipation. The use of an integrated serializer and an industry standard protocol simplifies integration of the ASIC into a complete camera system. The ASIC, called the High-Speed Image Preprocessor Targeted for Electron Readout, or HIPSTER, with a die area of 64.26 mm 2 , is packaged in a 480-ball grid array (BGA) and is fabricated in 180-nm CMOS technology. HIPSTER achieves typical differential nonlinearity (DNL) < 0.55 LSB, input-referred thermal noise of $114.5 μV-rms, and a bit error rate (BER) of better than 10-14. The power dissipation is 98 mW/channel.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Antipodal self-duality of square fishnet graphs

In strongly deformed planar 𝒩 = 4 super-Yang-Mills theory, or fishnet theory, a point-split single-trace correlation function of four dimension-𝑚 scalar operators is given by a single Feynman integral, which involves integrating over locations of a 𝑚 × 𝑚 grid of points. We show that for any integer 𝑚 this square fishnet graph is invariant under the combined action of a kinematic map and the antipode map of the Hopf algebra on multiple polylogarithms; i.e. it possesses an antipodal self-duality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Multiple Functional Bonds Integrated Interphases for Long Cycle Sodium-Ion Batteries

Sodium-ion batteries (SIBs) have garnered significant interest as one of the most promising energy suppliers for power grid energy storage. However, the poor electrode/electrolyte interfacial stability leads to continual electrolyte decomposition and transition metal dissolution, resulting in rapid performance degradation of SIBs. In this work, we propose a strategy integrating multiple functional bonds to regulate electrode/electrolyte interphase by triple-coupling of succinonitrile (SN), sodium hexafluorophosphate (NaPF 6 ) and fluorinated ethylene carbonate (FEC). Theoretical calculation and experiment results show that the solvation structure of Na + and ClO 4 – is effectively reconfigured by the solvated FEC, SN and PF 6 – in PC-based carbonate electrolyte. The newly developed electrolyte demonstrates increased Na + -FEC coordination, weakened interaction of Na + -PC and participation of SN and PF 6 – anions in solvation, resulting in the formation of a conformal interfacial layer comprising of sodium oxynitrides (NaN x O y ), sodium fluoride (NaF) and phosphorus oxide compounds (NaP x O y ). Consequently, a 3 Ah pouch full cell of hard carbon//NaNi 1/3 Fe 1/3 Mn 1/3 O 2 exhibits an excellent capacity retention of 90.4 % after 1000 cycles. Detailed postmortem analysis of interface chemistry is further illustrated by multiple characterization methods. Finally, this study provides a new avenue for developing electrolyte formulations with multiple functional bonds integrated interphases to significantly improve the long-term cycling stability of SIBs.

25 ENERGY STORAGE↗

DOE OSTI2020

The arbitrary-order virtual element method for linear elastodynamics models. Convergence, stability and dispersion-dissipation analysis.

We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under $\textit{h}$-refinement in both the energy and the $L^2$ norms, and optimal error estimates under $\textit{p}$-refinement in the energy norm. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally observed under p-refinement. Finally, we present a dispersion-dissipation analysis for both the semi-discrete and fully-discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion-dissipation properties.

97 MATHEMATICS AND COMPUTING↗