Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh data”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗
NASA NTRS2007

Simulating Space Capsule Water Landing with Explicit Finite Element Method

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Wang, John T.↗
NASA NTRS2011

NASA Soil Moisture Active Passive (SMAP) Mission Formulation

The Soil Moisture Active Passive (SMAP) Mission is one of the first Earth observation satellites being formulated by NASA in response to the 2007 National Research Council s Earth Science Decadal Survey [1]. SMAP s measurement objectives are high-resolution global measurements of near-surface soil moisture and its freeze-thaw state. These measurements would allow significantly improved estimates of water, energy and carbon transfers between the land and atmosphere. The soil moisture control of these fluxes is a key factor in the performance of atmospheric models used for weather forecasts and climate projections. Soil moisture measurements are also of great importance in assessing flooding and monitoring drought. Knowledge gained from SMAP s planned observations can help mitigate these natural hazards, resulting in potentially great economic and societal benefits. SMAP measurements would also yield high resolution spatial and temporal mapping of the frozen or thawed condition of the surface soil and vegetation. Observations of soil moisture and freeze/thaw timing over the boreal latitudes will contribute to reducing a major uncertainty in quantifying the global carbon balance and help resolve an apparent missing carbon sink over land. The SMAP mission would utilize an L-band radar and radiometer sharing a rotating 6-meter mesh reflector antenna (see Figure 1) [2]. The radar and radiometer instruments would be carried onboard a 3-axis stabilized spacecraft in a 680 km polar orbit with an 8-day repeating ground track. The instruments are planned to provide high-resolution and high-accuracy global maps of soil moisture at 10 km resolution and freeze/thaw at 3 km resolution, every two to three days (see Table 1 for a list of science data products). The mission is adopting a number of approaches to identify and mitigate potential terrestrial radio frequency interference (RFI). These approaches are being incorporated into the radiometer and radar flight hardware and ground processing designs.

Entekhabi, Dara↗
DOE OSTI2021

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4a Raw Data

Calculated formation energy of (0001 ̅) c-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (0001 ̅) c-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗
DOE OSTI2021

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4b Raw Data

Calculated formation energy of (1010 ̅) m-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (1010 ̅) m-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗
NASA NTRS2018

CFD Predictions of N+3 Cycle Emissions for a Three-Cup Gas-Turbine Combustor

The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a unique pre-filming type LDI-3 fuel injector for a three-cup, nineteen-element flametube configuration. All computations were performed with a consistent approach of mesh-generation, spray modeling, reduced finite-rate kinetics and turbulence-chemistry interaction, as developed for CFD analysis of single-element and multi-element LDI-3 designs with OpenNCC. Emissions and flowfield characteristics were predicted for a generic NASA N+3 engine cycle, with particular focus on the 7% and 30% ICAO power operating conditions. For both the conditions studied, the CFD analysis provided very good predictions for EINOx when compared with experimental data measured at NASA Glenn Research Center.

injectors↗
DOE OSTI2017

TORCHE - Toolbox for Reactor Cross-Flow Heat Exchangers

As interest within the scientific community for the development of small modular reactors (SMRs) grows, utilizing compact and cross-flow heat exchangers present many design advantages. Cross-flow configurations of tube bundles allow water, liquid metal, or a molten salt to be contained within the tube side or shell side in an efficient spatial arrangement. The goal of this project is to provide the user with an easy-to-use set of tools for cross-flow heat exchanger design. The tools are intended to require little to no geometric set-up like meshing in a CFD code in order to enable rapid preliminary design and optimization. With this set of simple Python scripts, the user can determine the necessary quantities and correctional factors in order to assess the heat exchanger characteristics. The two characteristics used in designing novel and efficient heat exchangers incorporated into this work currently are pressure drop and heat transfer coefficients. In support of these efforts, the literature has been scoured for sources of experimental data and proposed methods for calculating pressure drop and heat transfer in cross-flow tube banks. Several empirical correlations and methodologies for calculating pressure drop and heat transfer incorporated into the Toolbox for Reactor Cross-Flow Heat Exchangers (TORCHE). TORCHE is an expansion on the cross-flow pressure drop scripts reported by Haefner et al.. Presented in this summary is a synopsis of the empirical correlations and methodologies of pressure drops and heat transfer coefficient. A validation effort of these models is also presented to provide confidence in the implementation. For this validation of TORCHE, experimental data has been gathered and extracted from graphs in the literature. Error bars were not reported in the experimental data used in this summary and therefore is not included in the analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Multidimensional sum-up rounding for integer programming in optimal experimental design

Here we present a numerical method for approximating the solution of convex integer programs stemming from optimal experimental design. The statistical setup consists of a Bayesian framework for linear inverse problems for which the direct relationship is described by a discretized integral equation. Specifically, we aim to find the optimal sensor placement from a set of candidate locations where data are collected with measurement error. The convex objective function is a measure of the uncertainty, described here by the trace or log-determinant of the posterior covariance matrix, for the discretized linear inverse problem solution. The resulting convex integer program is relaxed, producing a lower bound. An upper bound is obtained by extending the sum-up rounding approach to multiple dimensions. For this extension, we analyze its accuracy as a function of the discretization mesh size for a rectangular domain. We show asymptotic optimality of the integer solution defining the upper bound for different experimental design criteria (A- and D-optimal), by proving the convergence to zero of the gap between the upper and lower bounds as the mesh size goes to zero. The technique is illustrated on a two-dimensional gravity surveying problem for both A-optimal and D-optimal sensor placement where our designs yield better results compared with a thresholding rounding approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

One-shot learning for solution operators of partial differential equations

Learning and solving governing equations of a physical system, represented by partial differential equations (PDEs), from data is a central challenge in many areas of science and engineering. Traditional numerical methods can be computationally expensive for complex systems and require complete governing equations. Existing data-driven machine learning methods require large datasets to learn a surrogate solution operator, which could be impractical. Here, we propose a solution operator learning method that requires only one PDE solution, i.e., one-shot learning, along with suitable initial and boundary conditions. Leveraging the locality of derivatives, we define a local solution operator in small local domains, train it using a neural network, and use it to predict solutions of new input functions via mesh-based fixed-point iteration or meshfree neural-network based approaches. We test our method on various PDEs, complex geometries, and a practical spatial infection spread application, demonstrating its effectiveness and generalization capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Signatures of primordial energy injection from axion strings

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultralight axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a subdominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants f a ≲ 10 12 GeV , up to model dependence on the ultraviolet completion, for axion masses m a ≲ 10 − 29 eV ; future CMB surveys could find striking evidence of axion strings with lower decay constants. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Communication-Avoiding 3D LU Factorization Algorithm for Sparse Matrices

We propose a new algorithm to improve the strong scalability of right-looking sparse LU factorization on distributed memory systems. Our 3D sparse LU algorithm uses a three-dimensional MPI process grid, aggressively exploits elimination tree parallelism and trades off increased memory for reduced per-process communication. We also analyze the asymptotic improvements for planar graphs (e.g., from 2D grid or mesh domains) and certain non-planar graphs (specifically for 3D grids and meshes). For planar graphs with n vertices, our algorithm reduces communication volume asymptotically in n by a factor of O{log n} and latency by a factor of O{log n}. For non-planar cases, our algorithm can reduce the per-process communication volume by 3× and latency by O{n^1/3} times. In all cases, the memory needed to achieve these gains is a constant factor. We implemented our algorithm by extending the 2D data structure used in superLU. Our new 3D code achieves speedups up to 27× for planar graphs and up to 3.3× for non-planar graphs over the baseline 2D superLU when run on 24,000 cores of a Cray XC30.

Sao, Piyush↗
DOE OSTI2020

Parallelization and Performance Portability in Hydrodynamics Codes

With the eve of Exascale computing, performance and portability are at the forefront of all scientific codes. Adding more cores and more energy to a system is no longer a sustainable way to achieve performance, and extra effort must now be made to improve performance in all areas of code and code development. Using hydrodynamic codes as a basis, this work explores numerous techniques to achieve performance in different ways. Adaptive mesh refinement (AMR) is a necessary technique to improve memory optimization in mesh-based simulations. However it is invasive and conventionally difficult to integrate into existing applications, so we present a new branch of AMR to create a smooth transition to these optimizations, which not only improves performance, but also greatly reduces developer effort. We introduce the concept of this improvement as Phantom-Cell AMR, and assess theoretically the improvements, as well as present an application of its use. Other work included involves and investigation into an efficient data structure that ensures optimal memory layout for cache performance, with a target of making codes performant and portable across all architectures. All of the work targets both performance and portability, not just on CPU hardware, but specifically across GPU architectures. Parallel performance is key to all of the methods presented, but the research makes a great effort to improve the portability of all applications to prepare for current high performance computing systems and those on the horizon.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Development of a Dynamically Configurable, Object-Oriented Framework for Distributed, Multi-modal Computational Aerospace Systems Simulation: Second Year Progress Report

Mesh generation has long been recognized as a bottleneck in the CFD process. While much research on automating the volume mesh generation process have been relatively successful,these methods rely on appropriate initial surface triangulation to work properly. Surface discretization has been one of the least automated steps in computational simulation due to its dependence on implicitly defined CAD surfaces and curves. Differences in CAD peometry engines manifest themselves in discrepancies in their interpretation of the same entities. This lack of "good" geometry causes significant problems for mesh generators, requiring users to "repair" the CAD geometry before mesh generation. The problem is exacerbated when CAD geometry is translated to other forms (e.g., IGES )which do not include important topological and construction information in addition to entity geometry. One technique to avoid these problems is to access the CAD geometry directly from the mesh generating software, rather than through files. By accessing the geometry model (not a discretized version) in its native environment, t h s a proach avoids translation to a format which can deplete the model of topological information. Our approach to enable models developed in the Denali software environment to directly access CAD geometry and functions is through an Application Programming Interface (API) known as CAPRI. CAPRI provides a layer of indirection through which CAD-specific data may be accessed by an application program using CAD-system neutral C and FORTRAN language function calls. CAPRI supports a general set of CAD operations such as truth testing, geometry construction and entity queries.

Afjeh, Abdollah A.↗
NASA NTRS1999

An Elasticity-Based Mesh Scheme Applied to the Computation of Unsteady Three-Dimensional Spoiler and Aeroelastic Problems

This paper presents a modification of the spring analogy scheme which uses axial linear spring stiffness with selective spring stiffening/relaxation. An alternate approach to solving the geometric conservation law is taken which eliminates the need for storage of metric Jacobians at previous time steps. Efficiency and verification are illustrated with several unsteady 2-D airfoil Euler computations. The method is next applied to the computation of the turbulent flow about a 2-D airfoil and wing with two and three- dimensional moving spoiler surfaces, and the results compared with Benchmark Active Controls Technology (BACT) experimental data. The aeroelastic response at low dynamic pressure of an airfoil to a single large scale oscillation of a spoiler surface is computed. This study confirms that it is possible to achieve accurate solutions with a very large time step for aeroelastic problems using the fluid solver and aeroelastic integrator as discussed in this paper.

Bartels, Robert E.↗
DOE OSTI2020

Scalable line and plane relaxation in a parallel structured multigrid solver

The efficient solution of sparse, linear systems that arise through the discretization of partial differential equations remains a key challenge for a range of high performance scientific simulations. One approach for reducing data movement and improving performance is by exposing and exploiting structure in a problem through the use of robust structured multilevel solvers. By choosing coarsening that preserves the structure of the problem, these methods maintain efficient structured computation and communication throughout the multigrid hierarchy. However, when coarsening is not permitted to be dependent on the operator, anisotropy must be addressed by the smoother — producing error compatible for coarse-grid correction with structured coarsening. Here, the components required in a scalable parallel structured solver are described with a focus on memory and communication efficiency of robust smoothers. While the implementation of communication and memory reduction techniques in smoothers integrated in a complete 3D solver present a significant engineering challenge, a novel approach is proposed that addresses these challenges systematically through a change to the solver’s execution model. Enabled by user-level threading paired with a set of data and communication abstractions, this approach permits seamless aggregation of communication in plane smoothers — directly reusing code for a 2D distributed multilevel cycle. Results show an effective reduction in communication costs for coarse-grid problems, and result in a speedup of 8.7x in smoothing routines shown in Fig. 12 using this approach. This produces a significant improvement to strong scalability while maintaining favorable weak scaling behavior. Finally, a parallel scaling study using a series of refined meshes is included that demonstrates the effectiveness of this approach in an application of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Handing off the outcome of binary neutron star mergers for accurate and long-term postmerger simulations

In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Neutronics Analysis of Heating and DPA of the Proton Beam Window and Core Vessel Nozzle Neutronics Heating Analysis for the Second Target Station

This reports details the analysis of the heating rates, dpa values and He production in the Proton Beam Window (PBW) assembly and the Core Vessel (CV) nozzle for the preliminary design for the Second Target Station (STS). This analysis focuses on the heating rates in the cooled beltline component of the CV. Neutronics calculations use a detailed unstructured mesh (UM) converted from a CAD model in SpaceClaim of the PBW and CV design of February 2024. The rest of the geometry is taken from the STS master model in constructive solid geometry (CSG). Heating is evaluated for four beam profiles, which are chosen as the most likely to produce the highest heating in the beltline component. The maximum energy deposition values in the beltline are 0.02-0.08 J/cc/pulse, depending on the beam profile. Detailed point cloud data of the full UM geometry is made available for mechanical analysis. The maximum dpa rate in the beltline steel is approximately 0.02 dpa/year, which corresponds to a lifetime of 3,000 years. These results are in good agreement with previous results for the CV. For the PBW itself, maximal heating, dpa, and He production rates are very similar to the previous analysis.

43 PARTICLE ACCELERATORS↗