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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

GPU accelerated fully space and time resolved numerical simulations of self-focusing laser beams in SBS-active media

A computer cluster equipped with Graphics Processing Units (GPUs) is used for simulating nonlinear optical wave packets undergoing Kerr self-focusing and stimulated Brillouin scattering in fused silica. We first recall the model equations in full (3+1) dimensions. These consist of two coupled nonlinear Schrödinger equations for counterpropagating optical beams closed with a source equation for light-induced acoustic waves seeded by thermal noise. Compared with simulations on a conventional cluster of Central Processing Units (CPUs), GPU-based computations allow us to use a significant (16 times) larger number of mesh points within similar computation times. Reciprocally, simulations employing the same number of mesh points are between 3 and 20 times faster on GPUs than on the same number of classical CPUs. Performance speedups close to 45 are reported for isolated functions evaluating, e.g., the optical nonlinearities. Since the field intensities may reach the ionization threshold of silica, the action of a defocusing electron plasma is also addressed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

High-fidelity numerical simulation of the dynamic beam equation

A high-fidelity finite difference approximation of the dynamic beam equation is derived. Different types of well-posed boundary conditions are analysed. The boundary closures are based on the summation-by-parts (SBP) framework and the boundary conditions are imposed using a penalty (SAT) technique, to guarantee linear stability. The resulting SBP–SAT approximation leads to fully explicit time integration. The accuracy and stability properties of the newly derived SBP–SAT approximations are demonstrated for both 1-D and 2-D problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗
DOE OSTI2015

Numerical simulation and analysis of the April 2013 Chicago floods

The weather event associated to record Chicago floods on April 2013 is investigated by using the Weather Research and Forecasting (WRF) model. Observations at Argonne National Laboratory and multi-sensor (weather radar and rain gauge) precipitation data from the National Weather Service were employed to evaluate the model’s performance. The WRF model captured the synoptic-scale atmospheric features well, but the simulated 24-h accumulated precipitation and short-period temporal evolution of precipitation over the heavy-rain region were less successful. To investigate the potential reasons for the model bias, four supplementary sensitivity experiments using various microphysics schemes and cumulus parameterizations were designed. Of the five tested parameterizations, the WRF Single-Moment 6-class (WSM6) graupel scheme and Kain-Fritsch (KF) cumulus parameterization outperformed the others, such as Grell-Dévényi (GD) cumulus parameterization, which underestimated the precipitation by 30–50% on a regional-average scale. Morrison microphysics and KF outperformed the others for the spatial patterns of 24-h accumulated precipitation. The spatial correlation between observation and Morrison-KF was 0.45, higher than those for other simulations. All of the simulations underestimated the precipitation over northeastern Illinois (especially at Argonne) during 0400–0800 UTC 18 April because of weak ascending motion or small moisture. In conclusion, all of the simulations except WSM6-GD also underestimated the precipitation during 1200–1600 UTC 18 April because of weak southerly flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

Numerical simulations of mounding and submerging flows of shear-thinning jets impinging in a container

Continuous jets of non-Newtonian fluids impinging on a fluid surface exhibit instabilities from jet buckling and coiling at low Reynolds numbers to delayed die swell, mounding, and air entrainment at higher Reynolds numbers. Filling containers with complex fluids is an important process for many industries, where the need for high throughput requires operating at high Reynolds numbers. In this regime, air entrainment can produce a visually unappealing product, causing a major quality control issue. Just prior to the onset of air entrainment, however, there exists an ideal filling regime which we term “planar filling,” as it is characterized by a relatively flat free surface that maintains its shape over time. In this paper, we create a steady-state, 2-D axisymmetric finite element model to study the transition from planar filling to the onset of air entrainment in a container filling process with generalized-Newtonian fluids. We use this model to explore the operating window for Newtonian and shear-thinning (or, more generally, deformation-rate-thinning) fluids, demonstrating that the flow behavior is characterized by a balance between inertial, viscous, and gravitational forces, as characterized by the Reynolds and Froude numbers. A scaling analysis suggests that the relevant parameters for calculating these dimensionless numbers are located where the jet impacts the liquid surface, and simulations show that the transition from planar filling to air entrainment often occurs when Re ~ O(10). Our study found that the bottom and side surfaces of the container drastically influence this transition to entrainment, stabilizing the flow.

42 ENGINEERING↗
DOE OSTI2018

Numerical simulations of depressurization-induced gas production from an interbedded turbidite gas hydrate-bearing sedimentary section in the offshore India: Site NGHP-02-16 (Area-B)

The recent National Gas Hydrate Program Expedition 02 (NGHP-02) identified the existence of gas hydrate-bearing sand reservoirs at a number of sites in the offshore of India including Site NGHP-02-16 in Area-B of the Krishna-Godavari Basin. The architecture of that gas hydrate accumulation is characterized by thin, gas hydrate-bearing, high quality sand layers interbedded with mud layers within a turbidite sedimentary interval. The lowest gas hydrate-bearing layer contacting a thinly-interbedded saline aquifer designates the base of the gas hydrate stability zone (BGHSZ). The proximity of the BGHSZ and the average temperature around 20 °C make the reservoir a favorable target for hydrate destabilization by means of the depressurization method. The results of the reservoir simulations indicate high gas production potential from this marine gas hydrate deposit with manageable concomitant water production using a well completion design that hydraulically isolates layers with water-saturated sands. Using a detailed geological input model, the predicted cumulative gas rates reach ~3.1 × 10 4 m 3 /day (~1.1 MMscf/day) after 90 days of continuous depressurization and demonstrate sustained production rates ~3.0 × 10 4 m 3 /day (~1.0 MMscf/day) after 5 years of production. The interbedded nature of this gas hydrate occurrence promotes the development of horizontal dissociation interfaces between gas hydrate-bearing sand and mud layers. As a result, non-uniform gas production along the horizontal interfaces becomes a primary determinant of reservoir performance. Simulation cases have been executed to determine the impact of the uncertainty in in situ reservoir permeability and the manner in which intrinsic permeability dynamically changes during dissociation in response to the imposed effective stress increase. Here, the cases based on low initial effective permeability and high sensitivity of compaction to stress result in the least favorable production predictions.

58 GEOSCIENCES↗
DOE OSTI2021

Numerical simulation study of impurity B transport during real-time B powder injection in EAST

To understand the beneficial effects of impurity boron (B) on the plasma, the transport of B species in both atomic and ionic states in the edge plasma and the ensuing evolution of the background plasma are simulated during real-time B powder injection by using two-dimensional fluid code SOLPS-ITER. The initial background plasma profiles for the simulation are reconstructed based on an EAST upper-single null discharge (shot No. 70601). With the injection rate of B atoms, 4.0 × 10 21 atoms/s, the simulation shows that effective charge number is less than 1.1 at the inner boundary R - R sep = -5.2 cm (background plasma density, 5.1 × 10 19 m -3 ) at the outer mid-plane indicating that the B impurity does not have much influence on the energy confinement. The present work finds that the impurity B exhibits very specific features in the edge plasma, mostly present in the narrow region just outside the separatrix. The particular distribution pattern can reduce the heat flux greatly, leading to the divertor detachment. The divertor detachment can be achieved easily under an injection rate of B powder in a reasonable range. This study suggests that the impurity B has a promising potential application in future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗