Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow simulation”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Meso- to microscale modeling of atmospheric stability effects on wind turbine wake behavior in complex terrain

Abstract. Terrain-induced flow phenomena modulate wind turbine performance and wake behavior in ways that are not adequately accounted for in typical wind turbine wake and wind plant design models. In this work, we simulate flow over two parallel ridges with a wind turbine on one of the ridges, focusing on conditions observed during the Perdigão field campaign in 2017. Two case studies are selected to be representative of typical flow conditions at the site, including the effects of atmospheric stability: a stable case where a mountain wave occurs (as in ∼ 50 % of the nights observed) and a convective case where a recirculation zone forms in the lee of the ridge with the turbine (as occurred over 50 % of the time with upstream winds normal to the ridgeline). We use the Weather Research and Forecasting Model (WRF), dynamically downscaled from the mesoscale (6.75 km resolution) to microscale large-eddy simulation (LES) at 10 m resolution, where a generalized actuator disk (GAD) wind turbine parameterization is used to simulate turbine wakes. We compare the WRF–LES–GAD model results to data from meteorological towers, lidars, and a tethered lifting system, showing good qualitative and quantitative agreement for both case studies. Significantly, the wind turbine wake shows different amounts of vertical deflection from the terrain and persistence downstream in the two stability regimes. In the stable case, the wake follows the terrain along with the mountain wave and deflects downwards by nearly 100 m below hub height at four rotor diameters downstream. In the convective case, the wake deflects above the recirculation zone over 40 m above hub height at the same downstream distance. Overall, the WRF–LES–GAD model is able to capture the observed behavior of the wind turbine wakes, demonstrating the model's ability to represent wakes over complex terrain for two distinct and representative atmospheric stability classes, and, potentially, to improve wind turbine siting and operation in hilly landscapes.

17 WIND ENERGY↗
DOE OSTI2024

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Microphysics in Simulations of Nonprecipitating Cumulus

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. Overall, the results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

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

Scalable Risk Assessment of Rare Events in Power Systems With Uncertain Wind Generation and Loads

Risk assessment of rare events has become increasingly important in power system planning and operation with the increasing integration of renewable energy and the presence of system uncertainties. However, quantifying the risk posed by rare events via the traditional method, i.e., Monte Carlo sampling (MCS), incurs substantial computational expense stemming from the vast ensemble of power flow simulations. To accelerate the assessment, this paper proposes a Deep Neural Network (DNN)-kernelized vector-valued Gaussian Process (VVGP) approach with excellent computational efficiency while maintaining high accuracy. Consequently, serving as a surrogate model for the power flow solver, the DNN-kernelized VVGP enables significantly faster but accurate risk assessment compared to the power flow solver. The developed surrogate model evaluates low-order N - k events that contain more than 90% instances by adeptly capturing the topological features while the high-order N - k events are assessed via a power flow solver, thereby striking a balance between computational efficiency and uncertainty quantification accuracy. Moreover, the model incorporates a Support Vector Machine (SVM) classifier to resample concerning low-probability tail events to counteract the biases potentially introduced during the DNN-kernelized VVGP evaluations. Simulations conducted on the modified IEEE 24-bus, 118-bus, and European 1354-bus systems demonstrate that the proposed method maintains the accuracy benchmark set by MCS while significantly reducing computational demands in large-scale power systems as compared to other state-of-the-art methods.

17 WIND ENERGY↗
DOE OSTI2014

Internal air flow analysis of a bladeless micro aerial vehicle hemisphere body using computational fluid dynamic

This paper explain the analysis of internal air flow velocity of a bladeless vertical takeoff and landing (VTOL) Micro Aerial Vehicle (MAV) hemisphere body. In mechanical design, before produce a prototype model, several analyses should be done to ensure the product's effectiveness and efficiency. There are two types of analysis method can be done in mechanical design; mathematical modeling and computational fluid dynamic. In this analysis, I used computational fluid dynamic (CFD) by using SolidWorks Flow Simulation software. The idea came through to overcome the problem of ordinary quadrotor UAV which has larger size due to using four rotors and the propellers are exposed to environment. The bladeless MAV body is designed to protect all electronic parts, which means it can be used in rainy condition. It also has been made to increase the thrust produced by the ducted propeller compare to exposed propeller. From the analysis result, the air flow velocity at the ducted area increased to twice the inlet air. This means that the duct contribute to the increasing of air velocity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Stress management for 3D through-silicon-via stacking technologies - The next frontier -

The status of the development of a Design-for-Stress simulation flow that captures the stress effects in packaged 3D-stacked Si products like integrated circuits (ICs) using advanced via-middle Through Si Via technology is outlined. The next set of challenges required to proliferate the methodology and to deploy it for making and dispositioning real Si product decisions are described here. These include the adoption and support of a Process Design Kit (PDK) that includes the relevant material properties, the development of stress simulation methodologies that operate at higher levels of abstraction in a design flow, and the development and adoption of suitable models required to make real product reliability decisions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Separating physical impacts from natural variability using piggybacking technique

Abstract. In a chaotic system, like moist convection, it is difficult to separate the impact of a physical process from effects of natural variability. This is because modifying even a small element of the system physics typically leads to a different system evolution and it is difficult to tell whether the difference comes from the physical impact or it merely represents a different flow realization. This paper discusses a relatively simple and computationally efficient modelling methodology that allows separation of the two. The methodology is referred to as the piggybacking approach. The idea is to use two sets of thermodynamic variables (the temperature, water vapor, and all aerosol, cloud, and precipitation variables) in a single cloud simulation. The two sets differ in a specific element of the physics, such as aerosol properties, microphysics parameterization, large-scale forcing, environmental profiles, etc. One thermodynamic set is coupled to the dynamics and drives the simulated flow, and the other set piggybacks the flow, that is, thermodynamic variables are carried by the flow but they do not affect it. By switching the two sets (i.e. the set driving the simulation becomes the piggybacking one, and vice versa), the impact on the cloud dynamics can be evaluated. This paper provides details of the method and reviews results of its application to such problems as the postulated deep convection invigoration in polluted environments, the impact of changes in environmental profiles (e.g., due to climate change) on convective dynamics, and the role of cloud-layer heterogeneities for shallow convective cloud field evolution. Prospects for applying piggybacking technique to other areas of atmospheric simulation (e.g., weather prediction or geoengineering) are also mentioned.

Grabowski, Wojciech W.↗
DOE OSTI2016

An overview of TOUGH-based geomechanics models

After the initial development of the first TOUGH-based geomechanics model 15 years ago based on linking TOUGH2 multiphase flow simulator to the FLAC3D geomechanics simulator, at least 15 additional TOUGH-based geomechanics models have appeared in the literature. This development has been fueled by a growing demand and interest for modeling coupled multiphase flow and geomechanical processes related to a number of geoengineering applications, such as in geologic CO 2 sequestration, enhanced geothermal systems, unconventional hydrocarbon production, and most recently, related to reservoir stimulation and injection-induced seismicity. This paper provides a short overview of these TOUGH-based geomechanics models, focusing on some of the most frequently applied to a diverse set of problems associated with geomechanics and its couplings to hydraulic, thermal and chemical processes.

58 GEOSCIENCES↗
DOE OSTI2023

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗
NASA NTRS1983

A chimera grid scheme

A mesh system composed of multiple overset body-conforming grids is described for adapting finite-difference procedures to complex aircraft configurations. In this so-called 'chimera mesh,' a major grid is generated about a main component of the configuration and overset minor grids are used to resolve all other features. Methods for connecting overset multiple grids and modifications of flow-simulation algorithms are discussed. Computational tests in two dimensions indicate that the use of multiple overset grids can simplify the task of grid generation without an adverse effect on flow-field algorithms and computer code complexity.

Steger, J. L.↗
NASA NTRS2011

Multidisciplinary Modeling Software for Analysis, Design, and Optimization of HRRLS Vehicles

The concept for Highly Reliable Reusable Launch Systems (HRRLS) under the NASA Hypersonics project is a two-stage-to-orbit, horizontal-take-off / horizontal-landing, (HTHL) architecture with an air-breathing first stage. The first stage vehicle is a slender body with an air-breathing propulsion system that is highly integrated with the airframe. The light weight slender body will deflect significantly during flight. This global deflection affects the flow over the vehicle and into the engine and thus the loads and moments on the vehicle. High-fidelity multi-disciplinary analyses that accounts for these fluid-structures-thermal interactions are required to accurately predict the vehicle loads and resultant response. These predictions of vehicle response to multi physics loads, calculated with fluid-structural-thermal interaction, are required in order to optimize the vehicle design over its full operating range. This contract with ResearchSouth addresses one of the primary objectives of the Vehicle Technology Integration (VTI) discipline: the development of high-fidelity multi-disciplinary analysis and optimization methods and tools for HRRLS vehicles. The primary goal of this effort is the development of an integrated software system that can be used for full-vehicle optimization. This goal was accomplished by: 1) integrating the master code, FEMAP, into the multidiscipline software network to direct the coupling to assure accurate fluid-structure-thermal interaction solutions; 2) loosely-coupling the Euler flow solver FEFLO to the available and proven aeroelasticity and large deformation (FEAP) code; 3) providing a coupled Euler-boundary layer capability for rapid viscous flow simulation; 4) developing and implementing improved Euler/RANS algorithms into the FEFLO CFD code to provide accurate shock capturing, skin friction, and heat-transfer predictions for HRRLS vehicles in hypersonic flow, 5) performing a Reynolds-averaged Navier-Stokes computation on an HRRLS configuration; 6) integrating the RANS solver with the FEAP code for coupled fluid-structure-thermal capability; and 7) integrating the existing NASA SRGULL propulsion flow path prediction software with the FEFLO software for quasi-3D propulsion flow path predictions, 8) improving and integrating into the network, an existing adjoint-based design optimization code.

Spradley, Lawrence W.↗
NASA NTRSDate not supplied

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

As computing technology continues to improve, simulating unsteady aerodynamic environments on launch vehicle geometries becomes more viable. For NASA’s Space Launch System, the flowfield is characterized by large regions of unsteady turbulent wall-bounded flow. Simulating these unsteady wall-bounded environments is of interest to the SLS Aerodynamics team. Wall-modeled LES methods are currently the best approach for simulating these environments in engineering applications. Less-expensive computationally than fully-resolved methods while maintaining many of the benefits in simulating unsteady environments. Before applying these methods blindly to large geometries, best practices can first be applied and understood on smaller canonical problems.

Bryson Frank↗
DOE OSTI2020

Six Degrees of Freedom (6DOF) Simulations of Supersonic Fragment Trajectories

In this report a process using existing technologies at Sandia National Laboratories (SNL) to simulate the six degrees-of-freedom (6DOF) trajectories of explosive fragments is described and tested. First, aerodynamic forces and moments as functions of orientation are computed using the SIERRA/Aero supersonic flow solver. The forces and moments are normalized and tabulated in a database. Second, this the aerodynamic coefficient database is imported into a 6DOF rigid body dynamics solver in order to compute the resulting trajectories. The supersonic flow simulations are tested for simple geometries and show good agreement with literature values. The simulation procedure is then demonstrated for an example fragment. The results of the example fragment indicate that the distance traveled in the early ight (from 2.5 km/s until decreasing down to 1 km/s) varies widely depending on the initial orientations. The fragment trajectory distribution and steady tumbling rate is explored. The study indicates that a 6DOF analysis will yield information about a spread of possible trajectories, while using an average drag coefficient can only represent the most likely trajectory.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1984

Metric-discontinuous zonal grid calculations using the Osher scheme

Computations on zonal grids - in particular, grids with metric discontinuities resulting from the interspersion of highly clustered regions with coarse regions - are possible using a fully conservative form of the Osher upwind scheme. These zonal grids can result from an abrupt clustering of points near solution discontinuities or near other flow features that require improved resolution. The zonal approach is shown to capture shocks with almost 'shock-fitting' quality but with minimal effort. Results for inviscid flow, including quasi-one-dimensional nozzle flow, supersonic flow over a cylinder, and blast-wave diffraction by a ramp, are presented. These calculations demonstrate the powerful capabilities of the Osher scheme used in conjunction with zonal grids in simulating flow fields with complex shock patterns.

Rai, M. M.↗
NASA NTRS2002

Launch Vehicle Sizing Benefits Utilizing Main Propulsion System Crossfeed and Project Status

To meet the goals for a next generation Reusable Launch Vehicle (RLV), a unique propulsion feed system concept was identified using crossfeed between the booster and orbiter stages that could reduce the Two-Stage-to-Orbit (TSTO) vehicle weight and Design, Development, Test and Evaluation (DDT&E) costs by approximately 25%, while increasing safety and reliability. The Main Propulsion System (MPS) crossfeed water demonstration test program addresses all activities required to reduce the risks for the MPS crossfeed system from a Technology Readiness Level (TRL) of 2 to 4 by the completion of testing and analysis by June 2003. During the initial period, that ended in March 2002, a subscale water flow test article was defined. Procurement of a subscale crossfeed check valve was initiated and the specifications for the various components were developed. The fluid transient and pressurization analytical models were developed separately and successfully integrated. The test matrix for the water flow test was developed to correlate the integrated model. A computational fluid dynamics (CFD) model of the crossfeed check valve was developed to assess flow disturbances and internal flow dynamics. Based on the results, the passive crossfeed system concept was very feasible and offered a safe system to be used in an RLV architecture. A water flow test article was designed to accommodate a wide range of flows simulating a number of different types of propellant systems. During the follow-on period, the crossfeed system model will be further refined, the test article will be completed, the water flow test will be performed, and finally the crossfeed system model will be correlated with the test data. This validated computer model will be used to predict the full-scale vehicle crossfeed system performance.

Chandler, Frank↗
DOE OSTI2015

Assessment of Saturation Changes Around Cavities and the Impact on Seismic Waves.

According to theory, seismic methods should be able to identify cavities. The high impedance contrast between the rock (fast velocity, high density) and the air (slow velocity, low density) should produce strong reflections and refractions off cavity walls. In practice, however, cavities are not always easily detected. Changes in pore saturation, fracturing, and stress near the wall of the cavity may decrease the impedance contrast and diffuse the rock/air interface. To determine the impact of these effects on seismic waves, we first aim to understand the hydrological processes in the medium surrounding the cavity, a tunnel in this case. We simulate flow through the vadose zone to provide insight into how a cavity alters saturation. Next we examine how this saturation affects seismic velocities. Finally, we use a wave propagation code to simulate seismic waves through the saturation models. Results are presented from hydrology models for three cases: the tunnel far above the water table, the tunnel just above the water table and within the capillary fringe, and the tunnel below the water table with the tunnel being drained. We also present the results of seismic wave propagation through these models to determine the impact these effects have on seismic waves.

58 GEOSCIENCES↗
DOE OSTI2023

MACHINE LEARNING-ENABLED PREDICTION OF TRANSIENT INJECTION MAP IN AUTOMOTIVE INJECTORS WITH UNCERTAINTY QUANTIFICATION

Accurate prediction of injection profiles is a critical aspect of linking injector operation with engine performance and emissions. However, highly resolved injector simulations can take one to two weeks of wall-clock time, which is incompatible with engine design cycles with desired turnaround times of less than a day. Hence, it is important to reduce the time-to-solution of the internal flow simulations by several orders of magnitude to make it compatible with engine simulations. This work demonstrates a data-driven approach for tackling the computational overhead of injector simulations, whereby the transient injection profiles are emulated for a side-oriented, single-hole diesel injector using a Bayesian machine-learning framework. First, an interpretable Bayesian learning strategy was employed to understand the effect of design parameters on the total void fraction field. Then, autoencoders are utilized for efficient dimensionality reduction of the flowfields. Gaussian process models are finally used to predict the spatiotemporal void fraction field at the injector exit for unknown operating conditions. The Gaussian process models produce principled uncertainty estimates associated with the emulated flowfields, which provide the engine designer with valuable information of where the data-driven predictions can be trusted in the design space. The Bayesian flowfield predictions are compared with the corresponding predictions from a deep neural network, which has been transfer-learned from static needle simulations from a previous work by the authors. The emulation framework can predict the void fraction field at the exit of the orifice within a few seconds, thus achieving a speed-up factor of up to 38 x 10(6) over the traditional simulation-based approach of generating transient injection maps.

machine learning↗
DOE OSTI2014

Sandia National Laboratories Environmental Fluid Dynamics Code V. 1 0.0 (Beta)

The DOE has funded Sandia National Labs (SNL) to develop an open-source modeling tool to guide the design and layout of marine hydrokinetic (MHK) arrays to maximize power production while minimizing environmental effects. This modeling framework simulates flows through and around MHK arrays while quantifying environmental responses. As an augmented version of US EPA's Environmental Fluid Dynamics Code (EFDC), SNL-EFDC includes: (1) a new module that simulates energy conversion (momentum withdrawal) by MHK devices with commensurate changes in the turbulent kinetic energy and its dissipation rate, (2) new, advanced sediment dynamics routines, and (3) augmented water quality modules.

Jones, Craig↗