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,207 records · Page 67

Updates to Simulation of a Single-Element Lean-Direct Injection Combustor Using Arbitary Polyhedral Meshes

This paper summarizes the procedures of (1) generating control volumes anchored at the nodes of a mesh; and (2) generating staggered control volumes via mesh reconstructions, in terms of either mesh realignment or mesh refinement, as well as presents sample results from their applications to the numerical solution of a single-element LDI combustor using a releasable edition of the National Combustion Code (NCC).

Power↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

An Object Oriented Extensible Architecture for Affordable Aerospace Propulsion Systems

Driven by a need to explore and develop propulsion systems that exceeded current computing capabilities, NASA Glenn embarked on a novel strategy leading to the development of an architecture that enables propulsion simulations never thought possible before. Full engine 3 Dimensional Computational Fluid Dynamic propulsion system simulations were deemed impossible due to the impracticality of the hardware and software computing systems required. However, with a software paradigm shift and an embracing of parallel and distributed processing, an architecture was designed to meet the needs of future propulsion system modeling. The author suggests that the architecture designed at the NASA Glenn Research Center for propulsion system modeling has potential for impacting the direction of development of affordable weapons systems currently under consideration by the Applied Vehicle Technology Panel (AVT). This paper discusses the salient features of the NPSS Architecture including its interface layer, object layer, implementation for accessing legacy codes, numerical zooming infrastructure and its computing layer. The computing layer focuses on the use and deployment of these propulsion simulations on parallel and distributed computing platforms which has been the focus of NASA Ames. Additional features of the object oriented architecture that support MultiDisciplinary (MD) Coupling, computer aided design (CAD) access and MD coupling objects will be discussed. Included will be a discussion of the successes, challenges and benefits of implementing this architecture.

Follen, Gregory J.↗

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗

Comparison Between Navier-Stokes and Thin-Layer Computations for Separated Supersonic Flow

In the numerical simulation of high Reynolds-number flow, one can frequently supply only enough grid points to resolve the viscous terms in a thin layer. As a consequence, a body-or stream-aligned coordinate system is frequently used and viscous terms in this direction are discarded. It is argued that these terms cannot be resolved and computational efficiency is gained by their neglect. Dropping the streamwise viscous terms in this manner has been termed the thin-layer approximation. The thin-layer concept is an old one, and similar viscous terms are dropped, for example, in parabolized Navier-Stokes schemes. However, such schemes also make additional assumptions so that the equations can be marched in space, and such a restriction is not usually imposed on a thin-layer model. The thin-layer approximation can be justified in much the same way as the boundary-layer approximation; it requires, therefore, a body-or stream-aligned coordinate and a high Reynolds number. Unlike the boundary-layer approximation, the same equations are used throughout, so there is no matching problem. Furthermore, the normal momentum equation is not simplified and the convection terms are not one-sided differenced for marching. Consequently, the thin-layer equations are numerically well behaved at separation and require no special treatment there. Nevertheless, the thin-layer approximation receives criticism. It has been suggested that the approximation is invalid at separation and, more recently, that it is inadequate for unsteady transonic flow. Although previous comparisons between the thin-layer and Navier-Stokes equations have been made, these comparisons have not been adequately documented.

Degani, David↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

Propagation of wide bandwidth signals through strongly turbulent ionized media

Analytic and numerical techniques are presented which directly address the problem of propagation of wide bandwidth signals through random ionized media. This work is applicable to the problems of satellite communication and space based radar observation through a disturbed ionospheric propagation channel that would result from a high altitude chemical release or nuclear detonation. An analytic solution is presented for the two-position, two-frequency mutual coherence function for spherical wave propagation in the strong scattering limit. This solution is used to derive simple expressions for the power impulse response function and to determine its relationship to the important parameters of decorrelation distance and coherence bandwidth which describe the disturbed propagation channel. Results for mean time delay and time delay jitter are presented and compared to direct simulation results and to other theoretical calculations. A numerical/analytical solution to the parabolic wave equation is presented in the form of a multiple phase-screen (MPS) propagation simulation. In this very general solution technique, the random medium is divided into a finite number of layers. The field fluctuation through each layer is obtained by replacing the layer by a centrally located thin phase-screen, whose statistical properties are determined from the statistics of the electron-density irregularities. The waveform then propagates from phase-screen to phase-screen via an exact solution to the Fresnel-Kirchnoff equation. For wide bandwidth waveforms, numerical solutions are obtained at a number of discrete frequencies centered about the carrier and then time-domain solutions are obtained by Fourier transform techniques. Detailed results are presented for a number of interesting cases including propagation of a 20 MHz bandwidth signal through a finite barium cloud at a carrier frequency of 100 MHz. One of the major uses of the MPS propagation simulation has been to provide realizations of the received signal after propagation through a disturbed channel. The MPS simulation obtains a general solution of the parabolic wave equation under both weak and strong scattering conditions. A second technique to directly obtain realizations of wide bandwidth waveforms is presented here. This technique is limited to the case of strong scattering but requires only a fraction of the computer resources needed for MPS signal generation. Detailed comparisons of the two signal generation techniques are presented.

Random Media↗

Extending mARC II Arc-Jet Test Duration via Design and Implementation of Model System Cooling Sleeve

The mARC II is a 30 kW arc-jet facility at NASA Ames Research Center used to generate high-enthalpy flows for low-cost thermal protection system (TPS) technology development. Sustained operation of downstream instrumentation and material samples is constrained by thermal loading transmitted through the arc-jet test environment, limiting achievable run times and experimental throughput. This work presents the design, integration, and validation of a cooling sleeve implemented on the sweep arm drive motor feedthrough to mitigate thermal accumulation during testing. The addition of the cooling sleeve is a simple, robust upgrade that translates directly into enhanced facility capability by supporting longer run durations, reduced turnaround time, and higher throughput.

numerical simulations↗

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗

A numerical simulation of barotropic instability. II Wave-wave interaction

A fully nonlinear numerical model of the point jet barotropic instability is used to test and confirm the hypothesis that the magnitude of the wave vorticity does not exceed the magnitude of the initial shear. This result arises directly from the local conservation of vorticity following a parcel and the fact that unstable waves are principally confined to the region where the zonal mean vorticity can be smoothed by the wave so as to eliminate the instability. Comparisons are made between fully nonlinear and quasi-linear models of the point jet instability and their tracer transport properties. Differences become particularly evident after wave saturation. The most important effect neglected by the wave-mean flow model appears to be the advection of wave vorticity by the most unstable mode. However, as equilibration of the instability proceeds, the globally averaged properties of both models are found to be similar.

Nielsen, J. E.↗

Application of computational physics within Northrop

An overview of Northrop programs in computational physics is presented. These programs depend on access to today's supercomputers, such as the Numerical Aerodynamical Simulator (NAS), and future growth on the continuing evolution of computational engines. Descriptions here are concentrated on the following areas: computational fluid dynamics (CFD), computational electromagnetics (CEM), computer architectures, and expert systems. Current efforts and future directions in these areas are presented. The impact of advances in the CFD area is described, and parallels are drawn to analagous developments in CEM. The relationship between advances in these areas and the development of advances (parallel) architectures and expert systems is also presented.

George, M. W.↗

Rainfall Morphology in Semi-Tropical Convergence Zones

Central Florida is the ideal test laboratory for studying convergence zone-induced convection. The region regularly experiences sea breeze fronts and rainfall-induced outflow boundaries. The focus of this study is the common yet poorly-studied convergence zone established by the interaction of the sea breeze front and an outflow boundary. Previous studies have investigated mechanisms primarily affecting storm initiation by such convergence zones. Few have focused on rainfall morphology yet these storms contribute a significant amount precipitation to the annual rainfall budget. Low-level convergence and mid-tropospheric moisture have both been shown to correlate with rainfall amounts in Florida. Using 2D and 3D numerical simulations, the roles of low-level convergence and mid-tropospheric moisture in rainfall evolution are examined. The results indicate that time-averaged, vertical moisture flux (VMF) at the sea breeze front/outflow convergence zone is directly and linearly proportional to initial condensation rates. This proportionality establishes a similar relationship between VMF and initial rainfall. Vertical moisture flux, which encompasses depth and magnitude of convergence, is better correlated to initial rainfall production than surface moisture convergence. This extends early observational studies which linked rainfall in Florida to surface moisture convergence. The amount and distribution of mid-tropospheric moisture determines how rainfall associated with secondary cells develop. Rainfall amount and efficiency varied significantly over an observable range of relative humidities in the 850- 500 mb layer even though rainfall evolution was similar during the initial or "first-cell" period. Rainfall variability was attributed to drier mid-tropospheric environments inhibiting secondary cell development through entrainment effects. Observationally, 850-500 mb moisture structure exhibits wider variability than lower level moisture, which is virtually always present in Florida. A likely consequence of the variability in 850-500 moisture is a stronger statistical correlation to rainfall, which observational studies have noted. The study indicates that vertical moisture flux forcing at convergence zones is critical in determining rainfall in the initial stage of development but plays a decreasing role in rainfall evolution as the system matures. The mid-tropospheric moisture (e.g. environment) plays an increasing role in rainfall evolution as the system matures. This suggests the need to improve measurements of magnitude/depth of convergence and mid-tropospheric moisture distribution. It also highlights the need for better parameterization of entrainment and vertical moisture distribution in larger-scale models.

Shepherd, J. Marshall↗

High-Resolution Simulation of Hurricane Bonnie (1998): The Organization of Vertical Motion - Part 1

Hurricanes are well known for their strong winds and heavy rainfall, particularly in the intense rainband (eyewall) surrounding the calmer eye of the storm. In some hurricanes, the rainfall is distributed evenly around the eye so that it has a donut shape on radar images. In other cases, the rainfall is concentrated on one side of the eyewall and nearly absent on the other side and is said to be asymmetric. This study examines how the vertical air motions that produce the rainfall are distributed within the eyewall of an asymmetric hurricane and the factors that cause this pattern of rainfall. We use a sophisticated numerical forecast model to simulate Hurricane Bonnie, which occurred in late August of 1998 during a special NASA field experiment designed to study hurricanes. The simulation results suggest that vertical wind shear (a rapid change in wind speed or direction with height) caused the asymmetric rainfall and vertical air motion patterns by tilting the hurricane vortex and favoring upward air motions in the direction of tilt. Although the rainfall in the hurricane eyewall may surround more than half of the eye, the updrafts that produce the rainfall are concentrated in very small-scale, intense updraft cores that occupy only about 10% of the eyewall area. The model simulation suggests that the timing and location of individual updraft cores are controlled by intense, small-scale vortices (regions of rapidly swirling flow) in the eyewall and that the updrafts form when the vortices encounter low-level air moving into the eyewall.

Braun, Scott A.↗

Over-Luminous Elliptical Galaxies

The first paper from our work has been completed and accepted for publication. Another paper presents a study of the ESO 30601 70 galaxy group, combining Chandra, XMM-Newton, and optical observations. We find that the system is a true fossil galaxy group - a group whose optical light is dominated by a single galaxy. The group X-ray emission is composed of a central, dense, cool core (10 kpc in radius) and an isothermal medium beyond the central 10 kpc. The region between 10 and 50 kpc (the cooling radius) has the same temperature as the gas from 50 to 400 kpc, although the gas cooling time between 10 and 50 kpc (2-6 Gyr) is shorter than the Hubble time. Thus, the ESO 3060170 group does not have a group-sized cooling core. We suggest that the group cooling core may have been heated by a central active galactic nucleus (AGN) outburst in the past and that the small, dense, cool core is the truncated relic of a previous cooling core. The Chandra observations also reveal a variety of X-ray features in the central region, including a finger, an edge-like feature, and a small tail, all aligned along a north-south axis, as are the galaxy light and group galaxy distribution. The proposed AGN outburst may cause gas to slosh around the center and produce these asymmetric features. The observed flat temperature profile to 1/3rvir is not consistent with the predicted temperature profile in recent numerical simulations. We compare the entropy profile of the ESO 3060170 group with those of three other groups and find a flatter relation than that predicted by simulations involving only shock heating, S approximately r approximately 0.85. This is direct evidence of the importance of non-gravitational processes in group centers. We derive the mass profiles within 1/3rvir and find that the ESO 3060170 group is the most massive fossil group known.

Forman, William↗

Multi-Scale Modeling of Low-Density Carbon-Phenolic Ablators

Protecting a spacecraft during atmospheric entry is one of highest risk factors that needs to be mitigated during design of a space exploration mission. At entry speeds from space, air turns into high-temperature plasma, and spacecraft Thermal Protection Systems (TPS) are needed to protect the vehicle payload. Modern successful material architectures of spacecraft shields use a porous carbon fiber substrate impregnated with phenolic as an ablator material. In the lecture, efforts to build a Predictive Material Modeling framework for porous ablators from micro-scale to macro-scale will be presented. Several numerical methods and techniques will be summarized that use voxelized images to compute geometrical properties of the porous substrate. These computed properties include porosity, specific surface area and tortuosity that are otherwise indirectly measured through experimental techniques. Direct simulation Monte Carlo (DSMC), a particle-based method for approximating the Boltzmann equation, is used to compute the permeability coefficient of the porous substrate based on its digitized representation. The method computes the flow within the microstructure, where the size of the pores may approach the mean-free-path of the flow. Finally, a high-fidelity model implemented in PATO (Porous-material Analysis Toolbox) is discussed, and some examples of ablative material response are presented, including for the first time 3D simulations of the full tiled heat shield for the Mars Science Laboratory (MSL) capsule.

Ablators↗

Analysis of High Contrast Imaging Techniques for Space Based Direct Planetary Imaging

We report on our ongoing investigations of a number of techniques for direct detection and imaging of Earth-like planets around nearby stellar sources. Herein, we give a quantitative analysis of these techniques and compare and contrast them via computer simulations. The techniques we will be reporting on are Bracewell Interferometry, Nisenson Apodized Square Aperture, and Coronagraphic masking techniques. We parameterize our results with respect to wavelength, aperture size, effects of mirror speckle, both mid- and high-spatial frequency, detector and photon noise as well pointing error. The recent numerous detections of Jupiter and Saturn like planets has driven a resurgence in research of space based high contrast imaging techniques for direct planetary imaging. Work is currently ongoing for concepts for NASA's Terrestrial Planet Finder mission and a number of study teams have been funded. The authors are members of one team.

Lyon, Richard G.↗

Numerical simulation of the geographical sources of water for Continental Scale Experiments (CSEs) Precipitation

There are several important research questions that the Global Energy and Water Cycle Experiment (GEWEX) is actively pursuing, namely: What is the intensity of the water cycle and how does it change? And what is the sustainability of water resources? Much of the research to address these questions is directed at understanding the atmospheric water cycle. In this paper, we have used a new diagnostic tool, called Water Vapor Tracers (WVTs), to quantify the how much precipitation originated as continental or oceanic evaporation. This shows how long water can remain in the atmosphere and how far it can travel. The model-simulated data are analyzed over regions of interest to the GEWEX community, specifically, their Continental Scale Experiments (CSEs) that are in place in the United States, Europe, Asia, Brazil, Africa and Canada. The paper presents quantitative data on how much each continent and ocean on Earth supplies water for each CSE. Furthermore, the analysis also shows the seasonal variation of the water sources. For example, in the United States, summertime precipitation is dominated by continental (land surface) sources of water, while wintertime precipitation is dominated by the Pacific Ocean sources of water. We also analyze the residence time of water in the atmosphere. The new diagnostic shows a longer residence time for water (9.2 days) than more traditional estimates (7.5 days). We emphasize that the results are based on model simulations and they depend on the model s veracity. However, there are many potential uses for the new diagnostic tool in understanding weather processes and large and small scales.

Bosilovich, Michael G.↗