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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Computational modelling and experimental tank testing of the multi float WaveSub under regular wave forcing

A submerged wave device generates energy from the relative motion of floating bodies. In WaveSub, three floats are joined to a reactor; each connected to a spring and generator. Electricity generated damps the orbital movements of the floats. The forces are non-linear and each float interacts with the others. Tuning to the wave climate is achieved by changing the line lengths, so there is a need to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design. to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design.

16 TIDAL AND WAVE POWER↗
NASA NTRS2002

Rocket Engine Oscillation Diagnostics

Rocket engine oscillating data can reveal many physical phenomena ranging from unsteady flow and acoustics to rotordynamics and structural dynamics. Because of this, engine diagnostics based on oscillation data should employ both signal analysis and physical modeling. This paper describes an approach to rocket engine oscillation diagnostics, types of problems encountered, and example problems solved. Determination of design guidelines and environments (or loads) from oscillating phenomena is required during initial stages of rocket engine design, while the additional tasks of health monitoring, incipient failure detection, and anomaly diagnostics occur during engine development and operation. Oscillations in rocket engines are typically related to flow driven acoustics, flow excited structures, or rotational forces. Additional sources of oscillatory energy are combustion and cavitation. Included in the example problems is a sampling of signal analysis tools employed in diagnostics. The rocket engine hardware includes combustion devices, valves, turbopumps, and ducts. Simple models of an oscillating fluid system or structure can be constructed to estimate pertinent dynamic parameters governing the unsteady behavior of engine systems or components. In the example problems it is shown that simple physical modeling when combined with signal analysis can be successfully employed to diagnose complex rocket engine oscillatory phenomena.

Nesman, Tom↗
DOE OSTI2020

A partial-volume correction for quantitative spectral x-ray radiography

The technique of spectral radiography has been shown to accurately discriminate a materials elemental composition by using known x-ray physics. This is achieved with a state-of-the-art pixelated spectrometer and a high-fidelity physical model. Similar to other radiographs, the system suffers from partial volume attenuation where the sample projection onto a pixel is only a fraction of the pixel area, leaving the rest unattenuated. This study updates the physical model to account for fractional area which presents new challenges to the existing algorithm. We present how this update is made and explore approaches to stabilize the algorithm through parameter scaling, constrained optimization, and hyperparameter selection.

Gillis, Wesley C.↗
DOE OSTI2019

Stripping off the implementation complexity of physics-based model predictive control for buildings via deep learning

Over the past decade, model predictive control (MPC) has been considered as the most promising solution for intelligent building operation. Despite extensive effort, transfer of this technology into practice is hampered by the need to obtain an accurate controller model with minimum effort, the need of expert knowledge to set it up, and the need of increased computational power and dedicated software to run it. A promising direction that tackles the last two problems was proposed by approximate explicit MPC where the optimal control policies are learned from MPC data via a suitable function approximator, e.g., a deep learning (DL) model. The main advantage of the proposed approach stems from simple evaluation at execution time leading to low computational footprints and easy deployment on embedded HW platforms. We present the energy savings potential of physics-based (also called 'white-box') MPC applied to an office building in Belgium. Moreover, we demonstrate how deep learning approximators can be used to cut the implementation and maintenance costs of MPC deployment without compromising performance. We also critically assess the presented approach by pointing out the major challenges and remaining open-research questions.

Drgona, Jan↗
DOE OSTI2015

The GENIE Neutrino Monte Carlo Generator: Physics and User Manual

GENIE is a suite of products for the experimental neutrino physics community. This suite includes i) a modern software framework for implementing neutrino event generators, a state-of-the-art comprehensive physics model and tools to support neutrino interaction simulation for realistic experimental setups (the Generator product), ii) extensive archives of neutrino, charged-lepton and hadron scattering data and software to produce a comprehensive set of data/MC comparisons (the Comparisons product), and iii) a generator tuning framework and fitting applications (the Tuning product). This book provides the definite guide for the GENIE Generator: It presents the software architecture and a detailed description of its physics model and official tunes. In addition, it provides a rich set of data/MC comparisons that characterise the physics performance of GENIE. Detailed step-by-step instructions on how to install and configure the Generator, run its applications and analyze its outputs are also included.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

COLLISIONS BETWEEN GRAVITY-DOMINATED BODIES. II. THE DIVERSITY OF IMPACT OUTCOMES DURING THE END STAGE OF PLANET FORMATION

Numerical simulations of the stochastic end stage of planet formation typically begin with a population of embryos and planetesimals that grow into planets by merging. We analyzed the impact parameters of collisions leading to the growth of terrestrial planets from recent N-body simulations that assumed perfect merging and calculated more realistic outcomes using a new analytic collision physics model. We find that collision outcomes are diverse and span all possible regimes: hit-and-run, merging, partial accretion, partial erosion, and catastrophic disruption. The primary outcomes of giant impacts between planetary embryos are approximately evenly split between partial accretion, graze-and-merge, and hit-and-run events. To explore the cumulative effects of more realistic collision outcomes, we modeled the growth of individual planets with a Monte Carlo technique using the distribution of impact parameters from N-body simulations. We find that fewer planets reached masses >0.7 M{sub Earth} using the collision physics model compared to simulations that assumed every collision results in perfect merging. For final planets with masses >0.7 M{sub Earth}, 60% are enriched in their core-to-mantle mass fraction by >10% compared to the initial embryo composition. Fragmentation during planet formation produces significant debris ({approx}15% of the final mass) and occurs primarily by erosion of the smaller body in partial accretion and hit-and-run events. In partial accretion events, the target body grows by preferentially accreting the iron core of the projectile and the escaping fragments are derived primarily from the silicate mantles of both bodies. Thus, the bulk composition of a planet can evolve via stochastic giant impacts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ageing Studies of Mega Battery Packs for Grid Storage Applications Using Physics Based Modeling

The lower Levelized Cost of Electricity (LCOE) from wind and solar photovoltaics has enabled for greater integration of variable energy resources with energy storage tech-nologies such as larger centralized lithium-ion battery megapacks (with 1MWh total energy and 500V rating) to provide utility-scale services to grid operators. The daily cycling of standalone lithium-ion grid storage will reduce the battery cells' capacity due to several degradation mechanisms. A comprehensive physics-based Python tm framework called Liionpack was developed to estimate these megapacks' remaining life and ageing. The study includes various degradation mechanisms coupled to the electrochemical-thermal model at the pack level. The effect of the inhomogeneities from cell-to-cell thermodynamic and kinetic properties for different working conditions including temperature and charge/discharge protocols on the ageing of a megapack are presented here.

Karra, Vikrant↗
DOE OSTI2014

The Wind Forecast Improvement Project (WFIP). A Public/Private Partnership for Improving Short Term Wind Energy Forecasts and Quantifying the Benefits of Utility Operations -- the Northern Study Area

This report contains the results from research aimed at improving short-range (0-6 hour) hub-height wind forecasts in the NOAA weather forecast models through additional data assimilation and model physics improvements for use in wind energy forecasting. Additional meteorological observing platforms including wind profilers, sodars, and surface stations were deployed for this study by NOAA and DOE, and additional meteorological data at or near wind turbine hub height were provided by South Dakota State University and WindLogics/NextEra Energy Resources over a large geographical area in the U.S. Northern Plains for assimilation into NOAA research weather forecast models. The resulting improvements in wind energy forecasts based on the research weather forecast models (with the additional data assimilation and model physics improvements) were examined in many different ways and compared with wind energy forecasts based on the current operational weather forecast models to quantify the forecast improvements important to power grid system operators and wind plant owners/operators participating in energy markets. Two operational weather forecast models (OP_RUC, OP_RAP) and two research weather forecast models (ESRL_RAP, HRRR) were used as the base wind forecasts for generating several different wind power forecasts for the NextEra Energy wind plants in the study area. Power forecasts were generated from the wind forecasts in a variety of ways, from very simple to quite sophisticated, as they might be used by a wide range of both general users and commercial wind energy forecast vendors. The error characteristics of each of these types of forecasts were examined and quantified using bulk error statistics for both the local wind plant and the system aggregate forecasts. The wind power forecast accuracy was also evaluated separately for high-impact wind energy ramp events. The overall bulk error statistics calculated over the first six hours of the forecasts at both the individual wind plant and at the system-wide aggregate level over the one year study period showed that the research weather model-based power forecasts (all types) had lower overall error rates than the current operational weather model-based power forecasts, both at the individual wind plant level and at the system aggregate level. The bulk error statistics of the various model-based power forecasts were also calculated by season and model runtime/forecast hour as power system operations are more sensitive to wind energy forecast errors during certain times of year and certain times of day. The results showed that there were significant differences in seasonal forecast errors between the various model-based power forecasts. The results from the analysis of the various wind power forecast errors by model runtime and forecast hour showed that the forecast errors were largest during the times of day that have increased significance to power system operators (the overnight hours and the morning/evening boundary layer transition periods), but the research weather model-based power forecasts showed improvement over the operational weather model-based power forecasts at these times.

17 WIND ENERGY↗
DOE OSTI2014

Temperature measuring analysis of the nuclear reactor fuel assembly

Study was based on rapid changes of measured temperature values from the thermocouple in the VVER 440 nuclear reactor fuel assembly. Task was to determine origin of fluctuations of the temperature values by experiments on physical model of the fuel assembly. During an experiment, heated water was circulating in the system and cold water inlet through central tube to record sensitivity of the temperature sensor. Two positions of the sensor was used. First, just above the central tube in the physical model fuel assembly axis and second at the position of the thermocouple in the VVER 440 nuclear reactor fuel assembly. Dependency of the temperature values on time are presented in the diagram form in the paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Neutral Current Elastic Scattering and the Strange Spin Structure of the Proton

Neutrinos can be used as a unique and informative probe of the structure within nucleons. This thesis presents the tools and methodology for studying the strange quark spin in the nucleon through neutral current elastic neutrino-proton scattering in the MicroBooNE experiment located at the Fermi National Accelerator Lab. An automated boosted decision tree based proton identification algorithm is used to select proton in liquid argon TPC data with a 70\% efficiency. After a set of protons are selected, a logistic regression model is used to select neutral current elastic proton interactions in MicroBooNE. Neutral current elastic proton interactions are selected with an 11\% efficiency and 30\% purity. The number of selected events measured in MicroBooNE data as a function of $Q^2$ is compared directly to the number of selected events in MicroBooNE simulation. An event reweighting scheme is used to vary the expected number of eventsin MicroBooNE simulation based on the values o f the physics parameters of interest. The likelihood for any given value of the physical parameters can then be easily calculated and the probability distributions of the physics parameters can be sampled using Markov Chain Monte Carlo. Future improvements to the detector physics models in MicroBooNE are required for a precise determination of the strange quark spin structure in the proton. With the current level of uncertainty on the detector physics models, we are only able to constrain the net strange quark spin, $\Delta s$ to the range $-1.8 < \Delta s < 3.8$ with 95\% confidence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1990

CFD validation for aerodynamic flows - Challenge for the '90's

The process of computational fluid dynamics (CFD) validation is described from two perspectives, numerical error validation and physical model validation. Errors associated with time and space discretization, grid refinement, numerical dissipation, and level of solution convergence, are all possible sources of numerical error. The second aspect of CFD code validation is associated with the governing equations and the physical models, e.g., chemistry and turbulence models. This type of validation requires comparisons with carefully conducted experiments and is the primary aspect of the validation process discussed in this paper. Examples of validation are shown for selected Reynolds-averaged Navier-Stokes codes for conditions from transonic through hypersonic. Since code applications are becoming more complex, it no longer suffices to use data from surface or integral quantities alone to provide the required validation. Flow field surveys and experimental boundary condition measurements are emerging as critical data that must be obtained for CFD code validation. Progress in the overall level of validation is shown to be improving. However, future challenges remain, some of which are enumerated.

Marvin, Joseph G.↗
DOE OSTI2018

Mechanism of Optical Charge Exchange of Magnetic Centers in BSO : Fe

High-resistance photosensitive crystals of Bi{sub 12}SiO{sub 20} (BSO) doped with iron ions were studied. X-ray diffraction analysis reveals the compression of a unit cell in a BSO : Fe crystal with increasing impurity concentration. Electron paramagnetic resonance demonstrates a decrease in the intensity of the EPR signal when the BSO: Fe crystal is exposed to light that generates photocarriers. It is found that the characteristic time of the EPR signal decrease is close to the value of the Maxwellian relaxation time measured with the help of the longitudinal electrooptical effect. A physical model of the mechanism of optical charge exchange of magnetic iron centers is discussed, based on the statement that the nature of the crystal bonds of the iron ion with ligands without structural modification of the crystal lattice changes during the photogeneration of carriers. A physical model is proposed, according to which a trivalent Fe{sup 3+} ion transforms into a divalent state of Fe{sup 2+} with a change in the total spin from 5/2 to 2. The compression of the unit cell with increasing iron ion concentration in the framework of the model under discussion is due to the transformation of atomic orbitals upon replacement of silicon ions by iron ions. The transformation process affects the cells unoccupied by iron, which is proved by the absence of a bifurcation of X-ray reflections and indicates the long-range nature of the intracrystalline interactions in sillenites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

SU-E-J-149: Secondary Emission Detection for Improved Proton Relative Stopping Power Identification

Purpose: This research investigates application of secondary prompt gamma (PG) emission spectra, resulting from nuclear reactions induced by protons, to characterize tissue composition along the particle path. The objective of utilizing the intensity of discrete high-energy peaks of PG is to improve the accuracy of relative stopping power (RSP) values available for proton therapy treatment planning on a patient specific basis and to reduce uncertainty in dose depth calculations. Methods: In this research, MCNP6 was used to simulate PG emission spectra generated from proton induced nuclear reactions in medium of varying composition of carbon, oxygen, calcium and nitrogen, the predominant elements found in human tissue. The relative peak intensities at discrete energies predicted by MCNP6 were compared to the corresponding atomic composition of the medium. Results: The results have shown a good general agreement with experimentally measured values reported by other investigators. Unexpected divergence from experimental spectra was noted in the peak intensities for some cases depending on the source of the cross-section data when using compiled proton table libraries vs. physics models built into MCNP6. While the use of proton cross-section libraries is generally recommended when available, these libraries lack data for several less abundant isotopes. This limits the range of their applicability and forces the simulations to rely on physics models for reactions with natural atomic compositions. Conclusion: Current end-of-range proton imaging provides an average RSP for the total estimated track length. The accurate identification of tissue composition along the incident particle path using PG detection and characterization allows for improved determination of the tissue RSP on the local level. While this would allow for more accurate depth calculations resulting in tighter treatment margins, precise understanding of proton beam behavior in tissue of various compositions is necessary requiring detailed simulations with a high degree of accuracy.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1979

Equilibrium velocities of a planetesimal population

The random velocities of planetesimal populations specified by maximum and minimum masses and a power-law number density are calculated iteratively. The calculation is based on two different physical models involving ratios of rates: (1) excitation of kinetic energy by gravitational perturbation and elastic collision equal to damping of kinetic energy by inelastic collisions; and (2) excitation of kinetic energy a ratio b (about 3 usually) to doubling of mass. Both physical models are developed approximately and then more precisely. The approximate model (2) agrees with results of Safronov (1972). Both precise models are applied to populations approximating those generated by the calculations of Greenberg et al. (1978), Cox and Lewis (1979), and Wetherill (1979).

Kaula, W. M.↗
NASA NTRS2001

Diagnosing Tropical Divergent Circulations From Satellite-Derived Diabatic Forcing and Variational Integral Constraints

Large-scale divergent circulations are part of the atmospheric dynamic response to diabatic heating from condensation, radiative processes, and surface energy fluxes. Vertical motion and the associated divergent wind is thus intimately tied to the hydrologic cycle and the global beat balance. Despite its importance, the divergent circulation is too small in comparison to the rotational flow to measure directly with any accuracy. Vertical motions are recovered diagnostically from reanalyses and, as such, are subject to shortcomings in model physics, numerics, and data availability. While reanalysis estimates of tropical divergent circulations are much improved over those from the Global Weather Experiment, there are still substantial differences between products from operational centers. This is because these circulations are still forced largely by model physics and only secondarily by observations. In order to produce a refined estimate of tropical divergence and its interannual variability we have used a number of remotely-sensed data sets along with variational constraints to improve upon reanalysis estimates. Among these are: precipitation from SSM/I and GPCP, TOVS Path-A vertical cloud distributions; ISCCP radiative cooling rates; TOA radiative fluxes from ERBS, surface radiative fluxes from the SRB project, and surface latent and sensible flux estimates from SSM/I. The TOVS Path-A data constrain the divergent outflow in precipitating regions to have the same vertical structure as observed cloudiness. Using integral constraints for moisture, heat, and mass balance, we retrieve consistent divergent wind flows. We examine the ability of this type analysis to capture regional details of ENSO related perturbations to the divergent wind and associated tropical energy balance. Precipitation from satellite is found to be the major constraint in supplying this horizontal structure. We also consider the ability of this analysis to quantify integrated (land vs. ocean) fluctuations in the global monsoonal flow and energy balance. This depends, of course, on the error characteristics of the data constraints. Some estimate of these properties and their relative importance is provided.

Robertson, Franklin R.↗
DOE OSTI2020

Wind ramp events validation in NWP forecast models during the second Wind Forecast Improvement Project (WFIP2) using the Ramp Tool and Metric (RT&M)

The second Wind Forecast Improvement Project (WFIP2) is a multi-agency field campaign held in the Columbia Gorge area (October 2015 - March 2017). The main goal of the project is to understand and improve the forecast skill of numerical weather prediction (NWP) models in complex terrain, particularly beneficial for the wind energy industry. This region is well-known for its excellent wind resource. One of the biggest challenges for wind power production is the accurate forecasting of wind ramp events (large changes of generated power over short periods of time). Poor forecasting of the ramps requires large and sudden adjustments in conventional power generation, ultimately increasing the costs of power. A Ramp Tool and Metric (RT&M) was developed during the first WFIP experiment, held in the U.S. Great Plains (September 2011 - August 2012). The RT&M was designed to explicitly measure the skill of NWP models at forecasting wind ramp events. Here we apply the RT&M to 80-m (turbine hub-height) wind speeds measured by 19 sodars and 3 lidars, and to forecasts from the High Resolution Rapid Refresh (HRRR), 3-km, and from the High Resolution Rapid Refresh Nest (HRRRNEST), 750-m horizontal grid spacing, models. The diurnal and seasonal distribution of ramp events are analyzed, finding a noticeable diurnal variability for spring and summer but less for fall and especially winter. Also, winter has fewer ramps compared to the other seasons. The model skill at forecasting ramp events, including the impact of the modification to the model physical parameterizations, was finally investigated.

Djalalova, Irina V.↗