Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical algorithms”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Ballistic Models for Radio Jets in Colliding Galaxies: 3C 278 (NGC 4782/4783)

We present here a general numerical simulation algorithm for ballistic radio jets. The evolution of the jets is determined by their response to mechanical forces (i.e., gravity and ram pressure). Using a previously derived collision model for the interacting pair of elliptical galaxies NGC 4782/4783, we have used this algorithm to simulate the specific two-sided jet morphology seen in the radio source associated with NGC 4782 (3C 278). This is the first time such simulations have been attempted for a galaxy pair whose relative orbit was determined independently from the jet modeling. The general morphological features of the jets in 3C 278 (e.g., position angles, lengths, curvature, and deflection angles) are well matched by the simulations, indicating that our model of the mechanical forces acting on the jets can indeed reproduce most of the basic details of the jet morphology and thereby provide an opportunity to model a combination of detailed optical, radio, and X-ray data for a wide variety of interacting radio-jet galaxies. Our models constrain the initial jet parameters (i.e., the ejection speed, direction, and starting time), the properties of the gaseous medium into which the jets are ejected, and the relative importance of gravitational deflection versus ram pressure bending in influencing the jet morphology. For 3C 278, we find that the effects of ram pressure dominate the structural evolution of the jets. In our best-fit simulation the jet ejection speed is of order 10(exp 4) km/s, to within 50%, suggesting that similar low-luminosity radio jets are also significantly subluminal. We also find that the jets are ejected within approximately 5 deg of the line of sight, the ISM in the nonjet galaxy has a much larger effect on the jet deflection than does the host galaxy, and the jet activity began just over 70 million years ago, roughly 50 million years before the pericenter passage of the two galaxies. Therefore, the onset of nuclear activity appears in 3C 278 to be related to the kinematically observed tidal shock that has been induced in the core of NGC 4782 as a result of its deeply penetrating collision with NGC 4783.

Borne, Kirk D.↗
NASA NTRS2003

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The May 2-7, 1998, Storm

Complete description of a self-consistent model for magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves, and back on waves, are considered self-consistently by solving both equations on a global magnetospheric scale under non steady-state conditions. In the paper by Khazanov et al. [2002] this self-consistent model has only been shortly outlined, and discussions of many the model related details have been omitted. For example, in present study for the first time a new algorithm for numerical finding of the resonant numbers for quasilinear wave-particle interaction is described, or it is demonstrated that in order to describe quasilinear interaction in a multi-ion thermal plasma correctly, both e and He(+) modes of electromagnetic ion cyclotron waves should be employed. The developed model is used to simulate the entire May 2-7, 1998 storm period. Trapped number fluxes of the ring current protons are calculated and presented along with their comparison with the data measured by the 3D hot plasma instrument Polar/HYDRA. Examining of the wave (MLT, L shell) distributions produced during the storm progress reveals an essential intensification of the wave emissions in about two days after main phase of storm. This result is well consistent with the earlier ground-based observations. Also the theoretical shapes and the occurrence rates for power spectral densities of electromagnetic ion cyclotron waves are studied. It is found that in about 2 days after the storm main phase on May 4, mainly non Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗
NASA NTRS1999

The NOAA-9 Earth Radiation Budget Experiment Wide Field-of-View Data Set

The Earth Radiation Budget Experiment (ERBE) consisted of wide field-of-view (WFOV) radiometers and scanning radiometers for measuring outgoing longwave radiation and solar radiation reflected from the Earth. These instruments were carried by the dedicated Earth Radiation Budget Satellite (ERBS) and by the NOAA-9 and -10 operational spacecraft. The WFOV radiometers provided data from which instantaneous fluxes at the top of the atmosphere (TOA) are computed by use of a numerical filter algorithm. Monthly mean fluxes over a 5-degree equal angle grid are computed from the instantaneous TOA fluxes. The WFOV radiometers aboard the NOAA-9 spacecraft operated from February 1985 through December 1992, at which time a failure of the shortwave radiometer ended the usable data after nearly 8 years. This paper examines the monthly mean products from that data set.

Bush, Kathryn A.↗
NASA NTRS2004

Fourth Computational Aeroacoustics (CAA) Workshop on Benchmark Problems

This publication contains the proceedings of the Fourth Computational Aeroacoustics (CAA) Workshop on Benchmark Problems. In this workshop, as in previous workshops, the problems were devised to gauge the technological advancement of computational techniques to calculate all aspects of sound generation and propagation in air directly from the fundamental governing equations. A variety of benchmark problems have been previously solved ranging from simple geometries with idealized acoustic conditions to test the accuracy and effectiveness of computational algorithms and numerical boundary conditions; to sound radiation from a duct; to gust interaction with a cascade of airfoils; to the sound generated by a separating, turbulent viscous flow. By solving these and similar problems, workshop participants have shown the technical progress from the basic challenges to accurate CAA calculations to the solution of CAA problems of increasing complexity and difficulty. The fourth CAA workshop emphasized the application of CAA methods to the solution of realistic problems. The workshop was held at the Ohio Aerospace Institute in Cleveland, Ohio, on October 20 to 22, 2003. At that time, workshop participants presented their solutions to problems in one or more of five categories. Their solutions are presented in this proceedings along with the comparisons of their solutions to the benchmark solutions or experimental data. The five categories for the benchmark problems were as follows: Category 1:Basic Methods. The numerical computation of sound is affected by, among other issues, the choice of grid used and by the boundary conditions. Category 2:Complex Geometry. The ability to compute the sound in the presence of complex geometric surfaces is important in practical applications of CAA. Category 3:Sound Generation by Interacting With a Gust. The practical application of CAA for computing noise generated by turbomachinery involves the modeling of the noise source mechanism as a vortical gust interacting with an airfoil. Category 4:Sound Transmission and Radiation. Category 5:Sound Generation in Viscous Problems. Sound is generated under certain conditions by a viscous flow as the flow passes an object or a cavity.

Dahl, Milo D.↗
NASA NTRS1997

Visual Computing Environment

The Visual Computing Environment (VCE) is a NASA Lewis Research Center project to develop a framework for intercomponent and multidisciplinary computational simulations. Many current engineering analysis codes simulate various aspects of aircraft engine operation. For example, existing computational fluid dynamics (CFD) codes can model the airflow through individual engine components such as the inlet, compressor, combustor, turbine, or nozzle. Currently, these codes are run in isolation, making intercomponent and complete system simulations very difficult to perform. In addition, management and utilization of these engineering codes for coupled component simulations is a complex, laborious task, requiring substantial experience and effort. To facilitate multicomponent aircraft engine analysis, the CFD Research Corporation (CFDRC) is developing the VCE system. This system, which is part of NASA's Numerical Propulsion Simulation System (NPSS) program, can couple various engineering disciplines, such as CFD, structural analysis, and thermal analysis. The objectives of VCE are to (1) develop a visual computing environment for controlling the execution of individual simulation codes that are running in parallel and are distributed on heterogeneous host machines in a networked environment, (2) develop numerical coupling algorithms for interchanging boundary conditions between codes with arbitrary grid matching and different levels of dimensionality, (3) provide a graphical interface for simulation setup and control, and (4) provide tools for online visualization and plotting. VCE was designed to provide a distributed, object-oriented environment. Mechanisms are provided for creating and manipulating objects, such as grids, boundary conditions, and solution data. This environment includes parallel virtual machine (PVM) for distributed processing. Users can interactively select and couple any set of codes that have been modified to run in a parallel distributed fashion on a cluster of heterogeneous workstations. A scripting facility allows users to dictate the sequence of events that make up the particular simulation.

Lawrence, Charles↗
NASA NTRS2008

AutoBayes Program Synthesis System Users Manual

Program synthesis is the systematic, automatic construction of efficient executable code from high-level declarative specifications. AutoBayes is a fully automatic program synthesis system for the statistical data analysis domain; in particular, it solves parameter estimation problems. It has seen many successful applications at NASA and is currently being used, for example, to analyze simulation results for Orion. The input to AutoBayes is a concise description of a data analysis problem composed of a parameterized statistical model and a goal that is a probability term involving parameters and input data. The output is optimized and fully documented C/C++ code computing the values for those parameters that maximize the probability term. AutoBayes can solve many subproblems symbolically rather than having to rely on numeric approximation algorithms, thus yielding effective, efficient, and compact code. Statistical analysis is faster and more reliable, because effort can be focused on model development and validation rather than manual development of solution algorithms and code.

Schumann, Johann↗
NASA NTRS2011

Ten Years of Cloud Properties from MODIS: Global Statistics and Use in Climate Model Evaluation

The NASA Moderate Resolution Imaging Spectroradiometer (MODIS), launched onboard the Terra and Aqua spacecrafts, began Earth observations on February 24, 2000 and June 24,2002, respectively. Among the algorithms developed and applied to this sensor, a suite of cloud products includes cloud masking/detection, cloud-top properties (temperature, pressure), and optical properties (optical thickness, effective particle radius, water path, and thermodynamic phase). All cloud algorithms underwent numerous changes and enhancements between for the latest Collection 5 production version; this process continues with the current Collection 6 development. We will show example MODIS Collection 5 cloud climatologies derived from global spatial . and temporal aggregations provided in the archived gridded Level-3 MODIS atmosphere team product (product names MOD08 and MYD08 for MODIS Terra and Aqua, respectively). Data sets in this Level-3 product include scalar statistics as well as 1- and 2-D histograms of many cloud properties, allowing for higher order information and correlation studies. In addition to these statistics, we will show trends and statistical significance in annual and seasonal means for a variety of the MODIS cloud properties, as well as the time required for detection given assumed trends. To assist in climate model evaluation, we have developed a MODIS cloud simulator with an accompanying netCDF file containing subsetted monthly Level-3 statistical data sets that correspond to the simulator output. Correlations of cloud properties with ENSO offer the potential to evaluate model cloud sensitivity; initial results will be discussed.

Platnick, Steven E.↗
NASA NTRS2003

Updated Integrated Mission Program

Integrated Mission Program (IMP) is a computer program for simulating spacecraft missions around the Earth, Moon, Mars, and/or other large bodies. IMP solves the differential equations of motion by use of a Runge-Kutta numerical-integration algorithm. Users control missions through selection from a large menu of events and maneuvers. Mission profiles, time lines, propellant requirements, feasibility analyses, and perturbation analyses can be computed quickly and accurately. A prior version of IMP, written in FORTRAN 77, was reported in Program Simulates Spacecraft Missions (MFS-28606), NASA Tech Briefs, Vol. 17, No. 4 (April 1993), page 60. The present version, written in double-precision Lahey FORTRAN 90, incorporates a number of improvements over the prior version. Some of the improvements modernize the code to take advantage of today's greater central-processing-unit speeds. Other improvements render the code more modular; provide additional input, output, and debugging capabilities; and add to the variety of maneuvers, events, and means of propulsion that can be simulated. The IMP user manuals (of which there are now ten, each addressing a different aspect of the code and its use) have been updated accordingly.

Dauro, Vincent A., Sr.↗
NASA NTRS2012

A Survey of Symplectic and Collocation Integration Methods for Orbit Propagation

Demands on numerical integration algorithms for astrodynamics applications continue to increase. Common methods, like explicit Runge-Kutta, meet the orbit propagation needs of most scenarios, but more specialized scenarios require new techniques to meet both computational efficiency and accuracy needs. This paper provides an extensive survey on the application of symplectic and collocation methods to astrodynamics. Both of these methods benefit from relatively recent theoretical developments, which improve their applicability to artificial satellite orbit propagation. This paper also details their implementation, with several tests demonstrating their advantages and disadvantages.

integration↗
NASA NTRS2013

Automated Loads Analysis System (ATLAS)

ATLAS is a generalized solution that can be used for launch vehicles. ATLAS is used to produce modal transient analysis and quasi-static analysis results (i.e., accelerations, displacements, and forces) for the payload math models on a specific Shuttle Transport System (STS) flight using the shuttle math model and associated forcing functions. This innovation solves the problem of coupling of payload math models into a shuttle math model. It performs a transient loads analysis simulating liftoff, landing, and all flight events between liftoff and landing. ATLAS utilizes efficient and numerically stable algorithms available in MSC/NASTRAN.

Gardner, Stephen↗
NASA NTRS2014

In-situ Microwave Brightness Temperature Variability from Ground-based Radiometer Measurements at Dome C in Antarctica Induced by Wind-formed Features

Space-borne microwave radiometers are among the most useful tools to study snow and to collect information on the Antarctic climate. They have several advantages over other remote sensing techniques: high sensitivity to snow properties of interest (temperature, grain size, density), subdaily coverage in the polar regions, and their observations are independent of cloud conditions and solar illumination. Thus, microwave radiometers are widely used to retrieve information over snow-covered regions. For the Antarctic Plateau, many studies presenting retrieval algorithms or numerical simulations have assumed, explicitly or not, that the subpixel-scale heterogeneity is negligible and that the retrieved properties were representative of whole pixels. In this presentation, we investigate the spatial variations of brightness temperature over arange of a few kilometers in the Dome C area (Antarctic Plateau).

Antarctica↗
NASA NTRS2017

Near Real-Time Probabilistic Damage Diagnosis Using Surrogate Modeling and High Performance Computing

This work investigates novel approaches to probabilistic damage diagnosis that utilize surrogate modeling and high performance computing (HPC) to achieve substantial computational speedup. Motivated by Digital Twin, a structural health management (SHM) paradigm that integrates vehicle-specific characteristics with continual in-situ damage diagnosis and prognosis, the methods studied herein yield near real-time damage assessments that could enable monitoring of a vehicle's health while it is operating (i.e. online SHM). High-fidelity modeling and uncertainty quantification (UQ), both critical to Digital Twin, are incorporated using finite element method simulations and Bayesian inference, respectively. The crux of the proposed Bayesian diagnosis methods, however, is the reformulation of the numerical sampling algorithms (e.g. Markov chain Monte Carlo) used to generate the resulting probabilistic damage estimates. To this end, three distinct methods are demonstrated for rapid sampling that utilize surrogate modeling and exploit various degrees of parallelism for leveraging HPC. The accuracy and computational efficiency of the methods are compared on the problem of strain-based crack identification in thin plates. While each approach has inherent problem-specific strengths and weaknesses, all approaches are shown to provide accurate probabilistic damage diagnoses and several orders of magnitude computational speedup relative to a baseline Bayesian diagnosis implementation.

Warner, James E.↗
NASA NTRS2018

A Higher-Order Space-Time Finite-Element Method for Moving-Body and Fluid-Structure Interaction Problems

We present a high-order finite-element method for moving body and fluid/structure interaction problems. Our solution strategy is based on a space-time discontinuous Galerkin (DG) spectral-element discretization which extends to arbitrary order of accuracy. The space-time DG discretization is a natural choice for moving body and fluid-structure interaction problems as moving surfaces are incorporated simply by considering curved space-time elements whose space-time faces align with the moving body. We present a discontinuous-Galerkin in time discretization for six-degree of motion modeling of rigid bodies, and a continuous-Galerkin discretization for equations of linear elasticity to generate curved space-time meshes. Numerical results for several simple 2D test cases are presented in order to verify the implementation of the different models. Finally we present a preliminary dynamic simulation of a parachute.

High-Order↗
NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗
NASA NTRS2022

A Transient Hydrodynamic Model of Screen Channel Liquid Acquisition Devices for In-Space Cryogenic Propellant Management

Screen channel liquid acquisition devices (LADs) will play a crucial role in future deep space travel. It is essential that vapor-free delivery of propellants during tank-to-tank transfer is ensured to maximize yield from storage tanks and prevent potential combustion instabilities. The screen channel LAD utilizes a fine screen wire mesh that can separate phases in a low Bond number (i.e. microgravity) environment using surface tension forces. This study presents the development and verification of a new model for transient screen compliance, one of the influential factors for screen channel LAD design. Screen compliance is crucial during LAD channel outflow transients because the slight deflection of the screen can provide needed mass to satisfy rapid outflow demands and reduce the pressure difference across the screen. The model is successfully verified against CFD simulations. In addition, the characteristic speed for the governing screen compliance equations is derived which allows for numerical stability criteria to be established. As shown in this study, the transient maximum pressure difference across the screen can greatly exceed the steady state maximum pressure difference across the screen in many cases.

Hydrodynamics Simulations↗
NASA NTRS1986

Parametric effects of CFL number and artificial smoothing on numerical solutions using implicit approximate factorization algorithm

An implicit approximate factorization algorithm is employed to quantify the parametric effects of Courant number and artificial smoothing on numerical solutions of the unsteady 3-D Euler equations for a windmilling propeller (low speed) flow field. The results show that propeller global or performance chracteristics vary strongly with Courant number and artificial dissipation parameters, though the variation is such less severe at high Courant numbers. Candidate sets of Courant number and dissipation parameters could result in parameter-dependent solutions. Parameter-independent numerical solutions can be obtained if low values of the dissipation parameter-time step ratio are used in the computations. Furthermore, it is realized that too much artificial damping can degrade numerical stability. Finally, it is demonstrated that highly resolved meshes may, in some cases, delay convergence, thereby suggesting some optimum cell size for a given flow solution. It is suspected that improper boundary treatment may account for the cell size constraint.

Daso, E. O.↗
NASA NTRS1979

Automated design using numerical optimization

Numerical optimization concepts are described with limited technical detail. The purpose is to provide the nonspecialist with sufficient information to judge the applicability of these methods to his particular design problem. The concepts are first described in physical terms to give a basic understanding of the iterative procedure employed by these methods. Next, the typical engineering task is presented and converted to a form ammenable to solution by numerical optimization. Basic algorithms for solving this problem are identified. Numerous applications are referenced, emphasizing the structural design discipline. The state of the art allows for the routine solution of nonlinear design problems of approximately 20 independent variables subject to 100 or more constraints. In many applications, much larger design problems may be solved. Selected references are provided which describe the methods and applications in more detail.

Vanderplaats, G. N.↗