Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphysics”

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

On the Interaction Between a Flow Instability and Inertial Particles in Supersonic Jet Plume-Surface Interaction

Extraterrestrial landings often require spacecraft to deploy supersonic jets directed at dusty planetary surfaces to ensure controlled descents. This interaction, termed plume-surface interaction (PSI), critically influences crater formation and the dynamics of ejected particles, which can obscure vision, damage onboard sensors, and impact nearby infrastructure with significant particle velocities. The physics of PSI remain poorly understood due to the complexities in simulating the involved multiphysics and experimentally replicating reduced atmospheric pressure conditions. This article discusses a collaborative effort to investigate PSI in reduced atmospheric pressure environments using a Mach 5.3 jet impinging on a granular bed. The results highlight observations of an azimuthal pattern in the ejecta distribution, characterized by alternating ejecta streaks of varying particle concentrations. This phenomenon suggests the presence of a low Reynolds number instability, a notable discovery given the supersonic nature of the jet.

PSI↗
NASA NTRS2024

Particle Impact Simulation and Ignition Prediction

An experimentally calibrated tool is needed to predict if a system is susceptible to failure by particle impact ignition (PI) based on use conditions, materials, and flow geometry. This tool will accelerate new components, evaluating existing hardware, and help disposition anomalies. - Conduct particle impact testing with in-situ diagnostics and complementary simulations on subset of key engineering materials (IN718, M400, 316L, 6061, Ti64, Zr) to develop a proof-of-concept predictive tool for assessing the risk of PI for idealized geometries (spherical particles) in realistic environments. - Assess particle/target interactions (coefficient of restitution, ignition, kindling) using instrumented particle impact rigs while systematically varying key parameters (materials, particle size, environment, target configuration). - Determine key field variables (temperature, strain, stress) in particle impacts using Multiphysics finite element and hydrocode simulations validated through comparison with experimental measurements and observations. - Synthesize experiments and simulations into constitutive models for PI that can be integrated with existing computational fluid dynamics (CFD) and Debris Transport Analysis (DTA) tools in future efforts.

Jonathan Tylka↗
NASA NTRSDate not supplied

Investigations of an Aeroelastic Optimization Benchmark Problem with MPhys and FUN3D

High-fidelity aeroelastic optimization capabilities are becoming more common in recent years. In this field, researchers typically present results on the design problems of greatest interest to them. The lack of consistency with respect to baseline configurations, design variables, objectives and constraints, etc., makes it difficult to compare approaches and creates an obstacle for the community to work together to advance the state of the art. The High-Fidelity Aeroelastic Optimization Benchmark Working Group has been established to create and study a series of benchmark aeroelastic optimization problems that the community can utilize to compare methods on the same optimization problem. This paper presents optimization results from NASA Langley Research Center for the first benchmark case of this series with analysis based on MPhys, a multiphysics library for OpenMDAO, and FUN3D, an unstructured computational fluid dynamics suite of tools. For the optimization case presented, it is demonstrated that a progression of increasing aerodynamic model fidelity reduces the wall time required to complete the optimization.

CFD↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2015

Validation Results for Core-Scale Oil Shale Pyrolysis

This report summarizes a study of oil shale pyrolysis at various scales and the subsequent development a model for in situ production of oil from oil shale. Oil shale from the Mahogany zone of the Green River formation was used in all experiments. Pyrolysis experiments were conducted at four scales, powdered samples (100 mesh) and core samples of 0.75”, 1” and 2.5” diameters. The batch, semibatch and continuous flow pyrolysis experiments were designed to study the effect of temperature (300°C to 500°C), heating rate (1°C/min to 10°C/min), pressure (ambient and 500 psig) and size of the sample on product formation. Comprehensive analyses were performed on reactants and products - liquid, gas and spent shale. These experimental studies were designed to understand the relevant coupled phenomena (reaction kinetics, heat transfer, mass transfer, thermodynamics) at multiple scales. A model for oil shale pyrolysis was developed in the COMSOL multiphysics platform. A general kinetic model was integrated with important physical and chemical phenomena that occur during pyrolysis. The secondary reactions of coking and cracking in the product phase were addressed. The multiscale experimental data generated and the models developed provide an understanding of the simultaneous effects of chemical kinetics, and heat and mass transfer on oil quality and yield. The comprehensive data collected in this study will help advance the move to large-scale in situ oil production from the pyrolysis of oil shale.

Staten, Josh↗
DOE OSTI2015

CASTING DEFECT MODELING IN AN INTEGRATED COMPUTATIONAL MATERIALS ENGINEERING APPROACH

To accelerate the introduction of new cast alloys, the simultaneous modeling and simulation of multiphysical phenomena needs to be considered in the design and optimization of mechanical properties of cast components. The required models related to casting defects, such as microporosity and hot tears, are reviewed. Three aluminum alloys are considered A356, 356 and 319. The data on calculated solidification shrinkage is presented and its effects on microporosity levels discussed. Examples are given for predicting microporosity defects and microstructure distribution for a plate casting. Models to predict fatigue life and yield stress are briefly highlighted here for the sake of completion and to illustrate how the length scales of the microstructure features as well as porosity defects are taken into account for modeling the mechanical properties. Thus, the data on casting defects, including microstructure features, is crucial for evaluating the final performance-related properties of the component. ACKNOWLEDGEMENTS This work was performed under a Cooperative Research and Development Agreement (CRADA) with the Nemak Inc., and Chrysler Co. for the project "High Performance Cast Aluminum Alloys for Next Generation Passenger Vehicle Engines. The author would also like to thank Amit Shyam for reviewing the paper and Andres Rodriguez of Nemak Inc. Research sponsored by the U. S. Department of Energy, Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Office, as part of the Propulsion Materials Program under contract DE-AC05-00OR22725 with UT-Battelle, LLC. Part of this research was conducted through the Oak Ridge National Laboratory's High Temperature Materials Laboratory User Program, which is sponsored by the U. S. Department of Energy, Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Program.

Sabau, Adrian S↗
DOE OSTI2015

Advanced Nodal P 3 /SP 3 Axial Transport Solvers for the MPACT 2D/1D Scheme

As part of its initiative to provide multiphysics simulations of nuclear reactor cores, the Consortium for Advanced Simulation of Light Water Reactors (CASL) is developing the Virtual Environment for Reactor Applications Core Simulator (VERA-CS). The MPACT code, which is the primary neutron transport solver of VERA-CS, employs the two-dimensional/one-dimensional (2D/1D) method to solve 3-dimensional neutron transport problems and provide sub-pin-level resolution of the power distribution. While 2D method of characteristics is used to solve for the transport effects within each plane, 1D-nodal methods are used axially. There have been extensive studies of the 2D/1D method with a variety nodal methods, and the P 3 /SP 3 solver has proved to be an effective method of providing higher-fidelity solutions while maintaining a low computational burden.The current implementation in MPACT wraps a one-node nodal expansion method (NEM) kernel for each moment, iterating between them and performing multiple sweeps to resolve flux distributions. However, it has been observed that this approach is more sensitive to convergence problems. This paper documents the theory and application two new nodal P 3 /SP 3 approaches to be used within the 2D/1D method in MPACT. These two approaches aim to provide enhanced stability compared with the pre-existing one-node approach. Results from the HY-NEM-SP 3 solver show that the accuracy is consistent with the one-node formulations and provides improved convergence for some problems; but the solver has issues with cases in thin planes. Although the 2N-SENM-SP 3 solver is still under development, it is intended to resolve the issues with HY-NEM-SP 3 but it will incur some additional computational burden by necessitating an additional 1D-CMFD-P 3 solver to generate the second moment cell-averaged scalar flux.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Assessment of Pellet-Clad Interaction Indicators in Watts Bar Unit 1 using the VERA Framework

The Consortium for Advanced Simulation of Light Water Reactors (CASL) aims to provide high-fidelity, multiphysics simulations of light water reactors (LWRs) by coupling a variety of codes within the Virtual Environment for Reactor Analysis (VERA). One of the primary goals of CASL is to predict local cladding failure through pellet-clad interaction (PCI). This capability is currently being pursued through several different approaches, such as with Tiamat, which is a simulation tool within VERA that more tightly couples the MPACT neutron transport solver, the CTF thermal hydraulics solver, and the MOOSE-based Bison-CASL fuel performance code. However, the process in this paper focuses on running fuel performance calculations with Bison-CASL to predict PCI using the multicycle output data from coupled neutron transport/thermal hydraulics simulations. In recent work within CASL, Watts Bar Unit 1 has been simulated over 12 cycles using the VERA core simulator capability based on MPACT and CTF. Using the output from these simulations, Bison-CASL results can be obtained without rerunning all 12 cycles, while providing some insight into PCI indicators. Multi-cycle Bison-CASL results are presented and compared against results from the FRAPCON fuel performance code. There are several quantities of interest in considering PCI and subsequent fuel rod failures, such as the clad hoop stress and maximum centerline fuel temperature, particularly as a function of time. Bison-CASL performs single-rod simulations using representative power and temperature distributions, providing high-resolution results for these and a number of other quantities. This will assist in identifying fuels rods as potential failure locations for use in further analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Determination of Experimental Fuel Rod Parameters using 3D Modelling of PCMI with MPS Defect

An in-reactor experiment is being designed in order to validate the pellet-cladding mechanical interaction (PCMI) behavior of the BISON fuel performance code. The experimental parameters for the test rod being placed in the Halden Research Reactor are being determined using BISON simulations. The 3D model includes a missing pellet surface (MPS) defect to generate large local cladding deformations, which should be measureable after typical burnup times. The BISON fuel performance code is being developed at Idaho National Laboratory (INL) and is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. BISON supports both 2D and 3D finite elements and solves the fully coupled equations for solid mechanics, heat conduction and species diffusion. A number of fuel performance effects are included using models for swelling, densification, creep, relocation and fission gas production & release. In addition, the mechanical and thermal contact between the fuel and cladding is explicitly modelled using a master-slave based contact algorithm. In order to accurately predict PCMI effects, the BISON code includes the relevant physics involved and provides a scalable and robust solution procedure. The depth of the proposed MPS defect is being varied in the BISON model to establish an optimum value for the experiment. The experiment will be interrupted approximately every 6 months to measure cladding radial deformation and provide data to validate BISON. The complete rodlet (~20 discrete pellets) is being simulated using a 180° half symmetry 3D model with MPS defects at two axial locations. In addition, annular pellets will be used at the top and bottom of the pellet stack to allow thermocouples within the rod to measure the fuel centerline temperature. Simulation results will be presented to illustrate the expected PCMI behavior and support the chosen experimental design parameters.

BISON↗
DOE OSTI2016

Analysis of Experimental Fuel Rod Parameters using 3D Modelling of PCMI with MPS Defect

An in-reactor experiment is being designed in order to validate the pellet-cladding mechanical interaction (PCMI) behavior of the BISON fuel performance code. The experimental parameters for the test rod being placed in the Halden Research Reactor are being determined using BISON simulations. The 3D model includes a missing pellet surface (MPS) defect to generate large local cladding deformations, which should be measureable after typical burnup times. The BISON fuel performance code is being developed at Idaho National Laboratory (INL) and is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. BISON supports both 2D and 3D finite elements and solves the fully coupled equations for solid mechanics, heat conduction and species diffusion. A number of fuel performance effects are included using models for swelling, densification, creep, relocation and fission gas production & release. In addition, the mechanical and thermal contact between the fuel and cladding is explicitly modelled using a master-slave based contact algorithm. In order to accurately predict PCMI effects, the BISON code includes the relevant physics involved and provides a scalable and robust solution procedure. The depth of the proposed MPS defect is being varied in the BISON model to establish an optimum value for the experiment. The experiment will be interrupted approximately every 6 months to measure cladding radial deformation and provide data to validate BISON. The complete rodlet (~20 discrete pellets) is being simulated using a 180° half symmetry 3D model with MPS defects at two axial locations. In addition, annular pellets will be used at the top and bottom of the pellet stack to allow thermocouples within the rod to measure the fuel centerline temperature. Simulation results will be presented to illustrate the expected PCMI behavior and support the chosen experimental design parameters.

BISON↗
DOE OSTI2016

STEADY STATE MODELING OF THE MINIMUM CRITICAL CORE OF THE TRANSIENT REACTOR TEST FACILITY

With the advent of next generation reactor systems and new fuel designs, the U.S. Department of Energy (DOE) has identified the need for the resumption of transient testing of nuclear fuels. The DOE has decided that the Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) is best suited for future testing. TREAT is a thermal neutron spectrum, air-cooled, nuclear test facility that is designed to test nuclear fuels in transient scenarios. These specific scenarios range from simple temperature transients to full fuel melt accidents. DOE has expressed a desire to develop a simulation capability that will accurately model the experiments before they are irradiated at the facility. It is the aim for this capability to have an emphasis on effective and safe operation while minimizing experimental time and cost. The multi physics platform MOOSE has been selected as the framework for this project. The goals for this work are to identify the fundamental neutronics properties of TREAT and to develop an accurate steady state model for future multiphysics transient simulations. In order to minimize computational cost, the effect of spatial homogenization and angular discretization are investigated. It was found that significant anisotropy is present in TREAT assemblies and to capture this effect, explicit modeling of cooling channels and inter-element gaps is necessary. For this modeling scheme, single element calculations at 293 K gave power distributions with a root mean square difference of 0.076% from those of reference SERPENT calculations. The minimum critical core configuration with identical gap and channel treatment at 293 K resulted in a root mean square, total core, radial power distribution 2.423% different than those of reference SERPENT solutions.

Critical↗
DOE OSTI2016

Least-Squares PN Formulation of the Transport Equation Using Self-Adjoint-Angular-Flux Consistent Boundary Conditions.

In this paper, we study the Least-Squares (LS) PN form of the transport equation compatible with voids in the context of Continuous Finite Element Methods (CFEM).We first deriveweakly imposed boundary conditions which make the LS weak formulation equivalent to the Self-Adjoint Angular Flux (SAAF) variational formulation with a void treatment, in the particular case of constant cross-sections and a uniform mesh. We then implement this method in Rattlesnake with the Multiphysics Object Oriented Simulation Environment (MOOSE) framework using a spherical harmonics (PN) expansion to discretize in angle. We test our implementation using the Method of Manufactured Solutions (MMS) and find the expected convergence behavior both in angle and space. Lastly, we investigate the impact of the global non-conservation of LS by comparing the method with SAAF on a heterogeneous test problem.

Least-Squares↗
DOE OSTI2017

Final Project Report - Revised Version

An over-arching goal of this research is to calibrate geophysical techniques for hydrate exploration, evaluation, and production monitoring. Extensive field data of hydrate-bearing sediments exist, but quantitative estimates of the amount and distribution of hydrates are difficult to determine. Thus, the primary project objectives were to relate seismic and acoustic velocities and attenuations to hydrate saturation and texture. The project aimed to collect seismic properties along with other measurements (e.g., complex resistivity, micro-focus x-ray computed tomography, etc.). The multiphysics dataset would enable researchers to understand not only the interaction between mineral surfaces and gas hydrates, but also how the hydrate formation method affects the hydrate-sediment system in terms of elastic properties.

03 NATURAL GAS↗
DOE OSTI2015

MODELLING OF FUEL BEHAVIOUR DURING LOSS-OF-COOLANT ACCIDENTS USING THE BISON CODE

This work presents recent developments to extend the BISON code to enable fuel performance analysis during LOCAs. This newly developed capability accounts for the main physical phenomena involved, as well as the interactions among them and with the global fuel rod thermo-mechanical analysis. Specifically, new multiphysics models are incorporated in the code to describe (1) transient fission gas behaviour, (2) rapid steam-cladding oxidation, (3) Zircaloy solid-solid phase transition, (4) hydrogen generation and transport through the cladding, and (5) Zircaloy high-temperature non-linear mechanical behaviour and failure. Basic model characteristics are described, and a demonstration BISON analysis of a LWR fuel rod undergoing a LOCA accident is presented. Also, as a first step of validation, the code with the new capability is applied to the simulation of experiments investigating cladding behaviour under LOCA conditions. The comparison of the results with the available experimental data of cladding failure due to burst is presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Least-Squares PN Formulation of the Transport Equation Using Self-Adjoint-Angular-Flux Consistent Boundary Conditions

In this paper, we study the Least-Squares (LS) PN form of the transport equation compatible with voids [1] in the context of Continuous Finite Element Methods (CFEM).We first deriveweakly imposed boundary conditions which make the LS weak formulation equivalent to the Self-Adjoint Angular Flux (SAAF) variational formulation with a void treatment [2], in the particular case of constant cross-sections and a uniform mesh. We then implement this method in Rattlesnake with the Multiphysics Object Oriented Simulation Environment (MOOSE) framework [3] using a spherical harmonics (PN) expansion to discretize in angle. We test our implementation using the Method of Manufactured Solutions (MMS) and find the expected convergence behavior both in angle and space. Lastly, we investigate the impact of the global non-conservation of LS by comparing the method with SAAF on a heterogeneous test problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Self-consistent hybrid neoclassical-magnetohydrodynamic simulations of axisymmetric plasmas

Neoclassical effects (e.g., conductivity reduction and bootstrap currents) have a profound impact on many magnetohydrodynamic (MHD) instabilities in toroidally-confined plasmas, including tearing modes, edge-localized modes, and resistive wall modes. High-fidelity simulations of such phenomena require a multiphysics code that self-consistently couples the kinetic and fluid models. We review a hybrid formulation from the recent literature AB that is appropriate for such studies. In particular, the formulation uses a set of time-dependent drift-kinetic equations (DKEs) to advance the non-Maxwellian part of the electron and ion distribution functions ( fNM ) with linearized Fokker-Planck-Landau collision operators. The form of the DKEs used were derived in a Chapman-Enskog-like fashion, ensuring that fNM carries no density, momentum, or temperature. Rather, these quantities are contained within the background Maxwellian and are evolved by a set of MHD equations which are closed by moments of fNM . We then present two DKE solvers based upon this formulation in axisymmetric toroidal geometries. The Neoclassical Ion-Electron Solver (NIES) solves the steady-state DKEs in the low-collisionality limit. Convergence and benchmark studies are discussed, providing a proof-of-principle that this new formulation can accurately reproduce results from the literature in the limit considered. We then present the DK4D code which evolves the finite-collisionality DKEs time-dependently. Computational methods used and successful benchmarks to other neoclassical models and codes are discussed. Furthermore, we couple DK4D to a reduced, transport-timescale MHD code. The resulting hybrid code is used to simulate the evolution of the current density in a large-aspect-ratio plasma in the presence of several different time-dependent pressure profiles. These simulations demonstrate the self-consistent, dynamic formation of the ohmic and bootstrap currents. In the slowly-evolving plasmas considered, these first-principle simulations are shown to verify a simpler, steady-state neoclassical model that is commonly used in transport codes. Future work will involve coupling DK4D to a spatially three-dimensional, extended MHD code, generalizing DK4D to nonaxisymmetric geometries, and simulating more quickly-evolving and realistic plasmas. The ultimate goal of this work is to perform self-consistent, hybrid simulations of complex tokamak instabilities and calculations of neoclassical toroidal viscosity. A J. J. Ramos, Phys. Plasmas 17 , 082502 (2010). B J. J. Ramos, Phys. Plasmas 18 , 102506 (2011).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗