Search NASA⌕ Search

SEARCH · Search NASA

Results for “mesh data structures”

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 343 records · Page 19

Molecular-dynamics study of amorphous SiO{sub 2} relaxation

Using Molecular-Dynamics simulation we observed the generation of amorphous SiO{sub 2} target from a randomly distributed Si and O atoms. We applied a sequence of annealing of the target with various temperature and quenching to room temperature. The relaxation time required by the system to form SiO{sub 4} tetrahedral mesh after a relatively long simulation time, is studied. The final amorphous target was analyzed using the radial distribution function method, which can be compared with the available theoretical and experimental data. We found that up to 70% of the target atoms form the tetrahedral SiO{sub 4} molecules. The number of formed tetrahedral increases following the growth function and the rate of SiO{sub 4} formation follows Arrhenius law, depends on the annealing temperature. The local structure of amorphous SiO{sub 2} after this treatment agrees well with those reported in some literatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MOSCATO Development and Integration in Fiscal Year 2024

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Investigation of vortex breakdown on a delta wing using Euler and Navier-Stokes equations

A numerical investigation of leading-edge vortex breakdown on a delta wing at high angles of attack is presented. The analysis has been restricted to low speed flows on a flat plate wing with sharp leading edges. Both Euler and Navier-Stokes (assuming fully laminar and fully turbulent flows) equations have been used in this study and the results are compared against experimental data. Predictions of vortex breakdown progression with angle of attack with both Euler and Navier-Stokes equations are shown to be consistent with the experimental data. However, the Navier-Stokes predictions show significant improvements in breakdown location at angles of attack where the vortex breakdown approaches the wing apex. The location of the primary vortex and the level of vorticity in the pre-breakdown regions in these flowfield solutions are affected very little by the viscous effects, even though the Navier-Stokes solutions exhibit viscous phenomena such as secondary and tertiary vortices which are not present in the Euler solutions. In the post-breakdown regions, however, the levels of vorticity in the primary vortex have increased differences between the Euler and Navier-Stokes solutions at comparable locations. Navier-Stokes indicates the presence of a secondary vortex even after the primary vortex is burst. The predicted trajectories of the primary vortex are in very good agreement with the test data, the laminar solutions providing the overall best comparison. The Euler shows a small displacement of the primary vortex, relative to experiment, due to the lack of secondary vortices. The turbulent Navier-Stokes, in general, fall between the Euler and laminar solutions. These findings are based on solutions from meshes that are not usually considered fine enough for resolving vortical flows. To further understand the vortex breakdown phenomenon and the fine details of the vortical flow structure, further calculations with the finer meshes and improved turbulence models are necessary.

S. Agrawal↗

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↗

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↗

Self-Excited, Multi-kHz Dynamics in a Linear, Semi-Bounded Detonation Channel

A linear, semi-bounded detonation channel was studied to characterize the injection, mixing, and ignition process in a two dimensional analog of a rotating detonation wave combustor (RDWC). The linear channel was developed to provide large optical accessibility in support of high-fidelity diagnostics. The canonical geometry was also chosen to support the generation of future numerical validation. A high degree of modularity was designed into the combustor to allow for testing multiple injector con figurations, and for reconfiguring the high-frequency instrumentation and optical window locations. Initial experiments have shown the presence of self-excited, multi-kHz combustion waves that bear the distinct characteristics of detonation waves. High-frequency pressure data details the dynamic response to detonation propagation in the combustor. The role of manifold-chamber coupling is investigated as the mechanism behind the formation of self-excited detonation waves in the combustor. We explore the formation of the waves using Schlieren imaging techniques and examine the structures present within the reaction zone and the interface between reactants and products. Companion computations were carried out using an in-house CFD solver GEMS(General Equations and Mesh Solver) to develop understanding of the structure and dynamics of the flow, flame, and their coupled interactions.

Schwinn, Kyle S.↗

Burner liner thermal/structural load modelling

The objective of this program is to develop a thermal data transfer computer program module for the Burner Liner Thermal-Structural Load Modelling Program. This will be accomplished by reviewing existing methodologies for thermal data transfer and selecting three heat transfer codes for application in this program, evaluating the selected codes to establish criteria for developing a computer program module to transfer thermal data from the heat transfer codes to selected stress analysis codes, developing the automated thermal load transfer module, and verifying and documenting the module. The overall objectives of this thermal transfer module are that it handle independent mesh configurations, perform the transfer in an accurate and efficient fashion and that the total system be flexible for future improvements.

Maffeo, R. J.↗

Numerical investigation of chemically reacting flows in ramjet dump combustors

The time-dependent Navier-Stokes equations, including second-order turbulence model, are numerically integrated by using four-stage Runge-Kutta scheme to predict the steady-state supersonic flow structures in ramjet dump combustors. The formulation is derived for reacting flows with finite-rate chemistry. In the present study, it is firstly attempted to assess the accuracy of existing high-order turbulence model in supersonic flows. The comparison shows reasonable agreement between calculated and measured data in terms of velocity distributions. It is indicated that a modified constant C-mu for calculating turbulent eddy viscosity is needed in the supersonic flow regime and the adaptive meshing is preferred to capture the recirculation zone. In the reacting flow calculation, the results from a test case of hydrogen and air combustion at premixed conditon show that the rearward facing step is able to increase flow residence time and stabilize the flame in supersonic flows.

Hsieh, Kwang-Chung↗

3D Navier-Stokes Time Accurate Solutions Using Multipartitioning Parallel Computation Methodology

A parallel CFD code solving 3D time accurate Navier-Stokes equations with multipartitioning parallel Methodology is being developed in collaboration with Ohio State University within the Air Vehicle Directorate, at Wright Patterson Air Force Base. The advantage of the multipartitioning parallel method is that the domain decomposition will not introduce domain boundaries for the implicit operators. A ring structure data communication is employed so that the implicit time accurate method can be implemented for multi-processors with the same accuracy as for the single processor. No sub-iteration is needed at the domain boundaries. The code has been validated for some typical unsteady flows, which include Coutte Flow, flow passing a cylinder. The code now is being employed for a large scale time accurate wall jet transient flow computation. 'ne preliminary results are promising. The mesh has been refined to capture more details of the flow field. The mesh refinement computation is in progress and would be difficult to successfully implement without the parallel computation techniques used. A modified version of the code with more efficient inversion of the diagonalized block matrix is currently being tested.

Zha, Ge-Cheng↗

Structural Analysis Using Computer Based Methods

The stiffness of a flex hose that will be used in the umbilical arms of the Space Launch Systems mobile launcher needed to be determined in order to properly qualify ground umbilical plate behavior during vehicle separation post T-0. This data is also necessary to properly size and design the motors used to retract the umbilical arms. Therefore an experiment was created to determine the stiffness of the hose. Before the test apparatus for the experiment could be built, the structure had to be analyzed to ensure it would not fail under given loading conditions. The design model was imported into the analysis software and optimized to decrease runtime while still providing accurate rest~lts and allow for seamless meshing. Areas exceeding the allowable stresses in the structure were located and modified before submitting the design for fabrication. In addition, a mock up of a deep space habitat and the support frame was designed and needed to be analyzed for structural integrity under different loading conditions. The load cases were provided by the customer and were applied to the structure after optimizing the geometry. Once again, weak points in the structure were located and recommended design changes were made to the customer and the process was repeated until the load conditions were met without exceeding the allowable stresses. After the stresses met the required factors of safety the designs were released for fabrication.

Dietz, Matthew R.↗

Computer models of the spacecraft wake

Until recently, computations of space plasma flow over a spacecraft have been unstable for ratios of spacecraft dimension to Debye length typical of the low Earth orbit environment. Calculations are presented of the spacecraft/environment interaction based on two computer codes, MACH and POLAR. MACH, an inside-out particle tracking code, was developed for the purpose of validating the physics of POLAR in regimes where these are no comprehensive theoretical or experimental results. While the spacecraft which can be treated by MACH are restricted to simple geometries, the methodology is more fundamental than POLAR. MACH generates self-consistent solutions within the context of quasisteady Vlasov plasma flow and achieves Debye ratios previously unobtainable. POLAR uses a three-dimensional finite-element representation of the vehicle in a staggered mesh. The plasma sheath is modeled by outside-in particle tracking. Solutions for the plasma flow, wake and vehicle charging are obtained by Vlasov-Poisson iteration; charge stabilization techniques make the results virtually insensitive to the Debye ratio. POLAR reproduces the Laframboise static plasma solutions for sperical probes and fits the Makita-Kuriki probe data for spheres in a flowing plasma in regions where comparisons are valid. POLAR and MACH solutions for the particle and electrostatic potential structure of the wake of a charged disk in a low-altitude flow are shown for Mach numbers 4, 5, and 8. New features of the solutions include ion focussing in the wake and a definitive determination of the sheath edge in the wake which shows that the sheath is not an equipotential.

Rubin, A. G.↗

Automated segmentation of soft X-ray tomography: Native cellular structure with submicron resolution at high-throughput for whole-cell quantitative imaging in yeast

Soft X-ray tomography (SXT) is an invaluable tool for quantitatively analyzing cellular structures at suboptical isotropic resolution. However, it has traditionally depended on manual segmentation, limiting its scalability for large datasets. Here, we leverage a deep learning-based autosegmentation pipeline to segment and label cellular structures in hundreds of cells across three Saccharomyces cerevisiae strains. This task-based pipeline uses manual iterative refinement to improve segmentation accuracy for key structures, including the cell body, nucleus, vacuole, and lipid droplets, enabling high-throughput and precise phenotypic analysis. Using this approach, we quantitatively compared the three-dimensional (3D) whole-cell morphometric characteristics of wild-type, VPH1-GFP, and vac14 strains, uncovering detailed strain-specific cell and organelle size and shape variations. We show the utility of SXT data for precise 3D curvature analysis of entire organelles and cells and detection of fine morphological features using surface meshes. Our approach facilitates comparative analyses with high spatial precision and statistical throughput, uncovering subtle morphological features at the single-cell and population level. This workflow significantly enhances our ability to characterize cell anatomy and supports scalable studies on the mesoscale, with applications in investigating cellular architecture, organelle biology, and genetic research across diverse biological contexts.

Chen, Jianhua [Lawrence Berkeley National Laborato↗

Using PAFEC as a preprocessor for COSMIC/NASTRAN

Programs for Automatic Finite Element Calculations (PAFEC) is a general purpose, three dimensional linear and nonlinear finite element program (ref. 1). PAFEC's features include free format input utilizing engineering keywords, powerful mesh generating facilities, sophisticated data base management procedures, and extensive data validation checks. Presented here is a description of a software interface that permits PAFEC to be used as a preprocessor for COSMIC/NASTRAN. This user friendly software, called PAFCOS, frees the stress analyst from the laborious and error prone procedure of creating and debugging a rigid format COSMIC/NASTRAN bulk data deck. By interactively creating and debugging a finite element model with PAFEC, thus taking full advantage of the free format engineering keyword oriented data structure of PAFEC, the amount of time spent during model generation can be drastically reduced. The PAFCOS software will automatically convert a PAFEC data structure into a COSMIC/NASTRAN bulk data deck. The capabilities and limitations of the PAFCOS software are fully discussed in the following report.

Gray, W. H.↗

Analytical and experimental study of vibrations in a gear transmission

An analytical simulation of the dynamics of a gear transmission system is presented and compared to experimental results from a gear noise test rig at the NASA Lewis Research Center. The analytical procedure developed couples the dynamic behaviors of the rotor-bearing-gear system with the response of the gearbox structure. The modal synthesis method is used in solving the overall dynamics of the system. Locally each rotor-gear stage is modeled as an individual rotor-bearing system using the matrix transfer technique. The dynamics of each individual rotor are coupled with other rotor stages through the nonlinear gear mesh forces and with the gearbox structure through bearing support systems. The modal characteristics of the gearbox structure are evaluated using the finite element procedure. A variable time steping integration routine is used to calculate the overall time transient behavior of the system in modal coordinates. The global dynamic behavior of the system is expressed in a generalized coordinate system. Transient and steady state vibrations of the gearbox system are presented in the time and frequency domains. The vibration characteristics of a simple single mesh gear noise test rig is modeled. The numerical simulations are compared to experimental data measured under typical operating conditions. The comparison of system natural frequencies, peak vibration amplitudes, and gear mesh frequencies are generally in good agreement.

Choy, F. K.↗

Towards district scale flood simulations using conventional and anisotropic porosity shallow water models with high-resolution topographic information

Current topographic survey technology provides high-resolution (HR) datasets for urban environments. Incorporating this HR information in models aiming to provide flood risk assessment is desirable because the flood wave propagation is depending on the urban topographic features, i.e. buildings, bridges and street networks. Conceptual, numerical and practical challenges arise from the application of shallow water models to HR urban flood modeling. For example, numerical challenges are occurrence of wet-dry fronts, geometric discontinuities in the urban environment and discontinuous solutions, i.e. shock waves. These challenges can be overcome by using a Godunov-type scheme. However, the computational cost of this type of schemes is high, such that HR two-dimensional shallow water simulations with practical relevance have to be run on supercomputers. The porous shallow water model is an alternative approach that aims to reduce computational cost by using a coarse resolution and accounting for unresolved processes by means of the porosity terms. Usually, a speedup between two and three orders of magnitude in comparison to HR simulations can be obtained. This study presents preliminary results of a practical test case concerning pluvial flooding in a district of the city of Nice, France, caused by the intense rainfall event on October 3rd, 2015. HR topography data set on a 1 m resolution is available for the district, whereby street features of infra-metric dimensions have been included. A reference solution is calculated by a HR shallow water model on a 1 m by 1 m structured computational grid. The porous shallow water model is run on a 10 m by 10 m grid and the influence of the drag source term is studied. The model calculations show a large deviation, which is caused by the poor meshing strategy of the porous shallow water (AP) model. The study also summarizes practical challenges that arise during the application of the AP and HR models to a large urban catchment. The main difficulty is to obtain a good mesh. In smaller scale investigations, the mesh is currently constructed by hand such that the cell edges align with buildings. This approach is not feasible for large scale urban catchments with a large number of buildings. Future steps that have to be taken, such as a strategy for automatic mesh generation, are reported on.

42 ENGINEERING↗

Validation and Simulation of ARES I Scale Model Acoustic Test -1- Pathfinder Development

The Ares I Scale Model Acoustics Test (ASMAT) is a series of live-fire tests of scaled rocket motors meant to simulate the conditions of the Ares I launch configuration. These tests have provided a well documented set of high fidelity measurements useful for validation including data taken over a range of test conditions and containing phenomena like Ignition Over-Pressure and water suppression of acoustics. To take advantage of this data, a digital representation of the ASMAT test setup has been constructed and test firings of the motor have been simulated using the Loci/CHEM computational fluid dynamics software. Within this first of a series of papers, results from ASMAT simulations with the rocket in a held down configuration and without water suppression have then been compared to acoustic data collected from similar live-fire tests to assess the accuracy of the simulations. Detailed evaluations of the mesh features, mesh length scales relative to acoustic signals, Courant-Friedrichs-Lewy numbers, and spatial residual sources have been performed to support this assessment. Results of acoustic comparisons have shown good correlation with the amplitude and temporal shape of pressure features and reasonable spectral accuracy up to approximately 1000 Hz. Major plume and acoustic features have been well captured including the plume shock structure, the igniter pulse transient, and the ignition overpressure. Finally, acoustic propagation patterns illustrated a previously unconsidered issue of tower placement inline with the high intensity overpressure propagation path.

Putnam, G. C.↗

Micro-CT analysis and mechanical properties of Ti spherical and polyhedral void composites made with saccharose as a space holder material

The paper describes the formation, morphology and mechanical properties of Ti void composites. The Ti void composites were made using 100 and 325 mesh Ti powder for solid scaffold formation. The spherical and polyhedral voids (pores) were formed using saccharose particles (table sugar) of different shapes. The Ti void composite morphology was investigated by microcomputed tomography and scanning electron microscopy. The Ti void composites of designed porosity of 50–70% were made. Compression test was applied for mechanical properties estimation. It has been found, that Ti void composites made from 100 mesh Ti and those having spherical pores have a higher strength and elastic modulus, i.e. for the designed porosity of 50% for 100 and 325 mesh Ti void composites, a compressive strength was 32.32 and 20.13 MPa, respectively. It has been shown that this is related to better sintering of the 100 mesh Ti powders compared with the 325 mesh Ti powders. A correlation between microcomputed tomography data and mechanical properties has also been shown. The Ti void composites, made with the use of saccharose as a space holder, described in this work should be a promising material for biomedical applications, where interconnected pores and good mechanical properties are required. - Highlights: • Ti scaffolds of the porosity of 50–70% were made. • Saccharose particles as space holder were applied. • The voids in the scaffolds were designed with spherical and polyhedral shape. • The scaffold structure was investigated by SEM and micro-CT. • Micro-CT data and mechanical properties of the Ti scaffold have been correlated.

36 MATERIALS SCIENCE↗

Parallel Block Structured Adaptive Mesh Refinement on Graphics Processing Units

Block-structured adaptive mesh refinement is a technique that can be used when solving partial differential equations to reduce the number of zones necessary to achieve the required accuracy in areas of interest. These areas (shock fronts, material interfaces, etc.) are recursively covered with finer mesh patches that are grouped into a hierarchy of refinement levels. Despite the potential for large savings in computational requirements and memory usage without a corresponding reduction in accuracy, AMR adds overhead in managing the mesh hierarchy, adding complex communication and data movement requirements to a simulation. In this paper, we describe the design and implementation of a native GPU-based AMR library, including: the classes used to manage data on a mesh patch, the routines used for transferring data between GPUs on different nodes, and the data-parallel operators developed to coarsen and refine mesh data. We validate the performance and accuracy of our implementation using three test problems and two architectures: an eight-node cluster, and over four thousand nodes of Oak Ridge National Laboratory’s Titan supercomputer. Our GPU-based AMR hydrodynamics code performs up to 4.87× faster than the CPU-based implementation, and has been scaled to over four thousand GPUs using a combination of MPI and CUDA.

97 MATHEMATICS AND COMPUTING↗