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,441 records · Page 80

Computerized Generation and Simulation of Meshing and Contact of New Type of Novikov-Wildhaber Helical Gears

In this report, we propose a new geometry for low-noise, increased-strength helical gears of the Novikov-Wildhaber type. Contact stresses are reduced as a result of their convex-concave gear tooth surfaces. The gear tooth surfaces are crowned in the profile direction to localize bearing contact and in the longitudinal direction to obtain a parabolic function of transmission errors. Such a function results in the reduction of noise and vibrations. Methods for the generation of the proposed gear tooth surfaces by grinding and hobbing are considered, and a tooth contact analysis (TCA) computer program to simulate meshing and contact is applied. The report also investigates the influence of misalignment on transmission errors and shift of bearing contact. Numerical examples to illustrate the developed approaches are proposed. The proposed geometry was patented by Ford/UIC (Serial Number 09-340-824, pending) on June 28, 1999.

Litvin, Faydor L.↗

Development of a Hybrid RANS/LES Method for Compressible Mixing Layer Simulations

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modem day aircraft and also those of hypersonic vehicles currently under development. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS/LES method on stretched, non-Cartesian grids. The hybrid RANS/LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two-dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Actual LES calculations, performed in three spatial directions, indicated an initial vortex shedding followed by rapid transition to turbulence, which is in agreement with experimental observations.

Georgiadis, Nicholas J.↗

Quantification of Ice Accretions for Icing Scaling Evaluations

The comparison of ice accretion characteristics is an integral part of aircraft icing research. It is often necessary to compare an ice accretion obtained from a flight test or numerical simulation to one produced in an icing wind tunnel or for validation of an icing scaling method. Traditionally, this has been accomplished by overlaying two-dimensional tracings of ice accretion shapes. This paper addresses the basic question of how to compare ice accretions using more quantitative methods. For simplicity, geometric characteristics of the ice accretions are used for the comparison. One method evaluated is a direct comparison of the percent differences of the geometric measurements. The second method inputs these measurements into a fuzzy inference system to obtain a single measure of the goodness of the comparison. The procedures are demonstrated by comparing ice shapes obtained in the Icing Research Tunnel at NASA Glenn Research Center during recent icing scaling tests. The results demonstrate that this type of analysis is useful in quantifying the similarity of ice accretion shapes and that the procedures should be further developed by expanding the analysis to additional icing data sets.

Ruff, Gary A.↗

Comparison of Meteoroid Flux Models for Near Earth Space

Over the last decade several new models for the sporadic interplanetary meteoroid flux have been developed. These include the Meteoroid Engineering Model (MEM), the Divine-Staubach model and the Interplanetary Meteoroid Engineering Model (IMEM). They typically cover mass ranges from 10-12 g (or lower) to 1 g and are applicable for model specific sun distance ranges between 0.2 A.U. and 10 A.U. Near 1 A.U. averaged fluxes (over direction and velocities) for all these models are tuned to the well established interplanetary model by Gr?n et. al. However, in many respects these models differ considerably. Examples are the velocity and directional distributions and the assumed meteoroid sources. In this paper flux predictions by the various models to Earth orbiting spacecraft are compared. Main differences are presented and analysed. The persisting differences even for near Earth space can be seen as surprising in view of the numerous ground based (optical, radar) and in-situ (captured IDPs, in-situ detectors and analysis of retrieved hardware) measurements and simulations. Remaining uncertainties and potential additional studies to overcome the existing model discrepancies are discussed.

Drolshagen, G.↗

Compiler-directed cache management in multiprocessors

The necessity of finding alternatives to hardware-based cache coherence strategies for large-scale multiprocessor systems is discussed. Three different software-based strategies sharing the same goals and general approach are presented. They consist of a simple invalidation approach, a fast selective invalidation scheme, and a version control scheme. The strategies are suitable for shared-memory multiprocessor systems with interconnection networks and a large number of processors. Results of trace-driven simulations conducted on numerical benchmark routines to compare the performance of the three schemes are presented.

Cheong, Hoichi↗

Structural Optimization for Reliability Using Nonlinear Goal Programming

This report details the development of a reliability based multi-objective design tool for solving structural optimization problems. Based on two different optimization techniques, namely sequential unconstrained minimization and nonlinear goal programming, the developed design method has the capability to take into account the effects of variability on the proposed design through a user specified reliability design criterion. In its sequential unconstrained minimization mode, the developed design tool uses a composite objective function, in conjunction with weight ordered design objectives, in order to take into account conflicting and multiple design criteria. Multiple design criteria of interest including structural weight, load induced stress and deflection, and mechanical reliability. The nonlinear goal programming mode, on the other hand, provides for a design method that eliminates the difficulty of having to define an objective function and constraints, while at the same time has the capability of handling rank ordered design objectives or goals. For simulation purposes the design of a pressure vessel cover plate was undertaken as a test bed for the newly developed design tool. The formulation of this structural optimization problem into sequential unconstrained minimization and goal programming form is presented. The resulting optimization problem was solved using: (i) the linear extended interior penalty function method algorithm; and (ii) Powell's conjugate directions method. Both single and multi-objective numerical test cases are included demonstrating the design tool's capabilities as it applies to this design problem.

El-Sayed, Mohamed E.↗

Experimental Results for Titan Aerobot Thermo-Mechanical Subsystem Development

This paper describes experimental results from a development program focused in maturing Titan aerobot technology in the areas of mechanical and thermal subsystems. Results from four key activities are described: first, a cryogenic balloon materials development program involving coupon and cylinder tests and culminating in the fabrication and testing of an inflated 4.6 m long prototype blimp at 93 K; second, a combined lab experiment and numerical simulation effort to assess potential problems resulting from radioisotope thermal generator waste heat generation near an inflated blimp; third, an aerial deployment and inflation development program consisting of laboratory and helicopter drop tests on a near full scale (11 m long) prototype blimp; and fourth, a proof of concept experiment demonstrating the viability of using a mechanically steerable high gain antenna on a floating blimp to perform direct to Earth telecommunications from Titan. The paper provides details on all of these successful activities and discusses their impact on the overall effort to produce mature systems technology for future Titan aerobot missions.

balloon↗

Wind Lidars for Aero-Assisted Entry, Descent, and Landing on Mars

NASA seeks to safely and consistently deliver 20 metric tons of payload to within 50m of the intended location on Mars. A study has been conducted to evaluate the utility of wind lidars to aid in aero-assisted entry, descent, and precision landing of vehicles carrying these payloads. Numerical simulation found that coherent-Doppler, infrared, aerosol-backscatter wind lidars on the Martian surface can measure winds nominally over a 15km radius hemisphere using eye-safe laser energies and optical apertures similar to those commercially available for terrestrial airport support. Such ground-based wind measurements around the target delivery site are useful for determining when to initiate atmospheric entry. Initial analyses of direct-Doppler, ultraviolet, molecular-backscatter wind lidars demonstrateforward-looking measurement of wind speed, atmospheric temperature, atmospheric density, and vehicle flight attitude beyond the vehicle boundary layer. These measurements would enable controlled flight of an aero-assisted cargo vehicle to the designated landing site.

Mars↗

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

MODIS Aerosol Observations used to Constrain Dust Distributions and Lifecycle in the NASA GEOS-5 Model

Approximately 240 Tg of mineral dust aerosol are transported annually from Saharan Africa to the Atlantic Ocean. Dust affects the Earth radiation budget, and plays direct (through scattering and absorption of radiation) and indirect (through modification of cloud properties and environment) roles in climate. Deposition of dust to the surface provides an important nutrient source to terrestrial and oceanic ecosystems. Dust is additionally a contributor to adverse air quality. Among the tools toward understanding the lifecycle and impacts of mineral dust aerosols are numerical models. Important constraints on these models come from quantitative satellite observations, like those from the space-based Moderate Resolution Imaging Spectroradiometer (MODIS). In particular, Kauhan et al. [2005] used MODIS aerosol observations to infer transport and deposition fluxes of Saharan dust over the Atlantic, Caribbean, and Amazonian basins. Those observations are used here to constrain the transport of dust and its interannual variability simulated in the NASA GEOS-5 general circulation model and data assimilation system. Significant uncertainty exists in the MODIS-derived fluxes, however, due to uncertainty in the wind fields provided by meteorological analyses in this region. That same uncertainty in the wind fields is manifest in our GEOS-5 simulations of dust distributions. Here we use MODIS observations to investigate the seasonality and location of the Saharan dust plume and explore through sensitivity analysis of our model the meteorological controls on the dust distribution, including dust direct radiative effects and sub-gridscale source and sink processes.

Colarco, P.↗

Scale disparity and spectral transfer in anisotropic numerical turbulence

To study the effect of cancellations within long-range interactions on local isotropy at the small scales, we calculate explicitly the degree of cancellation in distant interactions in the simulations of Yeung & Brasseur and Yeung, Brasseur & Wang using the single scale disparity parameter 's' developed by Zhou. In the simulations, initially isotropic simulated turbulence was subjected to coherent anisotropic forcing at the large scales and the smallest scales were found to become anisotropic as a consequence of direct large-small scale couplings. We find that the marginally distant interactions in the simulation do not cancel out under summation and that the development of small-scale anisotropy is indeed a direct consequence of the distant triadic group, as argued by Yeung, et. al. A reduction of anisotropy at later times occurs as a result of the isotropizing influences of more local energy-cascading triadic interactions. Nevertheless, the local-to-nonlocal triadic group persists as an isotropizing influence at later times. We find that, whereas long-range interactions, in general, contribute little to net energy transfer into or out of a high wavenumber shell k, the anisotropic transfer of component energy within the shell increases with increasing scale separations. These results are consistent with results by Zhou, and Brasseur & Wei, and suggest that the anisotropizing influences of long range interactions should persist to higher Reynolds numbers. The residual effect of the forced distant group in this low-Reynolds number simulation is found to be forward cascading, on average.

Zhou, YE↗

Numerical simulation of transition in a decelerating boundary layer

Transition in a decelerating flat-plate boundary layer is numerically simulated up to the beginning of three-dimensional breakdown, and the results are compared with an experiment. The adverse pressure gradient induced by deceleration increases the growth rate of disturbances and allows transition at lower Reynolds numbers. The primary instability is characterized by a wave packet, which undergoes three-dimensional distortion. Lambda vortices are locally observed, but they are not aligned with respect to the flow direction.

Yang, Kyung Soo↗

Formation of dwarf galaxies in tidal tails

The results are reported of numerical simulations of encounters between disk galaxies, each modeled with a central bulge, an exponential disk, and a spheroidal dark-matter halo. It is found that dwarf systems form in material drawn out during the encounter; these objects can capture large amounts of moderately enriched gas but retain little dark matter from their parents' haloes. They should therefore have lower mass-to-light ratios than galaxies formed directly by the collapse of primordial matter.

Barnes, Joshua E.↗

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗

The Interaction of the Solar Wind with Solar Probe Plus - 3D Hybrid Simulation. Report 2: The Study for the Distance 9.5Rs

Our paper is a 2.5D and 3D numerical plasma models of the interaction of the solar wind (SW) with the Solar Probe Plus spacecraft (SPPSC). These results should be interpreted as a basic plasma model for which the derived SW interaction with spacecraft (SC) could have consequences for both plasma wave and electron plasma measurements on board SC in the inner heliosphere. We observe an excitation of the low frequency Alfven and whistler type wave directed by the magnetic field with an amplitude of the electromagnetic field oscillation about of (0.015-0.06) V/m. The compression waves and the jumps in an electric field with an amplitude of about 1.5 V/m and (12-18) V/m were also observed. The observed strong electromagnetic perturbations may be a crucial point in the electromagnetic measurements, which were planned in future Solar Probe Plus mission.

Lipatov, Alexander S.↗

Wind-Optimal Trajectories for Multirotor eVTOL Aircraft on UAM Missions

Distributed electric propulsion powered electric vertical takeoff and landing aircraft are expected to enable urban air mobility. Low specific energy of onboard lithium-ion polymer batteries and wind conditions impose constraints on flight endurance. Therefore, from the safety and efficiency perspective, planning and flying minimum energy trajectories are important.The primary motivation for this paper is to determine if there is an operational benefit to flying wind-optimal trajectories for short flights (less than 60 miles) in the urban environment. This study employs wind-optimal trajectories for a NASA-proposed conceptual multirotor aircraft flight in the urban environment. The optimal control model presented includes a wind model for quantifying the effect of wind on the trajectory. The optimal control problem is numerically solved using the direct method. Energy consumption and flight duration flying wind-optimal trajectories between origins and destinations in a metropolitan area are compared with energy consumption and flight duration flying on the great-circle paths between the same origin and destination pairs to determine the operational benefit of wind-optimal routing for short flights.The flight duration results for different scenarios are validated using a simulation tool designed and developed by NASA for exploring advanced air traffic management concepts. In summary,this research study suggests that for short flights in an urban environment, the wind-optimal trajectories have a very slight (insignificant) operational advantage over the great-circle trajectories in terms of energy consumption and flight duration. As expected, headwinds and tailwinds along the route of flight have a significant impact on energy consumption and flight duration.

Urban Air Mobility↗

Numerical simulations of convectively generated stratospheric gravity waves

A 2D model of a mesoscale convective flow is used to simulate the excitation and vertical propagation of gravity waves. Data obtained show that, in the absence of storm-relative mean winds in the stratosphere, the primary mode of excitation of gravity waves is by mechanical forcing owing to oscillatory updrafts. The stratospheric response consists of waves whose periods match the primary periods of the forcing. Due to the tendency of the oscillating updrafts to propagate toward the rear of the storm, gravity wave propagation is limited primarily to the rearward direction. Results suggest that squall-line-generated gravity waves arise from mechanical forcing rather than thermal effects.

Fovell, R.↗

Identification of structural parameters from helicopter dynamic test data

A method is presented for obtaining the mass, stiffness, and damping parameters of a linear mathematical model, having fewer degrees of freedom than the structure it represents, directly from dynamic response measurements on the actual helicopter without a priori knowledge of the physical characteristics of the fuselage. The only input information required in the formulation is the approximate natural frequency of each mode and mobility data measured proximate to these frequencies with sinusoidal force excitation applied at only one point on the vehicle. The practicality and numerical soundness of the theoretical development was demonstrated through a computer simulation of an experimental program.

Giansante, N.↗