SEARCH · Search NASA

Results for “meshing framework”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Moments-based interface reconstruction, remap and advection

Here, we present a new moment-of-fluid (MOF 2 ) interface reconstruction method. It uses the zeroth, first, and second moments of the fragment of material inside a cell of the mesh to reconstruct a convex material polygon or a union of convex polygons that approximate the respective material fragment. The new method requires information about the material moments only for the cell under consideration. The MOF 2 method allows to exactly reproduce several convex shapes: corners, filaments, and some concave shapes: cell-complements to corners and filaments. Interface reconstruction is formulated as a local (for each cell), non-linear, equality constrained optimization problem, which does not require additional communication and allows for an efficient parallel implementation. We present an extensive set of test problems, both for interface reconstruction on a single cell, and for reconstruction of a variety of shapes on a variety of meshes. We describe how to perform two-material advection using the MOF 2 method and present the results for the classical advection tests. We also show the examples of material interface remapping needed in the framework of multi-material arbitrary Lagrangian-Eulerian methods, and give a brief description of a procedure that can be used to update the material moments on the Lagrangian stage of those methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Hydrothermal organic synthesis experiments

The serious scientific debate about spontaneous generation which raged for centuries reached a climax in the nineteenth century with the work of Spallanzani, Schwann, Tyndall, and Pasteur. These investigators demonstrated that spontaneous generation from dead organic matter does not occur. Although no aspects of these experiments addressed the issue of whether organic compounds could be synthesized abiotically, the impact of the experiments was great enough to cause many investigators to assume that life and its organic compounds were somehow fundamentally different than inorganic compounds. Meanwhile, other nineteenth-century investigators were showing that organic compounds could indeed be synthesized from inorganic compounds. In 1828 Friedrich Wohler synthesized urea in an attempt to form ammonium cyanate by heating a solution containing ammonia and cyanic acid. This experiment is generally recognized to be the first to bridge the artificial gap between organic and inorganic chemistry, but it also showed the usefulness of heat in organic synthesis. Not only does an increase in temperature enhance the rate of urea synthesis, but Walker and Hambly showed that equilibrium between urea and ammonium cyanate was attainable and reversible at 100 C. Wohler's synthesis of urea, and subsequent syntheses of organic compounds from inorganic compounds over the next several decades dealt serious blows to the 'vital force' concept which held that: (1) organic compounds owe their formation to the action of a special force in living organisms; and (2) forces which determine the behavior of inorganic compounds play no part in living systems. Nevertheless, such progress was overshadowed by Pasteur's refutation of spontaneous generation which nearly extinguished experimental investigations into the origins of life for several decades. Vitalism was dealt a deadly blow in the 1950's with Miller's famous spark-discharge experiments which were undertaken in the framework of the Oparin and Haldane hypotheses concerning the origin of life. These hypotheses were constructed on some basic assumptions which included a reduced atmosphere, and a low surface temperature for the early Earth. These ideas meshed well with the prevailing hypothesis of the 1940's and 50's that the Earth had formed through heterogeneous accretion of dust from a condensing solar nebula. Miller's experiments were extremely successful, and were followed by numerous other experiments by various investigators who employed a wide variety of energy sources for abiotic synthesis including spark discharges, ultra-violet radiation, heat, shock waves, plasmas, gamma rays, and other forms of energy. The conclusion reached from this body of work is that energy inputs can drive organic synthesis from a variety of inorganic starting materials.

Shock, Everett L.↗

NASA NTRS1999

A High-Order Method Using Unstructured Grids for the Aeroacoustic Analysis of Realistic Aircraft Configurations

A method for the prediction of acoustic scatter from complex geometries is presented. The discontinuous Galerkin method provides a framework for the development of a high-order method using unstructured grids. The method's compact form contributes to its accuracy and efficiency, and makes the method well suited for distributed memory parallel computing platforms. Mesh refinement studies are presented to validate the expected convergence properties of the method, and to establish the absolute levels of a error one can expect at a given level of resolution. For a two-dimensional shear layer instability wave and for three-dimensional wave propagation, the method is demonstrated to be insensitive to mesh smoothness. Simulations of scatter from a two-dimensional slat configuration and a three-dimensional blended-wing-body demonstrate the capability of the method to efficiently treat realistic geometries.

Atkins, Harold L.↗

NASA NTRS2004

Multidimensional Unstructured-Grid Liquid Rocket Engine Nozzle Performance and Heat Transfer Analysis

The objective of this study is to conduct a unified computational analysis for computing design parameters such as axial thrust, convective and radiative wall heat fluxes for regeneratively cooled liquid rocket engine nozzles, so as to develop a computational strategy for computing those parameters through parametric investigations. The computational methodology is based on a multidimensional, finite-volume, turbulent, chemically reacting, radiating, unstructured-grid, and pressure-based formulation, with grid refinement capabilities. Systematic parametric studies on effects of wall boundary conditions, combustion chemistry, radiation coupling, computational cell shape, and grid refinement were performed and assessed. Under the computational framework of this study, it is found that the computed axial thrust performance, flow features, and wall heat fluxes compared well with those of available data and calculations, using a strategy of structured-grid dominated mesh, finite-rate chemistry, and cooled wall boundary condition.

Wang, Ten-See↗

NASA NTRS2012

Integrated Design Engineering Analysis (IDEA) Environment Automated Generation of Structured CFD Grids using Topology Methods

This report documents the work performed from March 2010 to March 2012. The Integrated Design and Engineering Analysis (IDEA) environment is a collaborative environment based on an object-oriented, multidisciplinary, distributed framework using the Adaptive Modeling Language (AML) as a framework and supporting the configuration design and parametric CFD grid generation. This report will focus on describing the work in the area of parametric CFD grid generation using novel concepts for defining the interaction between the mesh topology and the geometry in such a way as to separate the mesh topology from the geometric topology while maintaining the link between the mesh topology and the actual geometry.

Kamhawi, Hilmi N.↗

DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2008

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

NASA NTRS2008

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

Saether, E.↗

NASA NTRS2014

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗

DOE OSTI2020

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗

DOE OSTI2017

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↗

DOE OSTI2019

Architecture Improvements and Retesting of the Microscopic Depletion Implementation in the MAMMOTH Reactor Physics Application

Nuclear fuel experiments conducted at the Transient Reactor Test Facility (TREAT) will be analyzed with the MAMMOTH reactor physics application, currently under development at Idaho National Laboratory. MAMMOTH natively couples with the BISON, RELAP-7, and Rattlesnake applications within the MOOSE framework. The implemented methods allow for simulating the irradiation of fuel in a nuclear reactor from the beginning of life in a nuclear reactor until it is placed in TREAT for fuel testing using the same finite element mesh throughout the analysis retaining a very high level of resolution and fidelity. The calculation of the isotopic distribution in fuel requires the solution to the decay and transmutation equations coupled to the neutron transport equation. This report focuses on three significant improvements to the MAMMOTH micro-depletion capability. First, the depletion architecture within MAMMOTH is modified to improve the flexibility and speed of burnup calculations. The results of various benchmark tests demonstrate that the refactoring of the previously implemented capability does not introduce any regression. The second improvement involves the expansion of the data functionality. MAMMOTH can now read the decay data files used by the Oak Ridge Isotope Generation (ORIGEN) code, thus enabling access to data of more than 1,600 isotopes. The third improvement is the capability to perform depletion at a constant power level which was verified via a comparison to ORIGEN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2013

Integration of the Total Lightning Jump Algorithm into Current Operational Warning Environment Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. In order to become a viable option for operational forecasters to incorporate into their severe storm monitoring process, the total lightning jump must be placed into the framework of several severe storm conceptual models (e.g., radar evolution, storm morphology) which forecasters have built through training and experience. Thus, one of the goals of this study is to examine and relate the lightning jump concept to often used radar parameters (e.g., dBZ vertical structure, VIL, MESH, MESO/shear) in the warning environment. Tying lightning trends and lightning jump occurrences to these radar based parameters will provide forecasters with an additional tool that they can use to build an accurate realtime depiction as to what is going on in a given environment. Furthermore, relating the lightning jump concept to these parameters could also increase confidence in a warning decision they have already made, help tip the scales on whether or not to warn on a given storm, or to draw the forecaster s attention to a particular storm that is rapidly developing. Furthermore the lightning information will add vital storm scale information in regions that are not well covered by radar, or when radar failures occur. The physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relation to updraft strength, updraft volume, precipitation -sized ice mass, etc.; however, very few have related the concept of the lightning jump and manifestation of severe weather to storm dynamics and microphysics using multi -Doppler and polarimetric radar techniques. Therefore, the second half of this study will combine the lightning jump algorithm and these radar techniques in order to place the lightning jump concept into a physical and dynamical framework. This analysis includes examining such parameters as mixed phase precipitation volume, charging zone, updraft strength and updraft volume. Such a study should provide increased understanding of and confidence in the strengths and limitations of the lightning jump algorithm in the storm warning process.

Schultz, Christopher J.↗

DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗

DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗

DOE OSTI2017

CPU Volume Rendering of Adaptive Mesh Refinement Data

Adaptive Mesh Refinement (AMR) methods are widespread in scientific computing, and visualizing the resulting data with efficient and accurate rendering methods can be vital for enabling interactive data exploration. In this work, we detail a comprehensive solution for directly volume rendering block-structured (Berger-Colella) AMR data in the OSPRay interactive CPU ray tracing framework. In particular, we contribute a general method for representing and traversing AMR data using a kd-tree structure, and four different reconstruction options, one of which in particular (the basis function approach) is novel compared to existing methods. Here, we demonstrate our system on two types of block-structured AMR data and compressed scalar field data, and show how it can be easily used in existing production-ready applications through a prototypical integration in the widely used visualization program ParaView.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗