SEARCH · Engineering Papers

Results for “Modeling and simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

DOE OSTI2016

HIGH-FIDELITY SIMULATION-DRIVEN MODEL DEVELOPMENT FOR COARSE-GRAINED COMPUTATIONAL FLUID DYNAMICS

Nuclear reactor safety analysis requires identifying various credible accident scenarios and determining their consequences. For a full-scale nuclear power plant system behavior, it is impossible to obtain sufficient experimental data for a broad range of risk-significant accident scenarios. In single-phase flow convective problems, Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) can provide us with high fidelity results when physical data are unavailable. However, these methods are computationally expensive and cannot be afforded for simulation of long transient scenarios in nuclear accidents despite extraordinary advances in high performance scientific computing over the past decades. The major issue is the inability to make the transient computation parallel, thus making number of time steps required in high-fidelity methods unaffordable for long transients. In this work, we propose to apply a high fidelity simulation-driven approach to model sub-grid scale (SGS) effect in Coarse Grained Computational Fluid Dynamics CG-CFD. This approach aims to develop a statistical surrogate model instead of the deterministic SGS model. We chose to start with a turbulent natural convection case with volumetric heating in a horizontal fluid layer with a rigid, insulated lower boundary and isothermal (cold) upper boundary. This scenario of unstable stratification is relevant to turbulent natural convection in a molten corium pool during a severe nuclear reactor accident, as well as in containment mixing and passive cooling. The presented approach demonstrates how to create a correction for the CG-CFD solution by modifying the energy balance equation. A global correction for the temperature equation proves to achieve a significant improvement to the prediction of steady state temperature distribution through the fluid layer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1994

Ice-atmosphere interactions in the central Arctic: Remotely-sensed and simulated ice concentration and motion

The response of the Beaufort Sea ice pack to the passage of strong low pressure systems is studied using SAR (synthetic aperture radar), AVHRR (advanced very high resolution radiometer), and SSM/I (special sensor microwave/imager) data combined with model simulations. The SAR, AVHRR, and modeled concentrations concur generally in showing a 1 to 5 percent decrease in ice fraction during the passage of the lows. SSM/I derived concentrations appear to underestimate overall concentration and overestimate the opening of the pack relative to the other data sets. Ice motion derived from SAR and AVHRR and model simulations using two ice rheologies agree reasonably well in direction and magnitude, although the cavitating fluid approximation tends to overestimate ice motion compared to a viscous plastic rheology. Ice motions from SAR, AVHRR, and buoy observations can be merged to yield uniform grids particularly well suited to synoptic scale studies and model validation.

Maslanik, J. A.↗

NASA NTRS1985

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗

NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Adverse Pressure Gradients due to Concave Streamline Curvature

Direct numerical simulations (DNS) of adverse-pressure-gradient turbulent boundary layers over a planar concave wall are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar concave wall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University (TAMU). The DNS was validated against the experimental results of TAMU for the same flow conditions and wall geometry. An analysis of the DNS datasets was also conducted to provide an assessment of the validity of Morkovin’s hypothesis and the strong Reynolds analog for turbulence subject to mechanical nonequilibrium. In addition to the DNS results, RANS predictions are obtained by using the Baldwin-Lomax (BL), Spalart-Allmaras (SA), and the k - w SST turbulence models. The comparisons between RANS and DNS showed little impact of an adverse pressure gradient on the accuracy of these models, at least up to an incompressible Clauser pressure gradient parameter of beta(sub inc) 1.22. While the Boussinesq assumption provided reasonable predictions for the Reynolds shear stress, it failed to adequately predict the normal components of the Reynolds stress.

turbulent boundary layers↗

DOE OSTI2018

On the Linearity of Local and Regional Temperature Changes from 1.5°C to 2°C of Global Warming

Because of the Paris Agreement it is imperative there is greater understanding of the consequences of limiting global warming to the target 1.5° and 2°C levels above preindustrial conditions. It is challenging to quantify changes across a small increment of global warming, so a pattern-scaling approach may be considered. Here we investigate the validity of such an approach by comprehensively examining how well local temperatures and warming trends in a 1.5°C world predict local temperatures at global warming of 2°C. Ensembles of transient coupled climate simulations from multiple models under different scenarios were compared and individual model responses were analyzed. For many places, the multimodel forced response of seasonal-average temperatures is approximately linear with global warming between 1.5° and 2°C. Yet, individual model results vary and large contributions from nonlinear changes in unforced variability or the forced response cannot be ruled out. In some regions, such as East Asia, models simulate substantially greater warming than is expected from linear scaling. Examining East Asia during boreal summer, we find that increased warming in the simulated 2°C world relative to scaling up from 1.5°C is related to reduced anthropogenic aerosol emissions. Our results indicate that, where forcings other than those due to greenhouse gas emissions change, the warming experienced in a 1.5°C world is a poor predictor for local climate at 2°C of global warming. In addition to the analysis of the linearity in the forced climate change signal, we find that natural variability remains a substantial contribution to uncertainty at these low-warming targets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗

NASA NTRS1991

Characterization of and sensor-model verification and control of the properties of PMR-13 during processing

This study presents an in situ sensor and a process-simulation model developed for monitoring and controlling the cure process of PMR-15. The time-temperature dependence of the buildup in the molecular network structure, extent of reaction, viscosity, flow, and consolidation during the cure of PMR-15 are discussed. The relationship of the time-temperature cycle used during imidization to the length of the endcapped chains formed is examined. The relationship of the time-temperature cure-processing cycle to the reaction kinetics, viscosity, flow, and consolidation during crosslinking is analyzed using frequency-dependent electromagnetic sensors and the Loos processing model. Application of the FDEMS sensing technique and the process-simulation model for quality assurance processing and automated on-line control of cure is discussed.

Kranbuehl, D.↗

DOE OSTI2016

Multiscale Model for Simulating Fission-Induced Recrystallization in U-Mo Alloy

U-Mo alloy is one of the most promising fuels for the future high performance research and test reactors due to its high uranium loading density and stable irradiation behavior. However, a major concern to its qualification is the fuel swelling at high fission density due to the formation of large fission-gas bubbles. The fission-induced recrystallization further expedites the formation and growth of gas bubbles. So far, the exact mechanism of recrystallization is still unclear. One possible driving force for recrystallization is the production of point defects and interstitial loops due to irradiation. Following this, we developed a multi-scale model by integrating the rate-theory and phase-field models to study microstructural evolutions in irradiated U-Mo alloy, i.e., fission-induced gas bubble evolution and recrystallization. Density function theory (DFT) calculations are utilized to predict the defect formation energies, diffusivities of defects, the interfacial energy of gas bubble, and the elastic constants of U-Mo alloy, etc. The dislocation density evolution is predicted by the rate-theory developed by Jeff. The grain subdivision occurs once the dislocation density is above a critical value. The grain growth coupled with gas bubbles evolution is modeled using a multi-phase phase-field model. The effects of fission rate and initial grain microstructure on the recrystallization process are systematically investigated. The simulated dislocation density, sub-grain size distribution, and volume fraction of the recrystallization grains under different fission density are comparable with experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age

Newton-Krylov solvers for ocean tracers have the potential to greatly decrease the computational costs of spinning up deep-ocean tracers, which can take several thousand model years to reach equilibrium with surface processes. One version of the algorithm uses offline tracer transport matrices to simulate an annual cycle of tracer concentrations and applies Newton’s method to find concentrations that are periodic in time. Here we present the impact of time-averaging the transport matrices on the equilibrium values of an ideal-age tracer. We compared annually-averaged, monthly-averaged, and 5-day-averaged transport matrices to an online simulation using the ocean component of the Community Earth System Model (CESM) with a nominal horizontal resolution of 1° × 1° and 60 vertical levels. We found that increasing the time resolution of the offline transport model reduced a low age bias from 12% for the annually-averaged transport matrices, to 4% for the monthly-averaged transport matrices, and to less than 2% for the transport matrices constructed from 5-day averages. The largest differences were in areas with strong seasonal changes in the circulation, such as the Northern Indian Ocean. As a result, for many applications the relatively small bias obtained using the offline model makes the offline approach attractive because it uses significantly less computer resources and is simpler to set up and run.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Engineering challenges for detectors at the ILC

Over the last years two proposals for experiments at the ILC have been developed, ILD and SID. Extensive R&D has been carried out around the world to develop the needed technologies. Furthermore a first round of engineering studies was made as part of the ILC TDR to understand the integration of these different sub-systems into coherent and integrated detector concepts. Among the key challenges for the sub detectors are the extreme low mass/low power requirements or the extreme channel densities needed in particle flow based detectors. Throughout these studies special care was taken to ensure that the engineering models and the simulation models, used in studies of the physics capabilities of the detectors, stay synchronized. In the near future, the models will need to be evolved to take the special requirements of the potential ILC site in Japan into account. Furthermore, the state of the integration of the detectors, and the future directions, will be discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Observations Indicate That Clouds Amplify Mechanisms of Southern Ocean Heat Uptake

Abstract The Southern Ocean has absorbed most of the excess heat associated with anthropogenic greenhouse gas emissions. Since Southern Ocean observations are sparse in certain regions and seasons, much of our knowledge of ocean heat uptake is based on climate model simulations. However, climate models still inadequately represent some properties of Southern Ocean clouds, and they have not identified the mechanisms by which clouds may affect Southern Ocean heat uptake (SOHU). Here, we use the ERA5 and JRA‐55 reanalyses to assess the influence of clouds and other atmospheric processes on SOHU from 1979 to 2020. We find that years with the highest SOHU between 45° and 65°S are dominated by ocean heat uptake anomalies during winter and spring, but not during summer or fall. Winter and spring cloud cover are up to 7% higher when SOHU is up to 5.5 W/m 2 higher than the climatological seasonal mean, with the largest increases in the South Pacific Ocean. Clouds also contain more liquid water. These changes in cloud properties increase downwelling longwave radiation, amplifying ocean heat uptake. Cloud changes are also concomitant with a more stable lower atmosphere, which suppresses turbulent heat fluxes out of the surface. Overall, we find that SOHU is likely not mediated by enhanced surface shortwave absorption over the observational time period. A better understanding of how atmospheric processes impact ocean heat uptake may help improve our understanding of ocean heat uptake mechanisms in the current generation of climate models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Increasing Hurricane Intensification Rate Near the US Atlantic Coast

Hurricanes often cause severe damage and loss of life, and storms that intensify close to the coast pose a particularly serious threat. While changes in hurricane intensification and environment have been examined at basin scales previously, near-coastal changes have not been adequately explored. In this study, we address this using a suite of observations and climate model simulations. Over the 40-year period of 1979–2018, the mean 24-hr hurricane intensification rate increased by ~1.2 kt 6-hr –1 near the US Atlantic coast. However, a significant increase in intensification did not occur near the Gulf coast over the same period. The enhanced hurricane intensification along the Atlantic coast is consistent with an increasingly favorable dynamic and thermodynamic environment there, which is well simulated by climate models over the historical period. Further, multi-model projections suggest a continued enhancement of the storm environment and hurricane intensification near the Atlantic coast in the future.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Thermal conductivity measurements of proton-heated warm dense aluminum

Thermal conductivity is one of the most crucial physical properties of matter when it comes to understanding heat transport, hydrodynamic evolution, and energy balance in systems ranging from astrophysical objects to fusion plasmas. In the warm dense matter regime, experimental data are very scarce so that many theoretical models remain untested. Here we present the first thermal conductivity measurements of aluminum at 0.5–2.7 g/cc and 2–10 eV, using a recently developed platform of differential heating. A temperature gradient is induced in a Au/Al dual-layer target by proton heating, and subsequent heat flow from the hotter Au to the Al rear surface is detected by two simultaneous time-resolved diagnostics. A systematic data set allows for constraining both thermal conductivity and equation-of-state models. Simulations using Purgatorio model or Sesame S27314 for Al thermal conductivity and LEOS for Au/Al release equation-of-state show good agreement with data after 15 ps. Discrepancy still exists at early time 0–15 ps, likely due to non-equilibrium conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Evolution of flux ropes in the magnetotail: A three-dimensional global hybrid simulation

Flux ropes in the Earth's magnetotail are widely believed to play a crucial role in energy transport during substorms and the generation of energetic particles. Previous kinetic simulations are limited to the local-scale regime, and thus cannot be used to study the structure associated with the geomagnetic field and the global-scale evolution of the flux ropes. Here, the evolution of flux ropes in the magnetotail under a steady southward interplanetary magnetic field are studied with a newly developed three-dimensional global hybrid simulation model for dynamics ranging from the ion Larmor radius to the global convection time scales. Magnetic reconnection with multiple X-lines is found to take place in the near-tail current sheet at geocentric solar magnetospheric distances x=−30R{sub E}∼−15R{sub E} around the equatorial plane (z=0). The magnetotail reconnection layer is turbulent, with a nonuniform structure and unsteady evolution, and exhibits properties of typical collisionless fast reconnection with the Hall effect. A number of small-scale flux ropes are generated through the multiple X-line reconnection. The diameter of the flux ropes is several R{sub E}, and the spatial scale of the flux ropes in the dawn-dusk direction is on the order of several R{sub E} and does not extend across the entire section of the magnetotail, contrary to previous models and MHD simulation results and showing the importance of the three-dimensional effects. The nonuniform and unsteady multiple X-line reconnection with particle kinetic effects leads to various kinds of flux rope evolution: The small-scale flux ropes propagate earthward or tailward after formation, and eventually merge into the near-Earth region or the mid-/distant-tail plasmoid, respectively. During the propagation, some of the flux ropes can be tilted in the geocentric solar magnetospheric (x,y) plane with respect to the y (dawn-dusk) axis. Coalescence between flux ropes is also observed. At the same time, the evolution of the flux ropes in the multiple X-line reconnection layer can also lead to the acceleration and heating of ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Light-Duty Vehicle Attribute Projections (Years 2015-2030)

This report describes the approach and results associated with projecting vehicle attributes for future light-duty vehicle classes and powertrains, from year 2015 to 2030. The focus is the California light-duty vehicle market and the projected attributes include fuel economy, acceleration, driving range (total and all-electric), and manufacturer suggested retail price from year 2015 to year 2030. Attributes are weighted by California vehicle sales, which are projected using a historically validated consumer choice model - the Automotive Deployment Option Projection Tool (ADOPT) - integrated with a modeling framework that simulates vehicle fuel economy, cost, and acceleration performance while optimizing vehicle components - the Future Automotive Systems Technology Simulator model (FASTSim). Both models were developed at the National Renewable Energy Laboratory and have been adapted to represent the future California light-duty vehicle market for the present report. The analysis includes several scenarios pertaining to electricity and, essentially, electric vehicle demand in California, as established by the Energy Commission. Results suggest that implementation of certain policies, such as the Corporate Average Fuel Economy standards, impacts vehicle attribute projections. The results also show that standards and policy targets are not exclusively met by changes in vehicle attributes, but also through shifts in market demand and sales for certain vehicle powertrains.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2019

Physics Models in Detector Simulation [PowerPoint]

This talk provides an overview of the detector simulation physics models needs of HEP experiments and is restricted to LHC experiments (ALICE, ATLAS, CMS, LHCb) and Belle II and will be focused on the Geant4 toolkit. Other experiments are discussed in other talks. This talk also discuss differences and commonalities among experiments and will identify opportunities for collaboration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Tidal Resource Gaps Analysis Technical Report

The tidal resource gaps project was created to address a growing body of evidence that models underpredict tidal current speeds compared to measurements at a number of the top-ranking tidal energy sites. In response, this project compared opportunistic tidal power density measurements from 16 tidal energy hot spot sites with estimates from resource assessment data to identify discrepancies. To improve the accuracy of resource estimates from model data, updated data from eight improved model simulations were obtained. Model improvements included grid refinement, domain coupling, and the use of unstructured or nested grids. New resource estimates were calculated from the updated models, and these data were used to update the tidal hot spots-a list of promising tidal energy sites around the United States.

16 TIDAL AND WAVE POWER↗

NASA NTRS1974

A study of two statistical methods as applied to shuttle solid rocket booster expenditures

The state probability technique and the Monte Carlo technique are applied to finding shuttle solid rocket booster expenditure statistics. For a given attrition rate per launch, the probable number of boosters needed for a given mission of 440 launches is calculated. Several cases are considered, including the elimination of the booster after a maximum of 20 consecutive launches. Also considered is the case where the booster is composed of replaceable components with independent attrition rates. A simple cost analysis is carried out to indicate the number of boosters to build initially, depending on booster costs. Two statistical methods were applied in the analysis: (1) state probability method which consists of defining an appropriate state space for the outcome of the random trials, and (2) model simulation method or the Monte Carlo technique. It was found that the model simulation method was easier to formulate while the state probability method required less computing time and was more accurate.

Perlmutter, M.↗