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 973 records · Page 54

Collision of Physics and Software in the Monte Carlo Application Toolkit (MCATK)

The topic is presented in a series of slides organized as follows: MCATK overview, development strategy, available algorithms, problem modeling (sources, geometry, data, tallies), parallelism, miscellaneous tools/features, example MCATK application, recent areas of research, and summary and future work. MCATK is a C++ component-based Monte Carlo neutron-gamma transport software library with continuous energy neutron and photon transport. Designed to build specialized applications and to provide new functionality in existing general-purpose Monte Carlo codes like MCNP, it reads ACE formatted nuclear data generated by NJOY. The motivation behind MCATK was to reduce costs. MCATK physics involves continuous energy neutron & gamma transport with multi-temperature treatment, static eigenvalue (k eff and α) algorithms, time-dependent algorithm, and fission chain algorithms. MCATK geometry includes mesh geometries and solid body geometries. MCATK provides verified, unit-test Monte Carlo components, flexibility in Monte Carlo application development, and numerous tools such as geometry and cross section plotters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

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 corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Importance of turbulence-chemistry interactions at low temperature engine conditions

The role of turbulence-chemistry interaction in autoignition and flame stabilization is investigated for spray flames at low temperature combustion (LTC) conditions by performing high-fidelity three-dimensional computational fluid dynamics (CFD) simulations. A recently developed Tabulated Flamelet Model (TFM) is coupled with a large eddy simulation (LES) framework and validated across a range of Engine Combustion Network (ECN) ambient temperature conditions for n-dodecane fuel. High resolution grids with 0.0625 mm minimum cell size and 25 million total cell count are implemented using adaptive mesh refinement over the spray and combustion regions. Simulations with these grids and multiple LES realizations, with a 103 species n-dodecane mechanism show good agreement with experimental data for all the ambient conditions investigated. This modeling approach with the computational cost advantage of tabulated chemistry is then extended towards understanding the auto-ignition and flame stabilization at an ambient temperature of 750 K. These low temperature conditions lead to substantially higher ignition delays and flame liftoff lengths, and significantly leaner combustion compared to conventional high temperature diesel combustion. These conditions also require the simulations to span significantly larger temporal and spatial dimensions thereby increasing the computational cost. The TFM approach is able to capture autoignition and flame liftoff length at the low temperature conditions. Significant differences with respect to mixing, species formation and flame stabilization are observed under low temperature compared to conventional diesel combustion. At higher ambient temperatures, formation of formaldehyde is observed in the rich region (phi > 1) followed by the formation of OH in the stoichiometric regions. Under low temperature conditions, formaldehyde is observed to form at leaner regions followed by the onset of OH formation in significantly lean regions of the flame. Qualitative differences between species formation and transient flame development for the high and low temperature conditions are presented. The two stage ignition process is further investigated by studying the species formation in mixture fraction space by solving 1D flamelet equations for different scalar dissipation rates and homogeneous reactor assumption. Results show that scalar dissipation causes these radicals to diffuse within the mixture fraction space. As a result, this significantly enhances ignition and plays a dominant role at such low temperature conditions which cannot be captured by the homogeneous reaction assumption based model.

33 ADVANCED PROPULSION SYSTEMS↗

A new class of finite element variational multiscale turbulence models for incompressible magnetohydrodynamics

New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and high Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ACTIVE GALACTIC NUCLEUS FEEDBACK AND ENTROPY INJECTION IN GALAXY CLUSTER CORES

We make the first estimate of non-gravitational energy profiles in galaxy cluster cores (and beyond) based on observational data. Comparing the observed entropy profiles within r{sub 500}, from the Representative XMM-Newton Cluster Structure Survey to simulated base entropy profiles without feedback from both adaptive mesh refinement (AMR) and smoothed particle hydrodynamic (SPH) non-radiative simulations, we estimate the amount of additional non-gravitational energy, E{sub ICM}, contained in the intracluster medium (ICM), as well as the total energy feedback, E{sub Feedback}, from active galactic nuclei (AGNs; the central AGNs in most cases) into the clusters. The total feedback energy scales with the mean spectroscopic temperature as E{sub Feedback}∝T{sub sp}{sup 2.52±0.08} and E{sub Feedback}∝T{sub sp}{sup 2.17±0.11} for the SPH and AMR baseline profiles. The mean non-gravitational energy per particle within r{sub 500} remaining in the ICM after energy lost during cooling is ε{sub ICM} = 2.8 ± 0.8 keV for the SPH theoretical relation and ε{sub ICM} = 1.7 ± 0.9 keV for the AMR theoretical relation. We use the NRAO/VLA Sky Survey source catalog to determine the radio luminosity, L{sub R} , at 1.4 GHz of the central source(s) of our sample. For T{sub sp} > 3 keV, the E{sub Feedback} correlates with L{sub R} , although with different normalization for cool-core and non-cool-core clusters. We show that AGNs could provide a significant portion of the feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of Load Balancing on Unstructured Adaptive Grid Computations for Distributed-Memory Multiprocessors

The computational requirements for an adaptive solution of unsteady problems change as the simulation progresses. This causes workload imbalance among processors on a parallel machine which, in turn, requires significant data movement at runtime. We present a new dynamic load-balancing framework, called JOVE, that balances the workload across all processors with a global view. Whenever the computational mesh is adapted, JOVE is activated to eliminate the load imbalance. JOVE has been implemented on an IBM SP2 distributed-memory machine in MPI for portability. Experimental results for two model meshes demonstrate that mesh adaption with load balancing gives more than a sixfold improvement over one without load balancing. We also show that JOVE gives a 24-fold speedup on 64 processors compared to sequential execution.

Biswas, Rupak↗

Impact of Load Balancing on Unstructured Adaptive Grid Computations for Distributed-Memory Multiprocessors

The computational requirements for an adaptive solution of unsteady problems change as the simulation progresses. This causes workload imbalance among processors on a parallel machine which, in turn, requires significant data movement at runtime. We present a new dynamic load-balancing framework, called JOVE, that balances the workload across all processors with a global view. Whenever the computational mesh is adapted, JOVE is activated to eliminate the load imbalance. JOVE has been implemented on an IBM SP2 distributed-memory machine in MPI for portability. Experimental results for two model meshes demonstrate that mesh adaption with load balancing gives more than a sixfold improvement over one without load balancing. We also show that JOVE gives a 24-fold speedup on 64 processors compared to sequential execution.

Sohn, Andrew↗

Computation of High Speed Jet Noise

The objective of this work was to use direct numerical simulation (DNS) techniques to study the physics of noise generation by a high-speed turbulent jet. A Mach 0.9, Reynolds number 3,600 jet was selected because of available experimental data. New numerical methods for generating disturbances at the nozzle and computing far-field sound were developed and reported in the course of this work. Over 25 million mesh points were used in the simulations which ran for over 50,000 timesteps and required over 50,000 processor hours on state-of-the-art parallel computer systems to complete. Figures show a visualization of the jet and sound field, a comparison of the mean flow development with the experiment, a directivity comparison with the experiment, and time spectrum comparison with the experiment. Agreement is seen to be excellent. These are fully document in the attached references. Full details of the work, detailed achievements and conclusions are discussed in appendices, which are copies of publications that resulted from this work. We have studied noise mechanisms in supersonic jets, the refraction of sound by turbulence in subsonic jets, and noise sources in conjunction with a DNS of a Mach 0.9 jet.

Freund, Jonathan B.↗

Method for improved prediction of bone fracture risk using bone mineral density in structural analysis

A non-invasive in-vivo method of analyzing a bone for fracture risk includes obtaining data from the bone such as by computed tomography or projection imaging which data represents a measure of bone material characteristics such as bone mineral density. The distribution of the bone material characteristics is used to generate a finite element method (FEM) mesh from which load capability of the bone can be determined. In determining load capability, the bone is mathematically compressed, and stress, strain force, force/area versus bone material characteristics are determined.

Cann, Christopher E.↗

Continuum Damage Mechanics Models for the Analysis of Progressive Failure in Open-Hole Tension Laminates

The performance of a state-of-the-art continuum damage mechanics model for interlaminar damage, coupled with a cohesive zone model for delamination is examined for failure prediction of quasi-isotropic open-hole tension laminates. Limitations of continuum representations of intra-ply damage and the effect of mesh orientation on the analysis predictions are discussed. It is shown that accurate prediction of matrix crack paths and stress redistribution after cracking requires a mesh aligned with the fiber orientation. Based on these results, an aligned mesh is proposed for analysis of the open-hole tension specimens consisting of different meshes within the individual plies, such that the element edges are aligned with the ply fiber direction. The modeling approach is assessed by comparison of analysis predictions to experimental data for specimen configurations in which failure is dominated by complex interactions between matrix cracks and delaminations. It is shown that the different failure mechanisms observed in the tests are well predicted. In addition, the modeling approach is demonstrated to predict proper trends in the effect of scaling on strength and failure mechanisms of quasi-isotropic open-hole tension laminates.

Song, Kyonchan↗

Assessment of the Fluid Network for the Effusion Holes in a RQL (Rich-burn, Quick-quench, Lean-burn) Gas Turbine Combustor

The transient processes involving lean blowout (LBO) in a RQL combustor such as the single cup combustor in the National Jet Fuels Combustion Program (NJFCP) referee rig are simulated with the reduced mechanisms based on Hybrid Chemistry for three different fuels, A2, C1 and C4. Experiments conducted in the NJFCP referee rig established a stable flame at a global equivalence ratio of 0.096 (the near LBO condition) as the starting condition for LBO tests performed with all A and C fuels. Simulations of the approach to LBO verify that a stable flame is established at the near-LBO condition for each fuel before reducing fuel flow rate (and thus equivalence ratio) in a stepwise manner. Main conclusions are: (1) For A2 fuel the equivalence ratio for LBO is 0.081 (2) For C1 fuel the equivalence ratio for LBO is 0.087 or 0.085depended on different step-down transition paths. (3) For C4 fuel the equivalence ratio for LBO is 0.087 or 0.084 depended on types of mesh or boundary conditions. (4) The fuel relative difference from LES on LBO for C1 fuel relative to A2 fuel is 7.40% or 4.94% versus 7.83% from the experimental data. (5) The fuel relative difference from LES on LBO for C4 fuel relative to A2 is 6.17% or 7.69% versus 5.21% from the experimental data

Lean blow out↗

VULCAN-CFD User Manual: Ver. 7.2.0

VULCAN-CFD offers a comprehensive set of capabilities to enable the simulation of continuum flowfields from subsonic to hypersonic conditions. The governing equations that are employed include allowances for both chemical and thermal nonequilibrium processes, coupled with a wide variety of turbulence models for both Reynolds-averaged and large eddy simulations. The software package can simulate two-dimensional, axisymmetric, or three-dimensional problems on structured multiblock meshes or fully unstructured meshes. A parabolic (i.e., space-marching) treatment can also be used for any subset of a structured mesh that can accommodate this solution strategy. The flow solver provides a significant level of geometric flexibility for structured grid simulations by allowing for arbitrary face-to-face C(0) continuous and non-C(0) continuous block interface connectivities. The unstructured grid paradigm allows for mixed element unstructured meshes that contain any combination of tetrahedral, prismatic, pyramidal, and hexahedral cell elements. The flow solver is also fully parallelized using MPI (Message Passing Interface) libraries in a data-parallel fashion, allowing for efficient simulations on modern High Performance Computing (HPC) systems. This document provides information related to the installation and execution of the VULCAN-CFD software package. A detailed description of the physical and numerical models available in the software are provided in the VULCAN-CFD Theory Manual.

VULCAN-CFD User Manual↗

Assessment of the Fluid Network for the Effusion Holes in a RQL (Rich-burn, Quick-quench, Lean-burn) Gas Turbine Combustor

The transient processes involving lean blowout (LBO) in a RQL combustor such as the single cup combustor in the National Jet Fuels Combustion Program (NJFCP) referee rig are simulated with the reduced mechanisms based on Hybrid Chemistry for three different fuels, A2, C1 and C4. Experiments conducted in the NJFCP referee rig established a stable flame at a global equivalence ratio of 0.096 (the near LBO condition) as the starting condition for LBO tests performed with all A and C fuels. Simulations of the approach to LBO verify that a stable flame is established at the near-LBO condition for each fuel before reducing fuel flow rate (and thus equivalence ratio) in a stepwise manner. Main conclusions are: (1) For A2 fuel the equivalence ratio for LBO is 0.081 (2) For C1 fuel the equivalence ratio for LBO is 0.087 or 0.085 depended on different step-down transition paths. (3) For C4 fuel the equivalence ratio for LBO is 0.087 or 0.084 depended on types of mesh or boundary conditions. (4) The fuel relative difference from LES on LBO for C1 fuel relative to A2 fuel is 7.40% or 4.94% versus 7.83% from the experimental data. (5) The fuel relative difference from LES on LBO for C4 fuel relative to A2 is 6.17% or 7.69% versus 5.21% from the experimental data.

Lean Blow out↗

Thermal behavior spiral bevel gears

An experimental and analytical study of the thermal behavior of spiral bevel gears is presented. Experimental data were taken using thermocoupled test hardware and an infrared microscope. Many operational parameters were varied to investigate their effects on the thermal behavior. The data taken were also used to validate the boundary conditions applied to the analytical model. A finite element-based solution sequence was developed. The three-dimensional model was developed based on the manufacturing process for these gears. Contact between the meshing gears was found using tooth contact analysis to describe the location, curvatures, orientations, and surface velocities. This information was then used in a three-dimensional Hertzian contact analysis to predict contact ellipse size and maximum pressure. From these results, an estimate of the heat flux magnitude and the location on the finite element model was made. The finite element model used time-averaged boundary conditions to permit the solution to attain steady state in a computationally efficient manner.Then time- and position-varying boundary conditions were applied to the model to analyze the cyclic heating and cooling due to the gears meshing and transferring heat to the surroundings, respectively. The model was run in this mode until the temperature behavior stabilized. The transient flash temperature on the surface was therefore described. The analysis can be used to predict the overall expected thermal behavior of spiral bevel gears. The experimental and analytical results were compared for this study and also with a limited number of other studies. The experimental and analytical results attained in the current study were basically within 10% of each other for the cases compared. The experimental comparison was for bulk thermocouple locations and data taken with an infrared microscope. The results of a limited number of other studies were compared with those obtained herein and predicted the same basic behavior.

Handschuh, Robert F.↗

Experiments and Simulation of a Borehole in Salt to Understand Heat, Brine, and Vapor Migration - 19192

Disposition of heat-generating nuclear waste (HGNW) remains a continuing technical and sociopolitical challenge. Numerous concepts for HGNW management have been proposed and examined internationally, including an extensive focus on geologic disposal. One proposed geologic material is salt because of its low permeability and viscoplastic deformation that causes self-repair of damage done to the salt by waste emplacement activities. Evaluating the safety and technical challenges of storing HGNW in a salt repository is an ongoing process involving experiments and supporting numerical simulation. Currently an experiment is underway at the Waste Isolation Pilot Plant (WIPP) to explore how the presence of a heat generating source affects phenomena such as brine migration, vapor transport, and mechanical changes to the bedded salt. A sub-horizontal heated borehole test is in progress in the underground at the WIPP that includes a centrally located 10.2 cm diameter borehole with an adjustable heater surrounded with smaller diameter boreholes instrumented with thermocouples. The central borehole contains an inflatable packer, heating block, brine sampler, and constantly flowing nitrogen gas circulation system. Air-injection tests performed in the central borehole provide pressure measurements that are used to constrain permeability of the system. The steady-state temperatures, as well as the rise and fall of temperature when the heater is cycled on and off, have been measured for up to 60 days. In the borehole, dry nitrogen gas circulation evaporates water and outflows to a desiccant container where water mass is measured daily during the experiment to quantify vapor removal. Although this test is a 'shake-down' for a planned second round of fresh borehole testing, we have gathered a rich dataset. These data allow us to build simulations using the Finite Element Heat and Mass transfer code (FEHM) to evaluate the experimental results, determine field-scale parameters, and identify code improvements to reproduce important physical processes that may not be accounted for at present. A 3-D numerical mesh, built using LaGrit software (lagrit.lanl.gov; Miller et al., 2007), includes increasing resolution around the central borehole. Modeling of the experiment allows for determining the local thermal conductivity and permeability of any damaged bedded salt around the borehole, where damage from drilling may change the permeability, porosity and saturation conditions, by parameter testing and inverse methods. Additionally, assumptions about brine and vapor flow and transport are being tested by comparing measured and simulated results. The combination of experimental data and model results provide additional data to help support the safety case for safe and effective HGNW disposition in bedded salt formations. Initial results from this experiment show that water flow into the borehole is within previous experimental results. Further, we have found that the design of the heater block is restricting energy to flow into the rock salt. Thus, this test has proven useful in design of the next generation experiment where infrared heating may be used to bypass issues caused by air gaps located around the current stainless steel block heater. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Benchmark Data for Evaluation of NASA Impedance Eduction Methods

Results achieved with three acoustic liner configurations are used to evaluate the effects of recent modifications to the NASA Normal Incidence Tube and Grazing Flow Impedance Tube. These include a calibration liner, a wire mesh liner, and a perforate liner. The effects of source type, source level, and mean flow Mach number on the impedances educed with these liners are explored. Existing models are used to predict the impedance for each test condition and are compared with educed impedance spectra to determine their validity. The results suggest that educed impedances are similar for either a stepped sine or swept sine source. The predicted models are acceptable for the perforate liner and are quite good for the calibration and wire mesh liners.

acoustic↗

Benchmark Data for Evaluation of NASA Impedance Eduction Methods

Results achieved with three acoustic liner configurations are used to evaluate the effects of recent modifications to the NASA Normal Incidence Tube and Grazing Flow Impedance Tube. These include a calibration liner, a wire mesh liner, and a perforate liner. The effects of source type, source level, and mean flow Mach number on the impedances educed with these liners are explored. Existing models are used to predict the impedance for each test condition, and are compared with educed impedance spectra to determine their validity. The results suggest that educed impedances are similar for either a stepped sine or swept sine source. The predicted models are acceptable for the perforate liner, and are quite good for the calibration and wire mesh liners.

acoustic↗

An Element-Based Concurrent Partitioner for Unstructured Finite Element Meshes

A concurrent partitioner for partitioning unstructured finite element meshes on distributed memory architectures is developed. The partitioner uses an element-based partitioning strategy. Its main advantage over the more conventional node-based partitioning strategy is its modular programming approach to the development of parallel applications. The partitioner first partitions element centroids using a recursive inertial bisection algorithm. Elements and nodes then migrate according to the partitioned centroids, using a data request communication template for unpredictable incoming messages. Our scalable implementation is contrasted to a non-scalable implementation which is a straightforward parallelization of a sequential partitioner.

concurrent partitioner finite element meshes distr↗