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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Diffusioosmotic flow reversals due to ion–ion electrostatic correlations

Existing theories of diffusioosmosis have neglected ion–ion electrostatic correlations, which are important in concentrated electrolytes. Here, we develop a mathematical model to numerically compute the diffusioosmotic mobilities of binary symmetric electrolytes across low to high concentrations in a charged parallel-plate channel. We use the modified Poisson equation to model the ion–ion electrostatic correlations and the Bikerman model to account for the finite size of ions. We report two key findings. First, ion–ion electrostatic correlations can cause a unique reversal in the direction of diffusioosmosis. Such a reversal is not captured by existing theories, occurs at ≈ 0.4 M for a monovalent electrolyte, and at a much lower concentration of ≈ 0.003 M for a divalent electrolyte in a channel with the same surface charge. This highlights that diffusioosmosis of a concentrated electrolyte can be qualitatively different from that of a dilute electrolyte, not just in its magnitude but also its direction. Second, we predict a separate diffusioosmotic flow reversal, which is not due to electrostatic correlations but the competition between the underlying chemiosmosis and electroosmosis. This reversal can be achieved by varying the magnitude of the channel surface charge without changing its sign. However, electrostatic correlations can radically change how this flow reversal depends on the channel surface charge and ion diffusivity between a concentrated and a dilute electrolyte. Furthermore, the mathematical model developed here can be used to design diffusioosmosis of dilute and concentrated electrolytes, which is central to applications such as species mixing and separation, enhanced oil recovery, and reverse electrodialysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Comprehensive Calibration Framework for the Northwest River Forecast Center

We present a comprehensive framework developed by the Northwest River Forecast Center for calibrating hydrologically diverse basins. The framework includes models for snow, soil moisture, routing, channel loss, and consumptive use. Data inputs include a wide range of open-access datasets for meteorology, land use, topography, and land cover. The framework uses conceptual hydrologic models to handle basins with various hydrologic regimes including rain-driven and snowmelt-dominated basins. We also develop a flexible automatic calibration system that can handle numerous unobservable model parameters in a computationally efficient manner. A single-basin automatic calibration run can typically be completed on a modern laptop in under 10 min. We found that model performance metrics for this new approach match the quality of the NWRFC's previous labor-intensive manual calibrations. The model performance also rivals that of a state-of-the-art deep learning model at a fraction of the computational cost. This framework presents a new standard for the quality of calibrations possible with lumped conceptual hydrologic models, combining careful data curation, an objective calibration framework, and expert local knowledge. In addition, we have made software packages available for the entire suite of National Weather Service River Forecast System models, including SAC-SMA, SNOW-17, and Lag-K. These modern interfaces are intended to increase accessibility and facilitate future research.

Forecasting↗
DOE OSTI2024

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗
DOE OSTI2023

Morrow1 Model - with Capillary Pressure

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Morrow1 model uses the Morrow1 relative permeability curve with capillary pressure assigned homogeneously

Capillary Pressure↗
DOE OSTI2023

Hydrostratigraphic Region 1 Model

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 1 model uses the HS1 relative permeability curves assigned heterogeneously by hydrostratigraphic unit. The model with capillary pressure applies the HS1 curves heterogeneously by hydrostratigraphic unit.

Capillary Pressure↗
DOE OSTI2023

Hydrostratigraphic Region 3 Model - with and without HSU Capillary Pressure

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 3 model uses the HS3 relative permeability curves assigned heterogeneously by hydrostratigraphic unit. The model with capillary pressure applies the HS3 curves heterogeneously by hydrostratigraphic unit.

Capillary Pressure↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Evaluating and modeling rate sensitivity in advanced reactor structural materials: 316H, Gr. 91, and A617

Advanced reactor systems operate at high temperatures to achieve better thermal efficiency. During transients from lower temperatures to high operating temperatures structural metals transition from rate insensitive behavior to rate dependent behavior. In the rate dependent regime creep and plastic deformation are not separable hence requiring a unified treatment with a viscoplastic model. Several bounding design methods in the ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B assume a decoupled treatment of creep and plasticity so establishing the boundary of the rate insensitive regimes for advanced reactor structural materials is necessary to set limits for the use of these bounding design methods. In the rate sensitive regime new bounding methods or fully inelastic analysis are required for the design of Division 5 Class A components. During temperature transients these components may operate in both the rate insensitive and rate sensitive regions, requiring computational models that can represent both types of behavior. This work establishes criteria for material rate sensitivity and develops a new composite method for modeling constitutive response in both the rate sensitive and rate insensitive regimes for three Division 5, Class A materials: 316H austenitic stainless steel, Grade 91 ferritic-martensitic steel, and Alloy 617, a nickel-based alloy. Both the criteria and new modeling strategy apply the concept of a normalized activation energy for thermally activated mechanisms developed by Kocks, Mecking, and co-workers. Specifically, this work quantifies the average rate sensitivity of the three materials from an experimental database collected from the literature. Additionally, this work develops rate sensitivity criteria and a complete model for each material that can transition between both deformation regimes and provides several example simulations demonstrating the accuracy of the models and the numerical performance of the composite modeling method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparison of Artificial Compressibility Methods for Coupled Laminar Fluid Flow and Heat Transfer

Artificial compressibility (AC) methods for incompressible flows perturb the mass continuity equation by the time derivative of the pressure field to achieve direct coupling between pressure and velocity without the need for a projection method. Several variations of the original AC method have been proposed (e.g., the entropically damped artificial compressibility [EDAC] method) that can, in some scenarios, improve on the original model. In this paper, we methodically compare the results of the original AC model and the EDAC model for both steady-state and transient problems with coupling between the momentum and energy equations. We demonstrate that the EDAC model and the AC model yield similar numerical results, and the EDAC model is 10%-20% faster than the AC model for a given problem.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
DOE OSTI2023

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Development of simplified model for injection rate prediction of diesel injectors during transient and steady operation

Determining the injection rate of the diesel injector during the transient short injection operations is a critical matter to enhance the combustion process by improving air/fuel mixture formation. In this regard, numerous injection models have been suggested in the literature to predict the injection rate accurately. However, such algorithms have not precisely predicted the rate of injection (ROI) in transient short injection events. In addition, most of those models require many input parameters to characterize hydraulic and mechanical subsystems in the injector, which can make them inconvenient to tune for different injector specifications, in other words not suitable for universal usage. Therefore, our study aims to develop a simple injection model that can predict injection rate precisely in the short and long injection durations considering relatively fewer input parameters. The model traces the pressure variation inside the sac based on the bulk modulus theory to predict nozzle exit velocity, thus the injection rate of the injector. In the algorithm, the bulk modulus of fuel (E) and discharge coefficients of the nozzle hole and needle seat (C d,H and C d,N ) have been considered to vary rather than set constant during the injection process since such flow parameters would have a critical impact on the flow development at the transient stage, thus affecting the injection rate of diesel injectors. Here with this modelling strategy, the accuracy of the new model on ROI prediction was confirmed with acceptable model error for both short and long injection durations by validating against experimental results. To understand the effect of the new model scheme on ROI prediction, ROI results with the new model have been discussed by comparing it with the ROI results predicted with constant flow parameters.

42 ENGINEERING↗
DOE OSTI2022

CFD-trained ANN Model for Approximating Near-occupant Condition in Real-time Simulations

The main drawback of Computational Fluid Dynamics (CFD) simulations has been the time and resource consuming nature which is not suitable for real-time applications. In this work, we first generated numerous CFD models of a given indoor space to obtain airspeed, temperature, and mean radiant temperature near an occupant as training data. Several artificial neural networks (ANN) models were trained using this CFD simulated data to approximate near real-time environmental conditions for a given occupant. This trained ANN model approach is a part of a real-time simulation of building operations using a combination of software and real hardware (HVAC equipment) approaches. The preliminary results suggest that the CFD- generated training data and the trained ANN model can accurately approximate such conditions in a real-time application, a method that has great potential in building simulation and building digital twin areas of research.

Zhang, Yun↗
DOE OSTI2020

Destructive reverse bias pinning in perovskite/silicon tandem solar modules caused by perovskite hysteresis under dynamic shading

In this work, we demonstrate how perovskite hysteresis can result in permanent reductions in power output in perovskite/silicon tandem modules—including irreversible hotspot-induced damage—from only brief periods of shading. We show that reverse bias events in which a perovskite cell is biased above a threshold voltage—which in this work we find to be as low as –1.1 V—produces a temporary reduction in power output that is of sufficient magnitude to keep the cell pinned in reverse bias after the shading event ends. As a hysteretic phenomena, this crucial failure mode may be overlooked by static models of perovskite-based solar cells. Higher reverse bias voltages exacerbate the temporary reduction in short-circuit photocurrent, which is also sensitive to the level of illumination under reverse bias. Numerical device modelling demonstrates that this effect is consistent with our understanding of perovskite hysteresis as a consequence of mobile ion-electron coupling controlling rates of non-radiative recombination over time. Measurements of the dynamic response of single-junction perovskite cells are extrapolated to two-terminal and four-terminal perovskite/silicon tandem module modelling. We validate these models with measurements from an equivalent electronic circuit that represents a two-terminal perovskite–silicon tandem mini module. Two module-level solutions are discussed that address this issue, which includes increasing the number of bypass diodes and choosing better suited silicon bottom cells with higher shunter resistance in two-terminal tandem modules.

14 SOLAR ENERGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Micro-scale Numerical Investigation of Internal and Interfacial Void Defects in Adhesive on Failure Behavior of Adhesively-Bonded Materials with Rough Surfaces

This paper studied the effects of air void defects on the failure behavior of adhesively-bonded materials under global shear via micro-scale computational modeling. The numerical results indicated that interfacial void defects can largely facilitate interfacial debonding of a weaker adherend/adhesive interface under shear. However, this is not the case for a stronger adherend/adhesive interface, showing the reduction on the nominal shear strength of an adhesive joint is mainly due to internal void defects. The reduced bonding performance due to voids can be improved by leveraging an appropriate surface roughness. This preliminary investigation is a first step to better understand the micro-mechanics of interfacial failure in the adhesion of a roughened/patterned adherend via surface modification(s) and an adhesive, and also shows the importance of minimizing interfacial void defects in particular at a weaker adherend/adhesive interface via different techniques.

Qiao, Yao↗