SEARCH · Engineering Papers

Results for “time step”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1981

A model of temporal variations in ozone density in the Martian atmosphere

Temporal variations of the ozone density profile in the Martian atmosphere at high latitudes are calculated for the course of a Martian year, taking into account seasonal and diurnal variations in temperature, water vapor and solar radiation. Calculations are based on a model including 35 neutral photochemical reactions, and vertical eddy diffusion using a time step of 12 min for the region from the surface to 240 km altitude. Results are found to be in better agreement with Mariner 9 observations of the time and magnitude of the seasonal maximum than previous model calculations. The diurnal variation is predicted to be small near the solstices, with the nighttime ozone density greater than the daytime and the magnitude of the difference dependent on season. Opposite temporal variations are predicted for ozone densities above and below about 25 km, and an ozone density maximum at 35-40 km is obtained. It is suggested that the effects of an aerosol layer may not be important in enhancing predicted ozone concentration, and may even decrease it.

Shimazaki, T.↗

NASA NTRS2002

Parametric Study of a YAV-8B Harrier in Ground Effect using Time-Dependent Navier-Stokes Computations

A process is described which enables the generation of 35 time-dependent viscous solutions for a YAV-8B Harrier in ground effect in one week. Overset grids are used to model the complex geometry of the Harrier aircraft and the interaction of its jets with the ground plane and low-speed ambient flow. The time required to complete this parametric study is drastically reduced through the use of process automation, modern computational platforms, and parallel computing. Moreover, a dual-time-stepping algorithm is described which improves solution robustness. Unsteady flow visualization and a frequency domain analysis are also used to identify and correlated key flow structures with the time variation of lift.

Pandya, Shishir↗

NASA NTRS1982

Flows over wings with leading-edge vortex separation

The unsteady cross flow analogy reduces the steady three dimensional separation flow problem into an unsteady two dimensional flow problem in which the section shape changes with time. The two dimensional VORSEP code is extended to the case of arbitrary body growth rates in order to generate the initial vortex structures for the three dimensional free vortex sheet (FVS) code. Automatic procedures to reduce the wing geometry definition to a set of cross flow plane sections corresponding to the locations of the time step solutions and to generate the effective source distribution on each cross flow section to represent the section normal growth across the following steps are incorporated in the VORSEP code. Also, the wake shedding model is improved by adopting a redistribution scheme which improves the stability of the free sheet development with time. The improved wake shedding model combined with the redistribution scheme alleviated the numerical instabilities associated with the vortex rollup.

Rao, B. M.↗

NASA NTRS1985

Efficient self-consistent viscous-inviscid solutions for unsteady transonic flow

An improved method is presented for coupling a boundary layer code with an unsteady inviscid transonic computer code in a quasi-steady fashion. At each fixed time step, the boundary layer and inviscid equations are successively solved until the process converges. An explicit coupling of the equations is described which greatly accelerates the convergence process. Computer times for converged viscous-inviscid solutions are about 1.8 times the comparable inviscid values. Comparison of the results obtained with experimental data on three airfoils are presented. These comparisons demonstrate that the explicitly coupled viscous-inviscid solutions can provide efficient predictions of pressure distributions and lift for unsteady two-dimensional transonic flow.

Howlett, J. T.↗

NASA NTRS1985

Efficient self-consistent viscous-inviscid solutions for unsteady transonic flow

Howlett, J. T.↗

NASA NTRS1997

A Multiblock Approach for Calculating Incompressible Fluid Flows on Unstructured Grids

A multiblock approach is presented for solving two-dimensional incompressible turbulent flows on unstructured grids. The artificial compressibility form of the governing equations is solved by a vertex-centered, finite-volume implicit scheme which uses a backward Euler time discretization. Point Gauss-Seidel relaxations are used to solve the linear system of equations at each time step. This work introduces a multiblock strategy to the solution procedure, which greatly improves the efficiency of the algorithm by significantly reducing the memory requirements while not increasing the CPU time. Results presented in this work shows that the current multiblock algorithm requires 70% less memory than the single block algorithm.

Sheng, Chunhua↗

NASA NTRS1991

Computational simulation of the creep-rupture process in filamentary composite materials

A computational simulation of the internal damage accumulation which causes the creep-rupture phenomenon in filamentary composite materials is developed. The creep-rupture process involves complex interactions between several damage mechanisms. A statistically-based computational simulation using a time-differencing approach is employed to model these progressive interactions. The finite element method is used to calculate the internal stresses. The fibers are modeled as a series of bar elements which are connected transversely by matrix elements. Flaws are distributed randomly throughout the elements in the model. Load is applied, and the properties of the individual elements are updated at the end of each time step as a function of the stress history. The simulation is continued until failure occurs. Several cases, with different initial flaw dispersions, are run to establish a statistical distribution of the time-to-failure. The calculations are performed on a supercomputer. The simulation results compare favorably with the results of creep-rupture experiments conducted at the Lawrence Livermore National Laboratory.

Slattery, Kerry T.↗

NASA NTRS1994

Calculations of Arctic ozone chemistry using objectively analyzed data in a 3-D CTM

A three-dimensional chemical transport model (CTM) (Kaminski, 1992) has been used to study the evolution of the Arctic ozone during the winter of 1992. The continuity equation has been solved using a spectral method with Rhomboidal 15 (R15) truncation and leap-frog time stepping. Six-hourly meteorological fields from the Canadian Meteorological Center global objective analysis routines run at T79 were degraded to the model resolution. In addition, they were interpolated to the model time grid and were used to drive the model from the surface to 10 mb. In the model, processing of Cl(x) occurred over Arctic latitudes but some of the initial products were still present by mid-January. Also, the large amounts of ClO formed in the model in early January were converted to ClNO3. The results suggest that the model resolution may be insufficient to resolve the details of the Arctic transport during this time period. In particular, the wind field does not move the ClO(x) 'cloud' to the south over Europe as seen in the MLS measurements.

Kaminski, J. W.↗

NASA NTRS2000

A Regional-Scale Assessment of Satellite Derived Precipitable Water Vapor Across The Amazon Basin

Atmospheric water vapor is widely recognized as a key climate variable, linking an assortment of poorly understood and complex processes. It is a major element of the hydrological cycle and provides a mechanism for energy exchange among many of the Earth system components. Reducing uncertainty in our current knowledge of water vapor and its role in the climate system requires accurate measurement, improved modeling techniques, and long-term prediction. Satellites have the potential to satisfy these criteria, as well as provide high resolution measurements that are not available from conventional sources. The focus of this paper is to examine the temporal and mesoscale variations of satellite derived precipitable water vapor (PW) across the Amazon Basin. This region is pivotal in the functioning of the global climate system through its abundant release of latent heat associated with heavy precipitation events. In addition, anthropogenic deforestation and biomass burning activities in recent decades are altering the conditions of the atmosphere, especially in the planetary boundary layer. A physical split-window (PSW) algorithm estimates PW using images from the GOES satellites along with the NCEP/NCAR Reanalysis data that provides the first guess information. Retrievals are made at a three-hourly time step during daylight hours in the Amazon Basin and surrounding areas for the months of June and October in 1988 (dry year) and 1995 (wet year). Spatially continuous fields are generated 5 times daily at 12Z, 15Z, 18Z, 21Z, and 00Z. These fields are then averaged to create monthly and 3 hourly monthly grids. Overall, the PSW estimates PW reasonable well in the Amazon with MAE ranging from 3.0 - 9.0 mm and MAE/observed mean around 20% in comparison to radiosonde observations. The distribution of PW generally mimics that of precipitation. Maximum values (42 - 52 mm) are located in the Northwest whereas minimum values (18 - 27 mm) are found along Brazil's East coast. Aside from the East coast, PW has a stronger north-south gradient than that of rainfall. As for the temporal variation of PW, June (1988 and 1995) experiences a peak about 1400 local time, corresponding to a maximum in air temperature. In contrast, October (1988 and 1995) experiences a maximum early in the day, 1100 local time, with a gradual decrease toward nighttime.

DeLiberty, Tracy↗

NASA NTRS2022

Unique Observational Constraints on the Seasonal and Longitudinal Variability of the Earth’s Planetary Albedo and Cloud Distribution Inferred From EPIC Measurements

Thorough comparison to observations is key to developing a credible climate model forecasting capability. Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar and emitted thermal radiation provide a unique observational perspective that permits a more reliable model/data comparison than is possible with the otherwise available satellite data. The uniqueness is in the DSCOVR satellite’s viewing geometry, which enables continuous viewing of the Earth’s sunlit hemisphere from its Lissajous orbit around the Lagrangian L1 point. The key instrument is the Earth Polychromatic Imaging Camera (EPIC), which views the Earth’s sunlit hemisphere with 1024-by-1024-pixel imagery in 10 narrow spectral bands from 317 to 780 nm, acquiring up to 22 high spatial resolution images per day. The additional feature is that the frequency of EPIC image acquisition is nearly identical to that of the climate GCM data generation scheme where climate data for the entire globe are ‘instantaneously’ calculated at 1-h radiation time-step intervals. Implementation of the SHS (Sunlit Hemisphere Sampling) EPIC-view geometry for the in-line GCM output data sampling establishes a precise self-consistency in the space-time data sampling between EPIC observational and GCM output data generation and sampling. The remaining problem is that the GCM generated data are radiative fluxes, while the EPIC measurements are backscatter-dependent radiances. Radiance to flux conversion is a complex problem with no simple way to convert GCM radiative fluxes into spectral radiances. The more expedient approach is to convert the EPIC spectral radiances into broadband radiances by MODIS/CERES-based regression relationships and then into solar radiative fluxes using the CERES angular distribution models. Averaging over the sunlit hemisphere suppresses the meteorological weather noise, but preserves the intra-seasonal larger scale variability. Longitudinal slicing by the Earth’s rotation permits a self-consistent model/data comparison of the longitudinal model/data differences in the variability of the reflected solar radiation. Ancillary EPIC Composite data provide additional cloud property information for climate model diagnostics. Comparison of EPIC-derived seasonal and longitudinal variability of the Earth’s planetary albedo with the GISS ModelE2 results shows systematic overestimate of cloud reflectivity over the Pacific Ocean with corresponding underestimates over continental land areas.

DSCOVR↗

NASA NTRS1994

On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow

In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Sweby, P. K.↗

NASA NTRS1981

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗

NASA NTRS2018

Progress Towards Integrating the Finite-Volume Cubed-Sphere (FV3) Dynamical Core Tangent Linear and Adjoint Models into JEDI

The Joint Effort for Data assimilation Integration (JEDI) -- led by the Joint Center for Satellite Data Assimilation (JCSDA) -- is an inter-organizational endeavor to develop a common framework for performing data assimilation. This extensive framework will ultimately provide solvers, observation operators, interpolation and model interfaces using object oriented modeling. Two partners involved in JEDI use or plan to use the Finite Volume Cubed-Sphere (FV3) dynamical core to produce weather forecasts; these are NASA's Global Modeling and Assimilation Office and NOAA's National Center for Environment Prediction. In this work we present an update on ongoing efforts to integrate the FV3 tangent linear and adjoint models into the prototype JEDI framework. We setup and run a simple cycled data assimilation experiment using 4DVAR on the cubed sphere grid and with the FV3 tangent linear and adjoint models. Development of the observation operators for JEDI is separately underway. Instead of using real observations a simplified set of simulated observations will be used. We discuss the steps required to bring the FV3 linearized model into the object oriented framework and consider what would be the computational requirements of running this configuration for an operational system. FV3 uses a small time-step to ensure that small scales are well resolved, however this presents design challenges when running 4DVAR with the adjoint. An approach to storing the FV3 model trajectory has been developed that maintains the flexibility of using automatic differentiation. We discuss how this approach is incorporated into the framework. Other important uses of adjoint models include computing observation impacts and singular vectors, we consider how these tools can be included in JEDI.

Holdaway, Daniel↗

NASA NTRS1995

Convergence acceleration of the Proteus computer code with multigrid methods

This report presents the results of a study to implement convergence acceleration techniques based on the multigrid concept in the two-dimensional and three-dimensional versions of the Proteus computer code. The first section presents a review of the relevant literature on the implementation of the multigrid methods in computer codes for compressible flow analysis. The next two sections present detailed stability analysis of numerical schemes for solving the Euler and Navier-Stokes equations, based on conventional von Neumann analysis and the bi-grid analysis, respectively. The next section presents details of the computational method used in the Proteus computer code. Finally, the multigrid implementation and applications to several two-dimensional and three-dimensional test problems are presented. The results of the present study show that the multigrid method always leads to a reduction in the number of iterations (or time steps) required for convergence. However, there is an overhead associated with the use of multigrid acceleration. The overhead is higher in 2-D problems than in 3-D problems, thus overall multigrid savings in CPU time are in general better in the latter. Savings of about 40-50 percent are typical in 3-D problems, but they are about 20-30 percent in large 2-D problems. The present multigrid method is applicable to steady-state problems and is therefore ineffective in problems with inherently unstable solutions.

Demuren, A. O.↗

NASA NTRS1995

Crack growth behavior of AISI-4340 steel during environmental exposure

AISI-4340 is observed to undergo stress corrosion cracking when subjected to a constant load during exposure to a 3.5% NaCl solution. Crack initiation, nucleation, and growth has been monitored as a function of time. Stepped regions consisting of fast and slow crack growth periods are shown to correspond to microstructural changes observed in the fracture surface of the steel. These regions of fast and slow crack rate variations with time show that the crack growth rates do not increase continuously with an increase in the stress intensity.

Giannuzzi, Lucille A.↗

NASA NTRS1983

Improvement of mixed time implicit-explicit algorithms for thermal analysis of structures

Computer implementation aspects and numerical evaluation of the recently introduced mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for a linear quadrilateral element is given herein for the methods introduced by Liu and co-workers. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗

NASA NTRS1984

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, H.↗

NASA NTRS1986

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, Hillel↗