Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Long-range prediction and the stratosphere

Abstract. Over recent years there have been concomitant advances in the development of stratosphere-resolving numerical models, our understanding of stratosphere–troposphere interaction, and the extension of long-range forecasts to explicitly include the stratosphere. These advances are now allowing for new and improved capability in long-range prediction. We present an overview of this development and show how the inclusion of the stratosphere in forecast systems aids monthly, seasonal, and annual-to-decadal climate predictions and multidecadal projections. We end with an outlook towards the future and identify areas of improvement that could further benefit these rapidly evolving predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Evaluation of ARM tethered-balloon system instrumentation for supercooled liquid water and distributed temperature sensing in mixed-phase Arctic clouds

A tethered-balloon system (TBS) has been developed and is being operated by Sandia National Laboratories (SNL) on behalf of the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) User Facility in order to collect in situ atmospheric measurements within mixed-phase Arctic clouds. Periodic tethered-balloon flights have been conducted since 2015 within restricted airspace at ARM's Advanced Mobile Facility 3 (AMF3) in Oliktok Point, Alaska, as part of the AALCO (Aerial Assessment of Liquid in Clouds at Oliktok), ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems), and POPEYE (Profiling at Oliktok Point to Enhance YOPP Experiments) field campaigns. The tethered-balloon system uses helium-filled 34 m 3 helikites and 79 and 104 m 3 aerostats to suspend instrumentation that is used to measure aerosol particle size distributions, temperature, horizontal wind, pressure, relative humidity, turbulence, and cloud particle properties and to calibrate ground-based remote sensing instruments. Supercooled liquid water content (SLWC) sondes using the vibrating-wire principle, developed by Anasphere Inc., were operated at Oliktok Point at multiple altitudes on the TBS within mixed-hase clouds for over 200 h. Sonde-collected SLWC data were compared with liquid water content derived from a microwave radiometer, Ka-band ARM zenith radar, and ceilometer at the AMF3, as well as liquid water content derived from AMF3 radiosonde flights. The in situ data collected by the Anasphere sensors were also compared with data collected simultaneously by an alternative SLWC sensor developed at the University of Reading, UK; both vibrating-wire instruments were typically observed to shed their ice quickly upon exiting the cloud or reaching maximum ice loading. Temperature sensing measurements distributed with fiber optic tethered balloons were also compared with AMF3 radiosonde temperature measurements. Combined, the results indicate that TBS-distributed temperature sensing and supercooled liquid water measurements are in reasonably good agreement with remote sensing and radiosonde-based measurements of both properties. From these measurements and sensor evaluations, tethered-balloon flights are shown to offer an effective method of collecting data to inform and constrain numerical models, calibrate and validate remote sensing instruments, and characterize the flight environment of unmanned aircraft, circumventing the difficulties of in-cloud unmanned aircraft flights such as limited flight time and in-flight icing.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

LIGHT-bgcArgo-1.0: using synthetic float capabilities in E3SMv2 to assess spatiotemporal variability in ocean physics and biogeochemistry

Since their advent over 2 decades ago, autonomous Argo floats have revolutionized the field of oceanography, and, more recently, the addition of biogeochemical and biological sensors to these floats has greatly improved our understanding of carbon, nutrient, and oxygen cycling in the ocean. While Argo floats offer unprecedented horizontal, vertical, and temporal coverage of the global ocean, uncertainties remain about whether Argo sampling frequency and density capture the true spatiotemporal variability in physical, biogeochemical, and biological properties. As the true distributions of, e.g., temperature or oxygen are unknown, these uncertainties remain difficult to address with Argo floats alone. Numerical models with synthetic observing systems offer one potential avenue to address these uncertainties. Here, we implement synthetic biogeochemical Argo floats into the Energy Exascale Earth System Model version 2 (E3SMv2), which build on the Lagrangian In Situ Global High-Performance Particle Tracking (LIGHT) module in E3SMv2 (E3SMv2-LIGHT-bgcArgo-1.0). Since the synthetic floats sample the model fields at model run time, the end user defines the sampling protocol ahead of any model simulation, including the number and distribution of synthetic floats to be deployed, their sampling frequency, and the prognostic or diagnostic model fields to be sampled. Using a 6-year proof-of-concept simulation, we illustrate the utility of the synthetic floats in different case studies. In particular, we quantify the impact of (i) sampling density on the float-derived detection of deep-ocean change in temperature or oxygen and on float-derived estimates of phytoplankton phenology, (ii) sampling frequency and sea-ice cover on float trajectory lengths and hence float-derived estimates of current velocities, and (iii) short-term variability in ecosystem stressors on estimates of their seasonal variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding the compound flood risk along the coast of the contiguous United States

Abstract. Compound flooding is a type of flood event caused by multiple flood drivers. The associated risk has usually been assessed using statistics-based analyses or hydrodynamics-based numerical models. This study proposes a compound flood (CF) risk assessment (CFRA) framework for coastal regions in the contiguous United States (CONUS). In this framework, a large-scale river model is coupled with a global ocean reanalysis dataset to (a) evaluate the CF exposure related to the coastal backwater effects on river basins, and (b) generate spatially distributed data for analyzing the CF hazard using a bivariate statistical model of river discharge and storm surge. The two kinds of risk are also combined to achieve a holistic understanding of the continental-scale CF risk. The estimated CF risk shows remarkable inter- and intra-basin variabilities along the CONUS coast with more variabilities in the CF hazard over the US west and Gulf coastal basins. Different risk assessment methods present significantly different patterns in a few key regions such as the San Francisco Bay area, the lower Mississippi River, and Puget Sound. Our results highlight the need to weigh different CF risk measures and avoid using single statistics-based or hydrodynamics-based CFRAs. Uncertainty sources in these CFRAs include the use of gauge observations, which cannot account for the flow physics or resolve the spatial variability of risks, and underestimations of the flood extremes and the dependence of CF drivers in large-scale models, highlighting the importance of understanding the CF risks for developing a more robust CFRA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The collapse of the Cordilleran–Laurentide ice saddle and early opening of the Mackenzie Valley, Northwest Territories, Canada, constrained by 10 Be exposure dating

Deglaciation of the northwestern Laurentide Ice Sheet in the central Mackenzie Valley opened the northern portion of the deglacial Ice-Free Corridor between the Laurentide and Cordilleran ice sheets and a drainage route to the Arctic Ocean. In addition, ice sheet saddle collapse in this section of the Laurentide Ice Sheet has been implicated as a mechanism for delivering substantial freshwater influx into the Arctic Ocean on centennial timescales. However, there is little empirical data to constrain the deglaciation chronology in the central Mackenzie Valley where the northern slopes of the ice saddle were located. Here, we present 30 new 10Be cosmogenic nuclide exposure dates across six sites, including two elevation transects, which constrain the timing and rate of thinning and retreat of the Laurentide Ice Sheet in the area. Our new 10Be dates indicate that the initial deglaciation of the eastern summits of the central Mackenzie Mountains began at ~15.8 ka (17.1–14.6 ka), ~1000 years earlier than in previous reconstructions. The main phase of ice saddle collapse occurred between ~14.9 and 13.6 ka, consistent with numerical modelling simulations, placing this event within the Bølling–Allerød interval (14.6–12.9 ka). Our new dates require a revision of ice margin retreat dynamics, with ice retreating more easterly rather than southward along the Mackenzie Valley. In addition, we quantify a total sea level rise contribution from the Cordilleran–Laurentide ice saddle region of ~11.2 m between 16 and 13 ka.

58 GEOSCIENCES↗
DOE OSTI2025

Subglacial discharge effects on basal melting of a rotating, idealized ice shelf

When subglacial meltwater is discharged into the ocean at the grounding line, it acts as a source of buoyancy, enhancing flow speeds along the ice base that result in higher basal melt rates. The effects of subglacial discharge have been well studied in the context of a Greenland-like, vertical calving front, where Earth's rotation can be neglected. Here we study these effects in the context of Antarctic ice shelves, where rotation is important. We use a numerical model to simulate ocean circulation and basal melting beneath an idealized three-dimensional ice shelf and vary the rate and distribution of subglacial discharge. For channelized discharge, we find that in the rotating case, total melt-flux anomaly increases with two-thirds power of the discharge, in contrast to existing non-rotating results for which the melt-flux anomaly increases with one-third power of the discharge. The higher melt-flux anomaly with discharge is attributed to a more extensive area of the ice-shelf base being exposed to direct high melting by the rising plume as it is deflected due to Earth's rotation and its path is prolonged. For distributed discharge, we find that in both the rotating and the non-rotating cases, the melt-flux anomaly increases with two-thirds power of the discharge. Furthermore, in the rotating case, the addition of channelized, subglacial discharge can produce either a higher or a lower ice-shelf basal melt-flux anomaly than the equivalent amount of distributed discharge, depending on its location along the grounding line relative to the directionality of the Coriolis force. This contrasts with previous results from non-rotating, vertical ice-cliff simulations, where distributed discharge was always found to be more efficient than channelized discharge at enhancing the terminus-averaged melt rate. The implication, based on our idealized simulations, is that melt-rate parameterizations attempting to include subglacial discharge effects that are not geometry- and rotation-aware may produce total melt-flux anomalies that are off by a factor of 2 or more.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗
DOE OSTI2025

Operational wind plants increase planetary boundary layer height: an observational study

As wind energy deployment grows, interactions between wind plants and the surrounding environment become more prevalent. The current investigation seeks to understand these interactions by characterizing the impact of wind plants on the planetary boundary layer height (PBLH), utilizing observations from the American WAKE experimeNt (AWAKEN) campaign. Given the ambiguity of the definition of PBLH under stable atmospheric conditions, where the impact of wind plants is expected to be strongest, a comparison of different methods for identifying PBLH is first conducted using data collected by multiple types of instruments. The Heffter method is selected as the thermodynamic method because it generates the most consistent results for the radiosonde and infrared spectrometer. A minimum vertical velocity variance method is used for a turbulence-based definition. Using both of these methods, the values of PBLH measured at spatially distributed sites are compared under a range of atmospheric conditions. Both methods show a clear increase in PBLH downstream of a wind plant for stable conditions. These impacts are strongest when the upstream PBLH is shallow (less than 0.25 km), with the thermodynamic method showing a PBLH increase of 35 %–39 % and the turbulence-based method showing a 143 % increase. At a site 20 km downstream of the wind plant, these effects are no longer observed, suggesting PBLH is not influenced by the wind plant at this distance. The results of this investigation show that wind plants can modify PBLH in their vicinity. As PBLH is a key parameter for numerical models, this insight into wind plant–atmosphere interactions can be used to validate and improve the accuracy of wind energy forecasts.

17 WIND ENERGY↗
DOE OSTI2023

Data from: Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

This data collection represents all UH8to2 model data used for the publication "Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model." This includes ocean and sea ice characteristics from the 2017-2020 UH8to2 model run. Further details about the model setup, forcing, and our interpretation of these data can be found in the publication.

58 GEOSCIENCES↗
DOE OSTI2021

Study of passively stable, fully detached divertor plasma regimes attained in innovative long-legged divertor configurations

Numerical modeling of divertor configurations with radially or vertically extended, tightly baffled, outer divertor legs has demonstrated the existence of a passively-stable fully detached divertor regime. In the simulations, long-legged divertors provide up to an order-of- magnitude increase in peak power handling capability compared to conventional divertors. The key physics for attaining the passively stable, fully detached regime in these simulations involves the interplay of strong convective plasma transport to the divertor leg outer sidewall, confinement of neutral gas in the divertor volume, geometric effects including a secondary X-point, and atomic radiation. New analysis shows that in this regime the detachment front location is set by the balance between the power entering the divertor leg and the losses to the walls of the divertor channel. Correspondingly, the maximum power that can be accommodated by the divertor, while still staying detached, increases with the poloidal length of the leg. The detached regime access window in terms of input power, density and impurity seeding concentration varies quantitatively depending on divertor geometry and modeling assumptions most specifically, cross-field transport to the side walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toroidal States of the 12 C Nucleus

Among the states of 12 C, there is an important subset of $K$ = 0 and $K$ = I planar intrinsic states in which the intrinsic motion of the nucleons are confined in the planar region defined by the three-alpha cluster or by their generated toroid. The intrinsic nuclear densities of these states are toroidal in nature. In this work we study these 12 C toroidal states from the generator-coordinate viewpoints in both the alpha cluster model and the toroidal shell model. Numerical solutions in the toroidal mean field approximation are examined to pave the way for future extensions and refinements.

12C Hoyle state↗
DOE OSTI2021

A Numerical Analysis of The Pool Boiling Process in Flooded Evaporators Using Different Refrigerants

The pool boiling process can be observed in several energy conversion processes including commercial and industrial refrigeration, industrial air-cooling operations, and power generation. The process becomes more involved when pool boiling in a tube bundle is considered. In the current study, a numerical model is developed to predict the key performance parameters of a flooded evaporator while considering a range of working fluids. A kettle reboiler configuration was considered, and a performance model was developed to account for boiling on individual tubes, merging of vapor bubbles, and movement under gravity. A volume of fluid (VOF) model was used to deal with the different phases in the simulation. Various fluids were considered in this study including HFE-7000, HFE-7300, and water. The trajectories of the bubbles were tracked, and the resulting information has been summarized in physical measurable quantities.

Zhang, Mingkan↗
DOE OSTI2022

The encasulation effect on thermal performance of micro-encapsulated phase change materials during energy absorption

Understanding the performance of the thermal storage materials is required to successfully meet the needs of next generation concentrating solar power projects. This study focuses on the thermal characteristics of the micro-encapsulated phase change material (MEPCM) with metallic core under high temperature. A numerical model was developed and verified to study the transient heat absorption process of MEPCMs made of tin, aluminum and copper core. Convective heat transfer rate, total energy absorption and volumetric energy density were presented. The effect of encapsulation on theses variables was analyzed under various boundary and initial conditions. Its impact on characteristic time variables was discussed. The method presented can be adopted to analyzing both manufacturing and application processes for MEPCMs of various compositions, during which melting and solidification happen. Given certain operating conditions, the results from the model can be referenced to obtain optimal thermal performance of these materials.

Benner, Jingru↗
DOE OSTI2022

Pu-238 Production Calculator and Estimation Tool

A MS Excel tool was created to help track Pu-238 production, Np-237 feedstock use, estimate production scenarios, and coordinate planning for shipments. As-run estimates were used as a basis for beginning and end of life isotopic concentrations. Multiple target types were created to demonstrate a proof of concept while actual results were not yet available The initial goals of this production tool were: • Track unirradiated, currently irradiated, and irradiated targets at the Advanced Test Reactor (ATR) • Estimate production yields from irradiated targets in various positions for each cycle • Provide past, current, and future year quantities of Np-237 and Pu-238 at ATR, amounts shipped to ATR and amounts shipped from ATR • Estimate production from current and future target designs based on data from numerical modeling (values will be provided by Idaho National Laboratory (INL)) • After proof of principle, additional work is planned that will further enhance the tool with: o Incorporate the High Flux Isotope Reactor (HFIR) reactor at Oak Ridge National Laboratory (ORNL) for irradiation o Incorporate irradiations and target shipments from ORNL to INL and back to track Np-237 transfers across the Department of Energy (DOE) Complex

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Design of Phase-Change Thermal Storage Device in a Heat Pump for Building Electric Peak Load Shaving: Preprint

Replacing carbon-intensive fossil fuel heating systems with electric heat pumps powered by renewables is a promising approach to decarbonize the building sector. However, one of the technical barriers of this approach is the large scale of heat demand, which will put excessive stress on the electricity grid. Integrating thermal energy storage (TES) into the heating systems can help alleviate this problem, by shifting thermal load and thus shaving peaks in the building electric load. Therefore, it is critical to understand how to design a thermal storage device in a heat pump for peak load shaving. In this study, we developed a numerical model for a cascaded vapor compression heat pump system integrating a phase change thermal storage device. This novel system can control the net thermal charging and discharging rate of the TES independently from the building's thermal load, which allows for precise control of electric power use. In the current study, we controlled the system to shave building electric peak load during a cold winter morning, and to charge TES during the relatively warm afternoon while still providing space heating. We used the model to evaluate the system performance with and without peak shaving. We also investigated the effect of PCM transition Tt on the peak reduction and electric energy saving potentials. The results show that peak shaving scheme effectively reduces the peak electric power consumption during a predefined discharge time window. When comparing to the no shaving case, for PCM Tt = 10 degrees C, the peak electric load reduction is 23.5%. When comparing to an air-source heat pump with back up electric heater, as Tt increases from 0 to 20 degrees C, the peak reduction increases from 46.1% to 50.9%. Integrating PCM Tt = 10 degrees C with peak shaving leads to the 45.5% of electric energy saving, which is the highest among the three transition temperatures.

electric energy saving↗
DOE OSTI2022

Development of New g-function Data for Simulating a Novel Shallow Bore Ground Heat Exchanger

Geothermal heat pump (GHP), which is also referred as ground source heat pump, is the most energy-efficient technology for space heating and cooling. However, the application of GHP is hindered by its high initial cost, of which approximately 30% is for the ground heat exchanger (GHE). In recent years, researchers have developed different types of shallow bore ground heat exchangers (SBGHE) intending to reduce the cost of drilling. The ability to predict the thermal response of an SBGHE is critical for sizing SBGHE. While g-functions have been commonly used for predicting the thermal response of conventional vertical bore ground heat exchangers (VBGHE), they cannot be directly used for predicting the performance of SBGHE because they did not account for the impacts of the seasonal variation of the soil temperature along with the depth of an SBGHE. In addition, an SBGHE has a larger thermal mass within the borehole than the VBGHE due to the larger borehole diameter of SBGHE. This study develops new g-function data pairs for a new design of SBGHE, which is named Underground Thermal Battery (UTB). Impacts of the seasonal variation of soil temperature on the thermal response of UTB were accounted for by superposing a time-dependent soil temperature onto the g-functions calculated with a numerical model that assumes constant undisturbed soil temperature. The TOUGH program was used to predict the thermal response of several configurations involving multiple UTBs. The results indicate that the proposed methodology is appropriate to generate g-functions for the UTB, and the g-function value of UTB is much lower than that of VBGHE in the time range of 15 min to 1 year due to the large thermal mass and convection heat transfer within the UTB.

Shi, Liang↗
DOE OSTI2022

Systems, methods and computer program products for collaborative agent control

Systems, methods and unmanned agents for collaboratively controlling agents in a collaborative network by one or more agents continuously simulating numeric models of one or more other agents in the network to dramatically reduce the computational bandwidth required between agents, and improve the quality of shared estimates of the agent locations as well as the locations and characteristics of other objects of interest, e.g. targets. Bandwidth is reduced by using the models to intelligently filter data before communicating.

Buerger, Stephen↗