Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Based Procedures”

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 865 records · Page 48

An Examination of Selected Datacom Options for the Near-Term Implementation of Trajectory Based Operations

A primary feature of the Next Generation Air Transportation System (NextGen) is trajectory based operations (TBO). Under TBO, aircraft flight plans are known to computer systems on the ground that aid in scheduling and separation. The Future Air Navigation System (FANS) was developed to support TBO, but relatively few aircraft in the US are FANSequipped. Thus, any near-term implementation must provide TBO procedures for non-FANS aircraft. Previous research has explored controller clearances, but any implementation must also provide procedures for aircraft requests. The work presented here aims to surface issues surrounding TBO communication procedures for non-FANS aircraft and for aircraft requesting deviations around weather. Three types of communication were explored: Voice, FANS, and ACARS,(Aircraft Communications Addressing and Reporting System). ACARS and FANS are datacom systems that differ in that FANS allows uplinked flight plans to be loaded into the Flight Management System (FMS), while ACARS delivers flight plans as text that must be entered manually via the Control Display Unit (CDU). Sixteen pilots (eight two-person flight decks) and four controllers participated in 32 20-minute scenarios that required the flight decks to navigate through convective weather as they approached their top of descents (TODs). Findings: The rate of non-conformance was higher than anticipated, with aircraft off path more than 20% of the time. Controllers did not differentiate between the ACARS and FANS datacom, and were mixed in their preference for Voice vs. datacom (ACARS and FANS). Pilots uniformly preferred Voice to datacom, particularly ACARS. Much of their dislike appears to result from the slow response times in the datacom conditions. As a result, participants frequently resorted to voice communication. These results imply that, before implementing TBO in environments where pilots make weather deviation requests, further research is needed to develop communication procedures that integrate voice and datacom.

trajectory based operations↗

Chapter 14: Dynamic Flux Analysis: An Experimental Approach of Fluxomics

Metabolic flux analysis represents an essential perspective to understand cellular physiology and offers quantitative information to guide pathway engineering. A valuable approach for experimental elucidation of metabolic flux is dynamic flux analysis, which estimates the relative or absolute flow rates through a series of metabolic intermediates in a given pathway. It is based on kinetic isotope labeling experiments, liquid chromatography-mass spectrometry (LC-MS), and computational analysis that relate kinetic isotope trajectories of metabolites to pathway activity. Herein, we illustrate the mathematic principles underlying the dynamic flux analysis and mainly focus on describing the experimental procedures for data generation. This protocol is exemplified using cyanobacterial metabolism as an example, for which reliable labeling data for central carbon metabolites can be acquired quantitatively. This protocol is applicable to other microbial systems as well and can be readily adapted to address different metabolic processes.

BASIC BIOLOGICAL SCIENCES↗

Inference for Nonparanormal Partial Correlation via Regularized Rank-Based Nodewise Regression

Abstract Partial correlation is a common tool in studying conditional dependence for Gaussian distributed data. However, partial correlation being zero may not be equivalent to conditional independence under non-Gaussian distributions. In this paper, we propose a statistical inference procedure for partial correlations under the high-dimensional nonparanormal (NPN) model where the observed data are normally distributed after certain monotone transformations. The NPN partial correlation is the partial correlation of the normal transformed data under the NPN model, which is a more general measure of conditional dependence. We estimate the NPN partial correlations by regularized nodewise regression based on the empirical ranks of the original data. A multiple testing procedure is proposed to identify the nonzero NPN partial correlations. The proposed method can be carried out by a simple coordinate descent algorithm for lasso optimization. It is easy-to-implement and computationally more efficient compared to the existing methods for estimating NPN graphical models. Theoretical results are developed to show the asymptotic normality of the proposed estimator and to justify the proposed multiple testing procedure. Numerical simulations and a case study on brain imaging data demonstrate the utility of the proposed procedure and evaluate its performance compared to the existing methods. Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database.

97 MATHEMATICS AND COMPUTING↗

Nonlinear panel flutter - A finite-element approach

A finite-element approach has been developed for computing nonlinear flutter characteristics of rectangular isotropic panels with stream-alined side edges, based on aerodynamic forces from supersonic two-dimensional quasi-steady aerodynamic theory. Stress distributions and panel oscillation frequencies were determined from the analysis. The finite-element formulation, solution procedure, and convergence characteristics are presented. Comparisons are made with linear flutter and large-amplitude vibration results and demonstrate that good accuracy is obtained. Non-linear flutter results are presented for effects of aerodynamic damping, length-width ratio, initial in-plane forces and boundary-support conditions. Comparisons with experimental results are also presented.

Mei, C.↗

Parameter identification methods for improving structural dynamic models

There is an increasing need to develop Parameter Identification methods for improving structural dynamic models, based on the inability of engineers to produce mathematical models which correlate with experimental data. This research explores the efficiency of combining Component Mode Synthesis (substructuring) methods with Parameter Identification procedures in order to improve analytical modeling of structural components and their connections. Improvements are computed in terms of physical stiffness and damping parameters in order that the physical characteristics of the model can be better understood. Connections involving both viscous and friction damping are investigated. Substructuring methods are utilized to reduce the complexity of the identification problem. Component and inter-component structural connection properties are evaluated and identified independently, thus simplifying the identification problem. It is shown that modal test data is effective for identifying modeling problems associated with structural components, and for determining the stiffness and damping properties of intercomponent connections. In general, Parameter Identification is improved when greater quantities of experimental data are available.

Lawrence, Charles↗

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗

Improved inhalation technology for setting safe exposure levels for workplace chemicals

Threshold Limit Values recommended as allowable air concentrations of a chemical in the workplace are often based upon a no-observable-effect-level (NOEL) determined by experimental inhalation studies using rodents. A 'safe level' for human exposure must then be estimated by the use of generalized safety factors in attempts to extrapolate from experimental rodents to man. The recent development of chemical-specific physiologically-based toxicokinetics makes use of measured physiological, biochemical, and metabolic parameters to construct a validated model that is able to 'scale-up' rodent response data to predict the behavior of the chemical in man. This procedure is made possible by recent advances in personal computer software and the emergence of appropriate biological data, and provides an analytical tool for much more reliable risk evaluation and airborne chemical exposure level setting for humans.

Stuart, Bruce O.↗

An RR Lyrae period shift in terms of the Fourier parameter Phi sub 31

The Fourier phase parameter Phi sub 31 has been determined for RRc stars in five globular clusters, NGC 6171, M5, M3, M53, and M15. The results indicate that the RRc stars in a given cluster show a sequence of Phi sub 31 increasing with period, and that the higher the cluster metallicity, the higher the sequence lies in a plot of Phi sub 31 with period. The Phi sub 31 values for the stars in NGC 6171 and M5 presented here are based on observations made with the University of Toronto 0.61 m telescope at Las Campanas, Chile, while those for M3, M53, and M15 are based on published data. A bootstrap procedure has been used to establish the uncertainties in the Fourier parameters. The physical significance of the relationship among Phi sub 31, period, and metallicity is not yet understood. It will need to be tested with hydrodynamic pulsation models computed with new opacities.

Clement, Christine M.↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

Vortical Flow Prediction Using an Adaptive Unstructured Grid Method

A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving practical vortical flow problems. The first test case concerns vortex flow over a simple 65' delta wing with different values of leading-edge bluntness, and the second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the windtunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.

Pirzadeh, Shahyar Z.↗

Vortical Flow Prediction Using an Adaptive Unstructured Grid Method

A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving practical vortical flow problems. The first test case concerns vortex flow over a simple 65deg delta wing with different values of leading-edge bluntness, and the second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the windtunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.

Pirzadeh, Shahyar Z.↗

A Geometry Based Infra-structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis (CFD, structural analysis, and etc.) are: Surface Generation - usually by employing a CAD system; Grid Generation - preparing the volume for the simulation; Flow Solver - producing the results at the specified operational point; and Post-processing Visualization - interactively attempting to understand the results For structural analysis, integrated systems can be obtained from a number of commercial vendors. For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. Specifically the problems with this procedure are: (1) File based. Information flows from one step to the next via data files with formats specified for that procedure. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. (3) One-Way communication. All information travels on from one phase to the next. Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive.

Haimes, Robert↗

A-Priori Tuning of Modified Magnussen Combustion Model

In the application of CFD to turbulent reacting flows, one of the main limitations to predictive accuracy is the chemistry model. Using a full or skeletal kinetics model may provide good predictive ability, however, at considerable computational cost. Adding the ability to account for the interaction between turbulence and chemistry improves the overall fidelity of a simulation but adds to this cost. An alternative is the use of simple models, such as the Magnussen model, which has negligible computational overhead, but lacks general predictive ability except for cases that can be tuned to the flow being solved. In this paper, a technique will be described that allows the tuning of the Magnussen model for an arbitrary fuel and flow geometry without the need to have experimental data for that particular case. The tuning is based on comparing the results of the Magnussen model and full finite-rate chemistry when applied to perfectly and partially stirred reactor simulations. In addition, a modification to the Magnussen model is proposed that allows the upper kinetic limit for the reaction rate to be set, giving better physical agreement with full kinetic mechanisms. This procedure allows a simple reacting model to be used in a predictive manner, and affords significant savings in computational costs for simulations.

Norris, A. T.↗

Systematically derived thermodynamic properties for alkane oxidation

Key combustion properties, such as ignition delay time, show strong sensitivity to the thermochemistry of the main species in the standard radical oxidation pathway at low temperatures (600-1000K). Significant uncertainties persist in current estimates of thermodynamic properties, particularly for the larger species of relevance to common practical fuels. In this article, we use advanced computational schemes to evaluate thermodynamic properties for the fuel, fuel radical, peroxy, hydroperoxy-alkyl radical, and hydroperoxide species for a set of 17 fuels containing up to 9 carbon atoms, and with various degrees of branching in the alkane backbone. The procedure, termed STAR-1D, combines conformer sampling to find the minimum geometry, B2PLYP-D3/cc-pVTZ harmonic frequency evaluations, and omega B97X-D/cc-pVTZ one-dimensional torsional mapping. It includes two physically-based scaling routines: a frequency-dependent scaling to B2PLYP-D3/cc-pVTZ anharmonic frequencies and a scaling of one-dimensional omega B97X-D/ccpVTZ torsional profiles to reproduce the product of the B2PLYP-D3/cc-pVTZ frequencies in the harmonic limit. Substantive comparisons with existing experimental databases, together with careful examinations of key theoretical assumptions, are used to explore the accuracy of the predictions. These computationally intensive explorations of 195 species were facilitated by automated thermochemistry software. Further, high accuracy 0 K heats of formation from a separate study are used with the STAR-1D computations to generate NASA polynomial representations. In addition to their intrinsic value, the present results also provide a reliable database for the optimization of group additivity or machine learning schemes for scaling to larger combustion systems. Towards this end a complementary extensive conformational analysis is carried out for the medium sized species and the thermodynamic properties of the lowest energy hydrogen-bonded and non-hydrogen bonded conformers are contrasted for larger hydroperoxy-alkyl radical species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗

Orbit determination and prediction, and computer programs

Orbit determination and prediction programs are needed to generate ephemerides for the satellite. Orbit determination is from tracking data consisting of angles only, and is based on a modified version of a method by R. E. Briggs and J. W. Slowey of the Smithsonian Institution. Trends in the data due to perturbations from a Keplerian orbit are removed before this process, and estimates of the orbital elements from individual passes are combined statistically to produce refined estimates. Ephemeris calculation is by a semi-analytic method in which deviations from a Keplerian orbit are obtained by integrating the perturbing forces. The programs to implement these procedures have been written for both the IBM 7090 and the IBM 1620 computers.

EPHEMERIS↗

A recurrence matrix solution for the dynamic response of aircraft in gusts

A systematic procedure developed for the calculation of the structural response of aircraft flying through a gust by use of difference equations in the solution of dynamic problems is first illustrated by means of a simple-damped-oscillator example. A detailed analysis is then given which leads to a recurrence matrix equation for the determination of the response of an airplane in a gust. The method takes into account wing bending and twisting deformations, fuselage deflection, vertical and pitching motion of the airplane, and some tail forces. The method is based on aerodynamic strip theory, but compressibility and three-dimensional aerodynamic effects can be taken into account approximately by means of over-all corrections. Either a sharp-edge gust or a gust of arbitrary shape in the spanwise or flight directions may be treated. In order to aid in the application of the method to any specific case, a suggested computational procedure is included. The possibilities of applying the method to a variety of transient aircraft problems, such as landing, are brought out. A brief review of matrix algebra, covering the extent to which it is used in the analysis, is also included. (author)

Houbolt, John C↗

High‐order cubature rules for tetrahedra

Summary In this paper, we construct new high‐order numerical integration schemes for tetrahedra, with positive weights and integration points that are in the interior of the domain. The construction of cubature rules is a challenging problem, which requires the solution of strongly nonlinear algebraic (moment) equations with side conditions given by affine inequality constraints. We present a robust algorithm based on a sequence of three modified Newton procedures to solve the constrained minimization problem. In the literature, numerical integration rules for the tetrahedron are available up to order p =15. We obtain integration rules for the tetrahedron from p =2 to p =20, which are computed using multiprecision arithmetic. For p ≤15, our approach provides integration rules that have the same or fewer number of integration points than existing rules; for p =16 to p =20, our rules are new. Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, rational, weakly singular and trigonometric test functions over tetrahedra with flat and curved faces. In all tests, improvements in accuracy is realized as p is increased, though in some cases nonmonotonic convergence is observed.

Jaśkowiec, Jan↗