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,207 records · Page 67

BIGHORN Computational Fluid Dynamics Theory, Methodology, and Code Verification & Validation Benchmark Problems

This document presents the theoretical background for a hybrid finite-element / finite-volume fluid flow solver, namely BIGHORN, based on the Multiphysics Object Oriented Simulation Environment (MOOSE) computational framework developed at the Idaho National Laboratory (INL). An overview of the numerical methods used in BIGHORN are discussed and followed by a presentation of the formulation details. The document begins with the governing equations for the compressible fluid flow, with an outline of the requisite constitutive relations. A second-order finite volume method used for solving the compressible fluid flow problems is presented next. A Pressure-Corrected Implicit Continuous-fluid Eulerian (PCICE) formulation for time integration is also presented. The multi-fluid formulation is being developed. Although multi-fluid is not fully-developed, BIGHORN has been designed to handle multi-fluid problems. Due to the flexibility in the underlying MOOSE framework, BIGHORN is quite extensible, and can accommodate both multi-species and multi-phase formulations. This document also presents a suite of verification & validation benchmark test problems for BIGHORN. The intent for this suite of problems is to provide baseline comparison data that demonstrates the performance of the BIGHORN solution methods on problems that vary in complexity from laminar to turbulent flows. Wherever possible, some form of solution verification has been attempted to identify sensitivities in the solution methods, and suggest best practices when using BIGHORN.

97 MATHEMATICS AND COMPUTING↗

NEAMS-IPL MOOSE Framework Activities

The Multiapp Picard iteration Milestone’s purpose was to support a framework level “tight-coupling” method within the hierarchical Multiapp’s execution scheme. This new solution scheme gives developers new choices for running multiphysics applications, particularly those with very strong nonlinear effects or those requiring coupling across disparate time or spatial scales. Figure 1 shows a typical Multiapp setup in MOOSE. Each node represents a separate simulation containing a separate equation system. MOOSE solves the equation system on each node in turn, in a user-controlled manner. Information can be aggregated or split and transferred from parent to child or child to parent as needed between solves. Performing a tightly coupled execution scheme using this method wasn’t possible in the original implementation. This is was due to the inability to back up to a previous state once a converged solution was accepted at a particular Multiapp level.

97 MATHEMATICS AND COMPUTING↗

SIERRA Code Coupling Module: Arpeggio User Manual Version 4.44

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Transforming Ordinary Buildings into Smart Buildings via Low-Cost, Self-Powering Wireless Sensors & Sensor Networks

The research objective of this project is to design and demonstrate a low-cost, compact, easy-to-deploy, maintenance-free sensor node technology, and a network of such sensors, which enable the monitoring of multiphysical parameters and can transform today’s ordinary buildings into smart buildings with environmental awareness. We develop the sensor node and network via engineering and integration of existing technologies, including high-efficiency mechanical energy harvesting, and ultralow-power integrated circuits (ICs) for sensing and wireless communication. Through integration and innovative power management via specifically designed low-power control circuits for wireless sensing applications, and tailoring energy-harvesting components to indoor applications, the target products will have smaller volume, higher efficiency, and much lower cost (in both manufacturing and maintenance) than the baseline technology. Our development and commercialization objective is to create prototypes for our target products under the CWRU-Intwine collaboration.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

MOOSE Framework Speed and Memory Optimizations

The Multiphysics Object Oriented Simulation Environment (MOOSE) framework has been developed to solve a variety of complex simulations ranging from groundwater flows to full reactor simulations. MOOSE excels at solving fully-implicit, fully-coupled, highly nonlinear models using the Jacobian Free Newton-Krylov (JFNK) method relying on years of expertise in several dependent packages including PETSc, libMesh, Hypre and many others. The modular design of MOOSE enables the rapid development of new physics models that can be coupled with relative ease within the framework for building up simulations of different levels of complexity. With the ongoing adoption of MOOSE to new application areas along with streamlining of existing simulations, the priorities must be adjusted to include code optimizations. This includes both execution speed and the memory footprint of MOOSE-based simulations. These two categories can be further broken down into ''micro-optimizations'' or source code optimizations and algorithmic optimizations where source code is restructured to accomplish the same tasks through a more efficient means. However, before any code optimizations can occur, one must have robust way of measuring these critical simulation properties to identify the best candidates for optimizations along with having a way to identify the effectiveness of proposed optimizations to the framework. This report discusses the preliminary work towards the implementation of a robust execution baseline of the framework and derivative applications under a variety of different simulations. In addition, several potential optimization target areas in both the MOOSE framework and dependent libraries have been identified. An implementation strategy and potential benefits are discussed.

97 MATHEMATICS AND COMPUTING↗

Transient Heating and Thermomechanical Stress Modeling of Ceramic HEPA Filters

The purpose of this report is to showcase an initial finite-element analysis model of a ceramic High-Efficiency Particulate (HEPA) Air filter design. Next generation HEPA filter assemblies are being developed at LLNL to withstand high-temperature fire scenarios by use of ceramics and advanced materials. The filters are meant for use in radiological and nuclear facilities, and are required to survive 500°C fires over an hour duration. During such conditions, however, collecting data under varying parameters can be challenging; therefore, a Finite Element Analysis model of the filter was conducted using COMSOL ® Multiphysics to analyze the effects of fire. Finite Element Analysis (FEA) modelling offers several opportunities: researchers can quickly and easily consider impacts of potential design changes, material selection, and flow characterization on filter performance. Specifically, this model provides stress references for the sealant at high temperatures. Modeling of full filter assemblies was deemed inefficient given the computational requirements, so a section of three tubes from the assembly was modeled. The model looked at the transient heating and thermomechanical stress development during a 500°C air flow at 6 CFM. Significant stresses were found at the ceramic-metal interfaces of the filter, and conservative temperature profiles at locations of interest were plotted. The model can be used for the development of sealants that minimize stresses at the ceramic-metal interface. Further work on the model would include the full filter assembly and consider heat losses to make more accurate predictions.

36 MATERIALS SCIENCE↗

FY17 Status Report on NEAMS Neutronics Activities

Under the U.S. DOE NEAMS program, the high-fidelity neutronics code system has been developed to support the multiphysics modeling and simulation capability named SHARP. The neutronics code system includes the high-fidelity neutronics code PROTEUS, the cross section library and preprocessing tools, the multigroup cross section generation code MC2-3, the in-house meshing generation tool, the perturbation and sensitivity analysis code PERSENT, and post-processing tools. The main objectives of the NEAMS neutronics activities in FY17 are to continue development of an advanced nodal solver in PROTEUS for use in nuclear reactor design and analysis projects, implement a simplified sub-channel based thermal-hydraulic (T/H) capability into PROTEUS to efficiently compute the thermal feedback, improve the performance of PROTEUS-MOCEX using numerical acceleration and code optimization, improve the cross section generation tools including MC2-3, and continue to perform verification and validation tests for PROTEUS.

47 OTHER INSTRUMENTATION↗

An Experimental and Modeling Study of NO x -CO Formation in High Hydrogen Content (HHC) Fuels Combustion in Gas Turbine Applications

Species concentration measurements specifically those associated with NO x can act as important validation targets for developing kinetic models to predict NO x emissions under syngas as well as natural gas combustion accurately. In this collaborative research effort that included both experimental measurements, model development and simulations a comprehensive kinetic model and a multiphysics computational fluid dynamics platform has been developed and validated against the experimental data available in the literature as well as those acquired under this project. The experimental data provide the necessary NO x and speciation data for conditions relevant to gas turbine operations but are not readily available in the literature. The comprehensive chemical kinetic model consists of CO/H 2 /NO x oxidation with the full implementation of NO x evolution pathways, including thermal, prompt, N 2 O and NNH paths. The experiments conducted included NO x perturbed oxidation of natural gas at elevated pressure in laminar flow reactor and syngas/air combustion in a McKenna Burner – Flow Tube setup. A wide range of equivalence ratio, operating pressure as well as H 2 /CO ratio (for syngas only) was investigated. Temperature and NO x concentrations were measured in the flame and post-combustion zone. Experiments were also conducted for seeded syngas where trace hydrocarbon was introduced. The proposed model has been extensively tested. Predictions from the model are compared against multiple experimental datasets over a wide range of venues and operating conditions. The experimental venues include shock tube, plug flow reactor, and stirred reactor experiments that cover pressures from 1 to 100 bar and equivalence ratios from 0.5 to 1.5. In general, the overall model predictions are in good agreement with global combustion targets, such as ignition delay time, as well as with more detailed measurements from flow reactors and stirred reactors. Simulations are conducted for a wide range of reacting mixtures (H 2 /O 2 /N 2 , CO/H 2 /O 2 and CO/H 2 O/O 2 /N 2 ) with initial NO and NO 2 perturbations to consider exhaust gas recirculation (EGR) conditions.

08 HYDROGEN↗

Multiscale Modeling of Deformation, Fracture and Failure of Fusion Materials and Structures Final Report

The main objective of the research performed at UCLA under grant number DE-FG02-03ER54708 was to develop a physical understanding of the underlying mechanisms that control high-temperature deformation and fracture of fusion structural materials. The key aspect is the synergistic interaction between thermo-mechanical damage phenomena (e.g. creep, fatigue, and fracture), and irradiation damage phenomena (helium and hydrogen bubbles, and irradiation-induced defects). Another thrust area of the research was on the development of constitutive equations for the mechanical behavior of irradiated structures in fusion components (mainly the First Wall and Plasma-Facing components), and in utilizing them in descriptions in large-scale Finite Element Modeling of fusion structures. A parallel experimental effort on thermo-mechanical and plasma-induced damage is designed to validate and confront the developed theoretical and computational predictions. The models have been experimentally validated through experimental tests in simulated plasma transients on specially-designed W samples. This grant started at UCLA on November 14, 2003 and has had successive continuations till November 14, 2017. Research performed under the grant was to develop a range of hierarchical models for the post-elastic deformation, fracture and failure of fusion structural materials. The developed multiscale modeling approach has been based on rigorous mathematical, physical, and computational methods at the forefront of computational materials science. At the fundamental and smallest length scale (nm-μm), we developed advanced rate theory and Monte Carlo approaches to model microstructure evolution, non-equilibrium phase transformations, and dislocation-defect interactions. Microscopic and mesoscopic models of radiation hardening, ductile-to-brittle-transition, post-elastic deformation, plastic instabilities, and fracture processes are based on Dislocation Dynamics (DD) and Grain Growth Dynamics (GGD). At the continuum level, we focused on modeling fracture and failure mechanisms for Virtual Integrated Testing (VISTA), and for understanding the limits of both materials-by-design versus structural component design through advancing the new approach of multiscale-multiphysics. The overall research has been strongly coupled with the U.S. national experimental program to ascertain and verify the range of investigated phenomena. In the following, we give a summary of research progress supported by this grant since its inception till November 14, 2017. Detailed research results are available in our publications, listed in the publications section of this report.

36 MATERIALS SCIENCE↗

Institutional Computing: Final Report Quantum Effects on Cosmology: Probing Physics Beyond the Standard Model with Big Bang Nucleosynthesis

The current one-year project allocation (w17 burst) supports the continuation of research performed in the two-year Institutional Computing allocation (w14 bigbangnucleosynthesis). The project has supported development and production runs resulting in several publications[1, 2, 3, 4] in peer-review journals and talks. Most signi cantly, we have recently achieved a signi cant improvement in code performance. This improvement was essential to the prospect of making further progress on this heretofore unsolved multiphysics problem that lies at the intersection of nuclear and particle theory and the kinetic theory of energy transport in a system with internal (quantum) degrees of freedom.

79 ASTRONOMY AND ASTROPHYSICS↗

A Templated Approach for Multi-Physics Modeling of Hybrid Energy Systems in Modelica

A prototypical hybrid energy system (HES) couples a primary thermal power generator (i.e., nuclear power plant) with one or more additional subsystems beyond the traditional balance of plant electricity generation system. The definition and architecture of an HES can be adapted based on the needs and opportunities of a given local market. For example, locations in need of potable water may be best served by coupling a desalination plant to the HES. A location near an oil refinery may have a need for emission-free hydrogen production. The flexible, multidomain capabilities of Modelica are being used to investigate the dynamics (e.g., thermal hydraulics and electrical generation/consumption) of such a hybrid system. This paper examines the simulation infrastructure created to enable the coupling of multiphysics subsystem models for HES studies. A demonstration of a tightly coupled nuclear hybrid energy system implemented using the Modelica based infrastructure is presented for two representative cases. An appendix is also included providing a step-by-step procedure for using the template-based infrastructure.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Integral Full Core Multi-Physics PWR Benchmark with Measured Data

In recent years, the importance of modeling and simulation has been highlighted extensively in the DOE research portfolio with concrete examples in nuclear engineering with the CASL and NEAMS programs. These research efforts and similar efforts worldwide aim at the development of high-fidelity multi-physics analysis tools for the simulation of current and next-generation nuclear power reactors. Like all analysis tools, verification and validation is essential to guarantee proper functioning of the software and methods employed. The current approach relies mainly on the validation of single physic phenomena (e.g. critical experiment, flow loops, etc.) and there is a lack of relevant multiphysics benchmark measurements that are necessary to validate high-fidelity methods being developed today. This work introduces a new multi-cycle full-core Pressurized Water Reactor (PWR) depletion benchmark based on two operational cycles of a commercial nuclear power plant that provides a detailed description of fuel assemblies, burnable absorbers, in-core fission detectors, core loading and re-loading patterns. This benchmark enables analysts to develop extremely detailed reactor core models that can be used for testing and validation of coupled neutron transport, thermal-hydraulics, and fuel isotopic depletion. The benchmark also provides measured reactor data for Hot Zero Power (HZP) physics tests, boron letdown curves, and three-dimensional in-core flux maps from 58 instrumented assemblies. The benchmark description is now available online and has been used by many groups. However, much work remains to be done on the quantification of uncertainties and modeling sensitivities. This work aims to address these deficiencies and make this benchmark a true non-proprietary international benchmark for the validation of high-fidelity tools. This report details the BEAVRS uncertainty quantification for the first two cycle of operations and serves as the final report of the project.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Manufacturing Induced Residual Stress Development and Deformation in PMDI-10 Structural Foam

In this report, we investigate how manufacturing conditions result in the warpage of moderate density PMDI polyurethane foam (12-50 lb/ft 3 ) when they are released from a mold. We have developed a multiphysics modeling framework to simulate the manufacturing process including resin injection, foaming and mold filling, gelation of the matrix, elevated cure, vitrification, cool down, and demolding. We have implemented this framework within the Sierra Mechanics Finite Element Code Suite. We couple Aria for flow, energy conservation, and foaming/curing kinetics with Adagio for the nonlinear viscoelastic solid response in a multi-staged simulation process flow. We calibrate a model for the PMDI-10S (10 lb/ft 3 free rise foam) through a suite of characterization data presented here to calibrate the solid cure behavior of the foam. The model is then used and compared to a benchmark experiment, the manufacturing and warpage over 1 year of a 10 cm by 10 cm by 2.5 cm foam "staple". This component features both slender and thick regions that warp considerably differently over time. Qualitative agreement between the model and the experiment is achieved but quantitative accuracy is not.

36 MATERIALS SCIENCE↗

NEAMS-IPL MOOSE Midyear Framework Activities

Multiphysics Object Oriented Simulation Environement (MOOSE) is a modular, pluggable framework for building complex simulations. The ability to add new objects with custom syntax is a core capability that makes MOOSE a powerful platform for coupling multiple applications together within a single environment. Documentation is not only a required component for a complete software quality assurance plan, but is critical to the overall success and wide-spread adoption of the platform by other researchers and institutions. Therefore, the notion of a modular, pluggable framework has been extended to the creation of documentation content. This report details the methodology for defining the NQA-1 compliant requirements in a manner that allows for the automatic generation of the Software Requirement Specification (SRS), Software Design Description (SDD), and Requirement Traceabiltity Matrix (RTM) from any repository revision.

97 MATHEMATICS AND COMPUTING↗

Workshop on Software Productivity for Extreme-scale Science, January 13-14, 2014, Hilton Hotel, Rockville, MD

While emerging extreme-scale computing systems provide unprecedented resources for scientific discovery, DOE’s scientific communities now face a crisis in extreme-scale scientific software productivity due to a powerful confluence of factors: ongoing disruptive changes in computer architectures and new frontiers in modeling, simulation, and analysis of complex multiscale and multiphysics phenomena. Research in mission-critical science and energy problems has become increasingly dependent on high-performance computing, but productivity in extreme-scale scientific application development has not kept pace. The principal issues are (1) a significant lag between extreme-scale hardware and algorithmic innovations and their effective use in applications, leading to poorly scaling codes; (2) obstacles in the combined use of independently developed software components; (3) lack of agile yet rigorous software engineering practices for high-performance computing; and (4) failure to consider the entire lifecycle of large scientific software efforts, leading to fragile, complicated application codes that are increasingly difficult to enhance.

97 MATHEMATICS AND COMPUTING↗

Modeling and Validation of Irradiation Damage in Ni-based Alloys for Long-Term LWR Applications

Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components, including Ni-Cr alloy components. Predictions of the LWR component integrity rely on the ductile-to-brittle transition temperature (DBTT) curve to determine the toughness of the alloy after a period of exposure to in-service reactor conditions. The macroscopic DBTT curve is governed by the microstructure evolution of the alloy, thus, a capability to predict the microstructure evolution in response to thermal ageing is being implemented into Grizzly. Within this report section, the implementation of the Avrami equation to estimate the phase fraction of the aged long range ordered precipitates is discussed. The initial evaluation of current crystal plasticity models within Grizzly to model these Ni-Cr alloys is also presented.

36 MATERIALS SCIENCE↗

Leveraging Intrinsic Principal Directions for Multifidelity Uncertainty Quantification

In this work we propose an approach for accelerating Uncertainty Quantification (UQ) analysis in the context of Multifidelity applications. In the presence of complex multiphysics applications, which often require a prohibitive computational cost for each evaluation, multifidelity UQ techniques try to accelerate the convergence of statistics by leveraging the in- formation collected from a larger number of a lower fidelity model realizations. However, at the-state-of-the-art, the performance of virtually all the multifidelity UQ techniques is related to the correlation between the high and low-fidelity models. In this work we proposed to design a multifidelity UQ framework based on the identification of independent important directions for each model. The main idea is that if the responses of each model can be represented in a common space, this latter can be shared to enhance the correlation when the samples are drawn with respect to it instead of the original variables. There are also two main additional advantages that follow from this approach. First, the models might be correlated even if their original parametrizations are chosen independently. Second, if the shared space between models has a lower dimensionality than the original spaces, the UQ analysis might benefit from a dimension reduction standpoint. In this work we designed this general framework and we also tested it on several test problems ranging from analytical functions for verification purpose, up to more challenging application problems as an aero-thermo-structural analysis and a scramjet flow analysis.

42 ENGINEERING↗

Report on Aging and Irradiation Theoretical Model Development for Alloy 690 in Grizzly

This report describes the current efforts to implement a model to predict the engineering scale behavior of irradiated and aged Ni 2 Cr model alloy and Alloy 690 into the Grizzly software. Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components. Engineering scale behavior of materials are strongly impacted by the underlying microstructure: the distribution and interaction of dislocations and defects within the individual grains dictates how the material will respond to loading. The evolution of these dislocations and defects determines how the engineering scale material will perform under different loading and environmental conditions. Mesoscale models can function as a bridge to connect microstructure evolution models to engineering models of the ductile-brittle transition temperature (DBTT) curve and fracture. The toughness of a material, which governs the engineering scale fracture behavior, is a function of both elastic deformation and plastic deformation. While the yield stress measurements capture the elastic contribution to a material’s toughness, information about the dislocation movement is required to understand the plastic component of toughness. Crystal plasticity models predict dislocation density evolution and thereby connect microstructure evolution to engineering scale behavior predictions. These models calculate the plastic strain within a metal as a function of dislocations and other crystal defects, including precipitates and interstitial loops, among others, by tracking dislocation movement and the interaction of these dislocations with the crystal defects. Within NiCr alloys, long-range ordered (LRO) precipitates form under aging conditions. These precipitates act as barriers to dislocation motion and are either sheared by dislocations or pin dislocations. Additional defects, such as stacking fault tetrahedra (SFTs) from irradiation damage, also impede dislocation motion in these alloys. A crystal plasticity mode for Ni 2 Cr model alloy and Alloy 690 is being currently implemented in the Grizzly code at INL. This report includes a description of the theoretical models incorporated into the Grizzly crystal plasticity code. These theoretical models include mechanisms for dislocation glade and twinning as well as dislocation-defect interaction relationship for both LRO precipitates and SFTs. We have selected these dislocation models based on the deformation mechanisms observed in experimental micropillar compression tests conducted on the Ni 2 Cr model alloy. The models are then extended to Alloy 690 with ongoing parameter fitting.

36 MATERIALS SCIENCE↗