Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

The Role of Numerical Simulation in the Design of Stimulation and Circulation Experiments for the EGS Collab Project

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additional plays an essential role in designing these meso-scale experiments. This paper describes specific calculations and numerical simulation approaches conducted in support of the designs for the stimulation and circulation experiments, and their monitoring, including those for predicting seismic energies during stimulation, stimulation pressure requirements, impacts of borehole orientation on fracture development and trajectories, influence of notch shapes on hydraulic fracture initiation, and pressure limits on fluid circulation to avoid fracture growth.

Numerical Simulation, EGS Collab, EGS, Sanford Und↗
DOE OSTI2015

Quantitative Determination of Grain Boundary Recombination Velocity in CdTe by Combination of Cathodoluminescence Measurements and Numerical Simulations

We developed a 2D numerical model simulating cathodoluminescence (CL) measurements in CdTe. Using this model we analyze how various material parameters impact the CL contrast and intensity observed in the measured signal, and determine if and when we can accurately determine the value of grain boundary recombination rate. In addition to grain boundary (GB) recombination, the grain size and its ratio to the carrier diffusion length impact the results of the measurement. Holding the grain interior and GB recombination rates constant, we find that as the grain size increases and becomes larger than the diffusion length, the observed CL contrast is larger. In a small grain size material the surface recombination lowers the overall intensity of the signal, but does not impact the observed contrast significantly. In a large grain size material, high surface recombination velocity can lower the observed contrast in a measurement. This model in combination with an experiment is used to quantify the grain boundary recombination velocity in polycrystalline CdTe before and after the CdCl2 treatment.

14 SOLAR ENERGY↗
DOE OSTI2017

Numerical simulation and dynamical response of a moored hydrokinetic turbine operating in the wake of an upstream turbine for control design

Herein numerical simulation of a downstream hydrokinetic turbine operating in the wake of an upstream turbine for feedback control design is presented. Wake effects from an upstream turbine are quantified in terms of wake velocity and amplified turbulence levels. These effects are integrated in an in-stream hydrokinetic turbine numerical simulation that utilizes a Blade Element Momentum approach with a dynamic wake inflow model. Simulations are carried out on a fixed turbine model to simulate operation in river or tidal channels with conventional foundations, as well as on a compliantly moored turbine model such as those designed to operate in open ocean currents.

13 HYDRO ENERGY↗
DOE OSTI2024

Direct Numerical Simulation of Involute Channel Turbulence

A direct numerical simulation (DNS) study was performed on turbulent flow in the high flux isotope reactor involute channel geometry to develop a numerical database and determine the differences compared with a flat parallel channel. The varying channel curvature along the walls was studied for differences in mean profiles. Parameters of interest include streamwise velocity, turbulent kinetic energy (TKE), and turbulence dissipation rate, as well as Reynolds stresses and turbulence transport terms. Profile sampling was carried out at 10 locations along the span of the involute. Additional DNS studies were performed on smaller domains of comparable curvature to the involute domain: a high curvature channel (high circular), a low curvature channel (low circular), and a flat channel (flat). Here, each of these four cases was compared against each other and to other DNS studies performed on parallel flows. The results indicate that the bulk involute channel flow does not differ significantly from a flat parallel channel flow and that the curvature of the walls does not significantly alter the mean flow parameters. However, the regions of the involute channel near the side walls exhibit relatively low magnitude twin recirculation structures driven toward the side walls from the centerline of the channel, which warrants further study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗
DOE OSTI2019

Direct numerical simulation and Reynolds-averaged Navier-Stokes modeling of the sudden viscous dissipation for multicomponent turbulence

Simulations of a turbulent multicomponent fluid mixture undergoing isotropic deformations are carried out to investigate the sudden viscous dissipation. This dissipative mechanism was originally demonstrated using simulations of an incompressible single-component fluid [S. Davidovits and N. J. Fisch, Phys. Rev. Lett. 116, 105004 (2016)]. By accounting for the convective and diffusive transfer of various species, the current work aims to increase the physical fidelity of previous simulations and their relevance to inertial confinement fusion applications. Direct numerical simulations of the compressed fluid show that the sudden viscous dissipation of turbulent kinetic energy is unchanged from the single-component scenario. More importantly, the simulations demonstrate that the mass fraction variance and covariance for the various species also exhibit a sudden viscous decay. Reynolds-averaged Navier-Stokes simulations were carried out using the $k - l$ model to assess its ability to reproduce the sudden viscous dissipation. Results show that the standard $k - l$ formulation does not capture the sudden decay of turbulent kinetic energy, mass-fraction variance, and mass-fraction covariance for simulations with various compression and expansion rates, or different exponents for the power-law model of viscosity. A new formulation of the $k - l$ model that is based on previous improvements to the $k - ε$ family of models is proposed, which leads to consistently good agreement with the direct numerical simulations for all the isotropic deformations under consideration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere

We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three-dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a four-dimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipole-dipole potential has two singularities which correspond to the following boundary conditions: when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects.

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

Planning Irreversible Electroporation in the Porcine Kidney: Are Numerical Simulations Reliable for Predicting Empiric Ablation Outcomes?

PurposeNumerical simulations are used for treatment planning in clinical applications of irreversible electroporation (IRE) to determine ablation size and shape. To assess the reliability of simulations for treatment planning, we compared simulation results with empiric outcomes of renal IRE using computed tomography (CT) and histology in an animal model.MethodsThe ablation size and shape for six different IRE parameter sets (70–90 pulses, 2,000–2,700 V, 70–100 µs) for monopolar and bipolar electrodes was simulated using a numerical model. Employing these treatment parameters, 35 CT-guided IRE ablations were created in both kidneys of six pigs and followed up with CT immediately and after 24 h. Histopathology was analyzed from postablation day 1.ResultsAblation zones on CT measured 81 ± 18 % (day 0, p ≤ 0.05) and 115 ± 18 % (day 1, p ≤ 0.09) of the simulated size for monopolar electrodes, and 190 ± 33 % (day 0, p ≤ 0.001) and 234 ± 12 % (day 1, p ≤ 0.0001) for bipolar electrodes. Histopathology indicated smaller ablation zones than simulated (71 ± 41 %, p ≤ 0.047) and measured on CT (47 ± 16 %, p ≤ 0.005) with complete ablation of kidney parenchyma within the central zone and incomplete ablation in the periphery.ConclusionBoth numerical simulations for planning renal IRE and CT measurements may overestimate the size of ablation compared to histology, and ablation effects may be incomplete in the periphery.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Scalable foundation models for numerical simulations on HPC platforms

In recent years, foundation models (FMs) have begun to reshape numerical simulations on high-performance computing (HPC) platforms. These large, pre-trained AI models enable rapid predictions across a broad range of physical domains, including Earth system modeling, fluid dynamics, materials science, as well as complex multi-modal simulations in aerospace engineering and fusion research. By training on diverse datasets, FMs learn intricate relationships and underlying physical behavior while also enabling the quantification of uncertainty in their predictions. This capability allows simulations that once required days of numerical calculation to be completed in minutes (FM inference), supporting real-time design optimization, uncertainty-aware decision making, and more comprehensive exploration of complex scenarios.

AI↗
DOE OSTI2018

Numerical simulation of runaway electrons: 3-D effects on synchrotron radiation and impurity-based runaway current dissipation

Numerical simulations of runaway electrons (REs) with a particular emphasis on orbit dependent effects in 3-D magnetic fields are presented. The simulations were performed using the recently developed Kinetic Orbit Runaway electron Code (KORC) that computes the full-orbit relativistic dynamics in prescribed electric and magnetic fields including radiation damping and collisions. The two main problems of interest are synchrotron radiation and impurity-based RE dissipation. Synchrotron radiation is studied in axisymmetric fields and in 3-D magnetic configurations exhibiting magnetic islands and stochasticity. For passing particles in axisymmetric fields, neglecting orbit effects might underestimate or overestimate the total radiation power depending on the direction of the radial shift of the drift orbits. For trapped particles, the spatial distribution of synchrotron radiation exhibits localized “hot” spots at the tips of the banana orbits. In general, the radiation power per particle for trapped particles is higher than the power emitted by passing particles. The spatial distribution of synchrotron radiation in stochastic magnetic fields, obtained using the MHD code NIMROD, is strongly influenced by the presence of magnetic islands. 3-D magnetic fields also introduce a toroidal dependence on the synchrotron spectra, and neglecting orbit effects underestimates the total radiation power. In the presence of magnetic islands, the radiation damping of trapped particles is larger than the radiation damping of passing particles. Results modeling synchrotron emission by RE in DIII-D quiescent plasmas are also presented. The computation uses EFIT reconstructed magnetic fields and RE energy distributions fitted to the experimental measurements. Qualitative agreement is observed between the numerical simulations and the experiments for simplified RE pitch angle distributions. However, it is noted that to achieve quantitative agreement, it is necessary to use pitch angle distributions that depart from simplified 2-D Fokker-Planck equilibria. Finally, using the guiding center orbit model (KORC-GC), a preliminary study of pellet mitigated discharges in DIII-D is presented. The dependence of RE energy decay and current dissipation on initial energy and ionization levels of neon impurities is studied. Furthermore, the computed decay rates are within the range of experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗
DOE OSTI2016

Estimation of heat transfer coefficients for biomass particles by direct numerical simulation using microstructured particle models in the Laminar regime

Here, direct numerical simulation of convective heat transfer from hot gas to isolated biomass particle models with realistic morphology and explicit microstructure was performed over a range of conditions with laminar flow of hot gas (500 degrees C). Steady-state results demonstrated that convective interfacial heat transfer is dependent on the wood species. The computed heat transfer coefficients were shown to vary between the pine and aspen models by nearly 20%. These differences are attributed to the species-specific variations in the exterior surface morphology of the biomass particles. We also quantify variations in heat transfer experienced by the particle when positioned in different orientations with respect to the direction of fluid flow. These results are compared to previously reported heat transfer coefficient correlations in the range of 0.1 < Pr < 1.5 and 10 < Re < 500. Comparison of these simulation results to correlations commonly used in the literature (Gunn, Ranz-Marshall, and Bird-Stewart-Lightfoot) shows that the Ranz-Marshall (sphere) correlation gave the closest h values to our steady-state simulations for both wood species, though no existing correlation was within 20% of both species at all conditions studied. In general, this work exemplifies the fact that all biomass feedstocks are not created equal, and that their species-specific characteristics must be appreciated in order to facilitate accurate simulations of conversion processes.

09 BIOMASS FUELS↗
DOE OSTI2021

Grid resolution requirement for resolving rare and high intensity wall-shear stress events in direct numerical simulations

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)]. Here, the present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Re τ = 180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Re τ = 10000 channel flow, one resolves about 90%–95% wall-shear stress events.

42 ENGINEERING↗
DOE OSTI2020

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗
DOE OSTI2014

Numerical simulations of epitaxial growth process in MOVPE reactor as a tool for design of modern semiconductors for high power electronics

In the present study numerical simulations of epitaxial growth of gallium nitride in Metal Organic Vapor Phase Epitaxy reactor AIX-200/4RF-S is addressed. Epitaxial growth means crystal growth that progresses while inheriting the laminar structure and the orientation of substrate crystals. One of the technological problems is to obtain homogeneous growth rate over the main deposit area. Since there are many agents influencing reaction on crystal area such as temperature, pressure, gas flow or reactor geometry, it is difficult to design optimal process. According to the fact that it's impossible to determine experimentally the exact distribution of heat and mass transfer inside the reactor during crystal growth, modeling is the only solution to understand the process precisely. Numerical simulations allow to understand the epitaxial process by calculation of heat and mass transfer distribution during growth of gallium nitride. Including chemical reactions in numerical model allows to calculate the growth rate of the substrate and estimate the optimal process conditions for obtaining the most homogeneous product.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Tomography based numerical simulation of the demagnetizing field in soft magnetic composites

The magneto-static behaviour of soft magnetic composites (SMCs) is investigated using tomography based direct numerical simulation. The microgeometry crucially affects the magnetic properties of the composite since a geometry dependent demagnetizing field is established inside the composite, which lowers the magnetic permeability. We determine the magnetic field information inside the SMC using direct numerical simulation of the magnetic field based on high resolution micro-computed tomography data of the SMC's microstructure as well as artificially generated data made of statistically homogeneous systems of identical fully penetrable spheres and prolate spheroids. Quasi-static electromagnetic behaviour and linear material response are assumed. The 3D magnetostatic Maxwell equations are solved using Whitney finite elements. Simulations show that clustering and percolation behaviour determine the demagnetizing factor of SMCs rather than the particle shape. The demagnetizing factor correlates with the slope of a 2-point probability function at its origin, which is related to the specific surface area of the SMC. Comparison with experimental results indicates that the relatively low permeability of SMCs cannot be explained by demagnetizing effects alone and suggests that the permeability of SMC particles has to be orders of magnitude smaller than the bulk permeability of the particle material.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Numerical simulation of tectonic plates motion and seismic process in Central Asia

An evolutionary approach proposed in [1, 2] combining the achievements of traditional macroscopic theory of solid mechanics and basic ideas of nonlinear dynamics is applied in a numerical simulation of present-day tectonic plates motion and seismic process in Central Asia. Relative values of strength parameters of rigid blocks with respect to the soft zones were characterized by the δ parameter that was varied in the numerical experiments within δ = 1.1–1.8 for different groups of the zonal-block divisibility. In general, the numerical simulations of tectonic block motion and accompanying seismic process in the model geomedium indicate that the numerical solutions of the solid mechanics equations characterize its deformation as a typical behavior of a nonlinear dynamic system under conditions of self-organized criticality.

58 GEOSCIENCES↗
DOE OSTI2016

Systematic parameter study of a nonlinear electromagnetic energy harvester with matched magnetic orientation: Numerical simulation and experimental investigation

Here, this paper reports the systematic parameter study of a tristable nonlinear electromagnetic energy harvester for ambient low-frequency vibration. Numerical simulations and experimental investigations are performed on the harvester which consists of a cantilever beam, a tip coil, two tip magnets and two external side magnets. The external side magnets are deployed symmetrically along a concave surface parallel to the trajectory of the cantilever tip with a controllable distance so that the magnetic orientation of the tip magnets are matched with that of the side magnets. Therefore, instead of the ternary position parameters (d, h, α), a binary parameters pair (d 0 , d) is used to characterize the position of the side magnets and the performance of the energy harvester. The magnetic force and magnetic field on the cantilever tip therefore depend on the relative distance in the tip displacement direction between the tip magnets and side magnets, but is independent of the position of the side magnets on the concave surface. The magnetic force (field)-distance relationship is measured experimentally and curve fitted to obtain explicit expressions, in order to characterize the magnetic force (field) when the side magnets are placed at varied positions along the concave surface. Numerical simulation is, then, performed to predict the electromagnetic voltage output and the bandwidth of the energy harvester. The simulation results coincided with the measured data. Significant broadband response is obtained experimentally and the maximum RMS power output is 40.2 mW at 0.45 g of excitation. Furthermore, the proposed structure showcasing the matched magnetic orientation is characterized by the binary parameters pair (d 0 , d) and the systematic parametric approach could contribute to the design and study of nonlinear broadband energy harvesters.

30 DIRECT ENERGY CONVERSION↗