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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Rings of the North Brazil Current: Their structure and behavior inferred from observations and a numerical simulation

Large anticyclonic rings are shed from the retroflecting North Brazil Current (NBC) near 8 deg N in the tropical western Atlantic. New subsurface velocity and temperature measurements within three such rings are presented here and are found to be consistent with previous in situ and remotely sensed NBC ring measurements. A high-resolution numerical model of the Atlantic Ocean forced by monthly wind stress and an imposed meridional overturning cell is found to shed NBC rings that approximate those observed. The model rings are more surface-intensified than those observed and somewhat smaller in diameter. Both observed and modeled NBC rings move northwestward along the coast of South America with a speed of 8-16 cm/s, considerably slower than predicted by analytical theories describing westward ring propagation. At least 2-3 rings per year separate from the NBC retroflection. Annually, 1-3 rings translate intact from their formation region near 50 deg W to the islands of the southeastern Caribbean, where they disintegrate after a lifetime of about 100 days. This study is the first to make use of subsurface temperature and velocity data to compute the volume of the anomalous ring core. NBC rings may be responsible for 3-4 Sv of direct mass transport across the equatorial-tropical gyre boundary or 20-25% of the total upper ocean cross-gyre transport required by the Atlantic meridional overturning cell. Translating NBC rings may contribute 20% of the total meridional heat transport by the ocean at this latitude.

Fratantoni, David M.↗
DOE OSTI2022

Airborne Lidar Measurements of XCO 2 in Synoptically Active Environment and Associated Comparisons With Numerical Simulations

Frontal boundaries have been shown to cause large changes in CO 2 mole-fractions, but clouds and the complex vertical structure of fronts make these gradients difficult to observe. It remains unclear how the column average CO 2 dry air mole-fraction (XCO 2 ) changes spatially across fronts, and how well airborne lidar observations, data assimilation systems, and numerical models without assimilation capture XCO 2 frontal contrasts (ΔXCO 2 , i.e., warm minus cold sector average of XCO 2 ). We demonstrated the potential of airborne Multifunctional Fiber Laser Lidar (MFLL) measurements in heterogeneous weather conditions (i.e., frontal environment) to investigate the ΔXCO 2 during four seasonal field campaigns of the Atmospheric Carbon and Transport-America (ACT-America) mission. Most frontal cases in summer (winter) reveal higher (lower) XCO 2 in the warm (cold) sector than in the cold (warm) sector. During the transitional seasons (spring and fall), no clear signal in ΔXCO 2 was observed. Intercomparison among the MFLL, assimilated fields from NASA's Global Modeling and Assimilation Office (GMAO), and simulations from the Weather Research and Forecasting-—Chemistry (WRF-Chem) showed that (a) all products had a similar sign of ΔXCO 2 though with different levels of agreement in ΔXCO 2 magnitudes among seasons; (b) ΔXCO 2 in summer decreases with altitude; and (c) significant challenges remain in observing and simulating XCO 2 frontal contrasts. A linear regression analyses between ΔXCO 2 for MFLL versus GMAO, and MFLL versus WRF-Chem for summer-2016 cases yielded a correlation coefficient of 0.95 and 0.88, respectively. The reported ΔXCO 2 variability among four seasons provide guidance to the spatial structures of XCO 2 transport errors in models and satellite measurements of XCO 2 in synoptically-active weather systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Discrete numerical simulations of torpedo anchor installation in granular soils

Torpedo anchors are a viable approach for mooring marine hydrokinetic (MHK) energy devices to the seafloor. These anchors can serve to maintain station and to provide the reaction force for an MHK device. The ability of the anchor to perform these duties is a strong function of its penetration depth during installation. This is a large-strain problem not amenable to typical continuum numerical approaches. In the current work, we propose that the discrete element method (DEM) is a more appropriate tool to investigate the shallow penetration of torpedo anchors in sands. The effects of anchor mass, impact velocity, and soil interparticle friction are considered in the DEM simulations. The relative maximum penetration depths for different penetration conditions are quantified and presented. Granular material response at the microscale during penetration are used to provide insight into system response. Energy dissipation in the assembly by both friction and collision at the particle scale are considered. Results show that anchor penetration increases approximately linearly with an increase in impact velocity or anchor weight. Penetration decreases with an increase in interparticle friction (i.e., soil strength). Observations of microscale behaviors and energy calculations are used to provide insight into overall system response.

58 GEOSCIENCES↗
DOE OSTI2016

A NUMERICAL SIMULATION OF COSMIC RAY MODULATION NEAR THE HELIOPAUSE. II. SOME PHYSICAL INSIGHTS

Cosmic ray (CR) transport near the heliopause (HP) is studied using a hybrid transport model, with the parameters constrained by observations from the Voyager 1 spacecraft. We simulate the CR radial flux along different directions in the heliosphere. There is no well-defined thin layer between the solar wind region and the interstellar region along the tail and polar directions of the heliosphere. By analyzing the radial flux curve along the direction of Voyager 2 , together with its trajectory information, the crossing time of the HP by Voyager 2 is predicted to be in 2017.14. We simulate the CR radial flux for different energy values along the direction of Voyager 1 . We find that there is only a modest modulation region of about 10 au wide beyond the HP, so that Voyager 1 observing the Local Interstellar Spectra is justified in numerical modeling. We analyze the heliospheric exit information of pseudo-particles in our stochastic numerical (time-backward) method, conjecturing that they represent the behavior of CR particles, and we find that pseudo-particles that have been traced from the nose region exit in the tail region. This implies that many CR particles diffuse directly from the heliospheric tail region to the nose region near the HP. In addition, when pseudo-particles were traced from the Local Interstellar Medium (LISM), it is found that their exit location (entrance for real particles) from the simulation domain is along the prescribed Interstellar Magnetic Field direction. This indicates that parallel diffusion dominates CR particle transport in the LISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Quasioptical modeling of wave beams with and without mode conversion. II. Numerical simulations of single-mode beams

This work continues a series of papers where we propose an algorithm for the quasi-optical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. Here, we present a corresponding code paraxial ray description and its test applications to single-mode beams in vacuum and also in inhomogeneous magnetized plasma. Here, the numerical results are compared with analytic formulas from Gaussian-beam optics and also with cold-plasma ray tracing, respectively. The quasi-optical simulations of mode-converting beams are reported in the third paper of this series.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗
NASA NTRS1992

Numerical simulations of drainage flows on Mars

Data collected by Viking Landers have shown that the meteorology of the near surface Martian environment is analogous to desertlike terrestrial conditions. Geological evidence such as dunes and frost streaks indicate that the surface wind is a potentially important factor in scouring of the martian landscape. In particular, the north polar basin shows erosional features that suggest katabatic wind convergence into broad valleys near the margin of the polar cap. The pattern of katabatic wind drainage off the north polar cap is similar to that observed on Earth over Antarctica or Greenland. The sensitivity is explored of Martian drainage flows to variations in terrain slope and diurnal heating using a numerical modeling approach. The model used is a 2-D sigma coordinate primitive equation system that has been used for simulations of Antarctic drainage flows. Prognostic equations include the flux forms of the horizontal scalar momentum equations, temperature, and continuity. Parameterization of both longwave (terrestrial) and shortwave (solar) radiation is included. Turbulent transfer of heat and momentum in the Martian atmosphere remains uncertain since relevant measurements are essentially nonexistent.

Parish, Thomas R.↗
NASA NTRS1991

Direct numerical simulation of transition and turbulence in a spatially evolving boundary layer

A high-order-accurate finite-difference approach to direct simulations of transition and turbulence in compressible flows is described. Attention is given to the high-free-stream disturbance case in which transition to turbulence occurs close to the leading edge. In effect, computation requirements are reduced. A method for numerically generating free-stream disturbances is presented.

Rai, Man M.↗
NASA NTRS2015

Numerical Simulations of Shock/Plume Interaction Using Structured Overset Grids

Computational simulations using structured overset grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting oblique shock/plume interaction effects to near-field sonic boom signatures. Standard second-order accurate as well as higher-resolution numerical discretizations are utilized and compared in the study. The numerical approach is compared with supersonic wind-tunnel data for three cases. The cases include an empty wind-tunnel at the operating conditions, an isolated shockgenerating diamond wedge within the tunnel, and a nozzle with diamond wedge configuration at five different nozzle pressure ratios. Solution sensitivity to numerical discretization is analyzed. Favorable comparisons between the computational results and experimental data of near-field pressure signatures are obtained. A simple prediction method for plume induced shock deflection is developed and results are compared with the CFD data.

Housman, Jeffrey A.↗
DOE OSTI2018

Numerical simulations of realistic grating compressors

Chirped pulse amplification (CPA), originally proposed by Strickland and Mourou, is a critical enabling technology for creating very short laser pulses with very high peak power. It consists of three optical components: a stretcher, an amplifier, and a compressor. An initial laser pulse of short duration and low energy content is first stretched out in time by the stretcher, resulting in a pulse with longer duration, low energy content, and very low peak power. The pulse is then amplified to increase its energy content without changing its duration in time. The amplified pulse is finally compressed in time in the compressor, resulting in a laser pulse with very short duration and very high peak power. The objective of this work is to develop numerical techniques that enable small-scale grating aberrations in the compressor to be analyzed through high performance computer simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Numerical simulation of transient moisture transfer into an electronic enclosure

Electronic systems are sometimes exposed to harsh environmental conditions of temperature and humidity. Moisture transfer into electronic enclosures and condensation can cause several problems such as corrosion and alteration in thermal stresses. It is therefore essential to study the local climate inside the enclosures to be able to protect the electronic systems. In this work, moisture transfer into a typical electronic enclosure is numerically studied using CFD. In order to reduce the CPU-time and make a way for subsequent factorial design analysis, a simplifying modification is applied in which the real 3D geometry is approximated by a 2D axial symmetry one. The results for 2D and 3D models were compared in order to calibrate the 2D representation. Furthermore, simulation results were compared with experimental data and good agreement was found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mid-height seismic isolation of equipment in nuclear power plants: Numerical simulations and design recommendations

Seismic isolation can mitigate earthquake demands on safety-class equipment and is being considered for application to some advanced nuclear reactors. For tall, slender vessels that could represent an advanced reactor, a steam generator or a heat exchanger, mid-height seismic isolation has been shown to be beneficial and practical. Here, this paper focuses on numerical modeling of a mid-height isolated tall, slender vessel with two primary goals 1) investigate and quantify the benefits of mid-height seismic isolation for a range of support structure stiffnesses, isolation systems, and seismic inputs, and 2) provide analysis and design recommendations for isolation of safety-class equipment. Results of response-history analysis were in good agreement with experimental measurements and demonstrated that mid-height isolation can substantially reduce seismic demands on tall, slender vessels for a range of support structure stiffnesses, isolation systems, and seismic inputs. Importantly, the reductions in horizontal spectral accelerations in the mid-height isolated vessel from the non-isolated condition were not affected by the stiffness of the support structure. Recommendations are made for analysis and design of isolated equipment and testing of isolators used for equipment protection in nuclear facilities, which are also applicable to other industries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1988

A numerical simulation of the full two-dimensional electrothermal de-icer pad

The ability to predict the time-temperature history of electrothermal de-icer pads is important in the subsequent design of improved and more efficient versions. These de-icer pads are installed near the surface of aircraft components, for the specific purpose of removing accreted ice. The proposed numerical model can incorporate the full 2-D geometry through a section of a region (i.e., section of an airfoil), that current 1-D numerical codes are unable to do. Thus, the effects of irregular layers, curvature, etc., can now be accounted for in the thermal transients. Each layer in the actual geometry is mapped via a body-fitted coordinate transformation into uniform, rectangular computational grids. The relevant heat transfer equations are transformed and discretized. To model the phase change that might occur in any accreted ice, in an enthalpy formulation the phase change equations are likewise transformed and discretized. The code developed was tested against numerous classical numerical solutions, as well as against experimental de-icing data on a UH1H rotor blade obtained from the NASA Lewis Research Center. The excellent comparisons obtained show that this code can be a useful tool in predicting the performance of current de-icer models, as well as in the designing of future models.

Masiulaniec, Konstanty C.↗
NASA NTRS1991

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.↗
NASA NTRS1994

Diagnostics of dust content in spiral galaxies: Numerical simulations of radiative transfer

In order to find the best observable diagnostics for the amount of internal extinction within spiral galaxies, we have constructed realistic models for disk galaxies with immersed dust layers. The radiative transfer including both scattering and absorption has been computed for a range of model galaxies in various orientations. Standard galaxy surface photometry techniques were then applied to the numerical data to illustrate how different observables such as total magnitude, color and luminosity distribution behave under given conditions of dust distribution. This work reveals a set of superior diagnostics for the dust in the disk. These include not only the integrated parameters, but also the apparent disk structural parameters, the amplitude of the asymmetry between the near and far sides of the galaxy as divided by the apparent major axis and their dependence on the orientation of the galaxy with respect to the observer. Combining the above diagnostics with our impressions of real galaxies, we arrive at the qualitative conclusion that galaxy disks are generally optically thin. Quantitative conclusions will appear in subsequent work.

Byun, Y. I.↗
NASA NTRS2003

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗
NASA NTRS1989

Numerical simulation of flow over a hypersonic aircraft using an explicit upwind PNS solver

A hypersonic flow field over a generic airplane configuration is simulated by solving the Parabolized Navier-Stokes (PNS) equations. The finite difference solution of the PNS equations is calculated using a noniterative space marching, explicit, upwind scheme recently developed by the authors. Special gridding techniques are used which allowed the sharp changes in surface geometry of the airplane configuration to be modelled without smoothing of corners. Comparisons of the PNS results to a solution of the Navier-Stokes equations demonstrates a good agreement of the numerical results in approximately 1/6 of the cpu time. This paper demonstrates that the explicit upwind algorithm for solving the PNS equations is an efficient method for simulating hypersonic flow fields about complete airplane configurations and should be considered as an alternative to solving the Navier-Stokes equations for flow fields where the PNS equations are valid.

Korte, John J.↗