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 307 records · Page 17

Coupling Exascale Multiphysics Applications: Methods and Lessons Learned

With the growing computational complexity of science and the complexity of new and emerging hardware, it is time to re-evaluate the traditional monolithic design of computational codes. One new paradigm is constructing larger scientific computational experiments from the coupling of multiple individual scientific applications, each targeting their own physics, characteristic lengths, and/or scales. We present a framework constructed by leveraging capabilities such as in-memory communications, workflow scheduling on HPC resources, and continuous performance monitoring. This code coupling capability is demonstrated by a fusion science scenario, where differences between the plasma at the edges and at the core of a device have different physical descriptions. This infrastructure not only enables the coupling of the physics components, but it also connects in situ or online analysis, compression, and visualization that accelerate the time between a run and the analysis of the science content. Results from runs on Titan and Cori are presented as a demonstration.

Choi, Jong Youl↗

Microwave sintering of complex shapes: From multiphysics simulation to improvements of process scalability

The microwave sintering homogeneity of large and complex shape specimens is analyzed. A new approach enabling the fabrication of complex shapes ceramics via 3D printing and microwave sintering is presented. The use of a dental microwave cavity is shown to enable a substantial level of densification of complex shape components while restricting the grain growth. The homogeneity of the processed samples during microwave sintering is studied by an electromagnetic-thermal-mechanical simulation. The realistic densification behavior, that phenomenologically takes into account the microwave effect, is included in the modeling framework. We observe that the simulation indicates the sharp correlation between the microwave field distribution in the cavity, the temperature profile, and the specimen's shape distortion.

36 MATERIALS SCIENCE↗

Coupled Multiphysics Modeling of Lithium-Ion Batteries for Automotive Crashworthiness Applications

Considerable advances have been made in battery safety models, but achieving predictive accuracy across a wide range of conditions continues to be challenging. Interactions between dynamically evolving mechanical, electrical, and thermal state variables make model prediction difficult during mechanical abuse scenarios. In this study, we develop a physics-based modeling approach that allows for choosing between different mechanical and electrochemical models depending on the required level of analysis. We demonstrate the use of this approach to connect cell-level abuse response to electrode-level and particle-level transport phenomena. A pseudo-two-dimensional model and simplified single-particle models are calibrated to electrical-thermal cycling data and applied to mechanically induced short-circuit scenarios to understand how the choice of electrochemical model affects the model prediction under abuse scenarios. These models are implemented using user-defined subroutines on ls-dyna finite element software and can be coupled with existing automotive crash safety models.

analysis and design of components↗

Multiphysics modeling of mesh piezoelectric atomizers

An ultrasonic clothes dryer was developed by researchers at Oak Ridge National Laboratory based on a novel approach of using high-frequency mechanical vibration instead of heat to extract moisture as a cold mist. This technology is based on direct mechanical coupling between mesh piezoelectric (PZT) transducers and wet fabric. The vibration introduces sufficient momentum to the droplets trapped in the fabric pores to atomize them and leave the garment in a cold state. In the vibrating transducer, deformation followed by the effects of boundary layer acoustic streaming results in ejection of the atomized droplets. The research presented bridges the vibration of a PZT mesh transducer to the induced acoustic field and to capillary-wave theory. Mathematical modeling studies free and forced vibrations of a mesh-like PZT structure, using the structural parameters identified by actuation testing in several case studies. Computational fluid–structure interaction modeling is performed to couple the vibrations of a PZT transducer with an in-contact droplet. The results obtained are used to investigate (1) the transverse deformation of the vibrating mesh transducer in contact with a droplet, (2) the resultant boundary layer acoustic streaming in the fluid surrounding the vibrating surface, and (3) the droplet deformation and fluid ejection. The physics of atomization are linked to the level of the near-wall droplet vibrations induced by the surface deformation of the transducer. Then the surface deformation is linked to the properties of the PZT mesh transducer and input actuation frequency and power.

Dupuis, Eric↗

Crack Growth Rate at Oxygen Electrode/Electrolyte Interface in Solid Oxide Electrolysis Cells Predicted by Experiment Coupled Multiphysics Modeling

Solid oxide electrolysis cell (SOEC) is a very efficient hydrogen production technology, but the cell degradation is a serious limiting factor for its long-term implementation. Oxygen electrode (OE) delamination is reported to be the critical degradation mechanism. In this study, we present a methodology to understand the delamination failure of the OE due to chemical stress in a better perspective. Several OE configurations were tested: baseline strontium-doped lanthanum cobalt iron oxide (LSCF) single layer design and tantalum-doped strontium cobalt oxide (SCT) - LSCF bilayer designs with different SCT loadings. An electro-chemo-mechanical model is developed to associate the electrochemical behavior of the cell with solid mechanics for calculating crack growth of the cell during long term test. The bilayer configuration with SCT 20 wt% has better performance as it survived in the long-term life test with the least crack length. This study implies that an additional nano-coating of SCT over the OE have improved the species transport and oxygen evolution with reduced chemical stress. As the operating current density decreases, it takes longer time for the cell to reach the delamination with the same critical crack length of 6.5 μ m (∼93% of the electrode/electrolyte interface length). Finally, it was concluded that chemical stress plays a significant role in interface delamination failure, however it may not be the only source of stresses at the interface.

08 HYDROGEN↗

Multiphysics Engineering Analysis for an Integrated Design of ITER Diagnostic First Wall and Diagnostic Shield Module Design

ITER diagnostic first walls (DFWs) and diagnostic shield modules (DSMs) inside the port plugs (PPs) are designed to protect diagnostic instrument and components from a harsh plasma environment and provide structural support while allowing for diagnostic access to the plasma. The design of DFWs and DSMs are driven by 1) plasma radiation and nuclear heating during normal operation 2) electromagnetic loads during plasma events and associate component structural responses. A multi-physics engineering analysis protocol for the design has been established at Princeton Plasma Physics Laboratory and it was used for the design of ITER DFWs and DSMs. The analyses were performed to address challenging design issues based on resultant stresses and deflections of the DFW-DSM-PP assembly for the main load cases. ITER Structural Design Criteria for In-Vessel Components (SDC-IC) required for design by analysis and three major issues driving the mechanical design of ITER DFWs are discussed. The general guidelines for the DSM design have been established as a result of design parametric studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics FISPACT-II and TENDL-2019 simulation: neutron-induced damage metrics

Nuclear interactions can be the source of atomic displacement, irradiation-induced defects and transmutation in structural materials. Such quantities are derived from, or can be correlated to, nuclear kinematic simulations of primary atomic energy distributions spectra and the quantification of the numbers of secondary defects produced per primary as a function of the available recoils, residual and emitted, energies. Recoil kinematics of neutral, residual, gas, proton, alpha particle emissions are now more rigorously treated based on recent, complete and enhanced nuclear data parsed in state-of-the-art processing tools. Defect production metrics are the starting point in the complex problem of correlating and simulating the behaviour of materials under irradiation, as experimental information is rare or scattered. Detailed, segregated primary knock-on-atom metrics are now becoming available as the starting point of further simulation processes of isolated and clustered defects in material lattices. This allows more materials, lattices, neutron incident energy ranges, and irradiation conditions to be explored with sufficient data to adequately cover both standard and novel applications and materials: the broader reactor applications landscape. The damage metrics of of materials are systematically explored under typical but different reactor's type environment. The inventory code FISPACT-II combined with the enhanced nuclear data forms of the TENDL-2019 libraries allow one to not only calculate dpa from mostly scattering events but to also properly predict gas production, nuclear heating and transmutation under the same conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiphysics simulation of supercritical CO 2 gasification for hydrogen production

A bench-scale supercritical CO 2 gasifier is being considered by the U.S. Department of Energy’s National Energy Technology Laboratory (NETL). A reacting computational fluid dynamics (CFD) model was developed using ANSYS/Fluent in order to investigate the impact of various operating parameters on syngas composition. The model was validated using available data in the literature. Simulations were conducted over a range of coal-slurry loadings, reaction temperatures, system pressures, and oxygen mass flowrates with the goal of optimizing hydrogen production for the considered reactor. The results from the simulations identified a positive correlation between slurry loading and hydrogen production, while reactor temperature and pressure had a limited impact on H 2 production at the conditions of interest. Additionally, the simulations indicated that oxygen mass flowrate has less of an impact on hydrogen production at higher slurry loadings.

08 HYDROGEN↗

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗

Development of Adaptive Model Refinement (AMoR) for Multiphysics and Multifidelity Problems

This project investigated the development and utilization of Adaptive Model Refinement (AMoR) for nuclear systems simulation applications. AMoR refers to utilization of several models of physical phenomena which differ in prediction fidelity. If the highest fidelity model is judged to always provide or exceeded the desired fidelity, than if one can determine the difference in a Quantity of Interest (QoI) between the highest fidelity model and lower fidelity models, one could utilize the fidelity model that would just provide the magnitude of the QoI desired. Assuming lower fidelity models require less computational resources, in this manner computational efficiency can be realized provided the QoI value can be accurately and efficiently evaluated. This work utilized Generalized Perturbation Theory (GPT) to evaluate the QoI, by convoluting the GPT solution with the residual of the highest fidelity model determined using the solution from lower fidelity models. Specifically, a reactor core neutronics problem and thermal-hydraulics problem were studied to develop and utilize AMoR. The highest fidelity neutronics model was based upon the 3D space-time, two-group, nodal diffusion equations as solved in the NESTLE computer code. Added to the NESTLE code was the ability to determine the time-dependent GPT neutron flux. The lower fidelity neutronics model was based upon the point kinetics equations along with utilization of a prolongation operator to determine the 3D space-time, two-group flux. The highest fidelity thermal-hydraulics model was based upon the space-time equations governing fluid flow in a closed channel around a heat generating fuel rod. The Homogenous Equilibrium Mixture (HEM) model was used for the fluid and Finite Difference Method was applied to both the coolant and fuel pin energy conservation equations. The lower fidelity thermal-hydraulic model was based upon the same equations as used for the highest fidelity model but now with coarse spatial meshing, corrected somewhat by employing effective fuel heat conduction values. The effectiveness of switching between the highest fidelity model and lower fidelity model as a function of time was assessed using the neutronics problem. Based upon work completed to date, one concludes that the time switching is effective in annealing out differences between the highest and lower fidelity solutions. The effectiveness of using a lower fidelity GPT solution, along with a prolongation operator, to estimate the QoI was also assessed. The utilization of a lower fidelity GPT solution was done in an attempt to avoid the high computational burden associated with solving for the highest fidelity GPT solution. Based upon work completed to date, one concludes that the lower fidelity adjoint solution is not sufficiently accurate with regard to estimating the QoI; however, a formulation has been revealed that may provide a path for addressing this shortcoming.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Full Core Multiphysics Simulation with Offline Mesh Deformation

In this report, building on previous reports issued in FY13 we describe our continued efforts to integrate thermal/hydraulics, neutronics, and structural mechanics modeling codes to perform coupled analysis of a representative fast sodium-cooled reactor core. The focus of the present report is a full core simulation with off-line mesh deformation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SHARP Multiphysics Tutorials

SHARP, developed under the NEAMS Reactor Product Line, is an advanced modeling and simulation toolkit for the analysis of advanced nuclear reactors. SHARP is comprised of three physics modules currently including neutronics, thermal hydraulics, and structural mechanics. SHARP empowers designers to produce accurate results for modeling physical phenomena that have been identified as important for nuclear reactor analysis. SHARP can use existing physics codes and take advantage of existing infrastructure capabilities in the MOAB framework and the coupling driver/solver library, the Coupled Physics Environment (CouPE), which utilizes the widely used, scalable PETSc library. This report aims at identifying the coupled-physics simulation capability of SHARP by introducing the demonstration example called sahex in advance of the SHARP release expected by Mar 2016. sahex consists of 6 fuel pins with cladding, 1 control rod, sodium coolant and an outer duct wall that encloses all the other components. This example is carefully chosen to demonstrate the proof of concept for solving more complex demonstration examples such as EBR II assembly and ABTR full core. The workflow of preparing the input files, running the case and analyzing the results is demonstrated in this report. Moreover, an extension of the sahex model called sahex_core, which adds six homogenized neighboring assemblies to the full heterogeneous sahex model, is presented to test homogenization capabilities in both Nek5000 and PROTEUS. Some primary information on the configuration and build aspects for the SHARP toolkit, which includes capability to auto-download dependencies and configure/install with optimal flags in an architecture-aware fashion, is also covered by this report. A step-by-step instruction is provided to help users to create their cases. Details on these processes will be provided in the SHARP user manual that will accompany the first release.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

(Project 13-5292) Correlating thermal and mechanical coupling based multiphysics behavior of nuclear materials through in-situ measurements

Irradiations and post characterization experiments were performed first on Zr samples. This step will help understand the effect of the 2.5% alloying elements on the behavior of Zircaloy-4 (PWR cladding material) when compared to pure Zr. Irradiation flux measurements and sample temperature calibrations were performed at different energies prior to the irradiation experiments. Irradiations were performed with two different energy regimes1: non-displacment energies and displacement energies. Time was also dedicated to optimize transmission electron microscopy (TEM) sample preparation conditions via electropolishing technique. This step is crucial to prepare TEM samples for the in-situ TEM/irradiation experiments (Year 2). In addition, Zircaloy-4 samples are being prepared for irradiation, and a setup is built by one of our collaborators (Dr. Mert Efe) to prepare ultrafine (UF) and nanocrystalline (NC) Zircaloy-4 samples for comparison with the commercial Zircaloy-4 samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗