Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Modeling”

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 451 records · Page 25

Aerothermal modeling program, phase 1

The physical modeling embodied in the computational fluid dynamics codes is discussed. The objectives were to identify shortcomings in the models and to provide a program plan to improve the quantitative accuracy. The physical models studied were for: turbulent mass and momentum transport, heat release, liquid fuel spray, and gaseous radiation. The approach adopted was to test the models against appropriate benchmark-quality test cases from experiments in the literature for the constituent flows that together make up the combustor real flow.

Sturgess, G. J.↗

Master-slave system with force feedback based on dynamics of virtual model

A master-slave system can extend manipulating and sensing capabilities of a human operator to a remote environment. But the master-slave system has two serious problems: one is the mechanically large impedance of the system; the other is the mechanical complexity of the slave for complex remote tasks. These two problems reduce the efficiency of the system. If the slave has local intelligence, it can help the human operator by using its good points like fast calculation and large memory. The authors suggest that the slave is a dextrous hand with many degrees of freedom able to manipulate an object of known shape. It is further suggested that the dimensions of the remote work space be shared by the human operator and the slave. The effect of the large impedance of the system can be reduced in a virtual model, a physical model constructed in a computer with physical parameters as if it were in the real world. A method to determine the damping parameter dynamically for the virtual model is proposed. Experimental results show that this virtual model is better than the virtual model with fixed damping.

Nojima, Shuji↗

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling↗

The GBS code for tokamak scrape-off layer simulations

We describe a new version of GBS, a 3D global, flux-driven plasma turbulence code to simulate the turbulent dynamics in the tokamak scrape-off layer (SOL), superseding the code presented by Ricci et al. (2012) [14]. The present work is driven by the objective of studying SOL turbulent dynamics in medium size tokamaks and beyond with a high-fidelity physics model. We emphasize an intertwining framework of improved physics models and the computational improvements that allow them. The model extensions include neutral atom physics, finite ion temperature, the addition of a closed field line region, and a non-Boussinesq treatment of the polarization drift. GBS has been completely refactored with the introduction of a 3-D Cartesian communicator and a scalable parallel multigrid solver. We report dramatically enhanced parallel scalability, with the possibility of treating electromagnetic fluctuations very efficiently. The method of manufactured solutions as a verification process has been carried out for this new code version, demonstrating the correct implementation of the physical model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting the performance uncertainty of a 1-MW pilot-scale carbon capture system after hierarchical laboratory-scale calibration and validation

A challenging problem in designing pilot-scale carbon capture systems is to predict, with uncertainty, the adsorber performance and capture efficiency under various operating conditions where no direct experimental data exist. Motivated by this challenge, we previously proposed a hierarchical framework in which relevant parameters of physical models were sequentially calibrated from different laboratory-scale carbon capture unit (C2U) experiments. Specifically, three models of increasing complexity were identified based on the fundamental physical and chemical processes of the sorbent-based carbon capture technology. Results from the corresponding laboratory experiments were used to statistically calibrate the physical model parameters while quantifying some of their inherent uncertainty. The parameter distributions obtained from laboratory-scale C2U calibration runs are used in this study to facilitate prediction at a larger scale where no corresponding experimental results are available. In this paper, we first describe the multiphase reactive flow model for a sorbent-based 1-MW carbon capture system then analyze results from an ensemble of simulations with the upscaled model. The simulation results are used to quantify uncertainty regarding the design’s predicted efficiency in carbon capture. In particular, we determine the minimum gas flow rate necessary to achieve 90% capture efficiency with 95% confidence.

bubbling bed↗

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗

Research reactor plate fuel multi-physics hotspot model for the HFIR

Ongoing development of a multi-physics hotspot model for nuclear research reactor plate fuel is presented with a focus on the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). The background of plate fueled research reactors in general and the HFIR in particular are discussed along with the need to convert from high enriched uranium (HEU) to low enriched uranium (LEU) fuel, which drives the need to improve the analysis capability of fuel thermal-fluid performance. HFIR HEU fuel plate fabrication is described to explain the origin of fuel segregations and cladding non-bonds that currently contribute to the hotspot flux peaking factors. Previous efforts at modeling hotspots caused by these defects are reviewed and areas for improvement associated with high resolution multi-physics simulation are identified. Modes of heat transfer through non-bonds, including thermal contact conductance, are examined. Conjugate heat transfer over hotspots is also explored to high resolution up to and through the wall boundary. The challenges associated with multi-physics verification and validation (VV), particularly the requirements beyond VV 20, are presented. An approach to Validation within the limits of existing data for conduction heat transfer through the fuel plate, and convection into the cooling fluid is set forth. Preliminary simulations are presented assessing hot spot modeling approaches for contact conductance, and fluid turbulence modeling within COMSOL. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Robust Statistical Approach for Determination of Graphite Nitridation Using Bayesian Model Comparison

A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at [Formula: see text].

Engineering↗

Bayesian seismic inversion based on rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies

We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to post-stack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rock-physics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Polarimetric signatures of sea ice. 1: Theoretical model

Physical, structral, and electromagnetic properties and interrelating processes in sea ice are used to develop a composite model for polarimetric backscattering signatures of sea ice. Physical properties of sea ice constituents such as ice, brine, air, and salt are presented in terms of their effects on electromagnetic wave interactions. Sea ice structure and geometry of scatterers are related to wave propagation, attenuation, and scattering. Temperature and salinity, which are determining factors for the thermodynamic phase distribution in sea ice, are consistently used to derive both effective permittivities and polarimetric scattering coefficients. Polarmetric signatures of sea ice depend on crystal sizes and brine volumes, which are affected by ice growth rates. Desalination by brine expulsion, drainage, or other mechanisms modifies wave penetration and scattering. Sea ice signatures are further complicated by surface conditions such as rough interfaces, hummocks, snow cover, brine skim, or slush layer. Based on the same set of geophysical parameters characterizing sea ice, a composite model is developed to calculate effective permittivities and backscattering covariance matrices at microwave frequencies to interpretation of sea ice polarimetric signatures.

Nghiem, S. V.↗

Physics-based Models, Machine Learning, and Experiment: Towards Understanding Complex Electrode Degradation

Degradation phenomena in Li-ion batteries are highly complex, coupled, and sensitive to use history and operating conditions. In this study, we show how tracking model parameters in continuum-level physics-based models, expedited by machine learning, can be useful in testing hypotheses for degradation mechanisms. An exemplary analysis using this approach is presented for a set of lithium trivanadate ( L i x V 3 O 8 ) cathodes cycled over a range of current rates. A simple cell revival process is combined with the parameter estimates over the course of cycling to extract valuable insights into cathode evolution and eliminate hypothesized degradation mechanisms for these cathodes. The presented approach is expected to be broadly applicable for degradation analysis of other electrodes.

25 ENERGY STORAGE↗

Deployment and performance of a Low-Energy-Threshold Skipper-CCD inside a nuclear reactor

Charge Coupled Devices (CCD) are being used for reactor neutrino experiments and have already demonstrated their potential in constraining new physics models. The prospect of a Skipper-CCD experiment looking for standard and beyond standard model (BSM) physics in a nuclear reactor has been evaluated for different benchmark scenarios. Here, we report the first installation of a 2-g Skipper-CCD inside the containment building of a 2 GW th nuclear power plant and analyze its performance throughout its first 18 months of operation. The sensor was successfully deployed at Atucha II, in Argentina, 12 meters away from the center of the reactor core. We discuss the challenges involved in the commissioning of the detector and present data acquired during reactor ON and reactor OFF periods, with the sensor functioning with a sub-electron readout noise of 0.17 e - . Based on an exposure of 56.8 g day reactor ON and two reactor OFF data sets with a total exposure of 118.1 g day we characterize the system and evaluate the sensitivity to CEvNS. We achieved a background rate of 33 kdru and a low threshold of 45 eV ee . The ongoing efforts to improve sensitivities to CEvNS and BSM interaction are also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗

Life Cycle Assessment of Vehicle Lightweighting: A Physics-Based Model To Estimate Use-Phase Fuel Consumption of Electrified Vehicles

Assessing the life-cycle benefits of vehicle lightweighting requires a quantitative description of mass-induced fuel consumption (MIF) and fuel reduction values (FRVs). We have extended our physics-based model of MIF and FRVs for internal combustion engine vehicles (ICEVs) to electrified vehicles (EVs) including hybrid electric vehicles (HEVs), plug-in hybrid electric vehicles (PHEVs), and battery electric vehicles (BEVs). We illustrate the utility of the model by calculating MIFs and FRVs for 37 EVs and 13 ICEVs. BEVs have much smaller MIF and FRVs, both in the range 0.04-0.07 L e /(100 km 100 kg), than those for ICEVs which are in the ranges 0.19-0.32 and 0.16-0.22 L/(100 km 100 kg), respectively. The MIF and FRVs for HEVs and PHEVs mostly lie between those for ICEVs and BEVs. Powertrain resizing increases the FRVs for ICEVs, HEVs and PHEVs. Lightweighting EVs is less effective in reducing greenhouse gas emissions than lightweighting ICEVs, however the benefits differ substantially for different vehicle models. The physics-based approach outlined here enables model specific assessments for ICEVs, HEVs, PHEVs, and BEVs required to determine the optimal strategy for maximizing the life-cycle benefits of lightweighting the light-duty vehicle fleet.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Development of a Physical-Topological Model for the Response of a High-Power Vertical DMOS Transistor to the Effect of Pulsed Gamma-Radiation

A method for developing models of semiconductor devices with two-dimensional nonuniform concentration profiles of donors and acceptors in working regions of a semiconductor device structure is for the first time proposed based on a complex of physical-topological modeling of charge carrier transport and process simulation of the forming processes of the device structure. The application of process simulation is due to the need to correctly define the parameters of the semiconductor device structure, which are used as initial data to calculate the electron transport according to the physical-topological model. The parameters of production processes of ion implantation, diffusion, and lithography, which were refined during process simulation, are determined for a high-power metal—oxide—semiconductor (MOS) transistor forming by the double diffusion method in accordance with known electrical characteristics and measured geometrical sizes of the structure. This results in two-dimensional distribution profiles of donors and acceptors in p–n junctions necessary to calculate the transistor breakdown under the effect of pulsed γ radiation. Breakdown processes are modeled with the help of the physical-topological model based on the Poisson and continuity equations as well as expressions for the diffusion and drift current densities in the transistor. Accounting for carriers formed at the instant of γ irradiation is implemented by the introduction of the dependence of the generation coefficient of electron–hole pairs on the radiation-dose power. The results of calculations correlate well with the experimental data, which makes it possible to give a conclusion regarding the adequacy of the proposed complex model.

36 MATERIALS SCIENCE↗

Warm neutral halos around molecular clouds. VI - Physical and chemical modeling

A combined physical and chemical modeling of the halos around molecular clouds is presented, with special emphasis on the H-to-H2 transition. On the basis of H I 21 cm observations, it is shown that the halos are extended. A physical model is employed in conjunction with a chemistry code to provide a self-consistent description of the gas. The radiative transfer code provides a check with H I, CO, and OH observations. It is concluded that the warm neutral halos are not gravitationally bound to the underlying molecular clouds and are isobaric. It is inferred from the observed extent of the H I envelopes and the large observed abundance of OH in them that the generally accepted rate for H2 information on grains is too large by a factor of two to three.

Andersson, B.-G.↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessment of fission product yields data needs in nuclear reactor applications

Studies on the build-up of fission products in fast reactors have been performed, with particular emphasis on the effects related to the physics of the nuclear fission process. Fission product yields, which are required for burn-up calculations, depend on the proton and neutron number of the target nucleus as well as on the incident neutron energy. Evaluated nuclear data on fission product yields are available for all relevant target nuclides in reactor applications. However, the description of their energy dependence in evaluated data is still rather rudimentary, which is due to the lack of experimental fast fission data and reliable physical models. Additionally, physics studies of evaluated JEFF-3.1.1 fission yields data have shown potential improvements, especially for various fast fission data sets of this evaluation. In recent years, important progress in the understanding of the fission process has been made, and advanced model codes are currently being developed. This paper deals with the semi-empirical approach to the description of the fission process, which is used in the GEF code being developed by K.-H. Schmidt and B. Jurado on behalf of the OECD Nuclear Energy Agency, and with results from the corresponding author's diploma thesis. An extended version of the GEF code, supporting the calculation of spectrum weighted fission product yields, has been developed. It has been applied to the calculation of fission product yields in the fission rate spectra of a MOX fuelled sodium-cooled fast reactor. Important results are compared to JEFF-3.1.1 data and discussed in this paper. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗