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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Beyond the Standard Model physics prospects at the Deep Underground Neutrino Experiment

The Deep Underground Neutrino Experiment (DUNE) is an international project for neutrino physics and proton-decay searches, currently in the design and planning stages. Once built, DUNE will consist of two detectors exposed to the world's most intense neutrino beam. The near detector at Fermilab will record neutrino interactions near the beginning of the beam line. The other, much larger, detector, comprising four 10-kton liquid argon time projection chambers (LArTPCs), will be installed at a depth of 1.5 km at the Sanford Underground Research Facility in South Dakota, about 1,300 km away from the neutrino source. The unique combination of the high-intensity neutrino beam, DUNE's high-resolution near detector system, and massive LArTPC far detector enables a variety of probes of Beyond-the-Standard-Model (BSM) physics, from the discovery of new particles (sterile neutrinos or dark matter), to precision tests of beyond the three-flavor mixing paradigm, non-standard neutrino interactions (NSIs), heavy neutral leptons, and the detailed study of rare processes (e.g., neutrino trident production). This article reviews these physics topics and discusses the prospects for their discovery at the DUNE experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Multi-physics modeling of grain growth during solidification in electron beam additive manufacturing of Inconel 718

While experimental work has shown promising results regarding control of additive manufacturing metal grain structure, the effects of processing parameters on the grain structure is difficult to understand and predict from experiment alone. To this end, a modeling framework is developed which sequentially couples a macro-scale, semi-analytic thermal model, and a meso-scale, cellular automata-based microstructure model. This framework is applied to electron beam additive manufacturing of Inconel 718 using a complex spot scan pattern. The model shows that, with the same scan pattern, variations in the spot time and electron-beam current produce thermal histories with significant spatial and temporal differences, which then produce complex solidification conditions from the interplay between molten pools in the same layer and subsequent layers, resulting in vastly different grain structures. Additionally, it is noted that the framework can significantly reduce the computational expenses for coupled thermal-metallurgical problems, and has the potential to be used for component level problems.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Multi-Higgs doublet models: physical parametrization, sum rules and unitarity bounds

If the scalar sector of the Standard Model is non-minimal, one might expect multiple generations of the hypercharge-1/2 scalar doublet analogous to the generational structure of the fermions. In this work, we examine the structure of a Higgs sector consisting of N Higgs doublets (where N ≥ 2). It is particularly convenient to work in the so-called charged Higgs basis, in which the neutral Higgs vacuum expectation value resides entirely in the first Higgs doublet, and the charged components of remaining N - 1 Higgs doublets are mass-eigenstate fields. We elucidate the interactions of the gauge bosons with the physical Higgs scalars and the Goldstone bosons and show that they are determined by an N × 2N matrix. This matrix depends on (N - 1)(2N - 1) real parameters that are associated with the mixing of the neutral Higgs fields in the charged Higgs basis. Among these parameters, N - 1 are unphysical (and can be removed by rephasing the physical charged Higgs fields), and the remaining 2(N - 1) 2 parameters are physical. We also demonstrate a particularly simple form for the cubic interaction and some of the quartic interactions of the Goldstone bosons with the physical Higgs scalars. These results are applied in the derivation of Higgs coupling sum rules and tree-level unitarity bounds that restrict the size of the quartic scalar couplings. In particular, new applications to three Higgs doublet models with an order-4 CP symmetry and with a Z 3 symmetry, respectively, are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Baryon isocurvature scenario in inflationary cosmology - A particle physics model and its astrophysical implications

A phenomenological model to produce isocurvature baryon-number fluctuations is proposed in the framework of inflationary cosmology. The resulting spectrum of density fluctuation is very different from the conventional Harrison-Zel'dovich shape. The model, with the parameters satisfying several requirements from particle physics and cosmology, provides an appropriate initial condition for the minimal baryon isocurvature scenario of galaxy formation discussed by Peebles.

Yokoyama, Jun'ichi↗
NASA NTRS2015

Physical Models of Layered Polar Firn Brightness Temperatures from 0.5 to 2 GHz

We investigate physical effects influencing 0.5-2 GHz brightness temperatures of layered polar firn to support the Ultra Wide Band Software Defined Radiometer (UWBRAD) experiment to be conducted in Greenland and in Antarctica. We find that because ice particle grain sizes are very small compared to the 0.5-2 GHz wavelengths, volume scattering effects are small. Variations in firn density over cm- to m-length scales, however, cause significant effects. Both incoherent and coherent models are used to examine these effects. Incoherent models include a 'cloud model' that neglects any reflections internal to the ice sheet, and the DMRT-ML and MEMLS radiative transfer codes that are publicly available. The coherent model is based on the layered medium implementation of the fluctuation dissipation theorem for thermal microwave radiation from a medium having a nonuniform temperature. Density profiles are modeled using a stochastic approach, and model predictions are averaged over a large number of realizations to take into account an averaging over the radiometer footprint. Density profiles are described by combining a smooth average density profile with a spatially correlated random process to model density fluctuations. It is shown that coherent model results after ensemble averaging depend on the correlation lengths of the vertical density fluctuations. If the correlation length is moderate or long compared with the wavelength (approximately 0.6x longer or greater for Gaussian correlation function without regard for layer thinning due to compaction), coherent and incoherent model results are similar (within approximately 1 K). However, when the correlation length is short compared to the wavelength, coherent model results are significantly different from the incoherent model by several tens of kelvins. For a 10-cm correlation length, the differences are significant between 0.5 and 1.1 GHz, and less for 1.1-2 GHz. Model results are shown to be able to match the v-pol SMOS data closely and predict the h-pol data for small observation angles.

layered polar firn emission↗
DOE OSTI2019

Standard Model Physics at the HL-LHC and HE-LHC

The successful operation of the Large Hadron Collider (LHC) and the excellent performance of the ATLAS, CMS, LHCb and ALICE detectors in Run-1 and Run-2 with $pp$ collisions at center-of-mass energies of 7, 8 and 13 TeV as well as the giant leap in precision calculations and modeling of fundamental interactions at hadron colliders have allowed an extraordinary breadth of physics studies including precision measurements of a variety physics processes. The LHC results have so far confirmed the validity of the Standard Model of particle physics up to unprecedented energy scales and with great precision in the sectors of strong and electroweak interactions as well as flavour physics, for instance in top quark physics. The upgrade of the LHC to a High Luminosity phase (HL-LHC) at 14 TeV center-of-mass energy with 3 ab$^{-1}$ of integrated luminosity will probe the Standard Model with even greater precision and will extend the sensitivity to possible anomalies in the Standard Model, thanks to a ten-fold larger data set, upgraded detectors and expected improvements in the theoretical understanding. This document summarises the physics reach of the HL-LHC in the realm of strong and electroweak interactions and top quark physics, and provides a glimpse of the potential of a possible further upgrade of the LHC to a 27 TeV $pp$ collider, the High-Energy LHC (HE-LHC), assumed to accumulate an integrated luminosity of 15 ab$^{-1}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

GeoThermalCloud framework for fusion of big data and multi-physics models in Nevada and Southwest New Mexico

Our GeoThermalCloud framework is designed to process geothermal datasets using a novel toolbox for unsupervised and physics-informed machine learning called SmartTensors. More information about GeoThermalCloud can be found at the GeoThermalCloud GitHub Repository. More information about SmartTensors can be found at the SmartTensors Github Repository and the SmartTensors page at LANL.gov. Links to these pages are included in this submission. GeoThermalCloud.jl is a repository containing all the data and codes required to demonstrate applications of machine learning methods for geothermal exploration. GeoThermalCloud.jl includes: - site data - simulation scripts - jupyter notebooks - intermediate results - code outputs - summary figures - readme markdown files GeoThermalCloud.jl showcases the machine learning analyses performed for the following geothermal sites: - Brady: geothermal exploration of the Brady geothermal site, Nevada - SWNM: geothermal exploration of the Southwest New Mexico (SWNM) region - GreatBasin: geothermal exploration of the Great Basin region, Nevada Reports, research papers, and presentations summarizing these machine learning analyses are also available and will be posted soon.

15 GEOTHERMAL ENERGY↗
NASA NTRSDate not supplied

Sampling the ENSO response to tropical volcanic eruptions in the GISS-E2-1 model: physical & biogeochemical ocean effects

Despite observational evidence that an El-Niño is likely to follow tropical volcanic eruptions (Predybaylo et al, 2017), modelling studies have so far reached no consensus on either the sign or physical mechanism of El Niño response to volcanism (Khodri et al, 2017) . A previous modelling study found that following the Pinatubo eruption, the tropical Pacific developed El Niño-like conditions that triggered a decrease in the upwelling of DIC and increased carbon uptake (Eddebbar et al., 2019). However, the carbon cycle response to volcanisms as it depends on ENSO phase remains unclear. Here we perform thirty ensemble control simulations using the GISS-E2-1-G model to sample the pre-industrial control ENSO cycle and thirty more to superimpose a Pinatubo-like volcanic eruption. We find that the model tends to favor the development of a weak El-Niño following the volcanic eruption only in the case of a volcano that erupts at the start of a neutral ENSO phase. In response to the weak El-Niño, some anomalies have the same sign as those that develop from a strong El-Niño, such as weaker winds, shallower mixed layer depths, less direct but more diffuse radiation at the surface of the ocean, less nitrate, and lower primary productivity. However, other anomalies show the opposite sign as those from a strong El-Niño, such as more precipitation, less surface DIC, and less surface ocean pCO 2 . In contrast to the neutral phase, the volcanic response during positive ENSO and negative ENSO phases are swamped by the model’s internal variability.

ENSO↗
NASA NTRS1989

A physical model for predicting bidirectional reflectances over bare soil

While most previous attempts to retrieve soil surface optical property characteristics have proceeded through a fitting of empirical functions to the data, an optimization technique is presently applied to a physically-based surface reflectance model developed for the study of planetary surfaces. This inversion procedure is shown to allow the direct estimation of the single-scattering coefficient, two parameters describing the 'hot spot' phenomenon, and two parameters describing the scattering phase function. A comparison of inversion technique results with both synthetic data and actual observations shows the model to be capable of predicting the observed bidirectional reflectances as well as directional-hemispherical reflectances; it can also build the complete radiance field over the upward hemisphere.

Pinty, Bernard↗
NASA NTRS1992

Physical model and experimental results of cathode erosion related to power supply ripple

This paper discusses the physical effects of power supply ripple on cathode erosion and cathode arc attachment in a water-cooled, 30 kW nitrogen arcjet. Experimental results are presented for 2 percent thoriated tungsten, which show that the long-term cathode erosion rate is a decreasing function of current ripple over the range 1-13 percent. Above this range, the cathode discharge becomes unstable, and the erosion rate rapidly increases. A qualitative model of this effect is given in terms of a magnetically induced radial motion of the arc column, and an overall increase in the cathode spot radius due to the higher peak current associated with higher ripple. The most important effect of power supply ripple is therefore shown to be its ability to collectively drive the cathode attachment away from the cathode center. This leads to an increase in the cathode attachment area, and a subsequent decrease in the cathode erosion rate.

Harris, W. J.↗
DOE OSTI2024

Multi-physics Modeling of Radiative Heat Transfer and Fluid Flow for the Reactor Cavity Cooling System

High-temperature gas-cooled reactors (HTGRs) are notable for their high thermal efficiency and potential for combined heat and power applications. These reactors are particularly appealing due to their advanced passive safety features. HTGRs utilize passive safety systems that function without requiring active components like pumps or compressors during emergencies. These reactor designs depend on a Reactor Cavity Cooling System (RCCS) to manage decay heat removal from the reactor pressure vessel (RPV) during accident conditions. The RCCS consists of vertical rectangular channels known as "risers" or riser ducts positioned around the RPV. These risers receive heat from the RPV through both convective and radiative heat transfer mechanisms. Understanding the interplay of multiple physical phenomena, such as fluid dynamics, heat transfer, and neutron interactions, is essential for the effective design and operation of nuclear reactors, particularly for systems like the RCCS. Multi-physics simulations provide a comprehensive approach to studying these interactions, offering detailed insights and enhancing accuracy. They are especially important in RCCS designs, where the interaction between radiative and convective heat transfer can significantly impact system performance. By leveraging multi-physics simulations, complex reactor behaviors can be modeled without compromising the fidelity of the underlying physical processes. This work aims to establish a robust methodology for coupling multiple physical processes in an air-cooled RCCS. By focusing on validating this multi-physics approach, the study involves designing test cases that simulate various conditions to verify the numerical models employed. The outcomes of this research will provide critical insights for accurately modeling and optimizing complex nuclear systems like the RCCS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

ExaWind: Exascale Predictive Wind Plant Flow Physics Modeling

The scientific goal of the ExaWind project is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. The primary application codes in the ExaWind environment are Nalu-Wind, an unstructured-grid computational fluid dynamics (CFD) code, AMR-Wind, a structured-grid CFD code, and OpenFAST, a whole-turbine simulation code. In this poster we present our progress towards simulating the ExaWind challenge problem, which is a predictive simulation of a wind farm with tens of megawatt-scale wind turbines dispersed over an area of 50 square kilometers.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2024

Centrifuge Physical Modeling of Underground Chemical Explosions

This report presents development of a new meso-scale geomechanics capability using an 8.84-m radius centrifuge for physically simulating underground chemical explosions and stress cage processes in scaled models.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Two-fluid physical modeling of superconducting resonators in the ARTEMIS framework

In this work, we implement a new London equation module for superconductivity in the GPU-enabled ARTEMIS framework, and couple it to a finite-difference time-domain solver for Maxwell's equations. We apply this two-fluid approach to model a superconducting coplanar waveguide (CPW) resonator. We validate our implementation by verifying that the theoretical skin depth and reflection coefficients can be obtained for several superconductive materials, with different London penetration depths, over a range of frequencies. Our convergence studies show that the algorithm is second-order accurate in both space and time, except at superconducting interfaces where the approach is spatially first-order. In our CPW simulations, we leverage the GPU scalability of our code to compare the two-fluid model to more traditional approaches that approximate superconducting behavior and demonstrate that superconducting physics can show comparable performance to the assumption of quasi-infinite conductivity as measured by the Q-factor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Diffstar: a fully parametric physical model for galaxy assembly history

We present Diffstar, a smooth parametric model for the in situ star formation history (SFH) of galaxies. The Diffstar model is distinct from traditional SFH models because it is parametrized directly in terms of basic features of galaxy formation physics. Diffstar includes ingredients for: the halo mass assembly history; the accretion of gas into the dark matter halo; the fraction of gas that is eventually transformed into stars, ϵ ms ; the time-scale over which this transformation occurs, τ cons ; and the possibility that some galaxies will experience a quenching event at time t q , and may subsequently experience rejuvenated star formation. We show that our model is sufficiently flexible to describe the average stellar mass histories of galaxies in both the IllustrisTNG (TNG) and UniverseMachine (UM) simulations with an accuracy of ~0.1 dex across most of cosmic time. We use Diffstar to compare TNG to UM in common physical terms, finding that: (i) star formation in UM is less efficient and burstier relative to TNG; (ii) UM galaxies have longer gas consumption time-scales, relative to TNG; (iii) rejuvenated star formation is ubiquitous in UM, whereas quenched TNG galaxies rarely experience sustained rejuvenation; and (iv) in both simulations, the distributions of ϵ ms , τ cons , and t q share a common characteristic dependence upon halo mass, and present significant correlations with halo assembly history. Finally, we conclude with a discussion of how Diffstar can be used in future applications to fit the SEDs of individual observed galaxies, as well as in forward-modelling applications that populate cosmological simulations with synthetic galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PREDICTING CME EJECTA AND SHEATH FRONT ARRIVAL AT L1 WITH A DATA-CONSTRAINED PHYSICAL MODEL

We present a method for predicting the arrival of a coronal mass ejection (CME) flux rope in situ, as well as the sheath of solar wind plasma accumulated ahead of the driver. For faster CMEs, the front of this sheath will be a shock. The method is based upon geometrical separate measurement of the CME ejecta and sheath. These measurements are used to constrain a drag-based model, improved by including both a height dependence and accurate de-projected velocities. We also constrain the geometry of the model to determine the error introduced as a function of the deviation of the CME nose from the Sun–Earth line. The CME standoff-distance in the heliosphere fit is also calculated, fit, and combined with the ejecta model to determine sheath arrival. Combining these factors allows us to create predictions for both fronts at the L1 point and compare them against observations. We demonstrate an ability to predict the sheath arrival with an average error of under 3.5 hr, with an rms error of about 1.58 hr. For the ejecta the error is less than 1.5 hr, with an rms error within 0.76 hr. We also discuss the physical implications of our model for CME expansion and density evolution. We show the power of our method with ideal data and demonstrate the practical implications of having a permanent L5 observer with space weather forecasting capabilities, while also discussing the limitations of the method that will have to be addressed in order to create a real-time forecasting tool.

79 ASTRONOMY AND ASTROPHYSICS↗