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 325 records · Page 18

High-Latitude Filtering in a Global Grid-Point Model Using Model Normal Modes

The aim of high-latitude filtering in the vicinity of the poles is to avoid the excessively short time steps imposed on an explicit time-differencing scheme by linear stability due to fast moving inertia-gravity waves near the poles. The model normal mode expansion toward the problem of high-latitude filtering in a global shallow water model using the same philosophy as that used by Daley for the problem for large timesteps in P.E. models with explicit time integration schemes was applied.

Takacs, L. L.↗

Calculation of viscous effects on transonic flow for oscillating airfoils and comparisons with experiment

A method is described for calculating unsteady transonic flow with viscous interaction by coupling a steady integral boundary-layer code with an unsteady, transonic, inviscid small-disturbance computer code in a quasi-steady fashion. Explicit coupling of the equations together with viscous -inviscid iterations at each time step yield converged solutions with computer times about double those required to obtain inviscid solutions. The accuracy and range of applicability of the method are investigated by applying it to four AGARD standard airfoils. The first-harmonic components of both the unsteady pressure distributions and the lift and moment coefficients have been calculated. Comparisons with inviscid calcualtions and experimental data are presented. The results demonstrate that accurate solutions for transonic flows with viscous effects can be obtained for flows involving moderate-strength shock waves.

Howlett, James T.↗

Holocinematographic velocimeter for measuring time-dependent, three-dimensional flows

Two simulatneous, orthogonal-axis holographic movies are made of tracer particles in a low-speed water tunnel to determine the time-dependent, three-dimensional velocity field. This instrument is called a Holocinematographic Velocimeter (HCV). The holographic movies are reduced to the velocity field with an automatic data reduction system. This permits the reduction of large numbers of holograms (time steps) in a reasonable amount of time. The current version of the HCV, built for proof-of-concept tests, uses low-frame rate holographic cameras and a prototype of a new type of water tunnel. This water tunnel is a unique low-disturbance facility which has minimal wall effects on the flow. This paper presents the first flow field examined by the HCV, the two-dimensional von Karman vortex street downstream of an unswept circular cylinder. Key factors in the HCV are flow speed, spatial and temporal resolution required, measurement volume, film transport speed, and laser pulse length. The interactions between these factors are discussed.

Beeler, George B.↗

Efficient computation of inviscid flow fields around complex configurations using a multiblock multigrid method

The applicability of a multigrid technique to block-structured, body-fitted meshes is examined focusing on three different strategies. In the first strategy data are exchanged between blocks in each stage of a five-stage Runge-Kutta time-stepping scheme which keeps a possible time lag between blocks to a minimum, but requires a large amount of I/O operations and storage. The second strategy is based on performing a complete Runge-Kutta cycle within a block before switching to the next. In the third strategy both a complete Runge-Kutta cycle and the residual evaluation for the restriction operator are done within a block, allowing a minimum of I/O and storage. The inviscid flow around a wing-body/engine-pylon configuration was computed on a mesh consisting of 11 computational blocks. It was found that both the first and the second strategies delivered converged results, but the third failed due to larger time lag between blocks.

Rossow, C.-C.↗

Predicting Accumulations of Ice on Aerodynamic Surfaces

LEWICE is a computer program that predicts the accumulation of ice on two-dimensional aerodynamic surfaces under conditions representative of the flight of an aircraft through an icing cloud. The software first calculates the airflow surrounding the body of interest, then uses the airflow to compute the trajectories of water droplets that impinge on the surface of the body. The droplet trajectories are also used to compute impingement limits and local collection efficiencies, which are used in subsequent ice-growth calculations and are also useful for designing systems to protect against icing. Next, the software predicts the shape of accumulating ice by modeling transfers of mass and energy in small control volumes. The foregoing computations are repeated over several computational time steps until the total icing exposure time is reached. Results of computations by LEWICE have been compared with an extensive database of measured ice shapes obtained from experiments, and have been shown to closely approximate those shapes under most conditions of interest to the aviation community.

Bidwell, Colin↗

Enhancing Discovery, Search, and Access of NASA Hydrological Data by Leveraging GEOSS

An ongoing NASA-funded project has removed a longstanding barrier to accessing NASA data (i.e., accessing archived time-step array data as point-time series) for selected variables of the North American and Global Land Data Assimilation Systems (NLDAS and GLDAS, respectively) and other EOSDIS (Earth Observing System Data Information System) data sets (e.g., precipitation, soil moisture). These time series (data rods) are pre-generated. Data rods Web services are accessible through the CUAHSI Hydrologic Information System (HIS) and the Goddard Earth Sciences Data and Information Services Center (GES DISC) but are not easily discoverable by users of other non-NASA data systems. The Global Earth Observation System of Systems (GEOSS) is a logical mechanism for providing access to the data rods. An ongoing GEOSS Water Services project aims to develop a distributed, global registry of water data, map, and modeling services cataloged using the standards and procedures of the Open Geospatial Consortium and the World Meteorological Organization. The ongoing data rods project has demonstrated the feasibility of leveraging the GEOSS infrastructure to help provide access to time series of model grid information or grids of information over a geographical domain for a particular time interval. A recently-begun, related NASA-funded ACCESS-GEOSS project expands on these prior efforts. Current work is focused on both improving the performance of the generation of on-the-fly (OTF) data rods and the Web interfaces from which users can easily discover, search, and access NASA data.

GEOSS↗

Application of Data Cubes for Improving Detection of Water Cycle Extreme Events

As part of an ongoing NASA-funded project to remove a longstanding barrier to accessing NASA data (i.e., accessing archived time-step array data as point-time series), for the hydrology and other point-time series-oriented communities, "data cubes" are created from which time series files (aka "data rods") are generated on-the-fly and made available as Web services from the Goddard Earth Sciences Data and Information Services Center (GES DISC). Data cubes are data as archived rearranged into spatio-temporal matrices, which allow for easy access to the data, both spatially and temporally. A data cube is a specific case of the general optimal strategy of reorganizing data to match the desired means of access. The gain from such reorganization is greater the larger the data set. As a use case of our project, we are leveraging existing software to explore the application of the data cubes concept to machine learning, for the purpose of detecting water cycle extreme events, a specific case of anomaly detection, requiring time series data. We investigate the use of support vector machines (SVM) for anomaly classification. We show an example of detection of water cycle extreme events, using data from the Tropical Rainfall Measuring Mission (TRMM).

water cycle extreme events↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗

Enhancement of Mutual Discovery, Search, and Access of Data for Users of NASA and GEOSS-Cataloged Data Systems

An ongoing NASA-funded Data Rods (time series) project has demonstrated the removal of a longstanding barrier to accessing NASA data (i.e., accessing archived time-step array data as point-time series) for selected variables of the North American and Global Land Data Assimilation Systems (NLDAS and GLDAS, respectively) and other NASA data sets. Data rods are pre-generated or generated on-the-fly (OTF), leveraging the NASA Simple Subset Wizard (SSW), a gateway to NASA data centers. Data rods Web services are accessible through the CUAHSI Hydrologic Information System (HIS) and the Goddard Earth Sciences Data and Information Services Center (GES DISC) but are not easily discoverable by users of other non-NASA data systems. An ongoing GEOSS Water Services project aims to develop a distributed, global registry of water data, map, and modeling services cataloged using the standards and procedures of the Open Geospatial Consortium and the World Meteorological Organization. Preliminary work has shown GEOSS can be leveraged to help provide access to data rods. A new NASA-funded project is extending this early work.

data access↗

A Time-Accurate Upwind Unstructured Finite Volume Method for Compressible Flow with Cure of Pathological Behaviors

A time-accurate, upwind, finite volume method for computing compressible flows on unstructured grids is presented. The method is second order accurate in space and time and yields high resolution in the presence of discontinuities. For efficiency, the Roe approximate Riemann solver with an entropy correction is employed. In the basic Euler/Navier-Stokes scheme, many concepts of high order upwind schemes are adopted: the surface flux integrals are carefully treated, a Cauchy-Kowalewski time-stepping scheme is used in the time-marching stage, and a multidimensional limiter is applied in the reconstruction stage. However even with these up-to-date improvements, the basic upwind scheme is still plagued by the so-called "pathological behaviors," e.g., the carbuncle phenomenon, the expansion shock, etc. A solution to these limitations is presented which uses a very simple dissipation model while still preserving second order accuracy. This scheme is referred to as the enhanced time-accurate upwind (ETAU) scheme in this paper. The unstructured grid capability renders flexibility for use in complex geometry; and the present ETAU Euler/Navier-Stokes scheme is capable of handling a broad spectrum of flow regimes from high supersonic to subsonic at very low Mach number, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics). Numerous examples are included to demonstrate the robustness of the methods.

Loh, Ching Y.↗

Study of Unsteady Flows with Concave Wall Effect

This paper presents computational fluid dynamic studies of the inlet turbulence and wall curvature effects on the flow steadiness at near wall surface locations in boundary layer flows. The time-stepping RANS numerical solver of the NASA Glenn-HT RANS code and a one-equation turbulence model, with a uniform inlet turbulence modeling level of the order of 10 percent of molecular viscosity, were used to perform the numerical computations. The approach was first calibrated for its predictabilities of friction factor, velocity, and temperature at near surface locations within a transitional boundary layer over concave wall. The approach was then used to predict the velocity and friction factor variations in a boundary layer recovering from concave curvature. As time iteration proceeded in the computations, the computed friction factors converged to their values from existing experiments. The computed friction factors, velocity, and static temperatures at near wall surface locations oscillated periodically in terms of time iteration steps and physical locations along the span-wise direction. At the upstream stations, the relationship among the normal and tangential velocities showed vortices effects on the velocity variations. Coherent vortices effect on the velocity components broke down at downstream stations. The computations also predicted the vortices effects on the velocity variations within a boundary layer flow developed along a concave wall surface with a downstream recovery flat wall surface. It was concluded that the computational approach might have the potential to analyze the flow steadiness in a turbine blade flow.

Wang, Chi R.↗

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool.

Chima, R. V.↗

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool. Previously announced in STAR as N84-22527

Chima, R. V.↗

Regenerative life support system research and concepts

Life support systems that involve recycling of atmospheres, water, food and waste are so complex that models incorporating all the interactions and relationships are vital to design, development, simulations, and ultimately to control of space qualified systems. During early modeling studies, FORTRAN and BASIC programs were used to obtain numerical comparisons of the performance of different regenerative concepts. Recently, models were made by combining existing capabilities with expert systems to establish an Intelligent Design Support Environment for simpliflying user interfaces and to address the need for the engineering aspects. Progress was also made toward modeling and evaluating the operational aspects of closed loop life support systems using Time-step and Dynamic simulations over a period of time. Example models are presented which show the status and potential of developed modeling techniques. For instance, closed loop systems involving algae systeMs for atmospheric purification and food supply augmentation, plus models employing high plants and solid waste electrolysis are described and results of initial evaluations are presented.

Source record↗

Modeling of Space Station power system components and their interactions

The authors demonstrate that the functional modeling approach is feasible in modeling large-scale spacecraft power systems such as the Space Station electric power system (EPS). The functional models for the solar array sequential shunt unit and the battery charging/discharging unit are presented. The functional modeling approach has also been applied to simulate other major EPS components. The usefulness of this approach in system-level studies has been demonstrated by incorporating these models into a power system model and simulating the system behavior under various conditions. The simulation results are consistent with those obtained experimentally. With larger time step size and reduced model complexity, the computation time required by the functional modeling approach is short. Functional models are effective in simulation studies of system-level issues, especially for large-scale power systems including many power electronic components.

Tam, Kwa-Sur↗

Solving modal equations of motion with initial conditions using MSC/NASTRAN DMAP. Part 1: Implementing exact mode superposition

Within the MSC/NASTRAN DMAP (Direct Matrix Abstraction Program) module TRD1, solving physical (coupled) or modal (uncoupled) transient equations of motion is performed using the Newmark-Beta or mode superposition algorithms, respectively. For equations of motion with initial conditions, only the Newmark-Beta integration routine has been available in MSC/NASTRAN solution sequences for solving physical systems and in custom DMAP sequences or alters for solving modal systems. In some cases, one difficulty with using the Newmark-Beta method is that the process of selecting suitable integration time steps for obtaining acceptable results is lengthy. In addition, when very small step sizes are required, a large amount of time can be spent integrating the equations of motion. For certain aerospace applications, a significant time savings can be realized when the equations of motion are solved using an exact integration routine instead of the Newmark-Beta numerical algorithm. In order to solve modal equations of motion with initial conditions and take advantage of efficiencies gained when using uncoupled solution algorithms (like that within TRD1), an exact mode superposition method using MSC/NASTRAN DMAP has been developed and successfully implemented as an enhancement to an existing coupled loads methodology at the NASA Lewis Research Center.

Abdallah, Ayman A.↗

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗

Computer Modeling of Electrostatic Aggregation of Granular Materials in Planetary and Astrophysical Settings

Electrostatic forces strongly influence the behavior of granular materials in both dispersed (cloud) systems and semi-packed systems. These forces can cause aggregation or dispersion of particles and are important in a variety of astrophysical and planetary settings. There are also many industrial and commercial settings where granular matter and electrostatics become partners for both good and bad. This partnership is important for human exploration on Mars where dust adheres to suits, machines, and habitats. Long-range Coulombic (electrostatic) forces, as opposed to contact-induced dipoles and van der Waals attractions, are generally regarded as resulting from net charge. We have proposed that in addition to net charge interactions, randomly distributed charge carriers on grains will result in a dipole moment regardless of any net charge. If grains are unconfined, or fluidized, they will rotate so that the dipole always induces attraction between grains. Aggregates are readily formed, and Coulombic polarity resulting from the dipole produces end-to-end stacking of grains to form filamentary aggregates. This has been demonstrated in USML experiments on Space Shuttle where microgravity facilitated the unmasking of static forces. It has also been demonstrated in a computer model using grains with charge carriers of both sign. Model results very closely resembled micro-g results with actual sand grains. Further computer modeling of the aggregation process has been conducted to improve our understanding of the aggregation process, and to provide a predictive tool for microgravity experiments slated for Space Station. These experiments will attempt to prove the dipole concept as outlined above. We have considerably enhanced the original computer model: refinements to the algorithm have improved the fidelity of grain behavior during grain contact, special attention has been paid to simulation time steps to enable establishment of a meaningful, quantitative time axis, and calibration of rounding accuracies have been conducted to test cumulative numerical influences in the model. The model has been run for larger grain populations, variable initial cloud densities, and we have introduced random net charging to individual grains, as well as a net charge to the cloud as a whole. The model uses 3 positive and 3 negative charges randomly distributed on each grain, with up to 160 grains contained within various size "boxes" that define the initial number densities in the clouds. Each charge represents localized charged region on a grain, but does not necessarily imply single quantized charge carriers. The Coulomb equations are then allowed to interact for each monopole: dipoles and any higher order charge coupling is a natural product of these "free" interactions over which the modeler exerts no influence. The charges are placed on surfaces of grains at random locations. A series of runs was conducted for neutral grains that had a perfect balance of negative and positive char carriers. Runs were also conducted with grains having additional fractional charges ranging between 0 and 1. By adding fractional charges of one sign, the model created grain populations in which all grains had excess charges the same sign, giving the cloud an overall net charge. This simulates clouds subjected to ionizing radiation (e. protoplanetary debris disk around a protosun), or any other process of charge biasing in a grain population (e.g., volcanic plumes). In another run series, random fractional charges of either sign were added to the grains so th some grains had a slight net positive charge while others had a slight net negative charge. This simulates triboelectrically-charged grain populations in which acquisition of an electron by one surface is at the expense creating a hole elsewhere. This dual sign charging was applied in two ways: in one case the cloud remained neutral by ensuring that all grain excess charges added to zero; in the other case, the cloud was permitted slight net char by not imposing a charge-balance condition. Additional information is contained in the original.

Marshall, J.↗