Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Microphysics in Simulations of Nonprecipitating Cumulus

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. Overall, the results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Numerical investigation of the flat band Bloch modes in a 2D photonic crystal with Dirac cones

A numerical method combining complex-k band calculations and absorbing boundary conditions for Bloch waves is presented. We use this method to study photonic crystals with Dirac cones. We demonstrate that the photonic crystal behaves as a zero-index medium when excited at normal incidence, but that the zero-index behavior is lost at oblique incidence due to excitation of modes on the flat band. We also investigate the formation of monomodal and multimodal cavity resonances inside the photonic crystals, and the physical origins of their different line-shape features.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

MIDA: A Multimodal Imaging-Based Detailed Anatomical Model of the Human Head and Neck

Computational modeling and simulations are increasingly being used to complement experimental testing for analysis of safety and efficacy of medical devices. Multiple voxel- and surface-based whole- and partial-body models have been proposed in the literature, typically with spatial resolution in the range of 1–2 mm and with 10–50 different tissue types resolved. We have developed a multimodal imaging-based detailed anatomical model of the human head and neck, named “MIDA”. The model was obtained by integrating three different magnetic resonance imaging (MRI) modalities, the parameters of which were tailored to enhance the signals of specific tissues: i) structural T1- and T2-weighted MRIs; a specific heavily T2-weighted MRI slab with high nerve contrast optimized to enhance the structures of the ear and eye; ii) magnetic resonance angiography (MRA) data to image the vasculature, and iii) diffusion tensor imaging (DTI) to obtain information on anisotropy and fiber orientation. The unique multimodal high-resolution approach allowed resolving 153 structures, including several distinct muscles, bones and skull layers, arteries and veins, nerves, as well as salivary glands. The model offers also a detailed characterization of eyes, ears, and deep brain structures. A special automatic atlas-based segmentation procedure was adopted to include a detailed map of the nuclei of the thalamus and midbrain into the head model. The suitability of the model to simulations involving different numerical methods, discretization approaches, as well as DTI-based tensorial electrical conductivity, was examined in a case-study, in which the electric field was generated by transcranial alternating current stimulation. The voxel- and the surface-based versions of the models are freely available to the scientific community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Autoclave treatment of the classical scrapie agent US No. 13-7 and experimental inoculation to susceptible VRQ/ARQ sheep via the oral route results in decreased transmission efficiency

Scrapie, a prion disease of sheep, is highly resistant to conventional deactivation. Numerous methods to deactivate scrapie have been tested in laboratory animal models, and adequate autoclave treatment can reduce or remove the infectivity of some classical scrapie strains depending on the heating parameters used. In this study, we autoclaved brain homogenate from a sheep with US scrapie strain 13–7 for 30 minutes at 121°C. Genetically susceptible VRQ/ARQ sheep were orally inoculated with 3 grams of the autoclaved brain homogenate. For comparison, a second group of sheep was inoculated with a non-autoclaved brain homogenate. Rectal biopsies were used to assess antemortem scrapie disease progression throughout the study. Five out of ten (5/10) sheep that received autoclaved inoculum ultimately developed scrapie after an experimental endpoint of 72 months. These sheep had a mean incubation period of 26.99 months. Two out of five (2/5) positive sheep had detectable PrP Sc in antemortem rectal biopsies, and two (2/5) other sheep had PrP Sc in postmortem rectal tissue. A single sheep (1/5) was positive for scrapie in the CNS, small intestine, and retropharyngeal lymph node but had negative rectal tissue. All of the sheep (10/10) that received non-autoclaved inoculum developed scrapie with a mean incubation period of 20.2 months and had positive rectal biopsies at the earliest timepoint (14.7 months post-inoculation). These results demonstrate that sheep are orally susceptible to US derived classical scrapie strain 13–7 after autoclave treatment at 121°C for 30 minutes. Differences in incubation periods and time interval to first positive rectal biopsies indicate a partial reduction in infectivity titers for the autoclaved inoculum group.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Simulating Fragmentation and Fluid-Induced Fracture in Disordered Media Using Random Finite-Element Meshes

Fracture and fragmentation are extremely nonlinear multiscale processes in which microscale damage mechanisms emerge at the macroscale as new fracture surfaces. Numerous numerical methods have been developed for simulating fracture initiation, propagation, and coalescence. In this paper, we present a computational approach for modeling pervasive fracture in quasi-brittle materials based on random close-packed Voronoi tessellations. Each Voronoi cell is formulated as a polyhedral finite element containing an arbitrary number of vertices and faces. Fracture surfaces are allowed to nucleate only at the intercell faces. Cohesive softening tractions are applied to new fracture surfaces in order to model the energy dissipated during fracture growth. The randomly seeded Voronoi cells provide a regularized discrete random network for representing fracture surfaces. The potential crack paths within the random network are viewed as instances of realizable crack paths within the continuum material. Mesh convergence of fracture simulations is viewed in a weak, or distributional, sense. The explicit facet representation of fractures within this approach is advantageous for modeling contact on new fracture surfaces and fluid flow within the evolving fracture network. Finally, applications of interest include fracture and fragmentation in quasi-brittle materials and geomechanical applications such as hydraulic fracturing, engineered geothermal systems, compressed-air energy storage, and carbon sequestration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Pressure Barrier Applicability to Polymer Flood Design

Effective oil displacement from a reservoir requires adequate and properly directed pressure gradients in areas of high oil saturation. If the polymer bank is too large or too viscous during a polymer flood, the pressure drop from the injection well to the polymer front may act as a pressure barrier by usurping most of the downstream driving force for oil displacement. Polymer injection pressures must be limited. The maximum allowable injection pressure is commonly constrained by caprock integrity, injection equipment, and/or regulations, even though fractures can be beneficial to polymer injectivity (and even sweep efficiency in some cases). This paper examines when the pressure-barrier concept limits the size and viscosity of the polymer bank during a polymer flood. Both analytical and numerical methods are used to address this issue. We examine the relevance of the pressure barrier concept for a wide variety of circumstances, including oil viscosities ranging from 10-cp (like at Daqing, China) to 1650-cp (like at Pelican Lake, Alberta), vertical wells (like at Tambaredjo, Suriname) versus horizontal wells (like at Milne Point, Alaska), single versus multiple layered reservoirs, permeability contrast, and with versus with crossflow between layers. We also examine the relation between the pressure-barrier concept and fractures and fracture extension during polymer injection. We demonstrate that in reservoirs with single layers, the pressure-barrier concept only limits the optimum viscosity of the injected polymer if the mobility of the polymer bank is less than the mobility of the displaced oil bank. The same is true for multi-zoned reservoirs with no crossflow between layers. Thus, for these cases, the optimum polymer viscosity is likely to be dictated by the mobility of the oil bank, unless other factors (like fracture extension) intervene. For multi-zoned reservoirs with free crossflow between layers, the situation is different. A compromise must be reached between injected polymer viscosity and the efficiency of oil recovery. The relevance of our findings is applied to operations for several existing polymer floods. This work is particularly relevant to viscous-oil reservoirs (like Pelican Lake and others) where the injected polymer viscosities are substantially lower than the oil viscosity

Wang, Dongmei↗
DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗
DOE OSTI2022

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Generative unfolding with distribution mapping

Machine learning enables unbinned, highly-differential cross section measurements. A recent idea uses generative models to morph a starting simulation into the unfolded data. We show how to extend two morphing techniques, Schrödinger Bridges and Direct Diffusion, in order to ensure that the models learn the correct conditional probabilities. This brings distribution mapping (DM) to a similar level of accuracy as the state-of-the-art conditional generative unfolding methods. Numerical results are presented with a standard benchmark dataset of single jet substructure as well as for a new dataset describing a 22-dimensional phase space of Z+2 -jets.

Butter, Anja↗
DOE OSTI2020

Toward a 3d Ising model with a weakly-coupled string theory dual

It has long been expected that the 3d Ising model can be thought of as a string theory, where one interprets the domain walls that separate up spins from down spins as two-dimensional string worldsheets. The usual Ising Hamiltonian measures the area of these domain walls. This theory has string coupling of unit magnitude. We add new local terms to the Ising Hamiltonian that further weight each spin configuration by a factor depending on the genus of the corresponding domain wall, resulting in a new 3d Ising model that has a tunable bare string coupling g_s g s . We use a combination of analytical and numerical methods to analyze the phase structure of this model as g_s g s is varied. We study statistical properties of the topology of worldsheets and discuss the prospects of using this new deformation at weak string coupling to find a worldsheet description of the 3d Ising transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Reduced-order model based feedback control of the Modified Hasegawa-Wakatani model

In this work, the development of model-based feedback control that stabilizes an unstable equilibrium is obtained for the Modified Hasegawa-Wakatani (MHW) equations, a classic model in plasma turbulence. First, a balanced truncation (a model reduction technique that has proven successful in flow control design problems) is applied to obtain a low dimensional model of the linearized MHW equation. Then a model based feedback controller is designed for the reduced order model using linear quadratic regulators (LQR). Finally, a linear quadratic gaussian (LQG) controller, which is more resistant to disturbances is deduced. The controller is applied on the non-reduced, nonlinear MHW equations to stabilize the equilibrium and suppress the transition to drift-wave induced turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Lattice Quantum Chromodynamics: SPI, mapping, site ordering, and QPX in Lattice QCD code on Mira (ALCF-2 Early Science Program Technical Report)

Lattice QCD is a numerical method to simulate QCD (Quantum Chromodynamics) including non-perturbative effects. Among other methods, lattice QCD is the only successful non-perturbative method that can be systematically improved from first principles. Lattice QCD plays an important role in High Energy Particle Physics (flavor physics, spectroscopy, and beyond the Standard Model physics) and Nuclear Physics (nucleon/nuclear spectrum & structure, quark-gluon plasma, and the QCD equation of state). Moreover, recently there are active developments in applying lattice gauge theory to quantum field theory in general, for example the theoretical description of graphene or cold atom systems. We define the QCD action on a finite and discrete ‘lattice’ like grid, so that we can compute the path integral (partition function) on a computer. A current lattice size is 96 3 ×192, and the corresponding number of degrees of freedom for the path integral is about 6×10 8 . Thus, it requires enormous computing resources. In this report, we describe our overall efforts to improve lattice QCD software on Mira, which include a new communication library using the Message Unit (MU) SPI (System Programming Interface), better node mapping strategy, lattice site reordering, and using the quad floating point (QPX) intrinsics.

Na, Heechang↗
DOE OSTI2013

Implementation of a Thermomechanical Model in Diablo for the Simulation of Selective Laser Melting

Modeling of selective laser melting requires consideration of both heat transfer and solid mechanics. The present work describes continuum modeling of SLM as envisioned for eventually supporting part-scale modeling of this fabrication process to determine end-state information such as residual stress and distortion. The determination of the evolving temperatures is dependent on the state of the material (powder or solid), the specified heating, the value of each of the constitutive parameters, and the configuration. Similarly, the current mechanical configuration is dependent on the temperatures, the powder-solid state, and the values of the constitutive parameters. A multi-mechanics formulation is required to properly describe such problems. The current report describes the problem formulation, numerical method, and constitutive parameters necessary to solve such a problem. Additionally, various verification and example problems are presented to illustrate the Diablo implementation.

42 ENGINEERING↗
DOE OSTI2013

Hierarchical Petascale Simulation Framework For Stress Corrosion Cracking

A number of major accomplishments resulted from the project. These include: • Data Structures, Algorithms, and Numerical Methods for Reactive Molecular Dynamics. We have developed a range of novel data structures, algorithms, and solvers (amortized ILU, Spike) for use with ReaxFF and charge equilibration. • Parallel Formulations of ReactiveMD (Purdue ReactiveMolecular Dynamics Package, PuReMD, PuReMD-GPU, and PG-PuReMD) for Messaging, GPU, and GPU Cluster Platforms. We have developed efficient serial, parallel (MPI), GPU (Cuda), and GPU Cluster (MPI/Cuda) implementations. Our implementations have been demonstrated to be significantly better than the state of the art, both in terms of performance and scalability. • Comprehensive Validation in the Context of Diverse Applications. We have demonstrated the use of our software in diverse systems, including silica-water, silicon-germanium nanorods, and as part of other projects, extended it to applications ranging from explosives (RDX) to lipid bilayers (biomembranes under oxidative stress). • Open Source Software Packages for Reactive Molecular Dynamics. All versions of our soft- ware have been released over the public domain. There are over 100 major research groups worldwide using our software. • Implementation into the Department of Energy LAMMPS Software Package. We have also integrated our software into the Department of Energy LAMMPS software package.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Localized Scale Coupling and New Educational Paradigms in Multiscale Mathematics and Science

We have experimentally demonstrated how microscale phenomena affect suspended particle behavior on the mesoscale, and how particle group behavior on the mesoscale influences the macroscale suspension behavior. Semi-analytical and numerical methods to treat flows on different scales have been developed, and a framework to combine these scale-dependent treatment has been described.

42 ENGINEERING↗
DOE OSTI2013

Modeling the Fracture of Ice Sheets on Parallel Computers

The objective of this project was to investigate the complex fracture of ice and understand its role within larger ice sheet simulations and global climate change. This objective was achieved by developing novel physics based models for ice, novel numerical tools to enable the modeling of the physics and by collaboration with the ice community experts. At the present time, ice fracture is not explicitly considered within ice sheet models due in part to large computational costs associated with the accurate modeling of this complex phenomena. However, fracture not only plays an extremely important role in regional behavior but also influences ice dynamics over much larger zones in ways that are currently not well understood. To this end, our research findings through this project offers significant advancement to the field and closes a large gap of knowledge in understanding and modeling the fracture of ice sheets in the polar regions. Thus, we believe that our objective has been achieved and our research accomplishments are significant. This is corroborated through a set of published papers, posters and presentations at technical conferences in the field. In particular significant progress has been made in the mechanics of ice, fracture of ice sheets and ice shelves in polar regions and sophisticated numerical methods that enable the solution of the physics in an efficient way.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Seismic Modeling Workshop: National Laboratory Overview/Discussion [Slides]

“full” waveform modeling refers to: • Use of numerical methods to accurately model the wave propagation in heterogeneous models. • Exploitation of the “full” waveform information or as much as it is physically possible • Iterative improvement of the tomographic images taking into account nonlinearity of the inversion.

58 GEOSCIENCES↗