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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Absorbing Boundary Conditions For Optical Pulses In Dispersive, Nonlinear Materials

This paper will present results in computational nonlinear optics. An algorithm will be described that provides absorbing boundary conditions for optical pulses in dispersive, nonlinear materials. A new numerical absorber at the boundaries has been developed that is responsive to the spectral content of the pulse. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of "light bullet" like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. Comparisons will be shown of calculations that use the standard boundary conditions and the new ones.

Goorjian, Peter M.↗
NASA NTRS2007

Software Tool for Computing Maximum Von Mises Stress

The maximum Van Mises stress and stress direction are of interest far analyzing launch accelerations such as with the Mass Acceleration Curves developed by JPL. Maximum launch stresses can be combined with appropriate load cases at consistent locations with resulting stress tensors. Maximum Van Mises stress is also of interest for understanding maximum operational loading such as traverse events. - For example, planetary traversing simulations may prescribe bounding acceleration values during traverse for a rover such as Mars Science Lab (MSL) in (X,Y,Z) of the rover. - Such accelerations can be really in any directions for many parts such as a mast or head mounted components which can be in numerous configurations and orientations when traversing a planet surface.

software↗
NASA NTRS2016

On the Assessment of Acoustic Scattering and Shielding by Time Domain Boundary Integral Equation Solutions

Based on the time domain boundary integral equation formulation of the linear convective wave equation, a computational tool dubbed Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) has recently been under development. The time domain approach has a distinct advantage that the solutions at all frequencies are obtained in a single computation. In this paper, the formulation of the integral equation, as well as its stabilization by the Burton-Miller type reformulation, is extended to cases of a constant mean flow in an arbitrary direction. In addition, a "Source Surface" is also introduced in the formulation that can be employed to encapsulate regions of noise sources and to facilitate coupling with CFD simulations. This is particularly useful for applications where the noise sources are not easily described by analytical source terms. Numerical examples are presented to assess the accuracy of the formulation, including a computation of noise shielding by a thin barrier motivated by recent Historical Baseline F31A31 open rotor noise shielding experiments. Furthermore, spatial resolution requirements of the time domain boundary element method are also assessed using point per wavelength metrics. It is found that, using only constant basis functions and high-order quadrature for surface integration, relative errors of less than 2% may be obtained when the surface spatial resolution is 5 points-per-wavelength (PPW) or 25 points-per-wavelength squared (PPW2).

Hu, Fang Q.↗
NASA NTRS2022

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

supermassive black holes↗
NASA NTRS1985

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗
NASA NTRS1994

Numerical Simulations of Light Bullets, Using the Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗
NASA NTRS1995

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that we currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Karr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗
NASA NTRS1994

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗
NASA NTRS2019

A new environment to simulate the dynamics in the close proximity of rubble-pile asteroids

This paper presents a new environment to simulate close-proximity dynamics around rubble-pile asteroids. The code provides methods for modeling the asteroid’s gravity field and surface through granular dynamics. It implements stateof-the-art techniques to model both gravity and contact interaction between particles: 1) mutual gravity as either direct N2 or Barnes-Hut GPU-parallel octree and 2) contact dynamics with a soft-body (force-based, smooth dynamics), hard-body (constraint-based, non-smooth dynamics), or hybrid (constraint-based with compliance and damping) approach. A very relevant feature of the code is its ability to handle complex-shaped rigid bodies and their full 6D motion. Examples of spacecraft close-proximity scenarios and their numerical simulations are shown.

Ferrari, Fabio↗
NASA NTRS1977

Development and application of a three dimensional numerical model for predicting pollutant and sediment transport using an Eulerian-Lagrangian marker particle technique

A computer coded Lagrangian marker particle in Eulerian finite difference cell solution to the three dimensional incompressible mass transport equation, Water Advective Particle in Cell Technique, WAPIC, was developed, verified against analytic solutions, and subsequently applied in the prediction of long term transport of a suspended sediment cloud resulting from an instantaneous dredge spoil release. Numerical results from WAPIC were verified against analytic solutions to the three dimensional incompressible mass transport equation for turbulent diffusion and advection of Gaussian dye releases in unbounded uniform and uniformly sheared uni-directional flow, and for steady-uniform plug channel flow. WAPIC was utilized to simulate an analytic solution for non-equilibrium sediment dropout from an initially vertically uniform particle distribution in one dimensional turbulent channel flow.

Pavish, D. L.↗
NASA NTRS1986

Zero-Lift Wave Drag of Complex Aircraft Configurations

WAVDRAG calculates supersonic zero-lift wave drag of complex aircraft configurations. Numerical model of aircraft used throughout design process from concept to manufacturing. Incorporates extended geometric input capabilities to permit use of more accurate mathematical model. Engineer defines aircraft components as fusiform or nonfusiform in terms of traditional parallel contours or nonintersecting contours in any direction. Laterally asymmetric configurations simulated. WAVDRAG calculates total drag and wave-drag coefficient of specified aircraft configuration.

Craidon, C., B.↗
NASA NTRS1992

Development of magnetohydrodynamic turbulence in coronal loops

It is proposed that a stationary spectrum of magnetohydrodynamic turbulence is generated in solar coronal loops by the interweaving of magnetic field lines driven by the turbulent velocity field in the convective region. Starting from the set of incompressible magnetohydrodynamic equations, two ordinary differential equations describing the stationary turbulent energy spectrum and the frequency of direct energy cascade for a given kinetic energy spectrum of the subphotospheric convective motions are derived. An estimate of the typical dissipation length scale, at which magnetic energy efficiently heats the plasma by Joule dissipation is obtained. The impact of recent results of two-dimensional magnetohydrodynamic numerical simulations on the coronal heating problem is also discussed.

Gomez, Daniel O.↗
NASA NTRS1993

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗
NASA NTRS2004

Reliability Quantification of the Flexure: A Critical Stirling Convertor Component

Uncertainties in the manufacturing, fabrication process, material behavior, loads, and boundary conditions results in the variation of the stresses and strains induced in the flexures and its fatigue life. Past experience and the test data at material coupon levels revealed a significant amount of scatter of the fatigue life. Owing to these facts, the design of the flexure, using conventional approaches based on safety factor or traditional reliability based on similar equipment considerations does not provide a direct measure of reliability. Additionally, it may not be feasible to run actual long term fatigue tests due to cost and time constraints. Therefore it is difficult to ascertain material fatigue strength limit. The objective of the paper is to present a methodology and quantified results of numerical simulation for the reliability of flexures used in the Stirling convertor for their structural performance. The proposed approach is based on application of finite element analysis method in combination with the random fatigue limit model, which includes uncertainties in material fatigue life. Additionally, sensitivity of fatigue life reliability to the design variables is quantified and its use to develop guidelines to improve design, manufacturing, quality control and inspection design process is described.

Shah, Ashwin R.↗
NASA NTRS2005

NASA's Global Precipitation Mission Ground Validation Segment

NASA is designing a Ground Validation Segment (GVS) as one of its contributions to the Global Precipitation Measurement (GPM) mission. The GPM GVS provides an independent means for evaluation, diagnosis, and ultimately improvement of the GPM spaceborne measurements and precipitation products. NASA's GPM GVS concept calls for a combination of direct observations executed within a Multidimensional Observing Volume (MOV) and model-based analyses executed by a Satellite Simulator Model (SSM). The MOV consists of ground-based instruments that measure local surface and atmospheric properties required for GPM validation. The SSM utilizes MOV measurements in a forward numerical model. The goal of the SSM forward modeling is calculation of the following properties: top-of-atmosphere microwave radiative quantities to within sensor noise of those measured by the GPM Core Satellite, precipitation quantities identical to those generated by the standard GPM precipitation retrieval algorithms, and quantitative/objective error estimates of both sets of quantities. At present, the GVS is in the early design stage and various scenarios have been generated to assess how it will be used in the GPM era. The GPM GVS will be operational in the year prior to the launch of the GPM core satellite, which has a launch date scheduled for December 2010.

GVS (GROUND VALIDATION SEGMENT)↗
NASA NTRS1991

Numerical simulation of spatially-evolving instability in plane channel flow

The spatial stability of plane channel flow is analyzed using a three-dimensional, time-dependent spectral/finite difference code (Danabasoglu et al., 1990) which integrates numerically the Navier-Stokes equations. The study centers on inflow disturbance amplitudes effects on the secondary instability. The resolution requirements along the spacewise direction, which become critical before the breakdown stage, are of particular interest. A direct comparison is made with the experiments of Nishioka et al. (1980).

Danabasoglu, G.↗
NASA NTRS2006

Numerical Simulations and an Experimental Investigation of a Finger Seal

Besides sealing, the other main goal of a successful finger seal design is to exhibit appropriate compliance to outside forces. The ability of the seal to ride or float along the rotor without rubbing or excessive heating is essential to the successful operation of the seal. The compliance of the finger must only occur in the radial plane; The seal needs to be as sturdy as possible in the axial direction. The compliant finger that moves radially outward with rotor growth and motion has to be able to ride the rotor back down as the rotor diameter recovers or the rotor moves "away". Thus there is an optimum stiffness for the finger.

Braun, Minel↗
NASA NTRS2012

Evaluation of Cloud Microphysics Simulated using a Meso-Scale Model Coupled with a Spectral Bin Microphysical Scheme through Comparison with Observation Data by Ship-Borne Doppler and Space-Borne W-Band Radars

Equivalent radar reflectivity factors (Ze) measured by W-band radars are directly compared with the corresponding values calculated from a three-dimensional non-hydrostatic meso-scale model coupled with a spectral-bin-microphysical (SBM) scheme for cloud. Three case studies are the objects of this research: one targets a part of ship-borne observation using 95 GHz Doppler radar over the Pacific Ocean near Japan in May 2001; other two are aimed at two short segments of space-borne observation by the cloud profiling radar on CloudSat in November 2006. The numerical weather prediction (NWP) simulations reproduce general features of vertical structures of Ze and Doppler velocity. A main problem in the reproducibility is an overestimation of Ze in ice cloud layers. A frequency analysis shows a strong correlation between ice water contents (IWC) and Ze in the simulation; this characteristic is similar to those shown in prior on-site studies. From comparing with the empirical correlations by the prior studies, the simulated Ze is overestimated than the corresponding values in the studies at the same IWC. Whereas the comparison of Doppler velocities suggests that large-size snowflakes are necessary for producing large velocities under the freezing level and hence rules out the possibility that an overestimation of snow size causes the overestimation of Ze. Based on the results of several sensitivity tests, we conclude that the source of the overestimation is a bias in the microphysical calculation of Ze or an overestimation of IWC. To identify the source of the problems needs further validation research with other follow-up observations.

Iguchi, T.↗