Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Inverse Load Calculation of Wind Turbine Support Structures - A Numerical Verification Using the Comprehensive Simulation Code FAST: Preprint (Revised)

Physically measuring the dynamic responses of wind turbine support structures enables the calculation of the applied loads using an inverse procedure. In this process, inverse means deriving the inputs/forces from the outputs/responses. This paper presents results of a numerical verification of such an inverse load calculation. For this verification, the comprehensive simulation code FAST is used. FAST accounts for the coupled dynamics of wind inflow, aerodynamics, elasticity and turbine controls. Simulations are run using a 5-MW onshore wind turbine model with a tubular tower. Both the applied loads due to the instantaneous wind field and the resulting system responses are known from the simulations. Using the system responses as inputs to the inverse calculation, the applied loads arecalculated, which in this case are the rotor thrust forces. These forces are compared to the rotor thrust forces known from the FAST simulations. The results of these comparisons are presented to assess the accuracy of the inverse calculation. To study the influences of turbine controls, load cases in normal operation between cut-in and rated wind speed, near rated wind speed and between ratedand cut-out wind speed are chosen. The presented study shows that the inverse load calculation is capable of computing very good estimates of the rotor thrust. The accuracy of the inverse calculation does not depend on the control activity of the wind turbine.

17 WIND ENERGY↗
DOE OSTI2022

Can reanalysis products outperform mesoscale numerical weather prediction models in modeling the wind resource in simple terrain?

Mesoscale numerical weather prediction (NWP) models are generally considered more accurate than reanalysis products in characterizing the wind resource at heights of interest for wind energy, given their finer spatial resolution and more comprehensive physics. However, advancements in the latest ERA-5 reanalysis product motivate an assessment on whether ERA-5 can model wind speeds as well as a state-of-the-art NWP model – the Weather Research and Forecasting (WRF) Model. We consider this research question for both simple terrain and offshore applications. Specifically, we compare wind profiles from ERA-5 and the preliminary WRF runs of the Wind Integration National Dataset (WIND) Toolkit Long-term Ensemble Dataset (WTK-LED) to those observed by lidars at a site in Oklahoma, United States, and in a United States Atlantic offshore wind energy area. We find that ERA-5 shows a significant negative bias (~-1ms-1) at both locations, with a larger bias at the land-based site. WTK-LED-predicted wind speed profiles show a limited negative bias (~-0.5ms-1) offshore and a slight positive bias (~+0.5ms-1) at the land-based site. On the other hand, we find that ERA-5 outperforms WTK-LED in terms of the centered root-mean-square error (cRMSE) and correlation coefficient, for both the land-based and offshore cases, in all atmospheric stability conditions. We find that WTK-LED's higher cRMSE is caused by its tendency to overpredict the amplitude of the wind speed diurnal cycle. At the land-based site, this is partially caused by wind plant wake effects not being accurately captured by WTK-LED.

17 WIND ENERGY↗
NASA NTRS1985

Numerical calculation of subsonic jets in crossflow with reduced numerical diffusion

A series of calculations are reported for two, subsonic jet in crossflow geometries. The parametric variation examined are the lateral spacing of a row of jets. The first series of calculations corresponds to a widely space jet geometry, S/D = 4, and the second series corresponds to closely spaced jets, S/D = 2. The calculations are done with alternate differencing schemes to illustrate the impact of numerical diffusion. The calculated jet trajectories agreed well with experimental data in the widely spaced jet geometry, but not in the closely spaced geometry.

Claus, R. W.↗
NASA NTRS1985

Numerical calculation of subsonic jets in crossflow with reduced numerical diffusion

A series of calculations are reported for two, subsonic jet in crossflow geometries. The parametric variation examined are the lateral spacing of a row of jets. The first series of calculations corresponds to a widely space jet geometry, S/D = 4, and the second series corresponds to closely spaced jets, S/D = 2. The calculations are done with alternate differencing schemes to illustrate the impact of numerical diffusion. The calculated jet trajectories agreed well with experimental data in the widely spaced jet geometry, but not in the closely spaced geometry.

Claus, R. W.↗
NASA NTRS1986

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗
DOE OSTI2016

Numerical simulations of depressurization-induced gas production from gas hydrates using 3-D heterogeneous models of L-Pad, Prudhoe Bay Unit, North Slope Alaska

Gas production potential is estimated using a three-dimensional reservoir model based on gas hydrate deposits located in the Prudhoe Bay region of the Alaska's North Slope. The model incorporates two hydrate-bearing sand units using detailed reservoir geological and structural information obtained from past and recent drilling programs. Geostatistical porosity models, conditioned to log data from 78 wells drilled in the vicinity of the Prudhoe Bay “L-Pad,” were developed, providing 3D heterogeneity in porosity, porosity-dependent hydrate saturation, and intrinsic permeability. The simulations utilize both vertical and inclined wellbores to induce depressurization of the reservoir at a constant bottom-hole pressure. The results show the superior performance of the inclined well design. Average gas production rates during the first five years were ~6.0 × 10 4 ST m 3 /day (~2.1 MMSCF/day) and ~2.7 × 10 4 ST m 3 /day (~1.0 MMSCF/day) for incline and vertical wells, respectively. After 30 years 5.3 × 10 8 - 5.7 × 10 8 ST m 3 (18.6–20.2 BSCF) and 6.2 × 10 8 - 6.4 × 10 8 ST m 3 (22.0–22.7 BSCF) gas were produced using the vertical and inclined well configurations, respectively. The analysis reveals that 2D reservoir models with homogeneous representations for porosity and hydrate saturation significantly underestimate production potential. The heterogeneity implemented in this work provides complex porous network and preferential pathways for mobile phase flow to a producing well. Consequently, no secondary hydrate formation around a wellbore is predicted in contrast to models that utilize uniform porosity and hydrate saturation distributions.

03 NATURAL GAS↗
DOE OSTI2013

Numerical investigation of electron self-injection in the nonlinear bubble regime

The process of electron self-injection in the nonlinear bubble wake generated by a short and intense laser pulse propagating in a uniform underdense plasma is studied by means of fully self-consistent particle-in-cell simulations and test-particle simulations. We consider a wake generated by a non-evolving laser driver traveling with a prescribed velocity, which then sets the structure and the velocity of the wake, so the injection dynamics is decoupled from driver evolution, but a realistic structure for the wakefield is retained. We show that a threshold for self-injection into a non-evolving bubble wake exists, and we characterize the dependence of the self-injection threshold on laser intensity, wake velocity, and plasma temperature for a range of parameters of interest for current and future laser-plasma accelerators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The impact of collisionality, FLR, and parallel closure effects on instabilities in the tokamak pedestal: Numerical studies with the NIMROD code

The extended-MHD NIMROD code [C. R. Sovinec and J. R. King, J. Comput. Phys. 229, 5803 (2010)] is verified against the ideal-MHD ELITE code [H. R. Wilson et al., Phys. Plasmas 9, 1277 (2002)] on a diverted tokamak discharge. When the NIMROD model complexity is increased incrementally, resistive and first-order finite-Larmour radius effects are destabilizing and stabilizing, respectively. Lastly, the full result is compared to local analytic calculations which are found to overpredict both the resistive destabilization and drift stabilization in comparison to the NIMROD computations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

NUMERICAL SIMULATIONS OF OPTICALLY THICK ACCRETION ONTO A BLACK HOLE. I. SPHERICAL CASE

Modeling the radiation generated by accreting matter is an important step toward realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al., have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but it is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Impact of Interface Recombination on Time Resolved Photoluminescence (TRPL) Decays in CdTe Solar Cells (Numerical Simulation Analysis): Preprint

Using Sentaurus Device Software, we analyze how bulk and interface recombination affect time-resolved photoluminescence (TRPL) decays in CdTe solar cells. This modeling analysis could improve the interpretation of TRPL data and increase the possibility of rapid defect characterization in thin-film solar cells. By illuminating the samples with photons of two different wavelengths, we try to deducethe spatial origin of the dominant recombination loss. Shorter-wavelength photons will be more affected by the interface recombination and drift compared to the longer ones. Using the two-wavelength TRPL characterization method, it may be possible to determine whether a specific change in deposition process has affected the properties of interface or the bulk of the absorber.

14 SOLAR ENERGY↗
DOE OSTI2017

Properties of Turbulence in the Reconnection Exhaust: Numerical Simulations Compared with Observations

The properties of the turbulence that develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. In this study, we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail, investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the dissipation mechanisms in the turbulent outflows similar to the ones observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.

79 ASTRONOMY AND ASTROPHYSICS↗