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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Comparison of remote sensing measurements with a two-scale polarimetric emission and scattering model for sea surfaces

Recently, it has been observed that the brightness temperatures of sea surfaces correlate with the azimuth angle of the ocean wind vector, indicating that ocean wind direction can not only be retrieved from the microwave backscatter, but also from the brightness temperature measurements. In this paper, comparison of the theoretical result calculated from a two-scale emission and scattering model with the Seasat and SSM/I model functions and aircraft measurements is presented and potential applications of microwave polarimetry are discussed. In our two-scale model, the modified reflectivities of large scale surfaces are calculated by extending the small perturbation method to the second order for small scale perturbation with anisotropic directional spectrum. It was found that the modified reflectivities derived from the second-order scattered field agree excellently well with the results obtained from a Monte Carlo simulation technique which numerically calculates the polarimetric reflectivities of one-dimensional random rough surfaces with a power-law spectrum. Without the second-order correction, the modified reflectivities of the rough surfaces are significantly over-estimated and sign errors are observed in the third and fourth Stokes parameters for thermal emissions. The surface spectrum parameters and two-scale cutoff are selected so that the calculated scattering and emission signatures agree with the reported model functions for Seasat and SSM/I. Subsequently, the polarimetric signatures of sea surfaces are illustrated to indicate the possibility of reducing the number of azimuthal looks required for spaceborne sensors in the remote sensing of ocean wind by using the polarimetric information. Furthermore, it is found that contrary to the dependence of backscattering coefficients on incidence angles, polarimetric brightness temperatures display a stronger wind direction dependence in the near nadir-looking direction than away from nadir. Finally, we discuss the potential of a single-look polarimetric radiometer for wind retrieval.

Yueh, S. H.↗
NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗
NASA NTRS1976

Turbulence and numerical wind tunnels

The feasibility and possible significance of aerodynamic simulations are considered. Despite certain limitations, there is reason to believe that many useful and important flows containing turbulence can be computed with an accuracy acceptable for use in aircraft designs. Attention is given to the equations for unsteady fluid flow, the direct calculation of unsteady viscous flows, the direct calculation of nonlinear equations, the computation of flows about bodies using turbulence modeling, future developments in numerical methods, and future developments in computers.

Lomax, H.↗
NASA NTRS1990

Two-dimensional simulation of a two-phase, regenerative pumped radiator loop utilizing direct contact heat transfer with phase change

An innovative pumped loop concept for 600 K space power system radiators utilizing direct contact heat transfer, which facilitates repeated startup/shutdown of the power system without complex and time-consuming coolant thawing during power startup, is under development. The heat transfer process with melting/freezing of Li in an NaK flow was studied through two-dimensional time-dependent numerical simulations to characterize and predict the Li/NaK radiator performance during startup (thawing) and shutdown (cold-trapping). Effects of system parameters and the criteria for the plugging domain are presented together with temperature distribution patterns in solid Li and subsequent melting surface profile variations in time.

Rhee, Hyop S.↗
NASA NTRS2019

Mechanism of Reconnection on Kinetic Scales Based on Magnetospheric Multiscale Mission Observations

We examine the role that ions and electrons play in reconnection using observations from the Magnetospheric Multiscale (MMS) mission on kinetic ion and electron scales, which are much shorter than magnetohydrodynamic scales. This study reports observations with unprecedented high resolution that MMS provides for magnetic field (7.8 ms) and plasma (30 ms for electrons and 150 ms for ions). We analyze and compare approaches to the magnetopause in 2016 November, to the electron diffusion region in the magnetotail in 2017 July followed by a current sheet crossing in 2018 July. Besides magnetic field reversals, changes in the direction of the flow velocity, and ion and electron heating, MMS observed large fluctuations in the electron flow speeds in the magnetotail. As expected from numerical simulations, we have verified that when the field lines and plasma become decoupled a large reconnecting electric field related to the Hall current (1–10 mV/m) is responsible for fast reconnection in the ion diffusion region. Although inertial accelerating forces remain moderate (1–2 mV/m), the electric fields resulting from the divergence of the full electron pressure tensor provide the main contribution to the generalized Ohm’s law at the neutral sheet (as large as 200 mV/m). In our view, this illustrates that when ions decouple electron physics dominates. The results obtained on kinetic scales may be useful for better understanding the physical mechanisms governing reconnection processes in various magnetized laboratory and space plasmas.

W. M. Macek↗
NASA NTRS1993

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗
NASA NTRS1990

Advances and trends in computational structures technology

The major goals of computational structures technology (CST) are outlined, and recent advances in CST are examined. These include computational material modeling, stochastic-based modeling, computational methods for articulated structural dynamics, strategies and numerical algorithms for new computing systems, multidisciplinary analysis and optimization. The role of CST in the future development of structures technology and the multidisciplinary design of future flight vehicles is addressed, and the future directions of CST research in the prediction of failures of structural components, the solution of large-scale structural problems, and quality assessment and control of numerical simulations are discussed.

Noor, A. K.↗
NASA NTRS2016

Computational Investigation of a Boundary-Layer Ingestion Propulsion System for the Common Research Model

This thesis will examine potential propulsive and aerodynamic benefits of integrating a boundary-layer ingestion (BLI) propulsion system with a typical commercial aircraft using the Common Research Model geometry and the NASA Tetrahedral Unstructured Software System (TetrUSS). The Numerical Propulsion System Simulation (NPSS) environment will be used to generate engine conditions for CFD analysis. Improvements to the BLI geometry will be made using the Constrained Direct Iterative Surface Curvature (CDISC) design method. Previous studies have shown reductions of up to 25% in terms of propulsive power required for cruise for other axisymmetric geometries using the BLI concept. An analysis of engine power requirements, drag, and lift coefficients using the baseline and BLI geometries coupled with the NPSS model are shown. Potential benefits of the BLI system relating to cruise propulsive power are quantified using a power balance method and a comparison to the baseline case is made. Iterations of the BLI geometric design are shown and any improvements between subsequent BLI designs presented. Simulations are conducted for a cruise flight condition of Mach 0.85 at an altitude of 38,500 feet and an angle of attack of 2deg for all geometries. A comparison between available wind tunnel data, previous computational results, and the original CRM model is presented for model verification purposes along with full results for BLI power savings. Results indicate a 14.3% reduction in engine power requirements at cruise for the BLI configuration over the baseline geometry. Minor shaping of the aft portion of the fuselage using CDISC has been shown to increase the benefit from boundary-layer ingestion further, resulting in a 15.6% reduction in power requirements for cruise as well as a drag reduction of eighteen counts over the baseline geometry.

Blumenthal, Brennan↗
NASA NTRS2025

A Comparison of Bank Control and Direct Force Control for Aerocapture at Uranus

We compare bank angle control and direct force control for aerocapture at Uranus. We utilize the Fully Numerical Predictor-corrector Aerocapture Guidance for both configurations and discuss modifications to the direct force control corrector to handle multiple roots. Using the Genesis flight mechanics simulation, we analyze both nominal and dispersed three-degree-of-freedom trajectories with perfect navigation and control. We present Monte Carlo results and discuss differences between the performance of the two control approaches. We find that increased bank angle rate and acceleration limits do not increase the performance of the bank control configuration. We analyze a variation with tighter angle of attack bounds and find that the tuning and performance of the two configurations converges, because direct force control is less able to modulate drag.

Sergio A Sandoval↗
NASA NTRS2015

Effects of Mesoscale Eddies in the Active Mixed Layer: Test of the Parametrisation in Eddy Resolving Simulations

In eddy resolving simulations, we test a mixed layer mesoscale parametrisation, developed recently by Canuto and Dubovikov [Ocean Model., 2011, 39, 200-207]. With no adjustable parameters, the parametrisation yields the horizontal and vertical mesoscale fluxes in terms of coarse-resolution fields and eddy kinetic energy (EKE). We compare terms of the parametrisation diagnosed from coarse-grained fields with the eddy mesoscale fluxes diagnosed directly from the high resolution model. An expression for the EKE in terms of mean fields has also been found to get a closed parametrisation in terms of the mean fields only. In 40 numerical experiments we simulated two types of flows: idealised flows driven by baroclinic instabilities only, and more realistic flows, driven by wind and surface fluxes as well as by inflow-outflow. The diagnosed quasi-instantaneous horizontal and vertical mesoscale buoyancy fluxes (averaged over 1-2 degrees and 10 days) demonstrate a strong scatter typical for turbulent flows, however, the fluxes are positively correlated with the parametrisation with higher (0.5-0.74) correlations at the experiments with larger baroclinic radius Rossby. After being averaged over 3-4 months, diffusivities diagnosed from the eddy resolving simulations are consistent with the parametrisation for a broad range of parameters. Diagnosed vertical mesoscale fluxes restratify mixed layer and are in a good agreement with the parametrisation unless vertical turbulent mixing in the upper layer becomes strong enough in comparison with mesoscale advection. In the latter case, numerical simulations demonstrate that the deviation of the fluxes from the parametrisation is controlled by dimensionless parameter estimating the ratio of vertical turbulent mixing term to mesoscale advection. An analysis using a modified omega-equation reveals that the effects of the vertical mixing of vorticity is responsible for the two-three fold amplification of vertical mesoscale flux. Possible physical mechanisms, responsible for the amplification of vertical mesoscale flux are discussed.

simulation↗
NASA NTRS2013

Numerical Simulations of the Partially-Ionized Gas in a 100-A LaB6 Hollow Cathode

Numerical simulations of a hollow cathode with a LaB6 emitter operating at 100 A have been performed for the first time using the 2-D Orificed Cathode (OrCa2D) code. Results for a variety of plasma properties are presented and compared with laboratory measurements. The large size of the device permits peak electron number densities in the cathode interior that are lower than those established in the NSTAR hollow cathode, which operates with a 7.3x lower discharge current and 3.2x lower mass flow rate. Also, despite the higher discharge current in the LaB6 cathode, the maximum electron current density is lower, by 4.2x, than that in the NSTAR cathode due to the larger orifice size. Simulations and direct measurements show that at 12 sccm of xenon flow the peak emitter temperature is in the range of 1594-1630 C. It is also found that the conditions for the excitement of current-driven streaming instabilities and ion-acoustic turbulence (IAT) are satisfied in this cathode, similarly to what was found in the past in its smaller counterparts like the NSTAR cathode. Based on numerical simulations, it has long been argued that these instabilities may be responsible for the anomalously large ion energies that have been measured in these discharges as well as for the enhancement of the plasma resistivity. Confirmation of the presence of IAT in this cathode is presented for the first time in a companion paper.

robotic exploration↗
NASA NTRS1993

Three-dimensional numerical simulation of near-surface flows over the Martian north polar cap

Measurements made by Viking Lander VL-2 (48 N) have shown that the near-surface wind and temperature regime on Mars displays striking similarities to terrestrial counterparts. The diurnal radiative cycle is responsible for establishment of a well-defined thermal circulation in which downslope (Katabatic) flows prevail during the nighttime hours and weak upslope (anabatic) conditions prevail during the daytime. Previous work has indicated that the slope flows are much like those found on Earth, particularly the Katabatic winds, which show striking similarities to drainage flows observed over Antarctica. The low-level wind regime appears to be an important factor in the scouring of the martian landscape. The north polar cap shows evidence of eolian features such as dunes, frost streaks, and grooves from Viking imagery. The direction of the prevailing wind can in cases be inferred from the eolian features. We examine the thermally induced flows that result from the radiative heating and cooling of the martian north polar region using a comprehensive three-dimensional atmospheric mesoscale numerical model. The same model has been used previously for simulation of Antarctic Katabatic winds. The model equations are written in terrain-following coordinates to allow for irregular terrain; prognostic equations include the flux forms of the horizontal momentum equations, temperature, continuity. A surface energy budget equation is also incorporated in which the surface temperature is determined. Explicit parameterization of both terrestrial (longwave) and solar (shortwave) radiation is included. Turbulent transfer of heat and momentum in the martian atmosphere is assumed to follow the similarity expressions in the surface boundary layer on Earth. The terrain heights for the martian north polar region have been obtained from the U.S. Geological Survey map and digitized onto a 57x57 grid with a spacing of 75 km. The resulting terrain map is shown in Fig. 1. The vertical grid consists of 15 levels (delta equals 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere is necessary to capture details of the boundary layer flows. The lowest level corresponds to a height of approximately 20 m above the ground, the seconds level 100 m.

Parish, Thomas R.↗
NASA NTRS1985

Middle atmosphere tides

The recent direction of theoretical and observational research in middle atmosphere tides is reviewed. The diurnal tide has received considerable emphasis, including observations of evanescent components at high latitudes, numerical simulations of nonmigrating propagating components due to longitudinally varying water vapor insolation absorption in the troposphere, and implications and origins of unsteadiness in tidal oscillations over periods of days. Recent analyses of Nimbus 7 LIMS temperatures and partial reflection drift measurements of neutral winds in the tropical mesosphere have also emphasized the diurnal propagating tide. A description of the Atmospheric Tides Middle Atmosphere Program (ATMAP) and an overview of ATMAP campaign results are also presented.

Forbes, J. M.↗
NASA NTRS1985

Samics on the IBM PC

An updated version of the Solar Array Manufacturing Industry Costing Standards (SAMICS) was presented. This version will run on the IBM PC-XT or compatibles, and embodies user friendly input screens and numerous help options. Use of this microcomputer version will still produce the old main frame direct cost driven analysis, although at some cost in turnaround time. A typical simulation will take about 4 hours. The program allows for unattended report printout. Also available is the IBM PC (or compatibles) version of the Improved Price Estimation Guidelines (IPEG). This program allows rapid analysis of process variables using SAMICS-generated coefficients.

Chamberlain, R.↗
NASA NTRS1988

Design for performance enhancement in feedback control systems with multiple saturating nonlinearities

A systematic control design methodology is introduced for multi-input/multi-output systems with multiple saturations. The methodology can be applied to stable and unstable open loop plants with magnitude and/or rate control saturations and to systems in which state limitations are desired. This new methodology is a substantial improvement over previous heuristic single-input/single-output approaches. The idea is to introduce a supervisor loop so that when the references and/or disturbances are sufficiently small, the control system operates linearly as designed. For signals large enough to cause saturations, the control law is modified in such a way to ensure stability and to preserve, to the extent possible, the behavior of the linear control design. Key benefits of this methodology are: the modified compensator never produces saturating control signals, integrators and/or slow dynamics in the compensator never windup, the directional properties of the controls are maintained, and the closed loop system has certain guaranteed stability properties. The advantages of the new design methodology are illustrated by numerous simulations, including the multivariable longitudinal control of modified models of the F-8 (stable) and F-16 (unstable) aircraft.

Kapasouris, Petros↗
NASA NTRS2008

Explaining Warm Coronal Loops

One of the great mysteries of coronal physics that has come to light in the last few years is the discovery that warn (- 1 INK) coronal loops are much denser than expected for quasi-static equilibrium. Both the excess densities and relatively long lifetimes of the loops can be explained with bundles of unresolved strands that are heated impulsively to very high temperatures. Since neighboring strands are at different stages of cooling, the composite loop bundle is multi-thermal, with the distribution of temperatures depending on the details of the "nanoflare storm." Emission hotter than 2 MK is predicted, but it is not clear that such emission is always observed. We consider two possible explanations for the existence of over-dense warm loops without corresponding hot emission: (1) loops are bundles of nanoflare heated strands, but a significant fraction of the nanoflare energy takes the form of a nonthermal electron beam rather then direct plasma heating; (2) loops are bundles of strands that undergo thermal nonequilibrium that results when steady heating is sufficiently concentrated near the footpoints. We present numerical hydro simulations of both of these possibilities and explore the observational consequences, including the production of hard X-ray emission and absorption by cool material in the corona.

Klimchuk, James A.↗
NASA NTRSDate not supplied

Comparing 3D and 2D CFD for Mars Helicopter Ingenuity Rotor Performance Prediction

Single and coaxial rotor performance simulations for the Mars Helicopter Ingenuity rotor are performed for representative Mars atmospheric conditions. Analyses are presented using both a high-fidelity 3D CFD model of the rotor and 2D CFD models of the airfoil sections for comprehensive analyses that use CAMRADII (Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics). When available, the airfoil performance calculations are generated using a numerical approach identical to that used in the high-fidelity 3D model, allowing for a direct comparison between the approaches. Experimental data from a validation campaign to explore higher thrust from an Ingenuity rotor is provided to substantiate a discussion on the simulation fidelity required for both coaxial and single rotor performance predictions. The data is in support of the Sample Recovery Helicopter (SRH) element that serves as the primary backup for tube retrieval as part of the Mars Sample Return (MSR) Campaign. Insights on modeling turbulence at low Reynolds numbers and its influence on the rotor figure of merit are discussed. Key rotor performance metrics are compared. A detailed investigation into differences between 2D and 3D rotor performance predictions, spanwise loading, and rotor stall behavior is included.

Mars Helicopter↗
NASA NTRS1993

Simulation of wave propagation in three-dimensional random media

Quantitative error analysis for simulation of wave propagation in three dimensional random media assuming narrow angular scattering are presented for the plane wave and spherical wave geometry. This includes the errors resulting from finite grid size, finite simulation dimensions, and the separation of the two-dimensional screens along the propagation direction. Simple error scalings are determined for power-law spectra of the random refractive index of the media. The effects of a finite inner scale are also considered. The spatial spectra of the intensity errors are calculated and compared to the spatial spectra of intensity. The numerical requirements for a simulation of given accuracy are determined for realizations of the field. The numerical requirements for accurate estimation of higher moments of the field are less stringent.

Coles, William A.↗