Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical methods and modeling”

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 55 records · Page 3

An infrared radiation routine for use in numerical atmospheric models

Previous methods for calculating radiative fluxes due to water vapor and CO2 absorption bands are extended to take into consideration the entire water vapor and CO2 bands, including e-type absorption in the window region and the overlapping of different absorptions. Cooling rate profiles in the water vapor bands for a tropical atmosphere were computed by a detailed line-by-line method and by a far-wing approximation method, and the error of both methods is less than 0.2 C/day. Cooling rate profiles in the 15 micron band including overlapping of CO2 absorption with water vapor were calculated by a method in which flux transmittance is computed by means of a linear expansion and the multiplication rule, and maximum errors of 0.3 C/day were found in comparison with the exact line-by-line method.

Chow, M.-D.↗

The present state and future directions of PDF methods

The objectives of the workshop are presented in viewgraph format, as is this entire article. The objectives are to discuss the present status and the future direction of various levels of engineering turbulence modeling related to Computational Fluid Dynamics (CFD) computations for propulsion; to assure that combustion is an essential part of propulsion; and to discuss Probability Density Function (PDF) methods for turbulent combustion. Essential to the integration of turbulent combustion models is the development of turbulent model, chemical kinetics, and numerical method. Some turbulent combustion models typically used in industry are the k-epsilon turbulent model, the equilibrium/mixing limited combustion, and the finite volume codes.

Pope, S. B.↗

Ship viscous flow: A report on the 1990 SSPA-IIHR Workshop

To assess the state of the art in ship viscous flow computation a Workshop was organized in 1990 by three organizations: SSPA Maritime Consulting AB, Chalmers University of Technology, and the Iowa Institute of Hydraulic Research. Two test cases were specified by the organizers and sent out to all interested research groups, which were asked to submit results in a prescribed format. In September 1990 a meeting was held at Chalmers University of Technology. All results had then been collected and presented in a common format, and the theories based on responses to a questionnaire sent out earlier. During the meeting, each research group was first given the opportunity to briefly introduce their method and results. Thereafter, a considerable time was spent on general discussions on the performance of the different methods considering the differences in the underlying theories. Specific items that were addressed were grid generation, governing equations, boundary conditions, turbulence modelling, and numerical methods. Practical aspects of the results, for instance from the point of view of propeller design, were also discussed. The Workshop Proceedings contain a description of the participating methods, and the results of both test cases. In the present paper, a summary of the Workshop and its results is presented.

Patel, Virendra C.↗

Toward Scientific Numerical Modeling

Ultimately, scientific numerical models need quantified output uncertainties so that modeling can evolve to better match reality. Documenting model input uncertainties and verifying that numerical models are translated into code correctly, however, are necessary first steps toward that goal. Without known input parameter uncertainties, model sensitivities are all one can determine, and without code verification, output uncertainties are simply not reliable. To address these two shortcomings, two proposals are offered: (1) an unobtrusive mechanism to document input parameter uncertainties in situ and (2) an adaptation of the Scientific Method to numerical model development and deployment. Because these two steps require changes in the computational simulation community to bear fruit, they are presented in terms of the Beckhard-Harris-Gleicher change model.

Kleb, Bil↗

Numerical modeling of spray combustion with an advanced VOF method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.

Chen, Yen-Sen↗

Modeling of growth, evaporation and sedimentation effects on transmission of visible and IR laser beams in artificial fogs

The dense polydisperse aerosol particles in a quiet chamber may spontaneously go through different microphysical processes including gravitational sedimentation, thermal coagulation, and growth or evaporation. In an earlier paper, we presented the results of a parametric study of the combined and separate effects of thermal coagulation and sedimentation on the time dependence of extinction of four visible and IR laser beams traversing an aerosol medium. As a continuation of this series of studies, the separate and combined effects of growth or evaporation and gravitational sedimentation on the time dependence of extinction of the same four visible and IR laser beams traversing in artificial fogs will be reported in this paper. The method of numerically modeling the change of water droplet size distribution with time due to growth/evaporation and the cutoff of larger aerosols due to gravitational sedimentation is described in detail. Factors governing the relative importance of these two processes are discussed. Results of this study show that the relative humidity or ambient temperature is a crucial parameter in determining the optical depth of the water droplet and aerosol media undergoing microphysical processes.

Yue, G. K.↗

Primitive numerical simulation of circular Couette flow - Carrousel wind tunnel nonturbulent solutions

The azimuthal-invariant, three-dimensional cylindrical, incompressible Navier-Stokes equations are solved to steady state for a finite-length, physically realistic model. The numerical method relies on an alternating-direction implicit scheme that is formally second-order accurate in space and first-order accurate in time. The equations are linearized and uncoupled by evaluating variable coefficients at the previous time iteration. Wall grid clustering is provided by a Roberts transformation in radial and axial directions. A vorticity-velocity formulation is found to be preferable to a vorticity-streamfunction approach. Subject to no-slip, Dirichlet boundary conditions, except for the inner cylinder rotation velocity (impulsive start-up) and zero-flow initial conditions, nonturbulent solutions are obtained for sub- and supercritical Reynolds numbers of 100 to 400 for a finite geometry where R(outer)/R(inner) = 1.5, H/R(inner) = 0.73, and H/Delta-R = 1.5. An axially-stretched model solution is shown to asymptotically approach the one-dimensional analytic Couette solution at the cylinder midheight. Flowfield change from laminar to Taylor-vortex flow is discussed as a function of Reynolds number. Three-dimensional velocities, vorticity, and streamfunction are presented via two-dimensional graphs and three-dimensional surface and contour plots.

Hasiuk, Jan↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

Improved numerical methods for turbulent viscous flows aerothermal modeling program, phase 2

The details of a study to develop accurate and efficient numerical schemes to predict complex flows are described. In this program, several discretization schemes were evaluated using simple test cases. This assessment led to the selection of three schemes for an in-depth evaluation based on two-dimensional flows. The scheme with the superior overall performance was incorporated in a computer program for three-dimensional flows. To improve the computational efficiency, the selected discretization scheme was combined with a direct solution approach in which the fluid flow equations are solved simultaneously rather than sequentially.

Karki, K. C.↗

Computation of viscous transonic flow about a lifting airfoil

The Reynolds averaged Navier-Stokes equations are solved numerically for the viscous transonic flow about a stationary NACA 64A010 airfoil in free air. This paper presents descriptions of the numerical method, turbulence models employed, and boundary conditions appropriate to simulation of free-air flight. Computed results are presented for the airfoil at a free-stream Mach number of 0.8, angles of attack of 0 and 2 deg, and a Reynolds number based on a chord of 4 x 10 to the 6th. For the lifting case, unsteady periodic motion was calculated along the aft portion of the airfoil and in its wake. Recent experimental results obtained by Johnson indicate periodicity aft of the shock closely approximates the computed frequency, but the amplitude of the disturbances was significantly less than the calculated amplitude.

Walitt, L.↗

A view toward future fluid dynamics computing

Advances in computational fluid dynamics are paced by simulation methodology and computer resources. Examples of three-dimensional fluid dynamic simulations are presented to illustrate recent developments in equation modeling and numerical methods and to point out the need for increased computer power. Electronic technology dictates that to fill this need, computers will be based on parallel processing principles. The identification of parallelism in three dimensions is illustrated by examining an implicit, approximate-factorization approach to the Navier-Stokes equations. Finally, two computer concepts aimed at satisfying the demands of the three-dimensional Reynolds averaged Navier-Stokes simulations are discussed.

Bailey, F. R.↗

A numerical study of transient, thermally-conductive solar wind

A numerical analysis of transient solar wind starting at the solar surface and arriving at 1 AU is performed by an implicit numerical method. The model hydrodynamic equations include thermal conduction terms for both steady and unsteady simulations. Simulation results show significant influence of thermal conduction on both steady and time-dependent solar wind. Higher thermal conduction results in higher solar wind speed, higher temperature, but lower plasma density at 1 AU. Higher base temperature at the solar surface gives lower plasma speed, lower temperature, but higher density at 1 AU. Higher base density, on the other hand, gives lower velocity, lower temperature, but higher density at 1 AU.

Han, S. M.↗

Dynamics of interacting galaxies

The dynamics of interacting galaxies as observed in the present epoch is reviewed. Topics discussed include numerical methods for modeling galactic interactions, the signatures of interacting galaxies due to tidal forces, and events that add mass to a galaxy. The review also covers major mergers between systems of comparable mass, forms of activity triggered or induced by galactic interactions, and galaxies' return to normality and related cosmological issues. Finally, some questions that yet have to be answered are examined.

Barnes, Joshua E.↗

NASA Tech Briefs, Februrary 2013

Topics covered include: Measurements of Ultra-Stable Oscillator (USO) Allan Deviations in Space; Gaseous Nitrogen Orifice Mass Flow Calculator; Validation of Proposed Metrics for Two-Body Abrasion Scratch Test Analysis Standards; Rover Low Gain Antenna Qualification for Deep Space Thermal Environments; Automated, Ultra-Sterile Solid Sample Handling and Analysis on a Chip; Measuring and Estimating Normalized Contrast in Infrared Flash Thermography; Spectrally and Radiometrically Stable, Wideband, Onboard Calibration Source; High-Reliability Waveguide Vacuum/Pressure Window; Methods of Fabricating Scintillators With Radioisotopes for Beta Battery Applications; Magnetic Shield for Adiabatic Demagnetization Refrigerators (ADR); CMOS-Compatible SOI MESFETS for Radiation-Hardened DC-to-DC Converters; Silicon Heat Pipe Array; Adaptive Phase Delay Generator; High-Temperature, Lightweight, Self-Healing Ceramic Composites for Aircraft Engine Applications; Treatment to Control Adhesion of Silicone-Based Elastomers; High-Temperature Adhesives for Thermally Stable Aero-Assist Technologies; Rockballer Sample Acquisition Tool; Rock Gripper for Sampling, Mobility, Anchoring, and Manipulation; Advanced Magnetic Materials Methods and Numerical Models for Fluidization in Microgravity and Hypogravity; Data Transfer for Multiple Sensor Networks Over a Broad Temperature Range; Using Combustion Synthesis to Reinforce Berms and Other Regolith Structures; Visible-Infrared Hyperspectral Image Projector; Three-Axis Attitude Estimation With a High-Bandwidth Angular Rate Sensor Change_Detection.m; AGATE: Adversarial Game Analysis for Tactical Evaluation; Ionospheric Simulation System for Satellite Observations and Global Assimilative; Modeling Experiments (ISOGAME); An Extensible, User- Modifiable Framework for Planning Activities; Mission Operations Center (MOC) - Precipitation Processing System (PPS) Interface Software System (MPISS); Automated 3D Damaged Cavity Model Builder for Lower Surface Acreage Tile on Orbiter; Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays; RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers; Powered Descent Guidance with General Thrust-Pointing Constraints; X-Ray Detection and Processing Models for Spacecraft Navigation and Timing; and Extreme Ionizing-Radiation-Resistant Bacterium

Source record↗

Infrared Dielectric Properties of Low-stress Silicon Nitride

Silicon nitride thin films play an important role in the realization of sensors, filters, and high-performance circuits. Estimates of the dielectric function in the far- and mid-IR regime are derived from the observed transmittance spectra for a commonly employed low-stress silicon nitride formulation. The experimental, modeling, and numerical methods used to extract the dielectric parameters with an accuracy of approximately 4% are presented.

dielectric parameters↗

Infrared Dielectric Properties of Low-Stress Silicon Oxide

Silicon oxide thin films play an important role in the realization of optical coatings and high-performance electrical circuits. Estimates of the dielectric function in the far- and mid-infrared regime are derived from the observed transmittance spectrum for a commonly employed low-stress silicon oxide formulation. The experimental, modeling, and numerical methods used to extract the dielectric function are presented.

Silicon oxide↗

A Numerical Method for Calculating Stellar Occultation Light Curves from an Arbitrary Atmospheric Model

We present a method for speeding up numerical calculations of a light curve for a stellar occultation by a planetary atmosphere with an arbitrary atmospheric model that has spherical symmetry. This improved speed makes least-squares fitting for model parameters practical. Our method takes as input several sets of values for the first two radial derivatives of the refractivity at different values of model parameters, and interpolates to obtain the light curve at intermediate values of one or more model parameters. It was developed for small occulting bodies such as Pluto and Triton, but is applicable to planets of all sizes. We also present the results of a series of tests showing that our method calculates light curves that are correct to an accuracy of 10(exp -4) of the unocculted stellar flux. The test benchmarks are (i) an atmosphere with a l/r dependence of temperature, which yields an analytic solution for the light curve, (ii) an atmosphere that produces an exponential refraction angle, and (iii) a small-planet isothermal model. With our method, least-squares fits to noiseless data also converge to values of parameters with fractional errors of no more than 10(exp -4), with the largest errors occurring in small planets. These errors are well below the precision of the best stellar occultation data available. Fits to noisy data had formal errors consistent with the level of synthetic noise added to the light curve. We conclude: (i) one should interpolate refractivity derivatives and then form light curves from the interpolated values, rather than interpolating the light curves themselves; (ii) for the most accuracy, one must specify the atmospheric model for radii many scale heights above half light; and (iii) for atmospheres with smoothly varying refractivity with altitude, light curves can be sampled as coarsely as two points per scale height.

Chamberlain, D. M.↗