SEARCH · Engineering Papers

Results for “direct 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

DOE OSTI2016

GAMMA-RAY BURSTS FROM MAGNETIC RECONNECTION: VARIABILITY AND ROBUSTNESS OF LIGHT CURVES

The dissipation mechanism that powers gamma-ray bursts (GRBs) remains uncertain almost half a century after their discovery. The two main competing mechanisms are the extensively studied internal shocks and the less studied magnetic reconnection. Here we consider GRB emission from magnetic reconnection accounting for the relativistic bulk motions that it produces in the jet's bulk rest frame. Far from the source the magnetic field is almost exactly normal to the radial direction, suggesting locally quasi-spherical thin reconnection layers between regions of oppositely directed magnetic field. We show that if the relativistic motions in the jet's frame are confined to such a quasi-spherical uniform layer, then the resulting GRB light curves are independent of their direction distribution within this layer. This renders previous results for a delta-function velocity-direction distribution applicable to a much more general class of reconnection models, which are suggested by numerical simulations. Such models that vary in their velocity-direction distribution differ mainly in the size of the bright region that contributes most of the observed flux at a given emission radius or observed time. The more sharply peaked this distribution, the smaller this bright region, and the stronger the light curve variability that may be induced by deviations from a uniform emission over the thin reconnection layer, which may be expected in a realistic GRB outflow. This is reflected both in the observed image at a given observed time and in the observer-frame emissivity map at a given emission radius, which are calculated here for three simple velocity-direction distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2011

Numerically Exact Computer Simulations of Light Scattering by Densely Packed, Random Particulate Media

Direct computer simulations of electromagnetic scattering by discrete random media have become an active area of research. In this progress review, we summarize and analyze our main results obtained by means of numerically exact computer solutions of the macroscopic Maxwell equations. We consider finite scattering volumes with size parameters in the range, composed of varying numbers of randomly distributed particles with different refractive indices. The main objective of our analysis is to examine whether all backscattering effects predicted by the low-density theory of coherent backscattering (CB) also take place in the case of densely packed media. Based on our extensive numerical data we arrive at the following conclusions: (i) all backscattering effects predicted by the asymptotic theory of CB can also take place in the case of densely packed media; (ii) in the case of very large particle packing density, scattering characteristics of discrete random media can exhibit behavior not predicted by the low-density theories of CB and radiative transfer; (iii) increasing the absorptivity of the constituent particles can either enhance or suppress typical manifestations of CB depending on the particle packing density and the real part of the refractive index. Our numerical data strongly suggest that spectacular backscattering effects identified in laboratory experiments and observed for a class of high-albedo Solar System objects are caused by CB.

Dlugach, Janna M.↗

NASA NTRS1996

Validation Tests of TASS for Application to 3-D Vortex Simulations

Direct analytical solutions can be useful in validating the Core formulation of numerical systems. In this document an exact analytical solution to the nonlinear Navier-Stokes equation is compared to the numerical results from the three-dimensional Terminal Area Simulation System (TASS). This exact solution, of which the derivation is included, is for Beltrami type flow. Direct comparison of TASS to the analytical Beltrami solution is then used in evaluating the accuracy of TASS.

Switzer, George F.↗

NASA NTRSDate not supplied

An Expanded Role for Python in Expediting System Simulation Development

The Python programming language has traditionally been used as a “scripting” language and is not generally recognized as a language for building system simulations, where C++, Java, and other compiled languages are typically used. While very powerful and gaining favor for a large set of programming tasks, Python is generally regarded as not having the speed to directly do the extensive numeric computations required in large-scale simulation. This paper aims to present an approach, supported by case studies, to re-think this assertion. Directly using Python in this new role has the potential to expedite an agile, spiral cycle of development of design-test-modify for systems simulation.

programming language↗

NASA NTRS2001

Impact and Collisional Processes in the Solar System

In the past year, we have successfully developed the techniques necessary to conduct impact experiments on ice at very low temperatures. We employ the method of embedding gauges within a target to measure the shock wave and material properties. This means that our data are not model dependent; we directly measure the essential parameters needed for numerical simulations of impact cratering. Since then we have developed a new method for temperature control of icy targets that ensures temperature equilibrium throughout a porous target. Graduate student, Sarah Stewart-Mukhopadhyay, is leading the work on ices and porous materials as the main thrust of her thesis research. Our previous work has focused on icy materials with no porosity, and we propose to extend our research to include porous ice and porous ice-silicate mixtures. There is little shockwave data for porous ice, and none of the data was acquired under conditions applicable to the outer solar system. The solid ice Hugoniot is only defined for initial temperatures above -20 C. Our program uniquely measures the properties of ice at temperatures directly applicable to the solar system. Previous experiments were conducted at ambient temperatures soon after removing the target from a cold environment, usually just below freezing, or in a room just below freezing. Since ice has an extremely complicated phase diagram, it is important to conduct experiments at lower temperatures to determine the true outcome of impacts in the outer solar system. This research is complementary to other programs on icy materials. Our work focuses on the inherent material properties by measuring the shock wave directly; this complements the macroscopic observations and immediately provides the parameters necessary to extend this research to the gravity regime. Our numerical simulations of impacts in porous ice under very low gravity conditions, such as found on comets, show that the final crater size and shape is very dependent on the dynamic strength of the material.

Ahrens, Thomas J.↗

DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗

DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗

DOE OSTI2020

On the peridynamic effective force state and multiphase constitutive correspondence principle

This article concerns modeling unsaturated deformable porous media as an equivalent single-phase and single-force state peridynamic material through the effective force state. The balance equations of linear momentum and mass of unsaturated porous media are presented by defining relevant peridynamic states. The energy balance of unsaturated porous media is utilized to derive the effective force state for the solid skeleton that is an energy conjugate to the nonlocal deformation state of the solid, and the suction force state. Through an energy equivalence, a multiphase constitutive correspondence principle is built between classical unsaturated poromechanics and peridynamic unsaturated poromechanics. The multiphase correspondence principle provides a means to incorporate advanced constitutive models in classical unsaturated porous theory directly into unsaturated peridynamic poromechanics. Finally, numerical simulations of localized failure in unsaturated porous media under different matric suctions are presented to demonstrate the feasibility of modeling the mechanical behavior of such three-phase materials as an equivalent single-phase peridynamic material through the effective force state concept.

42 ENGINEERING↗

DOE OSTI2021

Wind plants can impact long-term local atmospheric conditions

Abstract Long-term weather and climate observatories can be affected by the changing environments in their vicinity, such as the growth of urban areas or changing vegetation. Wind plants can also impact local atmospheric conditions through their wakes, characterized by reduced wind speed and increased turbulence. We explore the extent to which the wind plants near an atmospheric measurement site in the central United States have affected their long-term measurements. Both direct observations and mesoscale numerical weather prediction simulations demonstrate how the wind plants induce a wind deficit aloft, especially in stable conditions, and a wind speed acceleration near the surface, which extend $$\sim 30$$ ∼ 30 km downwind of the wind plant. Turbulence kinetic energy is significantly enhanced within the wind plant wake in stable conditions, with near-surface observations seeing an increase of more than 30% a few kilometers downwind of the plants.

17 WIND ENERGY↗

DOE OSTI2014

Breaking of relativistically intense longitudinal space charge waves: A description using Dawson sheet model

Spatio-temporal evolution of relativistically intense longitudinal space charge waves in a cold homogeneous plasma is studied analytically as well as numerically, as an initial value problem, using Dawson sheet model. It is found that, except for very special initial conditions which generates the well known longitudinal Akhiezer-Polovin mode, for all other initial conditions, the waves break through a novel mechanism called phase mixing at an amplitude well below the Akhiezer-Polovin limit. An immediate consequence of this is, that Akhiezer-Polovin waves break when subjected to arbitrarily small longitudinal perturbations. We demonstrate this by performing extensive numerical simulations. This result may be of direct relevance to ultrashort, ultraintense laser/beam pulse-plasma interaction experiments where relativistically intense waves are routinely excited.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1990

Efficient Monte Carlo simulation of rarefied flow in a small nozzle

The flow of carbon dioxide through a small conical nozzle is computed using two different numerical approaches. The first of these is the direct simulation Monte Carlo method in which the physical phenomena are modeled at the molecular level. The second technique solves the Navier-Stokes equations of fluid flow. The purpose of the paper is to compare the solutions generated by these two methods for the conditions prevailing in a small, low density nozzle. The comparison between the two sets of solutions reveals only small differences. However, the flow properties in the thick viscous layer close to the nozzle at the exit plane show important deviations which have implications for the subsequent determination of the plume backflow.

Boyd, Iain D.↗

NASA NTRS1989

Computational fluid dynamics at the Lewis Research Center: An overview

Lewis is a multidisciplinary Center with strong research and development programs in aeronautical and space propulsion, power, space communications, space experiments and materials. Computational fluid dynamics (CFD) is playing an important and growing role in most of these areas. Described here is how CFD is integrated into these programs and highlights elements of the CFD activities. Examples are presented of codes developed to predict flow fields in advanced propulsion systems and several of the code validation experiments are described. The CFD effort at Lewis ranges from basic research on new and improved algorithms through code development to the application of these codes to specific engineering problems. Because of the substantial improvement in CFD's predictive capability, its use at Lewis is on a steep growth path, spreading rapidly into new areas which had not traditionally taken advantage of the techniques of numerical simulation. Multidisciplinary codes and the future direction of CFD at Lewis are discussed.

Stubbs, Robert M.↗

NASA NTRS1989

Direct simulation of rarefied hypersonic flows

As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.

Moss, James N.↗

NASA NTRS1991

Shock wave-boundary layer interactions in rarefied gas flows

A numerical study is presented, using the direct simulation Monte Carlo (DSMC) method, of shock wave-boundary layer interactions in low density supersonic flows. Test cases include two-dimensional, axially-symmetric and three-dimensional flows. The effective displacement angle of the boundary layer is calculated for representative flat plate, wedge, and cone flows. The maximum pressure, shear stress, and heat transfer in the shock formation region is determined in each case. The two-dimensional reflection of an oblique shock wave from a flat plate is studied, as is the three-dimensional interaction of such a wave with a sidewall boundary layer.

Bird, G. A.↗

NASA NTRS2011

Impact of Aerosols on Convective Clouds and Precipitation

Aerosols are a critical factor in the atmospheric hydrological cycle and radiation budget. As a major reason for clouds to form and a significant attenuator of solar radiation, aerosols affect climate in several ways. Current research suggests that aerosol effects on clouds could further extend to precipitation, both through the formation of cloud particles and by exerting persistent radiative forcing on the climate system that disturbs dynamics. However, the various mechanisms behind these effects, in particular the ones connected to precipitation, are not yet well understood. The atmospheric and climate communities have long been working to gain a better grasp of these critical effects and hence to reduce the significant uncertainties in climate prediction resulting from such a lack of adequate knowledge. The central theme of this paper is to review past efforts and summarize our current understanding of the effect of aerosols on precipitation processes from theoretical analysis of microphysics, observational evidence, and a range of numerical model simulations. In addition, the discrepancy between results simulated by models, as well as that between simulations and observations will be presented. Specifically, this paper will address the following topics: (1) fundamental theories of aerosol effects on microphysics and precipitation processes, (2) observational evidence of the effect of aerosols on precipitation processes, (3) signatures of the aerosol impact on precipitation from large-scale analyses, (4) results from cloud-resolving model simulations, and (5) results from large-scale numerical model simulations. Finally, several future research directions on aerosol - precipitation interactions are suggested.

Tao, Wei-Kuo↗

NASA NTRS2012

Impact of Aerosols on Convective Clouds and Precipitation

Aerosols are a critical factor in the atmospheric hydrological cycle and radiation budget. As a major agent for clouds to form and a significant attenuator of solar radiation, aerosols affect climate in several ways. Current research suggests that aerosol effects on clouds could further extend to precipitation, both through the formation of cloud particles and by exerting persistent radiative forcing on the climate system that disturbs dynamics. However, the various mechanisms behind these effects, in particular the ones connected to precipitation, are not yet well understood. The atmospheric and climate communities have long been working to gain a better grasp of these critical effects and hence to reduce the significant uncertainties in climate prediction resulting from such a lack of adequate knowledge. Here we review past efforts and summarize our current understanding of the effect of aerosols on convective precipitation processes from theoretical analysis of microphysics, observational evidence, and a range of numerical model simulations. In addition, the discrepancy between results simulated by models, as well as that between simulations and observations, are presented. Specifically, this paper addresses the following topics: (1) fundamental theories of aerosol effects on microphysics and precipitation processes, (2) observational evidence of the effect of aerosols on precipitation processes, (3) signatures of the aerosol impact on precipitation from largescale analyses, (4) results from cloud-resolving model simulations, and (5) results from large-scale numerical model simulations. Finally, several future research directions for gaining a better understanding of aerosol--cloud-precipitation interactions are suggested.

Tao, Wei-Kuo↗

DOE OSTI2022

Data assimilation for burnup distribution of PWR with three-dimensional variational algorithm and artificial neutral network

In this paper, a data-assimilation method has been proposed and applied for the burnup distribution of PWR. The burnup distribution is significant to the safety and economy of the reactor, as it is essential for the fuel-reloading design and optimization. Due to the burnup distribution cannot be measured directly during the reactor operation, the numerical simulation is widely applied to determine the burnup distribution. However, there is a deviation between the numerical simulation and the actual core due to some unavoidable factors, such as component manufacturing deviation, uneven flow distribution and so on. These differences would induce the errors to the simulation values of power distributions and hence to the burnup distribution. To address this problem, a data-assimilation method for the burnup distribution has been proposed with the application of power-distribution measurements. In our research, the three-dimensional variational (3DVAR) algorithm was applied for burnup-distribution calibration and the artificial neutral network (ANN) was applied to establish the relation between power distribution and corresponding burnup distribution. As engineering verification, the proposed data-assimilation method has been applied to the CNP1000 PWR operated in China. The numerical results indicated that the burnup-distribution errors can be reduced notably, as the maximum value of relative errors for power distribution can be reduced from 5.35% to 3.96%. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1992

A domain-specific design architecture for composite material design and aircraft part redesign

Advanced composites have been targeted as a 'leapfrog' technology that would provide a unique global competitive position for U.S. industry. Composites are unique in the requirements for an integrated approach to designing, manufacturing, and marketing of products developed utilizing the new materials of construction. Numerous studies extending across the entire economic spectrum of the United States from aerospace to military to durable goods have identified composites as a 'key' technology. In general there have been two approaches to composite construction: build models of a given composite materials, then determine characteristics of the material via numerical simulation and empirical testing; and experience-directed construction of fabrication plans for building composites with given properties. The first route sets a goal to capture basic understanding of a device (the composite) by use of a rigorous mathematical model; the second attempts to capture the expertise about the process of fabricating a composite (to date) at a surface level typically expressed in a rule based system. From an AI perspective, these two research lines are attacking distinctly different problems, and both tracks have current limitations. The mathematical modeling approach has yielded a wealth of data but a large number of simplifying assumptions are needed to make numerical simulation tractable. Likewise, although surface level expertise about how to build a particular composite may yield important results, recent trends in the KBS area are towards augmenting surface level problem solving with deeper level knowledge. Many of the relative advantages of composites, e.g., the strength:weight ratio, is most prominent when the entire component is designed as a unitary piece. The bottleneck in undertaking such unitary design lies in the difficulty of the re-design task. Designing the fabrication protocols for a complex-shaped, thick section composite are currently very difficult. It is in fact this difficulty that our research will address.

Punch, W. F., III↗