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,405 records · Page 78

A two-dimensional model of plasma expansion in the ionosphere

A systematic parameter study is conducted of the motion of artificial plasma clouds across the geomagnetic field at ionospheric altitudes. The study is based on a two-dimensional numerical model. Taken into consideration are the effects solar cycle variations on the cloud evolution, as well as the effects of seasonal and geomagnetic activity conditions, and variations in the ionospheric release conditions. The effects of magnetospheric electric fields and thermospheric winds are also considered. The plasma cloud model and related assumptions are discussed, and numerical results are presented.

Ma, T.-Z.↗

Higher-order LaSDI: Reduced order modeling with multiple time derivatives

Solving complex partial differential equations (PDEs) is essential across scientific disciplines but often requires numerical models that can be prohibitively expensive in time-sensitive applications. Reduced-order models (ROMs) address this challenge by exploiting low-dimensional structure to create fast approximations. The Latent Space Dynamics Identification (LaSDI) framework has demonstrated success in learning ROMs for parameterized PDE families, but remains limited to first-order systems. Here, in this paper, we propose Higher-Order LaSDI (HLaSDI), which extends the LaSDI framework to PDEs with arbitrary order of time derivatives. This generalization significantly expands the applicability of LaSDI-based methods to systems previously outside their scope, including hyperbolic PDEs. We demonstrate HLaSDI’s accuracy and efficiency on several linear and nonlinear benchmark problems.

97 MATHEMATICS AND COMPUTING↗

A systematic method for selecting molecular descriptors as features when training models for predicting physiochemical properties

Machine learning has proven to be a powerful tool for accelerating biofuel development. Although numerous models are available to predict a range of properties using chemical descriptors, there is a trade-off between interpretability and performance. Neural networks provide predictive models with high accuracy at the expense of some interpretability, while simpler models such as linear regression often lack in accuracy. In addition to model architecture, feature selection is also critical for developing interpretable and accurate predictive models. We present a method for systematically selecting molecular descriptor features and developing interpretable machine learning models without sacrificing accuracy. Our method simplifies the process of selecting features by reducing feature multicollinearity and enables discoveries of new relationships between global properties and molecular descriptors. To demonstrate our approach, we developed models for predicting melting point, boiling point, flash point, yield sooting index, and net heat of combustion with the help of the Tree-based Pipeline Optimization Tool (TPOT). For training, we used publicly available experimental data for up to 8351 molecules. Our models accurately predict various molecular properties for organic molecules (mean absolute percent error (MAPE) ranges from 3.3% to 10.5%) and provide a set of features that are well-correlated to the property. This method enables researchers to explore sets of features that significantly contribute to the prediction of the property, offering new scientific insights. To help accelerate early stage biofuel research and development, we also integrated the data and models into a open-source, interactive web tool.

09 BIOMASS FUELS↗

Use of Remote Sensing and In-situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the Accelerated Climate Model for Energy. Final report

Representation of unresolved cloud processes has been one of the most challenging problems in numerical modeling of the atmosphere. During the last decades, considerable efforts have been spent to better represent the collective effects of those processes in terms of resolved scales (see Arakawa 2004 for a review). Still, the progress on this problem has been unacceptably slow as pointed out by Randall et al. (2003). The large part of the difficulty is inherent in the artificial separation of atmospheric processes into the explicitly simulated resolved scales and the highly parameterized subgrid scales. Arakawa (2004) and Arakawa et al. (2011) suggested two plausible approaches to eliminate the problem of artificial separation of atmospheric processes: “parameterize everything” and “resolve everything”.

58 GEOSCIENCES↗

Initialization of mesoscale models - The possible impact of remotely sensed data

Little or no improvement appears to have been achieved in mesoscale numerical prediction. This has been blamed on inherent unpredictability, insufficient spatial resolution, inadequate parameterization of important physical processes, and errors in the numerical discretization of the governing hydrodynamic equations. Attention is presently given to the more fundamental difficulty of forecasting when present conditions are not known with sufficient accuracy, and to the matching of remotely sensed sounding-capability satellite and all-weather Doppler radar data with mesoscale numerical models. This matching is not straightforward; simple and variational forms of four-dimensional assimilation, and Kalman filtering, are suggested approaches but their practical testing remains to be accomplished. The most severe problem faced is that of the temperature and humidity profiles' vertical resolution.

Gal-Chen, T.↗

Low-frequency sound propagation modeling over a locally-reacting boundary using the parabolic approximation

There is substantial interest in the analytical and numerical modeling of low-frequency, long-range atmospheric acoustic propagation. Ray-based models, because of frequency limitations, do not always give an adequate prediction of quantities such as sound pressure or intensity levels. However, the parabolic approximation method, widely used in ocean acoustics, and often more accurate than ray models for lower frequencies of interest, can be applied to acoustic propagation in the atmosphere. Modifications of an existing implicit finite-difference implementation for computing solutions to the parabolic approximation are discussed. A locally-reacting boundary is used together with a one-parameter impedance model. Intensity calculations are performed for a number of flow resistivity values in both quiescent and windy atmospheres. Variations in the value of this parameter are shown to have substantial effects on the spatial variation of the acoustic signal.

Robertson, J. S.↗

MAMI: Modeling of the Magnetosphere-Ionosphere-Atmosphere System

Major emphasis of the investigation was the development of theoretical and numerical models of the aurora and of high latitude ionospheric processes. In particular: (1) the NCAR TIGCM (Thermosphere-Ionosphere Global Circulation Model) was updated to include mid and low latitude electrodynamics (this version is called TIE-GCM); (2) the NCAR TIE-GCM was modified to include a more realistic representation of the aurora; (3) the UAF auroral electron transport model was modified to include local acceleration processes; (4) a local ionospheric and auroral model (AURORA) was developed to allow detailed studies of the aurora; (5) a proton-hydrogen transport code has been developed to model proton aurora; and (6) a theory for the production of suprathermal atoms and ions in the upper atmosphere was developed and applied to studies of atomic nitrogen transport and helium escape on open field lines. These models enable us to devise schemes for tile interpretation and quantitative analysis of data obtained by the POLAR spacecraft. Parameterizations were formulated and are available to the GGS investigators. The UVI and VIS teams have adopted these parameterizations and include them in their data analysis. We have developed software for the quantitative interpretation of UVI and VIS images. After the launch of the POLAR satellite we used the image data from UVI and VIS in combination with ground based data from SuperDARN and incoherent scatter radars, magnetometers, and in situ observations from DMSP and NOAA satellites to characterize the state of the ionosphere. A number of event studies have been carried out in cooperation with other CGS theory teams.

Source record↗

Observations and Modeling of the Transient General Circulation of the North Pacific Basin

Because of recent progress in satellite altimetry and numerical modeling and the accumulation and archiving of long records of hydrographic and meteorological variables, it is becoming feasible to describe and understand the transient general circulation of the ocean (i.e., variations with spatial scales larger than a few hundred kilometers and time scales of seasonal and longer-beyond the mesoscale). We have carried out various studies in investigation of the transient general circulation of the Pacific Ocean from a coordinated analysis of satellite altimeter data, historical hydrographic gauge data, scatterometer wind observations, reanalyzed operational wind fields, and a variety of ocean circulation models. Broadly stated, our goal was to achieve a phenomenological catalogue of different possible types of large-scale, low-frequency variability, as a context for understanding the observational record. The approach is to identify the simplest possible model from which particular observed phenomena can be isolated and understood dynamically and then to determine how well these dynamical processes are represented in more complex Oceanic General Circulation Models (OGCMs). Research results have been obtained on Rossby wave propagation and transformation, oceanic intrinsic low-frequency variability, effects of surface gravity waves, pacific data analyses, OGCM formulation and developments, and OGCM simulations of forced variability.

McWilliams, James C.↗

A two-dimensional analytical unit cell model for redox flow battery evaluation and optimization

Cell performance optimization is important for improving the overall system efficiency of a redox flow battery. To gain better insights into key controlling factors of system efficiency, this work first proposed a theoretical model for a unit flow battery cell by extending a two-dimensional analytic model to a full battery cell. Such a model is then used for cell performance optimization after validating it with experimental and numerical modeling data. With the model results, the activation, equilibrium, and pump energy losses are identified as the dominant sources of battery energy losses. Further, a guideline for reducing these sources is also proposed. Following the guideline, the mass transport coefficient is shown as a key control factor of the equilibrium energy loss and Coulombic efficiency (CE). Approaches are then proposed to improve CE and the overall system efficiency. The mass transport is also revealed as the mechanism of pump rate optimization where an optimal pump rate significantly reduces the equilibrium energy loss. With both low equilibrium and pump energy losses, an optimal electrode porosity or specific area design can further improve a battery's system efficiency based on an optimal porosity predicted by the present model. The model also demonstrates distinct behaviors and overestimation in the system efficiency when reduced to a zero-dimensional model with neglected mass transport resistance. With the new model, the guideline, and new insights, this work provides a reliable and efficient tool for the evaluation and optimization of redox flow battery design in practical applications.

25 ENERGY STORAGE↗

Understanding Differences in California Climate Projections Produced by Dynamical and Statistical Downscaling

Abstract We compare historical and end‐of‐century temperature and precipitation patterns over California from one dynamically downscaled simulation using the Weather Research and Forecast (WRF) model and two simulations statistically downscaled using Localized Constructed Analogs (LOCA). We uniquely separate causes of differences between dynamically and statistically based future climate projections into differences in historical climate (gridded observations versus regional climate model output) and differences in how these downscaling techniques explicitly handle future climate changes (numerical modeling versus analogs). In these methods, solutions between different downscaling techniques differ more in the future compared to the historical period. Changes projected by LOCA are insensitive to the choice of driving data. Only through dynamical downscaling can we simulate physically consistent regional springtime warming patterns across the Sierra Nevada, while the statistical simulations inherit an unphysical signal from their parent Global Climate Model (GCM) or gridded data. The results of our study clarify why these different techniques produce different outcomes and may also provide guidance on which downscaled products to use for certain impact analyses in California and perhaps other Mediterranean regimes.

54 ENVIRONMENTAL SCIENCES↗

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗

Modeling the Low-Pressure High-Voltage Branch of the Paschen Curve for Hydrogen and Deuterium

A physical and numerical model of the Townsend discharge in molecular hydrogen and deuterium has been developed to meet the needs of designing a plasma-based switching device for power grid applications. Here, the model allows to predict the low-pressure branch of the Paschen curve for applied voltage in the range of several hundred kilovolts. In the regime of interest, electrons are in a runaway state and ionization by ions and fast neutrals sustains the discharge. It was essential to correctly account for both gas-phase and surface interactions (electron emission and electron backscattering), especially in terms of their dependence on particle energy. The model yields results consistent with prior data obtained for lower voltage. The three-species (electrons, ions, and fast neutrals) model successfully captures the essential physics of the process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

pnnl/EZBattery

Cell performance optimization is important for improving the system efficiency of a redox flow battery. To gain better insights into key controlling factors Of system efficiency, this work first proposed a theoretical model for a unit cell by extending a two-dimensional analytic model to a full cell. The model is then used for cell performance optimization after validating it with experimental and numerical modeling data.

Bao, Jie↗

Dynamic Rupture Modeling of the 1999 Chi-Chi, Taiwan Earthquake Using 3DFinite Element Method. Sensitivity Analysis of Slip Rate Function to Model Parameters of Crustal Weak-Zone

The main objective of this study was the dynamic rupture modeling of the M7.3 1999 Chi-Chi, Taiwan earthquake. The purpose of the numerical modeling was two folded. First, using models of rupture dynamics, derive kinematic rupture characteristics, including spatial and temporal variations of the slip rate time history and rupture speed, in the long-period ground motion generation area (LMGA), strong motion generation areas (SMGAs) and background fault area for a thrust shallow rupture. Second, investigate the sensitivity of source time function characteristics in the LMGA and SMGA areas to model parameterization of low-velocity weak zone, such as thickness and shear wave velocity.

58 GEOSCIENCES↗

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The use of geomagnetic field models in magnetic surveys

The importance of global geomagnetic field models for the reduction of magnetic surveys is discussed. It is demonstrated that a numerical model with adequate secular variation correction, provides a suitable representation of the regional field. The limitations of the presently available models are reported, with emphasis on the International Geomagnetic Reference Field.

Regan, R. D.↗

Sympathetic reaction tests and analyses

The critical separation distances for explosive reactions of a solid rocket propellant were measured. Explosive reactions included low order explosion, low order detonation, and high order detonation. The effects of sample size, shape, damage and temperature on sympathetic reaction were determined experimentally. The sympathetic detonation of small cubes of solid rocket propellant was modelled numerically, using the Eulerian reactive hydrodynamic code 2DE with Forest Fire burn rates. The model was applied to cubes of 2.54 - 7.62 cm (1 - 3 in.), with agreement between calculated and experimental results.

Ricardson, D. E.↗