SEARCH · Engineering Papers

Results for “FLOW MODELS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Dynamic Simulation of a Sub-Critical Coal Fired Power Plant

In order to address the demanding operating conditions for remaining coal-fired power plants, a dynamic model and a suite of tools have been developed for studying load cycling and to find optimization opportunities. A sub-critical steam cycle power plant was modeled in a flow-sheet modeling tool, APROS™. The model represented the firing system, economizer, evaporator, superheat, and reheat systems. Four loads from 100% TMCR to 25% TMCR were calibrated and tested such that low-to-high cycling could be studied. The model was run through various load cycles; one of which is presented here. This modeling is a prototype for general use in developing cutting edge controls products and for maximizing economic, low-emissions, and efficient operation of the existing coal power fleet.

Braun, Timothy↗

DOE OSTI2022

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.

42 ENGINEERING↗

DOE OSTI2022

Transport and Retention of Particulate Organic Matter in Sand: Lab Experiments and Modelling

Flow-through reactor (FTR) experiments were conducted at three different downward vertical flow rates to study the transport and retention of particulate organic matter (Chlorella powder) in riverbed sediments. This data package includes data files (.csv) regarding effluent collected from the FTRs during the experiments (time, duration, volume of each sample collection; concentration of nonreactive bromide tracer; concentration of suspended Chlorella determined by absorbance measurements) and sediment (sand) slices collected from the FTRs after the experiments (mass of Chlorella retained in 1 cm depth intervals determined by loss-on-ignition method). More detailed descriptions of data are provided in the data dictionary, and detailed methods are available in the associated MSc thesis listed in ‘Related References.’ The data dictionary and file level meta data (FLMD) are included in the data package as both .csv and .xlsx files. This data package also includes timelapse videos (.mp4) of the experiments, and modelling scripts associated with the experiments. The modelling scripts are included as MATLAB ‘live scripts’ (.mlx, which can be run using MATLAB), and are also included as .pdf and .html files with output plots included (which can be viewed without a MATLAB license and installation). The data files and modelling scripts are organized into a .zip folder for each of the three FTR experiments (three flow rates).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Analytical model for quasi-linear flow response to resonant magnetic perturbation in resistive-inertial and viscous-resistive regimes

In this work, a quasi-linear model for plasma flow response to the resonant magnetic perturbation (RMP) in a tokamak has been rigorously developed in the resistive-inertial (RI) and viscousresistive (VR) regimes purely from the two-field reduced MHD model. Models for plasma response to RMP are commonly composed of equations for the resonant magnetic field response (i.e. the magnetic island) and the torque balance of plasma flow. However, in previous plasma response models, the magnetic island and the torque balance equations are often derived separately from reduced MHD and full MHD equations, respectively. By contrast, in this work we derive both the magnetic island response and the torque balance equations in a quasi-linear model for plasma flow response entirely from a set of two-field reduced MHD equations. Such a quasi-linear model can recover previous plasma flow response models within certain limits and approximations. Furthermore, the physical origins of quasi-linear forces and moments in the flow response equation are also accurately calculated and clarified self-consistently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗

DOE OSTI2022

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modelling of a Flow-Induced Oscillation, Two-Cylinder, Hydrokinetic Energy Converter Based on Experimental Data

The VIVACE Converter consists of cylindrical oscillators in tandem subjected to transverse flow-induced oscillations (FIOs) that can be improved by varying the system parameters for a given in-flow velocity: damping, stiffness, and in-flow center-to-center spacing. Compared to a single isolated cylinder, tandem cylinders can harness more hydrokinetic energy due to synergy in FIO. Experimental and numerical methods have been utilized to analyze the FIO and energy harnessing of VIVACE. A surrogate-based model of two tandem cylinders is developed to predict the power harvesting and corresponding efficiency by introducing a backpropagation neural network. It is then utilized to reduce excessive experimental or computational testing. The effects of spacing, damping, and stiffness on harvested power and efficiency of the established prediction-model are analyzed. At each selected flow velocity, optimization results of power harvesting using the prediction-model are calculated under different combinations of damping and stiffness. The main conclusions are: (1) The surrogate model, built on extensive experimental data for tandem cylinders, can predict the cylinder oscillatory response accurately. (2) Increasing the damping ratio range from 0–0.24 to 0–0.30 is beneficial for improving power efficiency, but has no significant effect on power harvesting. (3) In galloping, a spacing ratio of 1.57 has the highest optimal harnessed power and efficiency compared with other spacing values. (4) Two tandem cylinders can harness 2.01–4.67 times the optimal power of an isolated cylinder. In addition, the former can achieve 1.46–4.01 times the efficiency of the latter. (5) The surrogate model is an efficient predictive tool defining parameters of the Converter for improved energy acquisition.

backpropagation neural network↗

DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Active learning for the design of polycrystalline textures using conditional normalizing flows

Generative modeling has opened new avenues for solving previously intractable materials design problems. However, these new opportunities are accompanied by a drastic increase in the required amount of training data. This is in stark juxtaposition to the high expense and difficulty in curating such large materials datasets. In this work, we propose a novel framework for integrating generative models within an active learning loop. Further, this enables the training of generative models with datasets significantly smaller than what has previously been demonstrated, providing a direct route for their application in data constrained environments. The functionality of this framework is then demonstrated by addressing the challenge of designing polycrystalline textures associated with target anisotropic mechanical properties. The developed protocol exhibited a cost reduction between 14 to 18 times over a randomly sampled experimental design.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Physics-informed CoKriging model of a redox flow battery

Vanadium redox flow batteries (VRFBs) offer the capability to store large amounts of energy cheaply and efficiently, however, there is a need for fast and accurate models of the charge–discharge curve of a VRFB to potentially improve the battery capacity and performance. Here, we develop a multifidelity model for predicting the charge–discharge curve of a VRFB. In the multifidelity model, we use the Physics-informed CoKriging (CoPhIK) machine learning method that is trained on experimental data and constrained by the so-called “zero-dimensional” physics-based model. Here we demonstrate that the model shows good agreement with experimental results and significant improvements over existing zero-dimensional models. We show that the proposed model is robust as it is not sensitive to the input parameters in the zero-dimensional model. We also show that only a small amount of high-fidelity experimental datasets are needed for accurate predictions for the range of considered input parameters, which include current density, flow rate, and initial concentrations.

25 ENERGY STORAGE↗

DOE OSTI2022

Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This document details the application of the Plateau-to-River (P2R) Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the CA. The simulated flow field will support the simulation of fate and transport of contaminants for use in predicting dose as part of the CA. This document describes the changes made to model inputs in order to simulate the predictive flow field for the CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities

We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow

We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Assessment of reverse gun taylor cylinder experimental configuration

Experimental efforts for Taylor-anvil impact tests have often been limited to near room temperature. The ‘Reverse Gun’ method proposed by Gust in 1982 allows for the Taylor impact specimen to be uniformly heated without temperature losses before impact. Through the use of finite element analysis, we explore two topics in this work. First, we examine whether the reverse gun experimental configuration is comparable to the traditional Taylor-anvil setup. Second, we assess the accuracy of several commonly employed flow strength models in terms of their ability to predict the reverse gun experimental results which involve dynamic loading conditions and complex thermo-mechanical coupling. The reverse gun simulations are performed for tantalum targets at initial temperatures in the range 295 K to 1295 K and velocities from 135 m/s to 242 m/s. We show that with suitable care in the modeling of the preheated reverse gun experiments one can make valuable assessments of flow strength models. Given the conditions explored, these observations probe the thermal softening, strain hardening, and strain rate sensitivity of the material.

42 ENGINEERING↗

DOE OSTI2024

A Prefire Approach for Probabilistic Assessments of Postfire Debris‐Flow Inundation

Increases in wildfire activity and rainfall intensification are driving more postfire debris flows (PFDF) in many regions around the world. PFDFs are most common in the first postfire year and may even occur before a fire is fully controlled. This underscores the importance of assessing postfire hazards before a fire starts. Evaluation of PFDF hazards prior to fire can help strategize interventions lessening the negative effects of future fires. However, debris-flow runout and inundation analyses are not routine in PFDF hazard assessments, partially due to time constraints and substantial uncertainties in boundary conditions. Here, we propose a prefire PFDF inundation assessment framework using a debris-flow runout model based on the Herschel-Bulkley (HB) rheology (HEC-RAS v6.1). We constrain model inputs and parameters using Bayesian posterior analysis, rainfall-runoff simulations, and a debris-flow volume model. We use observations from recent PFDF incidents in northern Arizona, USA, to calibrate model components and then apply our prefire inundation assessment framework in a nearby unburned area. Specifically, we (a) identify yield stress as the most influential factor on inundation extent and arrival time in a HB model, (b) establish posterior distributions for model parameters suitable for forward modeling by leveraging uncertainties in field observations, and (c) implement a predictive forward analysis in an area that has not burned recently to evaluate PFDF inundation under several future fire scenarios. This study improves our ability to assess postfire debris-flow hazards before a fire begins and provides guidance for future applications of single-phase rheological models when assessing PFDF hazards.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A Prefire Approach for Probabilistic Assessments of Postfire Debris-Flow Inundation

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗