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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Effects of Convection during the Photodeposition of Polydiacetylene Thin Films

In this work, we describe a preliminary investigation of buoyancy-driven heat transfer during the growth of thin films from solution following exposure to ultraviolet (UV) light. Irradiation of the growth cell occurs at various directions relative to gravitational acceleration. Through numerical computations, the steady-state flow and temperature profiles are simulated during the course of light exposure. Light-induced polymerization accompanies a heat transfer process through a fairly complicated recirculating flow pattern. A scaling analysis shows that buoyancy-driven velocities only reduce by a factor of 10 for gravity levels as low as 10(exp -2) g(sub 0). Paley et al. observe what appears to be gravitationally sensitive particle development and inclusion in thin films using a photodeposition process. From this study, it is clear that production of homogeneous thin films would have to occur in the environment of a complicated flow pattern of recirculation with a nonuniform temperature distribution. Indeed, even when irradiation occurs from the top of the cell, the most stable stratified cell orientation, defects remain in our films due to the persistence of buoyancy-driven convection. To achieve homogeneity, minimal scattering centers, and possible molecular order, photodeposition of polymer films by UV light exposure must proceed in a reduced-convection environment. Fluid mechanics simulations are useful for establishing gravitational sensitivity to this recently discovered process (patent # 5,451,433) for preparing thin films having quite promising nonlinear optical characteristics.

Frazier, D. O.↗
NASA NTRS1997

Magnetic Damping of g-Jitter Driven Flows: 3-D Calculations

A 3-D numerical model is developed to represent the oscillating natural convection induced in a cylindrical cavity filled with Ga-doped germanium with and without the presence of an external magnetic field. The model is developed based on the penalty-finite element solution of the equations describing the transport of momentum, heat and solutal element as well as the electromagnetic field distribution in the melt pool. Automatic time step control is applied to help speed up the calculations. Numerical simulations are conducted to study the convection and magnetic damping effects as a function of frequency, directions and amplitudes of g-jitter and also the direction and magnitudes of the applied magnetic fields. The results show that the g-jitter driven flow is time dependent and exhibits a complex recirculating convection pattern in three dimensions and that an applied magnetic field can be employed to suppress this deleterious convective flow and both magnitude and orientation of the applied field are important in magnetic damping of the g-jitter induced convective flows.

Shang, D. Y.↗
NASA NTRS1997

Global MHD Magnetospheric Simulation of January 10, 1997 Encounter with Magnetic Cloud

We present a global MHD magnetospheric simulation of the encounter of the magnetic cloud with the terrestrial magnetosphere on January 10-11, 1997. The MHD simulation is driven by IMF and solar wind plasma measurements provided by Wind situated about 100 $R_E$ upstream from Earth. Field-aligned currents produced in the model are mapped down to the ionosphere and are directly compared to auroral images from the Polar UVI. Particular attention will be paid to the interval following the the initial shock wave arrival (around 0100 UT on the 1Oth) but preceding the passage of the magnetic cloud proper (commencing at about 0430 UT) for which there was continuous viewing of the entire auroral oval from apogee by Polar. This turbulent period is characterized by numerous dynamic pressure rises and dips and several northward and southward turnings of the IMF, all of which generate dynamic activity in the simulation that is reflected in the mapped field-aligned current patterns. As discussed by Brittnacher et al. in this session, the auroral morphology imaged by UVI during this period includes: shock wave-induced brightening of the oval followed by a pseudo-onset near midnight, several sun-aligned and curled arcs within the polar cap, an expanding polar cap cleared of arcs, and finally substorm onset at 0337 UT. These features will be directly compared to region 1 and 2 current systems as well as cusp currents in the simulation. Spacecraft magnetopause crossings are also predicted by the global simulation and will be compared to observed crossings, including numerous dayside crossings by Geotail on the 1Oth, and geosynchronous crossings early on January 11.

Elsen, R.↗
NASA NTRS1990

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗
NASA NTRS2003

Prospect of Using Numerical Dynamo Model for Prediction of Geomagnetic Secular Variation

Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Kuang, Weijia↗
NASA NTRS2019

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗
NASA NTRS1992

Experimental and numerical investigations of low-density nozzle and plume flows of nitrogen

New experimental data are used to show that the direct simulation Monte Carlo (DSMC) method provides an accurate description of a low-density flow in the nozzle and the near-field expansion of a small rocket for two slightly different experimental configurations. These results verify the DSMC method in an expansion flow for the first time. From a number of different gas/surface interaction models, it is found that fully diffuse reflection gives the best agreement with experiment.

Boyd, Iain D.↗
NASA NTRS2001

Formation of Electrostatic Potential Drops in the Auroral Zone

In order to examine the self-consistent formation of large-scale quasi-static parallel electric fields in the auroral zone on a micro/meso scale, a particle in cell simulation has been developed. The code resolves electron Debye length scales so that electron micro-processes are included and a variable grid scheme is used such that the overall length scale of the simulation is of the order of an Earth radii along the magnetic field. The simulation is electrostatic and includes the magnetic mirror force, as well as two types of plasmas, a cold dense ionospheric plasma and a warm tenuous magnetospheric plasma. In order to study the formation of parallel electric fields in the auroral zone, different magnetospheric ion and electron inflow boundary conditions are used to drive the system. It has been found that for conditions in the primary (upward) current region an upward directed quasi-static electric field can form across the system due to magnetic mirroring of the magnetospheric ions and electrons at different altitudes. For conditions in the return (downward) current region it is shown that a quasi-static parallel electric field in the opposite sense of that in the primary current region is formed, i.e., the parallel electric field is directed earthward. The conditions for how these different electric fields can be formed are discussed using satellite observations and numerical simulations.

Schriver, D.↗
NASA NTRS1979

A two dimensional model of viscous slot flow

The effects of viscosity on flows through slots are investigated. Using Briley's formulation of the Alternating Direction Implicit method together with the Thompson transformation, numerical solution of viscous incompressible flow through a slot is simulated, for Reynolds number 1000 and openness ratio one tenth. The problem is formulated in terms of stream function and vorticity. A wall iteration of Israeli type is employed to iteratively induce the no slip condition on velocity.

Cooke, C. H.↗
NASA NTRS1984

Nonlinear evolution of the resistive tearing mode

The nonlinear behavior of the tearing instability is investigated with the aid of numerical solutions of the resistive incompressible magnetohydrodynamic equations. This study is directed toward solar applications, such as flares. Simulations were carried out for values of the Lundquist number S from 100 to 1,000,000 and wavelength parameter alpha from 0.042 to 0.5. The influence of linear conditions, such as the constancy of psi, on nonlinear growth characteristics is investigated. For high S and low alpha, secondary-flow vortices, opposite in direction to the linear vortices, were found to generate a new magnetic island centered at the initial X point. The nonlinear spatial distributions of the physical variables differ greatly from the linear behavior, and growth slows considerably once nonlinear effects become important.

Steinolfson, R. S.↗
NASA NTRS1988

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

Chapman, Dean R.↗
NASA NTRS1989

Theoretical and numerical study of a three-dimensional turbulent boundary layer

A boundary layer created on an infinite flat plate by a time-dependent freestream velocity vector whose magnitude is independent of time but whose direction changes at a constant angular velocity is theoretically studied using Reynolds-number scaling laws and numerical simulations performed over a range of Reynolds numbers. Results obtained with a higher-order version of existing theories of the Ekman layer are shown to agree well with the numerical results at three Reynolds numbers. The present results can be extrapolated to the case of high Reynolds numbers. The Reynolds-averaged equations reduce to a one-dimensional steady problem, making possible the easy and accurate testing of turbulence models.

Spalart, Philippe R.↗
NASA NTRS1995

The tilts of corotating interaction regions at mid heliographic latitudes

Ulysses observations of corotating interaction regions (CIRs) at mid heliographic latitudes have shown that the flow downstream of the forward shock (or wave) on the leading edge of a CIR generally turns northward and westward, while the flow downstream of the reverse wave on the trailing edge generally turns southward and eastward. These systematic flow deflections are a natural consequence of large scale pressure gradients associated with the CIRs, and indicate that the forward waves tend to propagate toward and across the equator with increasing heliocentric distance, while the reverse waves tend to propagate toward the pole. Recent determinations of CIR shock normals using the Ulysses magnetic field data appear to confirm these plasma results (Burton, private communication). Numerical simulations indicate that these effects (which imply that CIRs are systematically tilted in the north-south direction at mid latitudes) are a natural consequence of the tilt of the solar magnetic dipole axis relative to the solar rotation axis. The present work utilizes a variety of techniques to analyze the flow deflections observed within CIRs from which we can infer the overall orientations of the CIRs and the speeds and directions of propagation of the waves. Where possible, the observations are quantitatively compared with the results of 3-dimensional MHD simulations.

Riley, Pete↗
NASA NTRS1995

Effects of Chemistry on Blunt-Body Wake Structure

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using direct simulation Monte Carlo (DSMC) and Navier-Stokes calculations. Particular emphasis is given to the effects of chemistry on the near-wake structure and on the surface quantities and the comparison of the DSMC results with the Navier-Stokes calculations. The flow conditions simulated are those experienced by a space vehicle at an altitude of 85 km and a velocity of 7 km/s during Earth entry. A steady vortex forms in the near wake for these freestream conditions for both chemically reactive and nonreactive air gas models. The size (axial length) of the vortex for the reactive air calculations is 25% larger than that of the nonreactive air calculations. The forebody surface quantities are less sensitive to the chemistry than the base surface quantities. The presence of the afterbody has no effect on the forebody flow structure or the surface quantities. The comparisons of DSMC and Navier-Stokes calculations show good agreement for the wake structure and the forebody surface quantities.

Dogra, Virendra K.↗
NASA NTRS2002

A Comparison of Trajectory Optimization Methods for the Impulsive Minimum Fuel Rendezvous Problem

In this paper we present a comparison of optimization approaches to the minimum fuel rendezvous problem. Both indirect and direct methods are compared for a variety of test cases. The indirect approach is based on primer vector theory. The direct approaches are implemented numerically and include Sequential Quadratic Programming (SQP), Quasi-Newton, Simplex, Genetic Algorithms, and Simulated Annealing. Each method is applied to a variety of test cases including, circular to circular coplanar orbits, LEO to GEO, and orbit phasing in highly elliptic orbits. We also compare different constrained optimization routines on complex orbit rendezvous problems with complicated, highly nonlinear constraints.

Hughes, Steven P.↗
NASA NTRS1985

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗
NASA NTRS1984

Simulation requirements for the Large Deployable Reflector (LDR)

Simulation tools for the large deployable reflector (LDR) are discussed. These tools are often the transfer function variety equations. However, transfer functions are inadequate to represent time-varying systems for multiple control systems with overlapping bandwidths characterized by multi-input, multi-output features. Frequency domain approaches are the useful design tools, but a full-up simulation is needed. Because of the need for a dedicated computer for high frequency multi degree of freedom components encountered, non-real time smulation is preferred. Large numerical analysis software programs are useful only to receive inputs and provide output to the next block, and should be kept out of the direct loop of simulation. The following blocks make up the simulation. The thermal model block is a classical heat transfer program. It is a non-steady state program. The quasistatic block deals with problems associated with rigid body control of reflector segments. The steady state block assembles data into equations of motion and dynamics. A differential raytrace is obtained to establish a change in wave aberrations. The observation scene is described. The focal plane module converts the photon intensity impinging on it into electron streams or into permanent film records.

Soosaar, K.↗
NASA NTRS2001

Cellular Spacing Selection During the Directional Solidification of Binary Alloys. A Numerical Approach

The evolution of cellular solid/liquid interfaces from an initially unstable planar front was studied by means of a two-dimensional computer simulation. The developed numerical model makes use of an interface tracking procedure and has the capability to describe the dynamics of the interface morphology based on local changes of the thermodynamic conditions. The fundamental physics of this formulation was validated against experimental microgravity results and the predictions of the analytical linear stability theory. The performed simulations revealed that in certain conditions, based on a competitive growth mechanism, an interface could become unstable to random perturbations of infinitesimal amplitude even at wavelengths smaller than the neutral wavelength, lambda(sub c), predicted by the linear stability theory. Furthermore, two main stages of spacing selection have been identified. In the first stage, at low perturbations amplitude, the selection mechanism is driven by the maximum growth rate of instabilities while in the second stage the selection is influenced by nonlinear phenomena caused by the interactions between the neighboring cells. Comparison of these predictions with other existing theories of pattern formation and experimental results will be discussed.

Catalina, Adrian V.↗