Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation 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 1,081 records · Page 60

Theoretical detection threshold of the proton-acoustic range verification technique

Purpose: Range verification in proton therapy using the proton-acoustic signal induced in the Bragg peak was investigated for typical clinical scenarios. The signal generation and detection processes were simulated in order to determine the signal-to-noise limits. Methods: An analytical model was used to calculate the dose distribution and local pressure rise (per proton) for beams of different energy (100 and 160 MeV) and spot widths (1, 5, and 10 mm) in a water phantom. In this method, the acoustic waves propagating from the Bragg peak were generated by the general 3D pressure wave equation implemented using a finite element method. Various beam pulse widths (0.1–10 μs) were simulated by convolving the acoustic waves with Gaussian kernels. A realistic PZT ultrasound transducer (5 cm diameter) was simulated with a Butterworth bandpass filter with consideration of random noise based on a model of thermal noise in the transducer. The signal-to-noise ratio on a per-proton basis was calculated, determining the minimum number of protons required to generate a detectable pulse. The maximum spatial resolution of the proton-acoustic imaging modality was also estimated from the signal spectrum. Results: The calculated noise in the transducer was 12–28 mPa, depending on the transducer central frequency (70–380 kHz). The minimum number of protons detectable by the technique was on the order of 3–30 × 10{sup 6} per pulse, with 30–800 mGy dose per pulse at the Bragg peak. Wider pulses produced signal with lower acoustic frequencies, with 10 μs pulses producing signals with frequency less than 100 kHz. Conclusions: The proton-acoustic process was simulated using a realistic model and the minimal detection limit was established for proton-acoustic range validation. These limits correspond to a best case scenario with a single large detector with no losses and detector thermal noise as the sensitivity limiting factor. Our study indicated practical proton-acoustic range verification may be feasible with approximately 5 × 10{sup 6} protons/pulse and beam current.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Large Scale GPU Accelerated PPMLR-MHD Simulations for Space Weather Forecast

PPMLR-MHD is a new magnetohydrodynamics (MHD) model used to simulate the interactions of the solar wind with the magnetosphere, which has been proved to be the key element of the space weather cause-and-effect chain process from the Sun to Earth. Compared to existing MHD methods, PPMLR-MHD achieves the advantage of high order spatial accuracy and low numerical dissipation. However, the accuracy comes at a cost. On one hand, this method requires more intensive computation. On the other hand, more boundary data is subject to be transferred during the process of simulation. In this work, we present a parallel hybrid solution of the PPMLR-MHD model implemented using the computing capabilities of both CPUs and GPUs. We demonstrate that our optimized implementation alleviates the data transfer overhead by using GPU Direct technology and can scale up to 151 processes and achieve significant performance gains by distributing the workload among the CPUs and GPUs on Titan at Oak Ridge National Laboratory. The performance results show that our implementation is fast enough to carry out highly accurate MHD simulations in real time.

Computer Science↗

Simulation of Gas Production from Multilayered Hydrate-Bearing Media with Fully Coupled Flow, Thermal, Chemical and Geomechanical Processes Using TOUGH + Millstone. Part 1: Numerical Modeling of Hydrates

TOUGH + Millstone has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH4-hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH + HYDRATE simulator, V2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the behavior of hydrate-bearing systems undergoing changes and includes the most recent advances in the description of the system properties, coupled seamlessly with (b) Millstone V1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of TOUGH + Millstone are demonstrated in the simulation and analysis of the system flow, thermal and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the first paper of this series, we discuss the physics underlying the T + H hydrate simulator, the constitutive relationships describing the physical, chemical (equilibrium and kinetic) and thermal processes, the states of the CH 4 + H 2 O system and the sources of critically important data, as well as the mathematical approaches used for the development of the of mass and energy balance equations and their solution. Finally, we provide verification examples of the hydrate code against numerical results from the simulation of laboratory and field experiments.

58 GEOSCIENCES↗

Watershed response of the Feather River Basin, California, United States of America to future climate changes

Increasing mean annual temperatures under climate change are expected to reduce seasonal snowpack, increase evapotranspiration (ET), and alter summer baseflow in headwater watersheds worldwide. Strong regional variability in hydrologic responses highlights the need for catchment-scale, physically based models to assess future flood risk and water availability. This study examines climate-driven changes in the hydrologic response of the Upper Feather River watershed in the Sierra Nevada Mountains, California. Four priority climate models and two representative concentration pathways (RCP4.5 and RCP8.5) are used to evaluate future hydrologic responses of the watershed. Hydrologic processes are simulated using an objectively calibrated Soil and Water Assessment Tool Plus model for a historical baseline (1986–2005) and a future period 2070–2099), driven by observed and projected precipitation and temperature. Results indicate a declining contribution of snowfall to annual precipitation, with peak snowfall and water yield shifting 1–3 months earlier. Long-term annual maximum flows are projected to increase considerably, whereas low-flow responses are mixed, with both increases and decreases projected by the end of the century. These findings highlight the need for adaptive watershed management to enhance flood protection, water storage, and drought resilience. Future water resource planning should also account for one-to-three-month shifts in peak water yield and surface runoff due to changes in snowmelt timing and a lower snowfall-to-rainfall ratio under climate change. The modeling framework and insights are transferable to other snow-dominated headwater watersheds experiencing climate-driven change.

Tigabu, T [UC Davis]↗

Investigating added value of regional climate modeling in North American winter storm track simulations

Extratropical Cyclone (EC) characteristics depend on a combination of large-scale factors and regional processes. However, the latter are considered to be poorly represented in global climate models (GCMs), partly because their resolution is too coarse. This work describes a framework using possibilities given by regional climate models (RCMs) to gain insight into storm activity during winter over North America (NA). Recent past climate period (1981–2005) is considered to assess EC activity over NA using the NCEP regional reanalysis (NARR) as a reference, along with the European reanalysis ERA-Interim (ERAI) and two CMIP5 GCMs used to drive the Canadian Regional Climate Model—version 5 (CRCM5) and the corresponding regional-scale simulations. While ERAI and GCM simulations show basic agreement with NARR in terms of climatological storm track patterns, detailed bias analyses show that, on the one hand, ERAI presents statistically significant positive biases in terms of EC genesis and therefore occurrence while capturing their intensity fairly well. On the other hand, GCMs present large negative intensity biases in the overall NA domain and particularly over NA eastern coast. In addition, storm occurrence over the northwestern topographic regions is highly overestimated. When the CRCM5 is driven by ERAI, no significant skill deterioration arises and, more importantly, all storm characteristics near areas with marked relief and over regions with large water masses are significantly improved with respect to ERAI. Conversely, in GCMdriven simulations, the added value contributed by CRCM5 is less prominent and systematic, except over western NA areas with high topography and over the Western Atlantic coastlines where the most frequent and intense ECs are located. Despite this significant added-value on seasonalmean characteristics, a caveat is raised on the RCM ability to handle storm temporal ‘seriality’, as a measure of their temporal variability at a given location. In fact, the driving models induce some significant footprints on the RCM skill to reproduce the intra-seasonal pattern of storm activity.

54 ENVIRONMENTAL SCIENCES↗

Modeling and Simulation of Porosity in Spray Deposition

Numerical models, consisting of interstitial porosity model and gas porosity model, have been developed to predict the porosity in the spray deposition process. The interstitial porosity model is based on the assumption of liquid feeding in the porous mushy layer of the deposit. The porosity fraction is related to the particle packing density and the amount of the feeding flow during solidification. The gas porosity model is based on the diffusion kinetics of gas pore formation. A numerical simulation of spray deposition for a billet was carried out by coupling the porosity model with the shape and heat-transfer models. The volume fraction and distribution of the porosity were investigated. As a result of the experimental verification, the simulation results showed good agreement with the experimental data. As applications, the models were used to investigate the influence of spray process parameters on porosity formation. The investigation demonstrated that the porosity is substantially influenced by the temperature of the impinging layer and the preheating of the substrate.

36 MATERIALS SCIENCE↗

Pilot Testing of a Modular Oxygen Production System Using Oxygen Binding Adsorbents (Final Report)

In this project, RTI and Air Liquide have focused on the development of innovative oxygen separation materials and technologies based on reversibly binding of oxygen to enable smaller and cheaper air separations. In alignment with that object the team has had multiple achievements to advance the technology including: Synthesized and characterized novel M-CoSalen (RTIO2Sorb-1) material which showed a dynamic oxygen adsorption capacity in excess of 1wt%. Developed synthesis routes for RTIO2Sorb-1 with commercially relevant techniques and produced batches of 0.25kg and over 4 kg of total synthesis. Synthesized and characterized extrudate forms of the RTIO2Sorb-1 with relevant crush strength and maintained dynamic oxygen adsorption capacity in near 1wt%. Developed techniques for forming structured beds in fiber shapes with conventional materials. Developed a model for O 2 sorption processes to analyze various bed configurations and cycle parameters. Studied classic O 2 solid sorbents reported in the literature to learn oxygen binding mechanisms. Design a 10 kg/day O 2 VPSA pilot system with 2 or 4-bed operation. Completed a techno-economic analysis based on the experimental results and modeling. The overall objective of the project was to design, fabricate, and test a modular O 2 production system and perform a techno-economic analysis (TEA) after testing to determine the cost-benefit of the advanced modular air separation system. The goals of this technology development project were to achieve a bed-size factor (BSF) of less than 600 lb-adsorbent/TPD O 2 (ton/day O 2 ) (as compared with the state-of-the-art BSF of 850), and O 2 -purity greater than 95% at a cost that is projected to be equivalent or lower than the current state of the art (SOTA), commercially available large-scale cryogenic air separation systems. To achieve these goals and objective, the project team executed on (1) oxygen binding adsorbent optimization and scale up, (2) adsorbent material formation process studies to form the adsorbent material into structured beds for rapid pressure swing adsorption (PSA) cycles with low pressure drop, fast mass transfer, and low attrition, (3) cycle development studies to optimize the PSA process, and (4) develop simulation tools for rapid cycle modeling and numerical evaluation/optimization.

20 FOSSIL-FUELED POWER PLANTS↗

RAY-RAMSES: a code for ray tracing on the fly in N-body simulations

We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

North Atlantic Oscillation response in GeoMIP experiments G6solar and G6sulfur: why detailed modelling is needed for understanding regional implications of solar radiation management

The realization of the difficulty of limiting global-mean temperatures to within 1.5 or 2.0°C above pre-industrial levels stipulated by the 21st Conference of Parties in Paris has led to increased interest in solar radiation management (SRM) techniques. Proposed SRM schemes aim to increase planetary albedo to reflect more sunlight back to space and induce a cooling that acts to partially offset global warming. Under the auspices of the Geoengineering Model Intercomparison Project, we have performed model experiments whereby global temperature under the high-forcing SSP5-8.5 scenario is reduced to follow that of the medium-forcing SSP2-4.5 scenario. Two different mechanisms to achieve this are employed: the first via a reduction in the solar constant (experiment G6solar) and the second via modelling injections of sulfur dioxide (experiment G6sulfur) which forms sulfate aerosol in the stratosphere. Results from two state-of-the-art coupled Earth system models (UKESM1 and CESM2-WACCM6) both show an impact on the North Atlantic Oscillation (NAO) in G6sulfur but not in G6solar. Both models show a persistent positive anomaly in the NAO during the Northern Hemisphere winter season in G6sulfur, suggesting an increase in zonal flow and an increase in North Atlantic storm track activity impacting the Eurasian continent and leading to high-latitude warming over Europe and Asia. These results are broadly consistent with previous findings which show similar impacts from stratospheric volcanic aerosol on the NAO and emphasize that detailed modelling of geoengineering processes is required if accurate impacts of SRM effects are to be simulated. Differences remain between the two models in predicting regional changes over the continental USA and Africa, suggesting that more models need to perform such simulations before attempting to draw any conclusions regarding potential continental-scale climate change under SRM.

36 MATERIALS SCIENCE↗

Integrated model construction for CO2-EOR monitoring via charged-wellbore casing controlled-source electromagnetics

We present a method for the integrated development of 3D electrical conductivity models for CO2-EOR monitoring with charged wellbore casing controlled-source electromagnetics (CWC-CSEM). The model is constructed through the multiphysics integration of reservoir simulation, seismic, borehole resistivity, and transient EM (TEM) data. The process is performed in two stages. First, a large background conductivity model is constructed from the combination of seismic, borehole, and TEM data. Second, a more detailed and dynamic component of the model is created within the production interval by converting the reservoir simulation parameters to conductivity through Archie’s equation. In this presentation, we demonstrate the complete workflow for building these 3D conductivity site models that can be updated throughout production using the Bell Creek oil field as an example. We then show application of the developed site model to simulating the predicted CWC-CSEM data as one step in the larger problem of reservoir imaging and monitoring of injected CO2 during enhanced oil recovery.

CO2, Sequestration, Monitoring, time-lapse, reserv↗

Fuel Fabrication and Single Stage Aqueous Process Modeling

The Material Protection, Accounting, and Control Technologies program utilizes modeling and simulation to assess Material Control and Accountability (MC&A) concerns for a variety of nuclear facilities. Single analyst tools allow for rapid design and evaluation of advanced approaches for new and existing nuclear facilities. A low enriched uranium (LEU) fuel conversion and fabrication facility simulator is developed to assist with MC&A for existing facilities. Measurements are added to the model (consistent with current best practices). Material balance calculations and statistical tests are also added to the model. In addition, scoping work is performed for developing a single stage aqueous reprocessing model. Preliminary results are presented and discussed, and next steps outlined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microstructure based modeling of friction stir welded joint of dissimilar metals using crystal plasticity

We report the friction stir welding (FSW) process shows promising results in joining dissimilar metals which are otherwise almost impossible to join using traditional welding techniques. Being a new technique, the deformation and the failure mechanism of the joints made by the FSW process needs to be investigated. In this work, a joint between AZ31 Mg alloy and DP590 steel is modeled using phenomenological crystal plasticity formulation on the mesoscale in the form of a representative volume element (RVE). The interface of the two materials is modeled using a cohesive zone model. A parametric study has been performed to understand the effect of grain size and interface fracture toughness as well as strength on the mechanical performance of the joint. It was found that the grain size of AZ31 Mg alloy, as well as DP590 steel, has little effect on the overall joint performance. On the other hand, interface fracture toughness and strength have a significant impact on the mechanical properties of the joint

36 MATERIALS SCIENCE↗

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗

Heat Pump Water Heater Using Wrapped-Tank Microchannel Condenser Coil for Charge Reduction – FY19 1st Quarter (Regular) Milestone Report: Development of Wrapped-Tank Condenser Model (Round Tube and Microchannel), Coupled with Water Tank Model

We recently grew the DOE/HPDM model library to enhance its heat pump water heater (HPWH) design and simulation capability. The new additions include wrapped-tank round tube coil, microchannel coil and a one-dimensional stratified water tank. A wrapped-tank condenser coil has strong coupling with a stratified water tank, which leads HPWH simulation a transient process. The new one-dimensional stratified water tank model is an improvement to the open-source EnergyPlus water tank model, by introducing a calibration factor to account for bulk mixing effect due to water draws and circulations, etc. The wrapped-tank condenser coil model, using a segment-to-segment modeling approach to simulate round tubes or microchannel ports. The heat exchanger models can accept arbitrary refrigerant types, e.g. HFCs, HFOs, natural refrigerants like CO 2 and propane, etc.

42 ENGINEERING↗

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↗

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Abstract. Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (tens of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) – Daily Surface Weather and Climatological Summaries (Daymet), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30 % of the watershed area). However, the 4 km-resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using the temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), and sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications for model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

Correlated prompt fission data in transport simulations

Detailed information on the fission process can be inferred from the observation, modeling and theoretical understanding of prompt fission neutron and γ-ray observables. Beyond simple average quantities, the study of distributions and correlations in prompt data, e.g., multiplicity-dependent neutron and γ-ray spectra, angular distributions of the emitted particles, n -n, n - γ, and γ - γ correlations, can place stringent constraints on fission models and parameters that would otherwise be free to be tuned separately to represent individual fission observables. The FREYA and CGMF codes have been developed to follow the sequential emissions of prompt neutrons and γ rays from the initial excited fission fragments produced right after scission. Both codes implement Monte Carlo techniques to sample initial fission fragment configurations in mass, charge and kinetic energy and sample probabilities of neutron and γ emission at each stage of the decay. This approach naturally leads to using simple but powerful statistical techniques to infer distributions and correlations among many observables and model parameters. The comparison of model calculations with experimental data provides a rich arena for testing various nuclear physics models such as those related to the nuclear structure and level densities of neutron-rich nuclei, the γ-ray strength functions of dipole and quadrupole transitions, the mechanism for dividing the excitation energy between the two nascent fragments near scission, and the mechanisms behind the production of angular momentum in the fragments, etc. Beyond the obvious interest from a fundamental physics point of view, such studies are also important for addressing data needs in various nuclear applications. The inclusion of the FREYA and CGMF codes into the MCNP6.2 and MCNPX - PoliMi transport codes, for instance, provides a new and powerful tool to simulate correlated fission events in neutron transport calculations important in nonproliferation, safeguards, nuclear energy, and defense programs. Here, this review provides an overview of the topic, starting from theoretical considerations of the fission process, with a focus on correlated signatures. It then explores the status of experimental correlated fission data and current efforts to address some of the known shortcomings. Numerical simulations employing the FREYA and CGMF codes are compared to experimental data for a wide range of correlated fission quantities. The inclusion of those codes into the MCNP6.2 and MCNPX - PoliMi transport codes is described and discussed in the context of relevant applications. The accuracy of the model predictions and their sensitivity to model assumptions and input parameters are discussed. Lastly, a series of important experimental and theoretical questions that remain unanswered are presented, suggesting a renewed effort to address these shortcomings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗