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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Development and numerical simulation of spherical microresonators based on SiO{sub 2} – GeO{sub 2} germanosilicate glasses for generation of optical frequency combs

We propose and theoretically investigate the possibility of using spherical, whispering gallery mode microresonators based on germanosilicate glasses with different GeO{sub 2} contents for generating optical frequency combs in the dissipative soliton regime under pumping at a wavelength of 1.55 or 2 μm. The dispersion and nonlinearity of microspheres of different radii are calculated and analysed, and their optimal characteristics and the expected parameters of the output radiation are determined. It is shown that the spectral widths of optical frequency combs formed in 0.8SiO{sub 2} – 0.2GeO{sub 2} and GeO{sub 2} glass microspheres pumped at λ = 1.55 and 2 μm can be ∼200 and ∼300 nm, respectively. In these cases, in addition to the dissipative soliton with a duration of ∼100 fs, the generation of dispersive waves is also observed. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Numerical simulation of dispersion and nonlinear characteristics of microstructured silica fibres with a thin suspended core in a wide range of their parameters

Dispersion and nonlinear characteristics of microstructured silica fibres with a thin suspended core surrounded by three, four or six air holes have been studied theoretically in the wavelength range 1 – 2 μm. It has been shown that, owing to strong fundamental mode confinement near the core, the Kerr nonlinearity coefficient can exceed the nonlinearity coefficient of standard telecom fibre SMF28e by two orders of magnitude. The large waveguide contribution allows for effective group velocity dispersion management. Estimates are presented that demonstrate the feasibility of using suspended core fibre exhibiting Kerr nonlinearity for generating non-classical light: a state with squeezed quantum fluctuations in one of the quadrature components of a cw laser signal at a wavelength near 1.55 μm. (paper)

36 MATERIALS SCIENCE↗
DOE OSTI2020

Numerical simulation of the beam self-cleaning process in a multimode graded-index fibre during propagation of a pump wave and a Stokes component

A model of coupled modes is proposed, which makes it possible to describe the propagation of a pump wave and a Stokes component in a multimode graded-index fibre. A decrease in the depression in the centre of the profile of the pump wave intensity distribution as a result of the influence of Kerr effects and random linear coupling of modes is demonstrated. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Numerical Simulation of Isothermal Flow Across Slant Five-Tube Bundle with Spectral Element Method Code Nek5000

The helical coil steam generator (HCSG) is a specific type of shell-and-tube heat exchanger known for having a higher heat transfer coefficient than many other designs. For this reason, they are considered in small modular reactor and high temperature reactor designs. Investigation of flow behavior in HCSGs is important for better design. In this paper we present our study for modeling the coolant flow in the primary side of the HCSG. We used Nek5000, an open source, high-order spectral element computational fluid dynamics code developed in Argonne National Laboratory. Nek5000 accepts only hexahedral mesh, which makes the meshing process for the complicated HCSG geometry very challenging. A tetrahedral-to-hexahedral meshing strategy was applied to bypass the geometric complexities. In this study large eddy simulation (LES) was performed at the Reynolds number of 9000 based on the inlet velocity and the tube diameter. The employed subgrid-scale model for LES relies on explicit filtering. First- and second-order statistics were compared to available experimental data. We report overall velocity and turbulent kinetic energy showed good agreement with particle image velocimetry data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2010

DETAILED NUMERICAL SIMULATIONS ON THE FORMATION OF PILLARS AROUND H II REGIONS

We study the structural evolution of turbulent molecular clouds under the influence of ionizing radiation emitted from a nearby massive star by performing a high-resolution parameter study with the iVINE code. The temperature is taken to be 10 K or 100 K, the mean number density is either 100 cm{sup -3} or 300 cm{sup -3}. Furthermore, the turbulence is varied between Mach 1.5 and Mach 12.5, the main driving scale of the turbulence is varied between 1 pc and 8 pc. We vary the ionizing flux by an order of magnitude, corresponding to allowing between 0.5% and 5% of the mass in the domain to be ionized immediately. In our simulations, the ionizing radiation enhances the initial turbulent density distribution and thus leads to the formation of pillar-like structures observed adjacent to H II regions in a natural way. Gravitational collapse occurs regularly at the tips of the structures. We find a clear correlation between the initial state of the turbulent cold cloud and the final morphology and physical properties of the structures formed. The most favorable regime for the formation of pillars is Mach 4-10. Structures and therefore stars only form if the initial density contrast between the high-density unionized gas and the gas that is going to be ionized is lower than the temperature contrast between the hot and the cold gas. The density of the resulting pillars is determined by a pressure equilibrium between the hot and the cold gas. A thorough analysis of the simulations shows that the complex kinematical and geometrical structure of the formed elongated filaments reflects that of observed pillars to an impressive level of detail. In addition, we find that the observed line-of-sight velocities allow for a distinct determination of different formation mechanisms. Comparing the current simulations to previous results and recent observations, we conclude that, e.g., the pillars of creation in M16 formed by the mechanism proposed here and not by the radiation driven implosion of pre-existing clumps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

TURBULENT MIXING AND LAYER FORMATION IN DOUBLE-DIFFUSIVE CONVECTION: THREE-DIMENSIONAL NUMERICAL SIMULATIONS AND THEORY

Double-diffusive convection, often referred to as semi-convection in astrophysics, occurs in thermally and compositionally stratified systems which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. This process has been given relatively little attention so far, and its properties remain poorly constrained. In this paper, we present and analyze a set of three-dimensional simulations of this phenomenon in a Cartesian domain under the Boussinesq approximation. We find that in some cases the double-diffusive convection saturates into a state of homogeneous turbulence, but with turbulent fluxes several orders of magnitude smaller than those expected from direct overturning convection. In other cases, the system rapidly and spontaneously develops closely packed thermo-compositional layers, which later successively merge until a single layer is left. We compare the output of our simulations with an existing theory of layer formation in the oceanographic context and find very good agreement between the model and our results. The thermal and compositional mixing rates increase significantly during layer formation and increase even further with each merger. We find that the heat flux through the staircase is a simple function of the layer height. We conclude by proposing a new approach to studying transport by double-diffusive convection in astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

TURBULENT LINEWIDTHS IN PROTOPLANETARY DISKS: PREDICTIONS FROM NUMERICAL SIMULATIONS

Submillimeter observations of protoplanetary disks now approach the acuity needed to measure the turbulent broadening of molecular lines. These measurements constrain disk angular momentum transport, and furnish evidence of the turbulent environment within which planetesimal formation takes place. We use local magnetohydrodynamic (MHD) simulations of the magnetorotational instability (MRI) to predict the distribution of turbulent velocities in low-mass protoplanetary disks, as a function of radius and height above the mid-plane. We model both ideal MHD disks and disks in which Ohmic dissipation results in a dead zone of suppressed turbulence near the mid-plane. Under ideal conditions, the disk mid-plane is characterized by a velocity distribution that peaks near v {approx_equal} 0.1c{sub s} (where c{sub s} is the local sound speed), while supersonic velocities are reached at z > 3H (where H is the vertical pressure scale height). Residual velocities of v Almost-Equal-To 10{sup -2} c{sub s} persist near the mid-plane in dead zones, while the surface layers remain active. Anisotropic variation of the linewidth with disk inclination is modest. We compare our MHD results to hydrodynamic simulations in which large-scale forcing is used to initiate similar turbulent velocities. We show that the qualitative trend of increasing v with height, seen in the MHD case, persists for forced turbulence and is likely a generic property of disk turbulence. Percentage level determinations of v at different heights within the disk, or spatially resolved observations that probe the inner disk containing the dead zone region, are therefore needed to test whether the MRI is responsible for protoplanetary disk turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE ROLE OF STREAMERS IN THE DEFLECTION OF CORONAL MASS EJECTIONS: COMPARISON BETWEEN STEREO THREE-DIMENSIONAL RECONSTRUCTIONS AND NUMERICAL SIMULATIONS

On 2009 September 21, a filament eruption and the associated coronal mass ejection (CME) were observed by the Solar Terrestrial Relations Observatory (STEREO) spacecraft. The CME originated from the southern hemisphere and showed a deflection of about 15 Degree-Sign toward the heliospheric current sheet (HCS) during the propagation in the COR1 field of view. The CME source region was near the central meridian, but no on-disk CME signatures could be seen from the Earth. The aim of this paper is to provide a physical explanation for the strong deflection of the CME observed on 2009 September 21. The two-sided view of the STEREO spacecraft allows us to reconstruct the three-dimensional travel path of the CME and the evolution of the CME source region. The observations are combined with a magnetohydrodynamic simulation, starting from a magnetic field configuration closely resembling the extrapolated potential field for that date. By applying localized shearing motions, a CME is initiated in the simulation, showing a similar non-radial evolution, structure, and velocity as the observed event. The CME gets deflected toward the current sheet of the larger northern helmet streamer due to an imbalance in the magnetic pressure and tension forces and finally gets into the streamer. This study shows that during solar minima, even CMEs originating from high latitude can be easily deflected toward the HCS, eventually resulting in geoeffective events. How rapidly they undergo this latitudinal migration depends on the strength of both the large-scale coronal magnetic field and the magnetic flux of the erupting filament.

79 ASTRONOMY AND ASTROPHYSICS↗