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

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Lightweight, Flexible Electromagnetic Shielding Composite Films Reinforced with Recycled Carbon Fibers and Carbon Nanofillers

Lightweight, flexible polymer composites are widely adopted in modern electronic devices and systems for high-efficiency electromagnetic interference (EMI) shielding. Reinforcing polymer composites with conductive fillers offers a promising alternative to conventional metal-based shielding material thanks to their low density, tunable properties, and flexibility. Herein, lightweight, flexible ultra-high molecular weight polyethylene composite films reinforced with recycled carbon fibers (rCFs) and carbon-based nanofillers, including graphene nanoplatelets (GNPs) and carbon nanotubes, are reported for EMI shielding applications. The incorporation of these nanofillers significantly reduces the required content of rCFs and improves processability while maintaining comparable shielding effectiveness. The addition of these nanofillers with rCFs enhances the EMI shielding effectiveness by up to 15 dB. Incorporating 1 wt% GNPs can replace 5 wt% rCFs and achieve comparable EMI shielding performance, while 5 wt% of either nanofillers can substitute for 10 wt% rCFs. Numerical modeling of electromagnetic wave transmission reveals that increasing nanofiller concentrations enhances both reflection and absorption losses, with absorption consistently dominating across all levels. Furthermore, this study not only provides insight into the synergistic contributions of rCFs and carbon nanofillers to shielding effectiveness but also paves the way for the design of sustainable, lightweight EMI shielding composite films for applications in electric vehicles, medical equipment, and portable electronics.

carbon nanofiller↗
DOE OSTI2021

A Fabrication Strategy for Reconfigurable Millimeter–Scale Metamaterials

Rather than depending on material composition to primarily dictate performance metrics, metamaterials can leverage geometry to achieve specific properties of interest. For example, reconfigurable metamaterials have enabled programmable shape transformations, tunable mechanical properties, and energy absorption. While several methods exist to fabricate such structures, they often place severe restrictions on manufacturing materials, or require significant manual assembly. Moreover, these arrays are typically composed of unit cells that are either macro–scale or micro–scale in dimension. Here, the fabrication gap is bridged, and laminate manufacturing is used to develop a method for designing reconfigurable metamaterials at the millimeter–scale, that is compatible with a wide range of materials, and that requires minimal manual assembly. In addition to showing the versatility of this fabrication method, how the use of laminate manufacturing affects the behavior of these multi–component arrays is also characterized. To this end, a numerical model that captures the deformations exhibited by the structures is developed, and an analytic model that predicts the strain of the structure under compressive stress is built. Overall, this approach can be leveraged to develop millimeter–scale metamaterials for applications that require reconfigurable materials, such as in the design of tunable acoustics, photonic waveguides, and electromagnetic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Understanding Recombination Pathways of CdSeTe/CdTe Devices Using Temperature Dependence of Fill Factor

Understanding the temperature dependence of photovoltaic parameters is critical for identifying performance limiting mechanisms in CdTe-based solar cells. In this study, we investigate the temperature dependence of fill factor (FF(T)) in CdSeTe/CdTe solar cells using a combination of experimental current density-voltage (JV) measurements and numerical simulations. The measured FF(T) exhibits a nonlinear response with temperature, peaking at moderate temperature values (10–30°C). Through numerical modeling, we reveal that this peaking behavior in FF(T) arises when the dominant recombination mechanism transitions from bulk to interface with increasing temperature. This transition is primarily driven by low hole mobility (0.6–1.2 cm2/V·s), while electron mobility, interface recombination velocities, and carrier lifetime also contribute to the effect. Our findings demonstrate that FF(T) analysis is a promising diagnostic tool for distinguishing recombination pathways and guiding the optimization of thin-film CdTe solar cells under real-world conditions.

CdTe photovoltaics↗
DOE OSTI2025

Subvoxel Control of Fiber Orientation via Multidirectional Shearing in 3D Printing

Anisotropy, the characteristic of materials exhibiting different properties based on their direction, is widespread in nature. Conventional manufacturing techniques often fall short in recreating such complex anisotropy. While 3D printing allows for precise fiber deposition in anisotropic composites, previous studies have only achieved bulk reorientation of fibers at the voxel level in two dimensions. This limits the replication of localized 3D anisotropies found in natural materials. To address this, a novel 3D printing technique is presented that enables subvoxel control of fiber orientation in all three directions using multidirectional shearing via nozzle rotation and inclination. The fiber orientation control in these experimental tests is driven by a numerical model, enabling a fully digital approach to program microstructures within a strand. This programmability is demonstrated through mechanical and thermal tests, illustrating a localized and controllable response to external stimuli. Achieving such complex anisotropy holds potential across several fields, including wearables and biomedical implants, lightweight composite structures, and energy storage technologies such as batteries and supercapacitors.

3D printing↗
DOE OSTI2022

Very small-scale, segregating-fluidized-bed experiments: A dataset for CFD-DEM validation and uncertainty quantification

IWe report fluidization experiments were conducted on a small scale and with a rapid response (short duration) to enable corresponding simulations at low-computational cost. Rise times are reported for four or fewer polyethylene particles (intruders) in an air-fluidized bed of ~5000 group D glass beads. Experimental inputs were completely characterized—particle properties, system dimensions and operating conditions—which is necessary for validating computational fluid mechanics (CFD)-discrete element method (DEM) including a comprehensive uncertainty quantification (UQ) analysis. Input uncertainties are reported as bounds or cumulative distribution functions of measured values. The staggering number of simulations required to complete a UQ analysis (~O[10 4 ] simulations corresponding to ~5 uncertain inputs) motivates this study. These segregating-bed experiments are designed to permit analogous CFD-DEM simulations to complete in less than a day on a single (~2.5 GHz) computational processor unit (CPU). Segregation times are reported for several operating conditions, intruder sizes, and initial configurations, providing a rich dataset for numerical model testing, validation and UQ.

42 ENGINEERING↗
DOE OSTI2023

Convectively coupled Rossby–Gravity waves in a field campaign: How they are captured in reanalysis products

Abstract Convectively coupled equatorial waves are a significant source of atmospheric variability in the tropics. Current numerical models continue to struggle in simulating the coupled diabatic heating fields that are responsible for the development and maintenance of these waves. This study investigates how the diabatic fields associated with Mixed Rossby–Gravity waves (MRGs) are represented in four reanalysis products by using a unique observational dataset from the TRMM‐KWAJEX (Tropical Rainfall Measuring Mission‐Kwajalein Experiment) field campaign. These reanalyses include ERA5, Japanese 55‐year Reanalysis (JRA‐55), Climate Forecast System Reanalysis (CFSR), and Modern‐Era Retrospective Analysis for Research and Applications (MERRA). We found that all four reanalyses captured the MRG structures in winds and temperature, and to a lesser degree in the humidity field except in the boundary layer. However, only the ERA5 and MERRA reanalyses captured the gradual rise and succession of the diabatic heating from boundary layer turbulence, shallow convection, cumulus congestus, and deep convection within the waves. ERA5 is the only product that also captured the gradual rise of the subgrid‐scale vertical transport of moist static energy. All reanalysis products underestimated the diabatic heating from cumulus congestus. Results provide observational basis on what aspects of MRG can be trusted and what cannot in the reanalysis products.

Wang, Xiaocong↗
DOE OSTI2023

Characterizing the seismic response of a molten salt nuclear reactor

Abstract One Generation IV nuclear reactor, which uses a fluoride salt as a coolant, graphite reflector blocks as a moderator, and circulating buoyant TRISO pebbles as fuel is at an advanced stage of development. To characterize the seismic behavior of components of this reactor, validate numerical models for analysis, and develop recommendations for design, a set of earthquake‐simulator experiments on a scaled model of the reactor vessel and its internals was executed on a six‐degree‐of‐freedom earthquake simulator. The model was seismically isolated at its base using two types of spherical sliding bearings. The scaled model involved representations of the prototype reactor vessel, core barrel, reflector blocks, coolant, and spherical fuel pebbles. The material and geometric properties of different test components were selected based on a dynamic similitude scaling analysis and an approximate length scale of 0.4. Four sets of three‐component earthquake motions were used as inputs for testing. Instrumentation on the test specimen recorded the dynamic responses of the outer vessel, core barrel, and reflector‐block assembly, the hydrodynamic responses (sloshing and hydrodynamic pressure) of the liquid coolant, pebble consolidation under earthquake shaking, and the behavior of the isolation systems. This paper describes the design of the experiments and presents key results from the tests. The dynamic responses of the outer vessel, core barrel, and the reflector blocks revealed that the components responded as a unit for the intense shaking used in the experiments. The sloshing response of the fluid in a thin annulus near the perimeter of the vessel was heavily damped. The change in the packing fraction of the pebble bed under repeated, intense 3D earthquake shaking was less than 3%. Seismically isolating the vessel substantially reduced demands on its internal components.

Engineering↗
DOE OSTI2021

Stream discharge and water quality data for East Fork Poplar Creek beginning 2012

This data note describes 15-min discharge and in situ water quality data at two locations along East Fork Poplar Creek in east Tennessee, USA. Data records include temperature, gauge height, water surface elevation above mean sea level, and volumetric discharge. Water quality measurements include temperature, specific conductance, pH, dissolved oxygen (percent saturation and concentration), turbidity, and less extensive fDOM data at one site. The data records begin in 2012 at one site and 2015 at the second site; monitoring at both sites is ongoing (as of 2021). The goal of this data collection is to improve understanding of watershed functions, hydrologic dynamics, and material flux. The data will contribute to site conceptual and numerical models, exposure and risk evaluation, remediation selection and design, and performance monitoring. The data are publicly available and can be accessed via unique url or DOI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Understanding ERE and iVOC Metrics for Graded CdSeTe Absorbers

PL-based external radiative efficiency (ERE) and implied open-circuit voltage (iVOC) metrics were introduced for thin-film solar absorbers to better understand the voltage deficit and diagnose losses in solar cells. Traditionally, elevated ERE and iVOC measurements are associated with diminished recombination within the solar device, a rationale heavily reliant on the assumption of a uniform bandgap and high carrier mobilities in the absorber. Recently, very low mobilities in CdSeTe absorbers (< 1 cm2/(V.s)) were measured using the light-induced transient grading technique. In this study, we use a detailed numerical model of iVOC to investigate the possible reasons of elevated iVOC in realistic CdSeTe absorbers with a graded Se profile. In particular, we examine how the bandgap nonuniformity and the reduced hole mobility in graded CdSeTe absorbers affect iVOC measurements. We show that high iVOC may result from inflated quasi-Fermi level splitting in the high-Se region in the front part of a CdSeTe absorber with slow hole transport. We reproduce the experimentally reported 360 mV increase in iVOC-VOC gap with reduced doping using a model with sub-1 cm2/(V.s) hole mobility in the high-Se region. Based on our results, we conclude that the iVOC metric (or ERE metric) should not be used as a sole metric of CdSeTe absorber quality. We discuss possible ways to extract useful information from the iVOC-VOC gap by supplementing the front-side illumination measurements with back-side illumination measurements.

14 SOLAR ENERGY↗
DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗
DOE OSTI2025

Multiphysics simulation of recent experiments on alkali‐silica reaction expansion in reinforced concrete members

Alkali‐silica reaction (ASR) is an important degradation process that causes volumetric expansion and damage in concrete, and is affected significantly by the local temperature, moisture and stress conditions that often vary across the regions of a structure. Numerical simulation is essential to predict the progression and effects of ASR on the performance of structures. Because of the interactions between thermal and moisture transport and mechanical deformation, it is important for numerical models to represent all these physical phenomena and the coupling between them. Simulations of ASR in reinforced concrete (RC) structures are further complicated by the need to capture interactions between concrete and embedded reinforcing bars. Here, this paper describes the implementation of a scalable, coupled‐physics ASR model for simulating RC structures and assesses the ability of that model to predict ASR‐induced expansion in recent laboratory tests on RC block and beam specimens. These laboratory tests and the simulation approach were selected because of their applicability to RC structural‐scale simulations. This validation study helps builds confidence the ability of this approach to model ASR expansion in large, complex RC structures, which is a current high‐priority need.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Developing a strategy for the national coordinated soil moisture monitoring network

Soil moisture is a critical land surface variable, affecting a wide variety of climatological, agricultural, and hydrological processes. Determining the current soil moisture status is possible via a variety of methods, including in situ monitoring, remote sensing, and numerical modeling. Although all of these approaches are rapidly evolving, there is no cohesive strategy or framework to integrate these diverse information sources to develop and disseminate coordinated national soil moisture products that will improve our ability to understand climate variability. The National Coordinated Soil Moisture Monitoring Network initiative has developed a national strategy for network coordination with NOAA's National Integrated Drought Information System. The strategy is currently in review within NOAA, and work is underway to implement the initial milestones of the strategy. This update reviews the goals and steps being taken to establish this national-scale coordination for soil moisture monitoring in the United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Changes in the Regional Water Cycle and Their Impact on Societies

ABSTRACT Changes in “blue water”, which is the total supply of fresh water available for human extraction over land, are quite closely related to changes in runoff or equivalently precipitation minus evaporation, . This article examines how climate change‐driven recent past and future changes in the regional water cycle relate to blue water availability and changes in human blue water demand. Although at the largest scales theoretical and numerical model predictions are in broad agreement with observations, at continental scales and below models predict large ranges of possible future and runoff especially at the scale of individual river catchments and for shorter timescale subseasonal floods and droughts. Nevertheless, it is expected that the occurrence and severity of floods will increase and that of droughts may increase, possibly compounded by human‐driven non‐climatic changes such as changes in land use, dam water impoundment, irrigation and extraction of groundwater. Contemporary assessments predict that increases in 21st century human water extraction in many highly‐populated regions are unlikely to be sustainable given projections of future . To reduce uncertainty in future predictions, there is an urgent need to improve modeling of atmospheric, land surface and human processes and how these components are coupled. This should be supported by maintaining the observing network and expanding it to improve measurements of land surface, oceanic and atmospheric variables. This includes the development of satellite observations stable over multiple decades and suitable for building reanalysis datasets appropriate for model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nonperturbative Topological Phenomena in QCD and Related Theories

Here, this book introduces a variety of aspects in nonperturbative Quantum Chromodynamics (QCD), focusing on the topological objects present in gauge theories. These objects, like magnetic monopoles, instantons, instanton-dyons, sphalerons, QCD flux tubes, etc, are first introduced individually and, later, treated collectively. As ensembles, they produce various phenomena that can be modeled numerically in lattice gauge theories and such collective effects, produced on the lattice, are extensively discussed in some chapters. In turn, the notion of duality, which is crucial in modern field/string theories, is elucidated by taking into consideration the electric-magnetic duality, the Poisson duality, and the AdS/CFT duality.This monograph is based on various lectures given by Edward Shuryak at Stony Brook during the last three decades and it is meant for advanced graduate students and young researchers in theoretical and mathematical physics who are willing to consolidate their knowledge in the topological phenomena encountered in fundamental QCD research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Joining of High Strength Low Ductility AA7055 by Friction Self-piercing Rivet

High strength aluminum (Al) alloy is one of higher specific strength materials for decarbonization in transportation industries. Because of low ductility at room temperature, conventional mechanical fastening such as self-piercing riveting produces cracks at the joint. In this work, we applied friction self-piercing riveting to join Al alloy (AA) 7055. No cracks were observed in the joints because of the improved local ductility of Al alloy by the generated frictional heat during joining step. Numerical modeling of joining process was applied to guide rivet geometry design and rivet material strength. Mechanical integrity of the AA7055 joints was assessed by lap shear tensile and cross-tension testing. Metallurgical characterizations revealed solid-state bonding formed not only between the rivet and surround Al materials, but also upper and lower Al sheets at the joint interface. Both solid-state bonding and mechanical interlocking between the flared rivet and bottom AA7055 sheet were the major joint mechanisms.

Lim, Yong Chae↗