Engineering Papers⌕ Search

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.

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

A Time Domain Analysis of Gust-Cascade Interaction Noise

The gust response of a 2 D cascade is studied by solving the full nonlinear Euler equations employing higher order accurate spatial differencing and time stepping techniques. The solutions exhibit the exponential decay of the two circumferential mode orders of the cutoff blade passing frequency (BPF) tone and propagation of one circumferential mode order at 2BPF, as would be expected for the flow configuration considered. Two frequency excitations indicate that the interaction between the frequencies and the self interaction contribute to the amplitude of the propagating mode.

Nallasamy, M.↗
NASA NTRS2003

NAS Parallel Benchmarks, Multi-Zone Versions

We describe an extension of the NAS Parallel Benchmarks (NPB) suite that involves solving the application benchmarks LU, BT and SP on collections of loosely coupled discretization meshes. The solutions on the meshes are updated independently, but after each time step they exchange boundary value information. This strategy, which is common among structured-mesh production flow solver codes in use at NASA Ames and elsewhere, provides relatively easily exploitable coarse-grain parallelism between meshes. Since the individual application benchmarks also allow fine-grain parallelism themselves, this NPB extension, named NPB Multi-Zone (NPB-MZ), is a good candidate for testing hybrid and multi-level parallelization tools and strategies.

vanderWijngaart, Rob F.↗
NASA NTRS2003

Research in Computational Astrobiology

We report on several projects in the field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. Research projects included modifying existing computer simulation codes to use efficient, multiple time step algorithms, statistical methods for analysis of astrophysical data via optimal partitioning methods, electronic structure calculations on water-nuclei acid complexes, incorporation of structural information into genomic sequence analysis methods and calculations of shock-induced formation of polycylic aromatic hydrocarbon compounds.

Chaban, Galina↗
NASA NTRS2003

Performance Characteristics of the Multi-Zone NAS Parallel Benchmarks

We describe a new suite of computational benchmarks that models applications featuring multiple levels of parallelism. Such parallelism is often available in realistic flow computations on systems of grids, but had not previously been captured in bench-marks. The new suite, named NPB Multi-Zone, is extended from the NAS Parallel Benchmarks suite, and involves solving the application benchmarks LU, BT and SP on collections of loosely coupled discretization meshes. The solutions on the meshes are updated independently, but after each time step they exchange boundary value information. This strategy provides relatively easily exploitable coarse-grain parallelism between meshes. Three reference implementations are available: one serial, one hybrid using the Message Passing Interface (MPI) and OpenMP, and another hybrid using a shared memory multi-level programming model (SMP+OpenMP). We examine the effectiveness of hybrid parallelization paradigms in these implementations on three different parallel computers. We also use an empirical formula to investigate the performance characteristics of the multi-zone benchmarks.

Jin, Haoqiang↗
NASA NTRS2003

Local and Remote Sources of Precipitation during the 1993 Midwestern United States Flood

During June and July of 1993, devastating and persistent rainfall caused substantial damage in the Midwest United States primarily focused on the Mississippi River Basin. Above noma1 springtime rains primed the region for flooding. The resulting wet soil was thought to be a major source of water for the heavy summer precipitation. On the other hand, northward moisture transport was extremely intense during the summer. In this study, we have implemented passive tracers to diagnose the geographic source of water for precipitation, called Water Vapor Tracers (WVT), into the NASA Data Assimilation Office (DAO) Data Assimilation System (DAS). The WVT diagnostics can produce detailed budgets of the geographical sources of water that precipitates in any region, and they consider physical and dynamical tendencies at the model time step. The 1993 summer has been reanalyzed with the new WVT diagnostics to better understand the local and remote source of water for precipitation during the flood. Because the DAS also uses observations of water in the atmospheric hydrologic budget, we will also evaluate the impact of the observations on the local and remote sources of water in the DAS.

Bosilovich, Michael↗
NASA NTRS2003

Numerical Simulation of High Drag Reduction in a Turbulent Channel Flow with Polymer Additives

The addition of small amounts of long chain polymer molecules to wall-bounded flows can lead to dramatic drag reduction. Although this phenomenon has been known for about fifty years, the action of the polymers and its effect on turbulent structures are still unclear. Detailed experiments have characterized two distinct regimes (Warholic et al. 1999), which are referred to as low drag reduction (LDR) and high drag reduction (HDR). The first regime exhibits similar statistical trends as Newtonian flow: the log-law region of the mean velocity profile remains parallel to that of the Newtonian ow but its lower bound moves away from the wall and the upward shift of the log-region is a function of drag reduction, DR. Although streamwise fluctuations are increased and transverse ones are reduced, the shape of the rms velocity profiles is not qualitatively modified. At higher drag reductions, of the order of 40-50%, the ow enters the HDR regime for which the slope of the log-law is dramatically augmented and the Reynolds shear stress is small (Warholic et al. 1999; Ptasinski et al. 2001). The drag reduction is eventually bounded by a maximum drag reduction (MDR) (Virk & Mickley 1970) which is a function of the Reynolds number. While several experiments report mean velocity profiles very close to the empirical profile of Virk & Mickley (1970) for MDR conditions, the observations regarding the structure of turbulence can differ significantly. For instance, Warholic et al. (1999) measured a near-zero Reynolds shear stress, whereas a recent experiment (Ptasinski et al. 2001) shows evidence of non-negligible Reynolds stress in their MDR flow. To the knowledge of the authors, only the LDR regime has been documented in numerical simulations (Sureshkumar et al. 1997; Dimitropoulos et al. 1998; Min et al. 2001; Dubief & Lele 2001; Sibilla & Baron 2002). This paper discusses the simulation of polymer drag reduced channel ow at HDR using the FENE-P (Finite Elastic non-linear extensibility-Peterlin) model which was used for the first LDR simulation by Sureshkumar et al. (1997). Flow and polymer parameters are close to realistic polymer drag reducing conditions. High drag reductions are achieved by using finite differences and a robust time stepping technique. A minimal channel flow is also used as a numerical experiment to investigate the effect of the outer region turbulent structures on the overall drag at HDR. The drag reducing action of the model is finally studied through the structure of energy transfers from the polymers to the velocity components. This investigation sheds some light on the details of polymer drag reduction.

Dubief, Yves↗
NASA NTRS2003

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗
NASA NTRS2004

Reliability of Complex Nonlinear Numerical Simulations

This work describes some of the procedure to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗
NASA NTRS2004

Discontinuous Galerkin Finite Element Method for Parabolic Problems

In this paper, we develop a time and its corresponding spatial discretization scheme, based upon the assumption of a certain weak singularity of parallel ut(t) parallel Lz(omega) = parallel ut parallel2, for the discontinuous Galerkin finite element method for one-dimensional parabolic problems. Optimal convergence rates in both time and spatial variables are obtained. A discussion of automatic time-step control method is also included.

Kaneko, Hideaki↗
NASA NTRS2004

Explanation Constraint Programming for Model-based Diagnosis of Engineered Systems

We can expect to see an increase in the deployment of unmanned air and land vehicles for autonomous exploration of space. In order to maintain autonomous control of such systems, it is essential to track the current state of the system. When the system includes safety-critical components, failures or faults in the system must be diagnosed as quickly as possible, and their effects compensated for so that control and safety are maintained under a variety of fault conditions. The Livingstone fault diagnosis and recovery kernel and its temporal extension L2 are examples of model-based reasoning engines for health management. Livingstone has been shown to be effective, it is in demand, and it is being further developed. It was part of the successful Remote Agent demonstration on Deep Space One in 1999. It has been and is being utilized by several projects involving groups from various NASA centers, including the In Situ Propellant Production (ISPP) simulation at Kennedy Space Center, the X-34 and X-37 experimental reusable launch vehicle missions, Techsat-21, and advanced life support projects. Model-based and consistency-based diagnostic systems like Livingstone work only with discrete and finite domain models. When quantitative and continuous behaviors are involved, these are abstracted to discrete form using some mapping. This mapping from the quantitative domain to the qualitative domain is sometimes very involved and requires the design of highly sophisticated and complex monitors. We propose a diagnostic methodology that deals directly with quantitative models and behaviors, thereby mitigating the need for these sophisticated mappings. Our work brings together ideas from model-based diagnosis systems like Livingstone and concurrent constraint programming concepts. The system uses explanations derived from the propagation of quantitative constraints to generate conflicts. Fast conflict generation algorithms are used to generate and maintain multiple candidates whose consistency can be tracked across multiple time steps.

Narasimhan, Sriram↗
NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗
NASA NTRS2004

Potential Predictability of U.S. Summer Climate with "Perfect" Soil Moisture

The potential predictability of surface-air temperature and precipitation over the United States continent was assessed for a GCM forced by observed sea surface temperatures and an estimate of observed ground soil moisture contents. The latter was obtained by substituting the GCM simulated precipitation, which is used to drive the GCM's land-surface component, with observed pentad-mean precipitation at each time step of the model's integration. With this substitution, the simulated soil moisture correlates well with an independent estimate of observed soil moisture in all seasons over the entire US continent. Significant enhancements on the predictability of surface-air temperature and precipitation were found in boreal late spring and summer over the US continent. Anomalous pattern correlations of precipitation and surface-air temperature over the US continent in the June-July-August season averaged for the 1979-2000 period increased from 0.01 and 0.06 for the GCM simulations without precipitation substitution to 0.23 and 0.3 1, respectively, for the simulations with precipitation substitution. Results provide an estimate for the limits of potential predictability if soil moisture variability is to be perfectly predicted. However, this estimate may be model dependent, and needs to be substantiated by other modeling groups.

Yang, Fanglin↗
NASA NTRS2004

Offline GCSS Intercomparison of Cloud-Radiation Interaction and Surface Fluxes

Simulations of deep tropical clouds by both cloud-resolving models (CRMs) and single-column models (SCMs) in the GEWEX Cloud System Study (GCSS) Working Group 4 (WG4; Precipitating Convective Cloud Systems), Case 2 (19-27 December 1992, TOGA-COARE IFA) have produced large differences in the mean heating and moistening rates (-1 to -5 K and -2 to 2 grams per kilogram respectively). Since the large-scale advective temperature and moisture "forcing" are prescribed for this case, a closer examination of two of the remaining external types of "forcing", namely radiative heating and air/sea hear and moisture transfer, are warranted. This paper examines the current radiation and surface flux of parameterizations used in the cloud models participating in the GCSS WG4, be executing the models "offline" for one time step (12 s) for a prescribed atmospheric state, then examining the surface and radiation fluxes from each model. The dynamic, thermodynamic, and microphysical fluids are provided by the GCE-derived model output for Case 2 during a period of very active deep convection (westerly wind burst). The surface and radiation fluxes produced from the models are then divided into prescribed convective, stratiform, and clear regions in order to examine the role that clouds play in the flux parameterizations. The results suggest that the differences between the models are attributed more to the surface flux parameterizations than the radiation schemes.

Tao, W.-K.↗
NASA NTRS1997

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗
NASA NTRS2004

Recent Enhancements to USM3D Unstructured Flow Solver for Unsteady Flows

The NASA USM3D unstructured flow solver is undergoing extensions to address dynamic flow problems in support of NASA and NAVAIR efforts to study the applicability of Computational Fluid Dynamics tools for the prediction of aircraft stability and control characteristics. The initial extensions reported herein include two second-order time stepping schemes, Detached-Eddy Simulation, and grid motion. This paper reports the initial code verification and validation assessment of the dynamic flow capabilities of USM3D. The cases considered are the classic inviscid shock-tube problem, low Reynolds number wake shedding from a NACA 0012 airfoil, high Reynolds number DES-based wake shedding from a 4-to-1 length-to-diameter cylinder, and forced pitch oscillation of a NACA 0012 airfoil with inviscid and turbulent flow.

Pandya, Mohagna J.↗
NASA NTRS1998

Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations

A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

Frink, N. T.↗
NASA NTRS1998

On Some Numerical Dissipation Schemes

Several schemes for introducing an artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar dissipation and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical solutions are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme. The coarse-grid accuracy for the original CUSP scheme is improved by modifying the limiter function used with the scheme, giving comparable accuracy to that obtained with the MATD scheme. The modifications reduce the background dissipation and provide control over the regions where the scheme can become first order.

Swanson, R. C.↗
NASA NTRS1997

Implicit/Multigrid Algorithms for Incompressible Turbulent Flows on Unstructured Grids

An implicit code for computing inviscid and viscous incompressible flows on unstructured grids is described. The foundation of the code is a backward Euler time discretization for which the linear system is approximately solved at each time step with either a point implicit method or a preconditioned Generalized Minimal Residual (GMRES) technique. For the GMRES calculations, several techniques are investigated for forming the matrix-vector product. Convergence acceleration is achieved through a multigrid scheme that uses non-nested coarse grids that are generated using a technique described in the present paper. Convergence characteristics are investigated and results are compared with an exact solution for the inviscid flow over a four-element airfoil. Viscous results, which are compared with experimental data, include the turbulent flow over a NACA 4412 airfoil, a three-element airfoil for which Mach number effects are investigated, and three-dimensional flow over a wing with a partial-span flap.

Anderson, W. Kyle↗