Engineering Papers⌕ Search

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.

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

At least 1,027 records · Page 57

Shock drift acceleration

This is a review of the fundamental physics of the interactions of charged particles treated individually while they interact with fast mode magnetohydrodynamic shocks. Numerical simulation and analytical theory are used to develop predictions of the expected characteristics of this process strong upstream anisotropies directed along the magnetic field and downstream anisotropies tending to be peaked more perpendicular to the field; relatively more of the enhancement of higher-energy particles occurring upstream; sensitive dependence on shock normal to magnetic field angle of the efficiency of energization. Observations which display all of the above characteristics are reviewed. Also discussed is the relationship of shock drift acceleration to the models for stochastic transport of charged particles in the vicinity of shocks. Extensions of this work in both the observational and theoretical approaches are discussed.

Armstrong, Thomas P.↗

On the noncoplanarity of the magnetic field within a fast collisionless shock

Within the magnetic ramp of fast collisionless plasma shocks observed with spacecraft instruments and simulated numerically, the magnetic field undergoes an excursion out of the plane of coplanarity. This rotation is consistently in the direction such that the electrostatic potential jump across the shock, as measured in the de Hoffman-Teller frame of the reference (HTF), is about 2-6 times smaller than the electrostatic potential jump measured in the normal incidence frame. The preferred direction is consistent with a basic whistler mode transition between the upstream and downstream orientations. The potential jump in the HTF is considerably smaller than the change in bulk flow energy across the shock, confirming the recent suggestion that magnetic forces contribute importantly to the slowing of the plasma in that frame. A further consequence is that suprathermal particles leaking back into the upstream region across the shock do not gain much energy from the cross-shock electric field.

Thomsen, M. F.↗

High-latitude dayside electric fields and currents during strong northward interplanetary magnetic field - Observations and model simulation

On July 23, 1983 the IMF turned strongly northward, becoming about 22 nT for several hours. Using a combined data set of ionospheric convection measurements made by the Sondre Stromfjord incoherent scatter radar and convection inferred from Greenland magnetometer measurements, the onset of the reconfiguration of the high-latitude ionospheric currents is found to occur about 3 min after the northward IMF encounters the magnetopause. The large-scale reconfiguration of currents, however, appears to evolve over a period of about 22 min. These observations and the results of numerical simulations indicate that the dayside polar-cap electric field observed during strong northward IMF is produced by a direct electrical current coupling with the solar wind.

Clauer, C. Robert↗

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang↗

Comparison of a 3-D CFD-DSMC Solution Methodology With a Wind Tunnel Experiment

A solution method for problems that contain both continuum and rarefied flow regions is presented. The methodology is applied to flow about the 3-D Mars Sample Return Orbiter (MSRO) that has a highly compressed forebody flow, a shear layer where the flow separates from a forebody lip, and a low density wake. Because blunt body flow fields contain such disparate regions, employing a single numerical technique to solve the entire 3-D flow field is often impractical, or the technique does not apply. Direct simulation Monte Carlo (DSMC) could be employed to solve the entire flow field; however, the technique requires inordinate computational resources for continuum and near-continuum regions, and is best suited for the wake region. Computational fluid dynamics (CFD) will solve the high-density forebody flow, but continuum assumptions do not apply in the rarefied wake region. The CFD-DSMC approach presented herein may be a suitable way to obtain a higher fidelity solution.

Glass, Christopher E.↗

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗

A Computational Methodology for Simulating Thermal Loss Testing of the Advanced Stirling Convertor

The U.S. Department of Energy (DOE) and Lockheed Martin Space Systems Company (LMSSC) have been developing the Advanced Stirling Radioisotope Generator (ASRG) for use as a power system for space science missions. This generator would use two highefficiency Advanced Stirling Convertors (ASCs), developed by Sunpower Inc. and NASA Glenn Research Center (GRC). The ASCs convert thermal energy from a radioisotope heat source into electricity. As part of ground testing of these ASCs, different operating conditions are used to simulate expected mission conditions. These conditions require achieving a particular operating frequency, hot end and cold end temperatures, and specified electrical power output for a given net heat input. In an effort to improve net heat input predictions, numerous tasks have been performed which provided a more accurate value for net heat input into the ASCs, including the use of multidimensional numerical models. Validation test hardware has also been used to provide a direct comparison of numerical results and validate the multi-dimensional numerical models used to predict convertor net heat input and efficiency. These validation tests were designed to simulate the temperature profile of an operating Stirling convertor and resulted in a measured net heat input of 244.4 W. The methodology was applied to the multi-dimensional numerical model which resulted in a net heat input of 240.3 W. The computational methodology resulted in a value of net heat input that was 1.7 percent less than that measured during laboratory testing. The resulting computational methodology and results are discussed.

Reid, Terry V.↗

The analysis of apparent optical depth profiles for interstellar absorption lines

Attention is given to the apparent optical depth method, a procedure for analyzing interstellar absorption lines. Observed absorption-line profiles are converted into profiles of apparent optical depth, and apparent column density per unit velocity. By comparing the latter for a given interstellar species having two or more absorption lines which differ in the product, the presence or absence of unresolved saturated structure in the profiles can be directly inferred. The method is illustrated using absorption-line data from the IUE satellite for the highly ionized gas toward HD 64760. Additional illustrations and a study of the method's accuracy are provided through a series of numerical simulations of multicomponent interstellar absorption situations. The method is compared to the standard growth curve method for deriving interstellar column densities. The principal value of the apparent optical depth method is that the absorption-line data are directly converted into a form that provides for direct scientific interpretations of the physical conditions in the interstellar absorbing medium as a function of velocity.

Savage, Blair D.↗

Dense Feature Tracking of Atmospheric Winds with Deep Optical Flow

Atmospheric winds are a key physical phenomenon impacting natural hazards, energy transport, ocean currents, large-scale circulation, and ecosystem fluxes. Observing winds is a complex process and presents a large gap in NASA’s Earth Observation System. Atmospheric motion vectors (AMVs) aim to fill this gap by making numerical estimates of cloud movement between sequences of multi-spectral satellite images, tracking clouds and water vapor. Recent imaging hardware and software advancements have enabled the use of numerical optical flow techniques to produce accurate and dense vector fields outperforming traditional methods. This work presents WindFlow as the first machine learning based system for feature tracking atmospheric motion using optical flow. Due to the lack of large-scale satellite-based observations, we leverage high-resolution numerical simulations from NASA's GEOS-5 Nature Run to perform supervised learning and transfer to satellite images. We demonstrate that our approach using deep learning based optical flow scales to ultra-high-resolution images of size 2881x5760 with less than 1 m/s bias and 2.5 m/s average error. Four network and learning architectures are compared and it is found that recurrent all-pairs field transforms (RAFT) produces the lowest errors on all metrics for wind speed and direction. Results on held out numerical outputs shows RAFT's good performance in each of the spatial, temporal, and physical dimensions. A comparison between WindFlow and an operational AMV product against rawinsonde observations show that RAFT transfers across simulations and thermal infrared satellite observations. This work shows that machine learning based optical flow is an efficient approach to generating robust feature tracking for AMVs consistently over large regions.

Atmospheric winds↗

Wave-current interaction near the Gulf Stream during the surface wave dynamics experiment

This paper presents a case study on the wave-current interaction near the local curvature of a Gulf Stream meander. The wave data were obtained from in situ measurements by a pitch-roll discus buoy during the Surface Wave Dynamics Experiment (SWADE) conducted off Wallops Island, Virginia, from October 1990 to March 1991. Owing to the advection of the Gulf Stream by the semidiurnal tide, the discus buoy was alternately located outside and inside the Gulf Stream. The directional wave measurements from the buoy show the changes in wave direction, wave energy, and directional spreading when waves encountered the current in the Gulf Stream meanders. A wave refraction model, using the ray-tracing method with an estimated Gulf Stream velocity field and meandering condition, was used to simulate wave refraction patterns and to estimate wave parameters at relative locations corresponding to buoy measurements. The numerical simulation shows that a focusing zone of wave rays was formed near the boundary and behind the crest of a simulated Gulf Stream meander. The focusing of wave rays causes changes in wave direction, increases in wave energy, and decreases in wave directional spreading, which are in good agreement with the results from the buoy measurements.

Wang, David W.↗

Numerical simulation of supersonic cone flow at high angle of attack

A conical symmetry assumption is applied to the full Navier-Stokes equations resulting in an equation set containing time and two coordinate directions as independent variables. The set is integrated by use of a finite difference technique for the particular case of sharp cones at incidence. Solutions are obtained and compared with experiment for auxiliary conditions corresponding to both laminar and turbulent flows. Closure for the turbulent flow case is provided by use of a scalor eddy viscosity model based on the mixing length hypothesis. Modifications to the eddy viscosity model were found which led to excellent surface pressure and surface flow direction agreement for turbulent flow at low supersonic Mach numbers. It is apparent that further work on the turbulence model is necessary for agreement at higher supersonic Mach numbers.

Mcrae, D. S.↗

Progress in Validation of Wind-US for Ramjet/Scramjet Combustion

Validation of the Wind-US flow solver against two sets of experimental data involving high-speed combustion is attempted. First, the well-known Burrows- Kurkov supersonic hydrogen-air combustion test case is simulated, and the sensitively of ignition location and combustion performance to key parameters is explored. Second, a numerical model is developed for simulation of an X-43B candidate, full-scale, JP-7-fueled, internal flowpath operating in ramjet mode. Numerical results using an ethylene-air chemical kinetics model are directly compared against previously existing pressure-distribution data along the entire flowpath, obtained in direct-connect testing conducted at NASA Langley Research Center. Comparison to derived quantities such as burn efficiency and thermal throat location are also made. Reasonable to excellent agreement with experimental data is demonstrated for key parameters in both simulation efforts. Additional Wind-US feature needed to improve simulation efforts are described herein, including maintaining stagnation conditions at inflow boundaries for multi-species flow. An open issue regarding the sensitivity of isolator unstart to key model parameters is briefly discussed.

Engblom, William A.↗

Computational fluid dynamics at NASA Ames and the numerical aerodynamic simulation program

Computers are playing an increasingly important role in the field of aerodynamics such as that they now serve as a major complement to wind tunnels in aerospace research and development. Factors pacing advances in computational aerodynamics are identified, including the amount of computational power required to take the next major step in the discipline. The four main areas of computational aerodynamics research at NASA Ames Research Center which are directed toward extending the state of the art are identified and discussed. Example results obtained from approximate forms of the governing equations are presented and discussed, both in the context of levels of computer power required and the degree to which they either further the frontiers of research or apply to programs of practical importance. Finally, the Numerical Aerodynamic Simulation Program--with its 1988 target of achieving a sustained computational rate of 1 billion floating-point operations per second--is discussed in terms of its goals, status, and its projected effect on the future of computational aerodynamics.

Peterson, V. L.↗

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

Short-time Lyapunov exponent analysis

A new technique for analyzing complicated fluid flows in numerical simulations has been successfully tested. The analysis uses short time Lyapunov exponent contributions and the associated Lyapunov perturbation fields. A direct simulation of the Taylor-Couette flow just past the onset of chaos demonstrated that this new technique marks important times during the system evolution and identifies the important flow features at those times. This new technique will now be applied to a 'minimal' turbulent channel.

Vastano, J. A.↗

Flash spectroscopic studies of the kinetics of the halorhodopsin photocycle

The photoreactions of halorhodopsin are complicated by the fact that the parent pigment and its photoproducts interact with chloride. Thus, in any photoreaction scheme at least four species have to be accounted for: HR565 and HR578 Cl-, as well as HR640 and HR520 Cl-. A photocycle scheme proposed earlier places the two main photointermediates of halorhodopsin, HR520 Cl- and HR640, into a single photocycle, with a chloride-dependent equilibrium between them [Oesterhelt, D., Hegemann, P., & Tittor, J. (1985) EMBO J. 4, 2351-2356]. This scheme, with the additional feature of direct photoproduction of HR640 from HR565, was tested in this work by using numerical solutions of the appropriate differential equations to simulate flash-induced absorption changes at 500 nm (production of HR520 Cl-) and at 660 nm (production of HR640). The time scale of the simulation was ms following the flash. Comparison of the simulated curves with experimental traces yielded a unique set of three rate constants. The proposed photocycle scheme and these rate constants predict well the shapes and amplitudes of flash traces at various chloride concentrations. It appears from the photocycle scheme, and the numerical values of rate constants, that chloride is bound with high affinity to the parent halorhodopsin molecule, but with much lower affinity to its main photointermediate. This may be the consequence of the fact that in the parent halorhodopsin in the retinal configuration is all-trans, but in the two photointermediates it is 13-cis.

Non-NASA Center↗

Digital Lunar Exploration Sites (DLES)

It has been almost 50 years since humans last set foot on the Moon. With NASA’s Artemis program, the United States and its international and commercial partners are embarking on a new endeavor to explore the lunar surface. Before we return, we will have simulated every aspect of these future missions. Many of these models and simulations (M&S) will rely on well-known and commonly-used technologies, some of which trace their origins back to the Apollo program. However, M&S has advanced significantly, as have the underlying computational capabilities. As a result, we are able to model many more aspects of the Artemis vehicles and support systems with significantly improved detail and confidence. Accurately and efficiently modeling the lunar environment will be critical to simulating the Artemis elements and mission activities. This includes characterizing and modeling lunar topography, smaller craters, exposed surface rocks, lunar regolith, surface lighting, and ambient thermal environment. These are all necessary for understanding fundamental behaviors and performance of vehicles and support systems in the lunar environment and are often determining factors in the selection of exploration sites and defining mission profiles. The Astromaterials Research and Exploration Sciences (ARES) team in the Exploration Integration and Science Directorate at NASA’s Johnson Space Center (JSC) and the NASA Exploration Systems Simulations (NExSyS) team in the Simulation and Graphics Branch at JSC are developing and maintaining the Digital Lunar Exploration Sites (DLES) data, documentation, and software packages. This data is being fed directly into a diverse collection of graphics and simulation environments, where it is used to construct the closest known truth for numerous potential Lunar South Pole landing sites. DLES is available to projects across NASA and particularly the Artemis program to support coordinated digital representations of the lunar environment. This paper describes the fundamental need for DLES, the science data sets that are going into DLES, some of the processes used to integrate this data into DLES products, the basic products that constitute DLES, and some examples of DLES in use.

Lunar↗

Implementing a VLBI Time Delay Model for Earth-orbiting Satellites: Partial derivatives and Verification

This document describes the partial derivatives for commonly estimated parameters including antenna positions, satellite position, and satellite velocity from a near-field VLBI delay model for Earth-orbiting satellites. This model was presented in the Journal of Geodesy by Jaron & Nothnagel (2019). In this context, we present a streamlined version of the near-field VLBI delay. From this simplified model, we deduce a delay rate expression and calculate partial derivatives, maintaining only the terms with significant impact on the computed derivatives’ magnitude. To verify the accuracy of the simplified model and the partial derivatives computed from it, we have created a simple simulation in Matlab of an Earth-orbiting satellite at the altitude of a typical Global Navigation Satellite Systems (GNSS) satellite. From this simulation, we compare the simplified and original delays, and we verify the magnitude and direction of the partial derivatives against the numerically computed derivatives from the original delay model. The partial derivatives and simplified VLBI delay model detailed here are implemented in Fortran in the open-source library VTD.

Geodesy↗