Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational physics”

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 253 records · Page 14

An Introduction of Numerical Methods for Unconstrained Nonlinear Equations

Many applications in computational physics are nonlinear problems. These nonlinear problems include quantitative object reconstruction problems, inverse transport problems, nuclear weapon design problems, and particle tracking on unstructured meshes. The numerical algorithms used to solve nonlinear equations should be robust, reliable, and efficient. Numerical methods for unconstrained nonlinear equations are introduced in the slides.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Correlation of Injury Simulation with Clinical Assessment of Traumatic Brain Injury

This report contains a summary of our efforts to correlate head injury simulations predicting intracranial fluid cavitation with clinical assessments of brain injury from blunt impact to the head. Magnetic resonance imaging (MRI) data, collected on traumatic brain injury (TBI) subjects by researchers at the MIND Institute of New Mexico, was acquired for the current work. Specific blunt impact TBI case histories were selected from the TBI data for further study and possible correlation with simulation. Both group and single-subject case histories were examined. We found one single-subject case that was particularly suited for correlation with simulation. Diffusion tensor image (DTI) analysis of the TBI subject identified white matter regions within the brain displaying reductions in fractional anisotropy (FA), an indicator of local damage to the white matter axonal structures. Analysis of functional magnetic resonance image (fMRI) data collected on this individual identified localized regions of the brain displaying hypoactivity, another indicator of brain injury. We conducted high fidelity simulations of head impact experienced by the TBI subject using the Sandia head-neck-torso model and the shock physics computer code CTH. Intracranial fluid cavitation predictions were compared with maps of DTI fractional anisotropy and fMRI hypoactivity to assess whether a possible correlation exists. The ultimate goal of this work is to assess whether one can correlate simulation predictions of intracranial fluid cavitation with the brain injured sites identified by the fMRI and DTI analyses. The outcome of this effort is described in this report.

59 BASIC BIOLOGICAL SCIENCES↗

Remote Sensing of Dynamic Events [Slides]

Research into using advanced methods for remote sensing applications. A particular interest is in remote sensing of dynamic events involving high explosives (HE) and estimating parameters (e.g., temperature) that characterize the fireball. However, the research is interested in machine learning (ML) for infrared remote sensing in general. The motivation is to improve remote sensing methods of dynamic events to help with validation of computational physics codes.

47 OTHER INSTRUMENTATION↗

On the Development of Godunov Methods and Their Application to Gasdynamics

We are interested in developing the most accurate numeral models to simulate Gadynamics of explosions. But we must first recognize that there is no unique way to discretize the partial differential equations (PDEs) of gas dynamics, and also, there is no unique method to integrate the discretized forms of these equations. This has spawned a variety of numerical methods (religions*) that attempt to solve this conundrum, for example: Artificial Viscosity method created by the American mathematician John von Neumann, Flux-Corrected-Transport (FCT) method invented by Dr. Jay Boris, founder of the Laboratory for Computational Physics (LCP) at NRL, WENO method developed by Professor Stan Osher of UCLA and Professor Bjorn Enquist of U. Texas, and Godunov method by Academician Sergei Konstantinovich Godunov, Professor of the Sobolev Institute of Mathematics, RAS, in Novosibirsk to name just a few.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-Resolution Characterization of the Coupling Effects of Molten Salts, High Temperature and Irradiation on Intergranular Fracture

This project focused on providing a fundamental physico-chemical understanding of the coupling mechanisms of corrosion- and radiation-induced degradation at material-salt interfaces in Ni-based alloys operating in emulated Molten Salt Reactor(MSR) environments through the use of a unique suite of aging experiments, in-situ nanoscale characterization experiments on these materials, and multi-physics computational models. The technical basis and capabilities described in this report bring us a step closer to accelerate the deployment of MSRs by closing knowledge gaps related to materials degradation in harsh environments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Correlation of Injury Simulation with Clinical Assessment of Traumatic Brain Injury

This report contains a summary of our efforts to correlate head injury simulations predicting intracranial fluid cavitation with clinical assessments of brain injury from blunt impact to the head. Magnetic resonance imaging (MRI) data, collected on traumatic brain injury (TBI) subjects by researchers at the MIND Institute of New Mexico, was acquired for the current work. Specific blunt impact TBI case histories were selected from the TBI data for further study and possible correlation with simulation. Both group and single-subject case histories were examined. We found one single-subject case that was particularly suited for correlation with simulation. Diffusion tensor image (DTI) analysis of the TBI subject identified white matter regions within the brain displaying reductions in fractional anisotropy (FA), an indicator of local damage to the white matter axonal structures. Analysis of functional magnetic resonance image (fMRI) data collected on this individual identified localized regions of the brain displaying hypoactivity, another indicator of brain injury. We conducted high fidelity simulations of head impact experienced by the TBI subject using the Sandia head-neck-torso model and the shock physics computer code CTH. Intracranial fluid cavitation predictions were compared with maps of DTI fractional anisotropy and fMRI hypoactivity to assess whether a possible correlation exists. The ultimate goal of this work is to assess whether one can correlate simulation predictions of intracranial fluid cavitation with the brain injured sites identified by the fMRI and DTI analyses. The outcome of this effort is described in this report.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of Transient Coupled Chemical-Thermal-Fluid Multiphysics Simulation for Unsealed, Vented Aluminum-clad Spent Nuclear Fuel Storage Canister

As the second step towards modeling three-dimensional (3D) multi-physics computational fluid dynamics (CFD) of the unsealed containers in the IFSF facility, the previous methodology used for developing the sealed DOE standard canister model was utilized here. This 3D CFD model considers unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a larger network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than a helium backfilled canister. Similar to the sealed system, the set of full radiolytic chemical reactions of bulk moist air were established, and then reduced to a more manageable set of equations for the computational efficiency of the 3D CFD simulation. Due to the similar size of the unsealed canisters to the sealed DOE canisters, it is assumed that the chemical and thermal response times to changes in the conditions remains similar such that the same time-stepping algorithm can be utilized. This report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to INL INTEC CPP-603 over a long period of time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Transient Coupled Chemical-Thermal-Fluid Field Simulation for Sealed Aluminum-clad Spent Nuclear Fuel Storage Canister

As the first step toward developing three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model for unsealed and vented canister storage system, a 3D CFD model coupled with bulk gas radiolysis reactions was developed first for sealed DOE standard canisters filled with inert gas and trace amount of air and water. The workflow for constructing canister-scale 3D CFD models and coupling with gas phase radiolysis reactions were established, which can be readily extended to unsealed, vented canister storage system. This interim milestone report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields and hydrogen gas concentrations within sealed DOE standard canisters over long period of time. In addition, a mesh refinement test was done to show that increasing the models mesh refinement had negligible impact upon the temperature profiles.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling of Reactor Design and Optimization for Scale-Up of the Catalyxx Process for Ethanol Conversion to Higher Alcohol Biofuels

This report summarizes the results of a collaborative efforts between Oak Ridge National Laboratory (ORNL) and Catalyxx Inc. to investigate scale-up of Catalyxx’s Ethanol upgrading to higher alcohols process. The study was funded by the U.S. Department of Energy (DOE) Bioenergy Technologies Office (BETO) under CRADA (Cooperative Research and Development Agreement) No: NFE-20-08396. The project is part of the Direct Funding Opportunity (DFO) for Computational Science to Enable Bioenergy program which utilized computational toolsets developed by the Consortium for Computational Physics and Chemistry, a multi-laboratory consortium in BETO. The report here summarizes a packed-bed reactor modeling effort spanning the range from lab to industrial scale (from 4 gram to 5-ton catalyst beds), and examining reactor design, process optimization strategies, and suggested design and operating conditions for Catalyxx’s ethanol upgrading plants. The results in this report have been shared in monthly steering meetings and presentations are available in the shared data house owned by Catalyxx. The modeling effort helped to define optimum operation conditions for maximum alcohol selectivity and yield: i.e., temperature control scenarios ranging from adiabatic to isothermal, feed rate, pressure, and inlet H 2 /Ethanol ratio. The modeling results were verified at lab-(4 gram) and pre-pilot (4 kg) scales and has been used to evaluate a 5-ton packed-bed reactor and identify operating conditions to maximize the butanol yield. Special focus was given to understanding mass-transfer effects in the pre-pilot and pilot-scale reactors, over the domain of flow rate, pressure, feed composition, pellet size, shape, porosity, bed voidage, and reactor dimensions (i.e., length/diameter). Modeling was also used to evaluate innovative reactor design concepts such as water removal to improve alcohol selectivity and yield, and a reactor with an additional side inlet to facilitate quenching. These concepts were thoroughly explored, and potential benefits were disclosed. The results in this report are summarized and described qualitatively to protect the IP rights of Catalyxx. The details have been shared with the Catalyxx team in the regular steering meetings. At the end of the project, Catalyxx Inc. announced a successful demonstration of pilot scale operation in Seville, Spain.

02 PETROLEUM↗

DFO Computational Modeling Project - Catalytic Upgrading of Bio-based Furfural to 1,5-Pentanediol: A New Renewable Monomer for the Coatings Industry (Final Report)

This report summarizes the results of a collaborative effort between Oak Ridge National Laboratory (ORNL) and Pyran™ Inc. to utilize modeling capabilities developed by the CCPC (Consortium for Computational Physics and Chemistry) to assist Pyran in scaling up its proprietary process for thermocatalytic conversion of furfural to 1,5 pentanediol (PDO). Pyran’s bio-based PDO is a direct replacement for petroleum-based PDO and 1,6 hexanediol (HDO) currently used in polymers, additives, coatings, adhesives and sealants. The current market is in excess of $\$$1B/yr. According to Pyran, their production process results in 95+% reduction in fossil CO 2 at lower production costs compared to the existing petroleum-based routes for PDO & HDO. Other chemicals based on intermediates from this process have current markets in excess of $\$$10B/yr.

36 MATERIALS SCIENCE↗

Catalytic Upgrading of Bio-based Furfural to 1,5-Pentanediol: A New Renewable Monomer for the Coatings Industry

This report summarizes the results of a collaborative effort between Oak Ridge National Laboratory (ORNL) and Pyran™ Inc. to utilize modeling capabilities developed by the CCPC (Consortium for Computational Physics and Chemistry) to assist Pyran in scaling up its proprietary process for thermocatalytic conversion of furfural to 1,5 pentanediol (PDO). Pyran’s bio-based PDO is a direct replacement for petroleum-based PDO and 1,6 hexanediol (HDO) currently used in polymers, additives, coatings, adhesives and sealants. The current market is in excess of $1B/yr. According to Pyran, their production process results in 95+% reduction in fossil CO 2 at lower production costs compared to the existing petroleum-based routes for PDO & HDO. Other chemicals based on intermediates from this process have current markets in excess of $10B/yr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Test and Validate Distributed Coaxial Cable Sensors for in situ Condition Monitoring of Coal-Fired Boiler Tubes

This project aims to test, validate, and advance the technology readiness level (from TRL5 to TRL7) of a novel low-cost distributed stainless-steel/ceramic coaxial cable sensing (SSC-CCS) technology for in situ monitoring of the boiler tube temperature in existing coal-fired power plants. The novel SSC-CCS sensing technology and associated condition-based monitoring (CBM) software to be demonstrated in this project will lead to an improved understanding of the boiler tube failure mechanisms and a prognostic system to improve the overall performance, reliability, and flexibility of the nation’s coal-fired power plant fleet. A boiler tube monitoring system with distributed coaxial cable temperature sensors and a sensor acquisition system was constructed. The high-temperature coaxial cable sensor with a length of 1.3m was made by using a quartz tube (1mm inner diameter (ID) and 6mm outer diameter (OD)) to concentrically separate a 304 stainless-steel (SS) rod (1mm OD) and SS tube (7.94mm OD and 6.16mm ID). The sensor acquisition system includes a vector network analyzer (VNA), a radio frequency (RF) power amplifier, multiple switches and a USB hub. The distributed stainless-steel quartz coaxial cable sensor (SSQ-CCS) had a linear response to temperature with a resolution uncertainty of σ = 0.77℃. To withstand the harsh conditions of 3,300 steam pressures and 800℃ high temperatures, the sensor was shielded by a protective tube made of the same material as the boiler tube. The protection tube had an OD of 1.5 inches and a thickness of 0.25 inches. In the laboratory tests, the sensor showed good sensitivity and fast response. The drift was bounded between +0.33% and -0.67% during a test at 600℃ for 350 hours, indicating good stability of the sensor. A field test was conducted where four sensors were welded on four superheat tubes (SH-Ts) at a coal-fired power station over 400 days. Conventional thermocouples were welded to the superheater tubes alongside the coaxial cable sensors for the purpose of comparison. Two sensors were capable of distributed sensing, with three multiplexed sensing sections. The other two sensors were single section. During the 400-day test period, the power plant experienced startups and shutdowns. At the steady state operations, the temperature of the boiler tube is about 600℃ (1112°F). The sensors recorded the entire coal-firing processes (start-up, steady state, and shut-down) and the glitch event. A GSM modem and a Watchdog were added to the system to ensure reliable data recording. The GSM modem sent daily messages to plant managers and Clemson team to inform the status of the sensor system. If the system was not normally working, the Watchdog would reboot the system automatically. The new coaxial cable based distributed sensing technology has been proven to be successful in both laboratory and field tests. A comprehensive four-stage multi-physics computational framework has been developed to assist the design, optimization, installation, and operation of SSQ-CCS. With the consideration of various operation conditions, we predict the distributions of flue gas temperatures within coal-fired boilers, the temperature correlation between the boiler tube and SSQ-CCS, and the safety of SSQ-CCS. A conditional-based monitoring system is implemented as well. The computational framework developed in this work can guide the future operation of coal-fired plants and other power plants for the safety prediction of boiler operations.

01 COAL, LIGNITE, AND PEAT↗

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗

Applications of numerical codes to space plasma problems

Solar wind, earth's bowshock, and magnetospheric convection and substorms were investigated. Topics discussed include computational physics, multifluid codes, ionospheric irregularities, and modeling laser plasmas.

Northrop, T. G.↗

Time-dependent computational studies of flames in microgravity

The research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of the NASA Microgravity Science and Applications Program is described. The primary focus was on investigating fundamental questions concerning the propagation and extinction of premixed flames in Earth gravity and in microgravity environments. The approach was to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. A combination of one-dimensional and two-dimensional simulations was used to investigate the effects of curvature and dilution on ignition and propagation of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the mechanism leading to cellular instability, and to study the effects of gravity on the transition to cellular structure. A flame in a microgravity environment can be extinguished without external losses, and the mechanism leading to cellular structure is not preferential diffusion but a thermo-diffusive instability. The simulations have also lead to a better understanding of the interactions between buoyancy forces and the processes leading to thermo-diffusive instability.

Oran, Elaine S.↗

Laboratory studies, analysis, and interpretation of the spectra of hydrocarbons present in planetary atmospheres including cyanoacetylene, acetylene, propane, and ethane

Combining broadband Fourier transform spectrometers (FTS) from the McMath facility at NSO and from NRC in Ottawa and narrow band TDL data from the laboratories with computational physics techniques has produced a broad range of results for the study of planetary atmospheres. Motivation for the effort flows from the Voyager/IRIS observations and the needs of Voyager analysis for laboratory results. In addition, anticipation of the Cassini mission adds incentive to pursue studies of observed and potentially observable constituents of planetary atmospheres. Current studies include cyanoacetylene, acetylene, propane, and ethane. Particular attention is devoted to cyanoacetylen (H3CN) which is observed in the atmosphere of Titan. The results of a high resolution infrared laboratory study of the line positions of the 663, 449, and 22.5/cm fundamental bands are presented. Line position, reproducible to better than 5 MHz for the first two bands, are available for infrared astrophysical searches. Intensity and broadening studies are in progress. Acetylene is a nearly ubiquitous atmospheric constituent of the outer planets and Titan due to the nature of methane photochemistry. Results of ambient temperature absolute intensity measurements are presented for the fundamental and two two-quantum hotband in the 730/cm region. Low temperature hotband intensity and linewidth measurements are planned.

Blass, William E.↗