Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “implicit time integration”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Golden Ratio Versus Pi as Random Sequence Sources for Monte Carlo Integration

We discuss here the relative merits of these numbers as possible random sequence sources. The quality of these sequences is not judged directly based on the outcome of all known tests for the randomness of a sequence. Instead, it is determined implicitly by the accuracy of the Monte Carlo integration in a statistical sense. Since our main motive of using a random sequence is to solve real world problems, it is more desirable if we compare the quality of the sequences based on their performances for these problems in terms of quality/accuracy of the output. We also compare these sources against those generated by a popular pseudo-random generator, viz., the Matlab rand and the quasi-random generator ha/ton both in terms of error and time complexity. Our study demonstrates that consecutive blocks of digits of each of these numbers produce a good random sequence source. It is observed that randomly chosen blocks of digits do not have any remarkable advantage over consecutive blocks for the accuracy of the Monte Carlo integration. Also, it reveals that pi is a better source of a random sequence than theta when the accuracy of the integration is concerned.

Sen, S. K.↗
NASA NTRS1984

Development of a dynamically adaptive grid method for multidimensional problems

An approach to solution adaptive grid generation for use with finite difference techniques, previously demonstrated on model problems in one space dimension, has been extended to multidimensional problems. The method is based on the popular elliptic steady grid generators, but is 'dynamically' adaptive in the sense that a grid is maintained at all times satisfying the steady grid law driven by a solution-dependent source term. Testing has been carried out on Burgers' equation in one and two space dimensions. Results appear encouraging both for inviscid wave propagation cases and viscous boundary layer cases, suggesting that application to practical flow problems is now possible. In the course of the work, obstacles relating to grid correction, smoothing of the solution, and elliptic equation solvers have been largely overcome. Concern remains, however, about grid skewness, boundary layer resolution and the need for implicit integration methods. Also, the method in 3-D is expected to be very demanding of computer resources.

Holcomb, J. E.↗
NASA NTRS1990

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗
NASA NTRS2010

Thermal-Acoustic Analysis of a Metallic Integrated Thermal Protection System Structure

A study is undertaken to investigate the response of a representative integrated thermal protection system structure under combined thermal, aerodynamic pressure, and acoustic loadings. A two-step procedure is offered and consists of a heat transfer analysis followed by a nonlinear dynamic analysis under a combined loading environment. Both analyses are carried out in physical degrees-of-freedom using implicit and explicit solution techniques available in the Abaqus commercial finite-element code. The initial study is conducted on a reduced-size structure to keep the computational effort contained while validating the procedure and exploring the effects of individual loadings. An analysis of a full size integrated thermal protection system structure, which is of ultimate interest, is subsequently presented. The procedure is demonstrated to be a viable approach for analysis of spacecraft and hypersonic vehicle structures under a typical mission cycle with combined loadings characterized by largely different time-scales.

Behnke, Marlana N.↗
NASA NTRS1988

Unsteady transonic airfoil computation using the integral solution of full-potential equation

The shock-capturing integral-equation scheme developed by Kandil and Hu (1987) for the analysis of steady transonic flow over airfoils is extended to the unsteady case. The full potential formulation of the governing equations is reviewed; the solution method is outlined; and results for a NACA 0012 airfoil in forced pitching oscillation at Mach 0.755 are presented in extensive graphs and briefly characterized. The present technique is shown to require significantly less computation time than finite-difference or finite-volume methods, and to give shock-motion predictions in good agreement with those of an implicit finite-volume Euler solver; the surface-pressure peaks are slightly underpredicted.

Kandil, O. A.↗
NASA NTRS2006

Precision of Sensitivity in the Design Optimization of Indeterminate Structures

Design sensitivity is central to most optimization methods. The analytical sensitivity expression for an indeterminate structural design optimization problem can be factored into a simple determinate term and a complicated indeterminate component. Sensitivity can be approximated by retaining only the determinate term and setting the indeterminate factor to zero. The optimum solution is reached with the approximate sensitivity. The central processing unit (CPU) time to solution is substantially reduced. The benefit that accrues from using the approximate sensitivity is quantified by solving a set of problems in a controlled environment. Each problem is solved twice: first using the closed-form sensitivity expression, then using the approximation. The problem solutions use the CometBoards testbed as the optimization tool with the integrated force method as the analyzer. The modification that may be required, to use the stiffener method as the analysis tool in optimization, is discussed. The design optimization problem of an indeterminate structure contains many dependent constraints because of the implicit relationship between stresses, as well as the relationship between the stresses and displacements. The design optimization process can become problematic because the implicit relationship reduces the rank of the sensitivity matrix. The proposed approximation restores the full rank and enhances the robustness of the design optimization method.

Patnaik, Surya N.↗
NASA NTRS2010

Software Model Checking of ARINC-653 Flight Code with MCP

The ARINC-653 standard defines a common interface for Integrated Modular Avionics (IMA) code. In particular, ARINC-653 Part 1 specifies a process- and partition-management API that is analogous to POSIX threads, but with certain extensions and restrictions intended to support the implementation of high reliability flight code. MCP is a software model checker, developed at NASA Ames, that provides capabilities for model checking C and C++ source code. In this paper, we present recent work aimed at implementing extensions to MCP that support ARINC-653, and we discuss the challenges and opportunities that consequentially arise. Providing support for ARINC-653 s time and space partitioning is nontrivial, though there are implicit benefits for partial order reduction possible as a consequence of the API s strict interprocess communication policy.

Thompson, Sarah J.↗
NASA NTRS2011

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗
NASA NTRS1992

Diabatic dynamic initialization

A generalized dynamical adjustment procedure was applied to a diabatic model to produce balanced initial conditions. Namely, backward adiabatic model integration is followed by forward diabatic model integration, with a high frequency (low pass) filter in the form of the Euler-backward time differencing scheme being applied throughout the whole integration. As a result of the application of such a diabatic dynamic initialization procedure within the Goddard Laboratory for Atmospheres (GLA) 4-D data assimilation system, the following properties of forecasts from initialized fields are achieved right from the beginning of the usual forecast integration: (1) the forecast tendencies (and fields) are free of any noise due to imbalance in initial conditions; (2) the shocks related to an initial imbalance between model physics and dynamics, and especially the substantial initial imbalance of precipitation and evaporation fields, or the initial spin-up effect, are practically removed. Diabatic dynamic initialization has been compared with implicit nonlinear normal mode initialization, and found to be superior in removing the initial spin-up effect, and also in improving the tropical structure. The diabatic dynamic initialization procedure was successfully tested for the GLA system with the use of all conventional data and the GLA satellite data retrievals. It allows a smooth data insertion without any shocks or imbalances, which is highly desirable for efficient functioning of 4-D data assimilation systems. The developed initialization procedure is computationally efficient and in principle easily applicable to different forecast models.

Fox-Rabinovitz, Michael S.↗
NASA NTRS1990

Numerical computation of viscous flow about unconventional airfoil shapes

A new two-dimensional computer code was developed to analyze the viscous flow around unconventional airfoils at various Mach numbers and angles of attack. The Navier-Stokes equations are solved using an implicit, upwind, finite-volume scheme. Both laminar and turbulent flows can be computed. A new nonequilibrium turbulence closure model was developed for computing turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present code was evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the present code in conjunction with the new nonequilibrium turbulence model gives excellent results.

Ahmed, S.↗
NASA NTRS2010

Trajectory Optimization: OTIS 4

The latest release of the Optimal Trajectories by Implicit Simulation (OTIS4) allows users to simulate and optimize aerospace vehicle trajectories. With OTIS4, one can seamlessly generate optimal trajectories and parametric vehicle designs simultaneously. New features also allow OTIS4 to solve non-aerospace continuous time optimal control problems. The inputs and outputs of OTIS4 have been updated extensively from previous versions. Inputs now make use of objectoriented constructs, including one called a metastring. Metastrings use a greatly improved calculator and common nomenclature to reduce the user s workload. They allow for more flexibility in specifying vehicle physical models, boundary conditions, and path constraints. The OTIS4 calculator supports common mathematical functions, Boolean operations, and conditional statements. This allows users to define their own variables for use as outputs, constraints, or objective functions. The user-defined outputs can directly interface with other programs, such as spreadsheets, plotting packages, and visualization programs. Internally, OTIS4 has more explicit and implicit integration procedures, including high-order collocation methods, the pseudo-spectral method, and several variations of multiple shooting. Users may switch easily between the various methods. Several unique numerical techniques such as automated variable scaling and implicit integration grid refinement, support the integration methods. OTIS4 is also significantly more user friendly than previous versions. The installation process is nearly identical on various platforms, including Microsoft Windows, Apple OS X, and Linux operating systems. Cross-platform scripts also help make the execution of OTIS and post-processing of data easier. OTIS4 is supplied free by NASA and is subject to ITAR (International Traffic in Arms Regulations) restrictions. Users must have a Fortran compiler, and a Python interpreter is highly recommended.

Riehl, John P.↗
NASA NTRS1997

A comparative modeling of supernova 1993J

The light curve of Supernova 1993J is calculated using two computational radiation transport approaches. The two approaches are represented by the computer codes STELLA and EDDINGTON. The emphasis is on the shock breakout and the photometry in the U, B and V bands during the first 120 days of the supernova. The STELLA model includes implicit hydrodynamics and is able to model early supernova evolution before the expansion is homologous. The STELLA model employs multi-group photonics and is able to follow the radiation as it decouples from the matter. The EDDINGTON code uses an algorithm for integrating the transport equation which assumes homologous expansion and uses a finer frequency resolution. The agreement between the two codes is considered to be satisfactory only in the case where compatible physical assumptions are made concerning the opacity. The assumptions are justified. The continuum spectrum for SN 1993J is predicted near the shock breakout to be superior to that predicted by standard single energy group hydrocodes. The uncertainties involved in current time dependent models of supernova light curves are discussed.

Blinnikov, Sergei↗
NASA NTRS1995

Diabatic forcing and intialization with assimilation of cloud water and rainwater in a forecast model

In this study, diabatic forcing, and liquid water assimilation techniques are tested in a semi-implicit hydrostatic regional forecast model containing explicit representations of grid-scale cloud water and rainwater. Diabatic forcing, in conjunction with diabatic contributions in the initialization, is found to help the forecast retain the diabatic signal found in the liquid water or heating rate data, consequently reducing the spinup time associated with grid-scale precipitation processes. Both observational Special Sensor Microwave/Imager (SSM/I) and model-generated data are used. A physical retrieval method incorporating SSM/I radiance data is utilized to estimate the 3D distribution of precipitating storms. In the retrieval method the relationship between precipitation distributions and upwelling microwave radiances is parameterized, based upon cloud ensemble-radiative model simulations. Regression formulae relating vertically integrated liquid and ice-phase precipitation amounts to latent heating rates are also derived from the cloud ensemble simulations. Thus, retrieved SSM/I precipitation structures can be used in conjunction with the regression-formulas to infer the 3D distribution of latent heating rates. These heating rates are used directly in the forecast model to help initiate Tropical Storm Emily (21 September 1987). The 14-h forecast of Emily's development yields atmospheric precipitation water contents that compare favorably with coincident SSM/I estimates.

Raymond, William H.↗
NASA NTRS1991

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Tamma, Kumar K.↗
NASA NTRS1972

Ephemeris of a highly eccentric orbit - Explorer 28.

An ephemeris has been obtained for Explorer 28 (IMP 3) which agrees well with 2 years of radio observations and with SAO observations a year later. This ephemeris is generated over the 3 year lifetime by a numerical integration method utilizing a set of initial conditions at launch and without requiring further differential correction. Because highly eccentric orbits are difficult to compute with acceptable accuracy and because a long continuous arc has been obtained which compares with actual data to a known precision, this ephemeris may be used as a standard for computing highly eccentric orbits in the earth-moon system. Orbit improvement was used to obtain the initial conditions which generated the ephemeris. This improvement was based on correcting the energy by adjusting the semimajor axis to match computed times of perigee passage with the observed. This procedure may generate errors in semimajor axis to compensate for model errors in the energy; however, this compensation error is also implicit in orbit determination itself.

Lowrey, B. E.↗
NASA NTRSDate not supplied

Full Stall Simulations of a Redesigned Ventilation Fan for the ISS

The concept of a stall is studied rigorously in the aerospace industry. From a design standpoint, instabilities such as stall are undesirable– operation in the stall regime has a tremendous impact on aerodynamic performance as well as structural integrity. In extreme cases, operating in stall conditions can cause failure. Further, stall can cause a loss of lift on aircraft wings or a loss of thrust in aircraft engines. In any case, stall continues to be a topic of interest in the aerospace industry. The present work aims to analyze the stall characteristics of a ventilation fan that was recently designed for the International Space Station (ISS). Although the ventilation fan has a rotor-stator design, this paper considers a rotor-only configuration. The FUN3D Computational Fluid Dynamic (CFD) solver developed by NASA Langley Research Center was used to simulate the operational characteristics of the ventilation fan. FUN3D solves the Unsteady Reynolds-Averaged Naiver-Stokes (URANS) equations using implicit time marching and a dynamic overset grid. The FUN3D solver was originally written for exterior flow fields; however, this work represents an extension of the FUN3D solver to turbomachinery or interior flow fields. The FUN3D results for the rotor-only ventilation fan accurately captured the operational characteristics inherent to compressors– the results shared similar performance trends when compared to the experimental results for the rotor-stator case. A peak adiabatic efficiency of 96% occurred at a MFR of 105.2 CFM, the minimum aerodynamically stable point. A computationally stable stall occurred at a mass flow rate of 43.8 CFM where the adiabatic efficiency dropped to 69%.

CFD↗
NASA NTRS1985

A hybrid numerical technique for predicting the aerodynamic and acoustic fields of advanced turboprops

A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.

Homicz, G. F.↗
NASA NTRS1998

Thermal Infrared Remote Sensing for Analysis of Landscape Ecological Processes: Methods and Applications

Thermal Infrared (TIR) remote sensing data can provide important measurements of surface energy fluxes and temperatures, which are integral to understanding landscape processes and responses. One example of this is the successful application of TIR remote sensing data to estimate evapotranspiration and soil moisture, where results from a number of studies suggest that satellite-based measurements from TIR remote sensing data can lead to more accurate regional-scale estimates of daily evapotranspiration. With further refinement in analytical techniques and models, the use of TIR data from airborne and satellite sensors could be very useful for parameterizing surface moisture conditions and developing better simulations of landscape energy exchange over a variety of conditions and space and time scales. Thus, TIR remote sensing data can significantly contribute to the observation, measurement, and analysis of energy balance characteristics (i.e., the fluxes and redistribution of thermal energy within and across the land surface) as an implicit and important aspect of landscape dynamics and landscape functioning. The application of TIR remote sensing data in landscape ecological studies has been limited, however, for several fundamental reasons that relate primarily to the perceived difficulty in use and availability of these data by the landscape ecology community, and from the fragmentation of references on TIR remote sensing throughout the scientific literature. It is our purpose here to provide evidence from work that has employed TIR remote sensing for analysis of landscape characteristics to illustrate how these data can provide important data for the improved measurement of landscape energy response and energy flux relationships. We examine the direct or indirect use of TIR remote sensing data to analyze landscape biophysical characteristics, thereby offering some insight on how these data can be used more robustly to further the understanding and modeling of landscape ecological processes.

Quattrochi, Dale A.↗