Search NASA⌕ Search

SEARCH · Search NASA

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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 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 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 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 OSTI2025

cymyc: $\underline{C}$alabi-$\underline{Y}$au $\underline{M}$etrics, $\underline{Y}$ukawas, and $\underline{C}$urvature

We introduce cymyc, a high-performance Python library for numerical investigation of the geometry of a large class of string compactification manifolds and their associated moduli spaces. We develop a well-defined geometric ansatz to numerically model tensor fields of arbitrary degree on a large class of Calabi-Yau manifolds. cymyc includes a machine learning component which incorporates this ansatz to model tensor fields of interest on these spaces by finding an approximate solution to the system of partial differential equations they should satisfy.

differential and algebraic geometry↗
DOE OSTI2026

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗
DOE OSTI2025

Hydrodynamic characterization of the coastal pioneer array ocean observing system

Ocean observation buoys require relatively small amounts of power, yet traditionally necessitate costly resupply trips for battery replacement. With the offshore location of the buoys and small power requirements, wave energy may be an effective solution for providing consistent and reliable power to support the buoy instrumentation. The US National Science Foundation Ocean Observatories Initiative (OOI) includes arrays of point absorber-like buoy systems used for ocean observation that have been deployed at multiple locations including the Southern Mid-Atlantic Bight. A study is currently underway to design a pitch resonator wave energy converter to supplement existing renewable energy generation for powering observation instrumentation. This paper details field measurements from surface moorings of the OOI Coastal Pioneer Array, which informs the subsequent development of a numerical model for the moored observation system. The model is developed in Wave Energy Converter Simulator (WEC-Sim), which leverages the Simscape multibody solver within the MATLAB/Simulink framework and linear potential flow theory to simulate the hydrodynamic interactions and multibody dynamics in 6 degrees of freedom. Multiple tuning variables are considered to produce a model for the system that matches well with empirical data (about 8% error). In conclusion, the WEC-Sim model will serve as a platform for integrating the pitch resonator wave energy converter concept and deployment preparation (detailed design including power take-off and control systems, response evaluation, etc.).

hydrodynamic modeling↗
DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗
DOE OSTI2025

Investigation of ventilation-coupled high energy density sensible thermal energy storage

Low-cost energy storage will play an important role in supporting the decarbonization of the energy sector. Here, a novel approach to thermal energy storage for buildings is proposed, in which a tank of antifreeze liquid can be used to heat outdoor air in mechanical ventilation systems. If heated to a high enough temperature, the fluid can potentially undergo temperature swings of nearly 100°C during discharge in cold climates. Because the energy density of sensible storage systems scales proportionally to the system temperature change, this concept has the potential to offer higher energy densities than other liquid-based sensible storage devices used in building applications. A one-dimensional numerical model of the storage system was developed and experimentally validated using 30 wt% (wt%) potassium acetate as an aqueous antifreeze solution. Using a commercial hot water tank, energy densities of 47.0 kW-hour per cubic meter (kWh/m 3 ) were demonstrated in a laboratory setting. Material energy densities of approximately 74.7 (kWh/m 3 ) were measured. Design improvements may boost this energy storage density even further. Because of the system's relative simplicity, ventilation-coupled sensible storage has the potential to be an easily deployable, low-cost energy storage solution for building systems.

25 ENERGY STORAGE↗
DOE OSTI2024

Observations of airflow around a supertall curved building and its impact on temperature and humidity in Houston's urban center

Characterization of realistically shaped skyscrapers embedded in non-uniform neighbourhoods experiencing intricate weather patterns remains inadequately investigated. Aiming to close this gap, the Center for Multiscale Applied Sensing team deployed its mobile observatory in the street canyons around the curved Wells Fargo Plaza skyscraper in downtown Houston, TX. Three deployments allowed airflow observations under different inflow wind and thermodynamic stability conditions. Doppler lidar measurements reveal that when inflow hits the curved wall of the skyscraper, perpendicular canyons experience similar vortex configurations creating two windward and two leeward circulations. Windsond measurements support that buoyancy within the deep street canyons can generate thermal updrafts as strong as 2 m s -1 which is sufficient to overturn the mechanical downwash under gentle wind conditions. Canyons experiencing venting during the daytime were observed to be more thermodynamically stable at night while thermodynamically stable canyons during the day were observed to be more thermodynamically unstable at night owing to the accumulation of heat near street-level. Fourier decomposition of the vertical velocity measurements shows that in all cases flow exhibited high Reynolds numbers and was composed of turbulent eddies of predominantly 6 and 15 min periods. Here, this observational dataset provides insights to assess wind load, pedestrian comfort, urban air mobility, and natural ventilation and may be used as a benchmark for numerical model and wind tunnel studies attempting to represent realistically complex urban conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
DOE OSTI2025

Development of fiber orientation in injection molding: Comparison of glass fiber, carbon fiber and their hybrid composites

Fiber orientation distribution (FOD) in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for glass fiber (GF), carbon fiber (CF), and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varied with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical framework for integrated additive manufacturing-compression molding (AM-CM) of thermoplastic composites

Additive manufacturing-compression molding (AM-CM) has emerged as a transformative technology in advanced composite manufacturing. Additive manufacturing (AM) offers high design flexibility and the ability to produce complex geometries with precisely aligned fibers in the preferred orientation. Compression molding (CM) enhances composite materials by providing excellent dimensional stability, reduced porosity, high production rates, and a smooth surface finish. Despite these advantages, extensive integrated analysis is required to optimize processing conditions for improved fiber orientation distribution (FOD) and porosity control. Here, this study develops a comprehensive numerical model to simulate the AM-CM manufacturing process. The model isolates the effects of both the AM and CM phases while also capturing their integration. Additionally, it accounts for heat transfer, temperature-dependent viscosity, and fiber orientation in the extruded fiber-filled polymer, accurately representing material behavior during processing. This approach enables the analysis of interactions between deposited beads of complex strand shapes and their interface regions after full compression. Moreover, the model predicts key parameters such as polymer flowability, fiber orientation, and temperature evolution in AM-CM parts. By optimizing processing conditions, it facilitates a controlled and predictable microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Geologic hydrogen: From natural occurrences to anthropogenic generation – A review of fundamentals, potential, challenges and prospects

Growing demand for hydrogen is exposing the environmental and economic limits of reforming-based and carbon-managed supply chains, while the scale-up of electrolytic capacity remains capital-constrained. Geologic hydrogen, defined as molecular H₂ generated and stored within the Earth's crust offers a complementary, potentially lower-cost resource, yet exploration is still ad hoc. This review (1) revisits a global inventory of confirmed hydrogen seeps and subsurface occurrences; (2) analyzes the controlling reactions, migration pathways, and trapping conditions governing these occurrences; (3) proposes a process-based geologic hydrogen system concept analogous to, yet distinct from, the petroleum system; and (4) evaluates potential geologic hydrogen systems within the United States as a representative case study. Here, we contrast natural systems powered by serpentinization, mantle degassing or radiolysis with anthropogenic systems that stimulate the same reactions or convert in-situ hydrocarbons. Stable hydrogen accumulations require generation rates that outpace combined physical, chemical and microbial losses; the Bourakébougou field (Mali) exemplifies a self-recharging, free-gas reservoir sustained by meteoric-water serpentinization beneath an efficient caprock. Prospective geologic hydrogen resources are likely to occur in regions where iron-rich lithologies, deep-seated faults, and low-permeability sealing formations coexist. Applying this principle, we highlight three promising hydrogen play types in U.S. geological terrains: ophiolite belts (Appalachian and Californian regions), the Midcontinent Rift and the Lake Superior banded‑iron formations. Multiphysics numerical models and positive-unlabeled machine-learning workflows help to accelerate play screening and de-risk future production; yet, reaction kinetics, stimulation strategies, and full techno-economic and life-cycle assessments remain pivotal knowledge gaps.

Anthropogenic hydrogen generation↗
DOE OSTI2025

Uniform secondary current distribution at disk electrodes under Tafel kinetics enabled by concentric current-shielding rings

Non-uniform secondary current distribution at rotating disk electrodes (RDE) is a common problem when using resistive electrolyte media or large applied currents. In a recent publication, we have shown that auxiliary electrodes such as the ring of a rotating ring disk electrode (RRDE) can help suppress current non-uniformities at the disk enabling reliable electroanalytical measurements. However, this previous work considered linear kinetics where current distribution non-uniformities were moderate. In the present contribution, we consider current distribution non-uniformities encountered under Tafel kinetics. We show, for the case of 2,5-dihydroxy-1,4-benzoquinone (DHBQ) reduction, that optimally chosen ring conditions serve to provide effective shielding at the disk edge rendering the overall disk current distribution to be uniform. Numerical modeling and scaling analysis (using the Wagner number) are presented to aid a user in determining the optimal ring current density for achieving uniform disk current distribution under Tafel kinetics. This approach is especially useful when studying soluble-soluble redox transitions for which, unlike deposit distribution in electrodeposition, the current distribution non-uniformity is not visually apparent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗