Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical methods and modeling”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A Numerical Investigation of Turbulent Flow in Noncircular Ducts

Since turbulent duct flows play an important role in engineering, continuous efforts to investigate this problem have been made. Until recently, these investigations were limited mostly to experiments and various semi-analytic methods. Owing to the improvement of turbulence modeling, the prevailing method of prediction is now mainly numerical. The majority of these studies deals, however, with turbulent flows in circular or square ducts and only limited information is available for straight noncircular ducts. In view of this situation, we propose to conduct a numerical investigation of turbulent flow in a class of ducts, whose cross sections vary from a circle to a near square. Turbulent flow in a noncircular duct is characterized by the presence of secondary flow for which a more refined turbulence model than the k-epsilon equations is required. In order to show that the calculated results are credible, various modes of verification were used to examine the results for a selected configuration, including an accuracy check by a scaling law and observing the decay of secondary flow as the cross section changes from a noncircular shape to a circle. After this was done, computations were performed for other configurations and with different Reynolds numbers from which wall shear stresses and friction factors are plotted.

Kao, Hsiao C.↗
DOE OSTI2024

Implementing Ordinary Differential Equation Solvers in Rust Programming Language for Modeling Vehicle Powertrain Systems: Preprint

Efficient and accurate ordinary differential equation (ODE) solvers are necessary for powertrain and vehicle dynamics modeling. However, current commercial ODE solvers can be financially prohibitive, leading to a need for accessible, effective, open-source ODE solvers designed for powertrain modeling. Rust is a compiled programming language that has the potential to be used for fast and easy-to-use powertrain models, given its exceptional computational performance, robust package ecosystem, and short time required for modelers to become proficient. However, of the three commonly used (>3,000 downloads) packages in Rust with ODE solver capabilities, only one has more than four numerical methods implemented, and none are designed specifically for modeling physical systems. Therefore, the goal of the Differential Equation System Solver (DESS) was to implement accurate ODE solvers in Rust designed for the component-based problems often seen in powertrain modeling. DESS is a text-based software package that provides a flexible framework for building and solving systems of ODEs. This allows DESS to be included as a dependency for automotive powertrain models that require a variety of solvers and solver configurations. Seven explicit ODE solver methods have been implemented in DESS: Euler’s, Heun’s, midpoint, Ralston’s, classic Runge-Kutta, Bogacki-Shampine, and Cash-Karp. These represent five fixed-step methods and two adaptive-step methods. This paper shows that the solver implementations increase accuracy and computational efficiency compared to Euler's method when modeling a system of three thermal masses in Rust. DESS also includes features designed for modeling component-based physical systems. Users can define relationships between nodes in their system, which the package then translates into a system of equations, leading to simpler and more intuitive code. In the case of a three-thermal-mass system, the user can specify node thermal properties (e.g., thermal capacitance), how nodes are interconnected, and thermal conductance between nodes rather than providing a system of equations. The core contribution from this work is an open-source, text-based Rust package with ODE solvers for automotive powertrain modeling to support cost-free, fast, and accurate simulation.

ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1986

Reynolds stress closure in jet flows using wave models

A wave model was developed for the two dimensional shear layer. This configuration is being used as a test case for the closure schemes. Numerical methods are under development to solve the nonseparable Rayleigh equation. A model problem is being used to assist in the algorithm development. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique are being examined.

Morris, P. J.↗
NASA NTRS2015

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗
DOE OSTI2025

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗
NASA NTRS2017

Probabilistic Forecasting of Surface Ozone with a Novel Statistical Approach

The recent change in the Environmental Protection Agency's surface ozone regulation, lowering the surface ozone daily maximum 8-h average (MDA8) exceedance threshold from 75 to 70 ppbv, poses significant challenges to U.S. air quality (AQ) forecasters responsible for ozone MDA8 forecasts. The forecasters, supplied by only a few AQ model products, end up relying heavily on self-developed tools. To help U.S. AQ forecasters, this study explores a surface ozone MDA8 forecasting tool that is based solely on statistical methods and standard meteorological variables from the numerical weather prediction (NWP) models. The model combines the self-organizing map (SOM), which is a clustering technique, with a step wise weighted quadratic regression using meteorological variables as predictors for ozone MDA8. The SOM method identifies different weather regimes, to distinguish between various modes of ozone variability, and groups them according to similarity. In this way, when a regression is developed for a specific regime, data from the other regimes are also used, with weights that are based on their similarity to this specific regime. This approach, regression in SOM (REGiS), yields a distinct model for each regime taking into account both the training cases for that regime and other similar training cases. To produce probabilistic MDA8 ozone forecasts, REGiS weighs and combines all of the developed regression models on the basis of the weather patterns predicted by an NWP model. REGiS is evaluated over the San Joaquin Valley in California and the northeastern plains of Colorado. The results suggest that the model performs best when trained and adjusted separately for an individual AQ station and its corresponding meteorological site.

Ozone variables↗
NASA NTRS1988

Computer modeling of convection in the liquid metal during solidification

Numerical methods are used to study thermally/gravitationally induced convection during solidification. The primary object is aimed at understanding the flow field in the liquid metal pool due to convection. An Al-0.5Fe weight percent alloy is assumed to fill a small scale vertical rectangular enclosure model for this study. The Alternating Direction Implicit method and Successive Overrelaxation method are employed to solve these equations, which include the energy balance equation coupled with the mass and momentum equations. Both transient and steady state conditions under 1 g and low gravity (0.001 g) are investigated. The solutions for the two-dimensional isotherms and streamlines under 1 g conditions are graphically presented. The quantitative results show that the convection which occurs in the liquid metal pool is very small and is likely to have a negligible effect on the solute redistribution ahead of the solid-liquid interface.

Lott, R. L., Jr.↗
NASA NTRS1996

Dynamic Models for LES of Turbulent Front Propagation With a Spectral Method

Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Im, H. G.↗
NASA NTRS1981

Computational methods of robust controller design for aerodynamic flutter suppression

The development of Riccati iteration, a tool for the design and analysis of linear control systems is examined. First, Riccati iteration is applied to the problem of pole placement and order reduction in two-time scale control systems. Order reduction, yielding a good approximation to the original system, is demonstrated using a 16th order linear model of a turbofan engine. Next, a numerical method for solving the Riccati equation is presented and demonstrated for a set of eighth order random examples. A literature review of robust controller design methods follows which includes a number of methods for reducing the trajectory and performance index sensitivity in linear regulators. Lastly, robust controller design for large parameter variations is discussed.

Anderson, L. R.↗
NASA NTRS2000

Modeling the Compression of Merged Compact Toroids by Multiple Plasma Jets

A fusion propulsion scheme has been proposed that makes use of the merging of a spherical distribution of plasma jets to dynamically form a gaseous liner. The gaseous liner is used to implode a magnetized target to produce the fusion reaction in a standoff manner. In this paper, the merging of the plasma jets to form the gaseous liner is investigated numerically. The Los Alamos SPHINX code, based on the smoothed particle hydrodynamics method is used to model the interaction of the jets. 2-D and 3-D simulations have been performed to study the characteristics of the resulting flow when these jets collide. The results show that the jets merge to form a plasma liner that converge radially which may be used to compress the central plasma to fusion conditions. Details of the computational model and the SPH numerical methods will be presented together with the numerical results.

Thio, Y. C. Francis↗
NASA NTRS1991

A critical comparison of two-equation turbulence models

Several two-equation models were proposed and tested against benchmark flows by various researchers. For each study, different numerical methods or codes were used to obtain the results which were reported to be an improvement over other models. However, these comparisons may be overshadowed by the different numerical schemes used to obtain the results. With this in mind, several existing two-equation turbulence models, including k-epsilon, k-tau, k-omega, and q-omega models, are implemented into a common flow solver code for near wall turbulent flows. The quality of each model is based on several criteria, including robustness and accuracy of predicting the turbulent quantities.

Lang, N. J.↗
NASA NTRS1994

Center for modeling of turbulence and transition: Research briefs, 1993

This research brief contains the progress reports of the research staff of the Center for Modeling of Turbulence and Transition (CMOTT) from June 1992 to July 1993. It is also an annual report to the Institute for Computational Mechanics in Propulsion located at Ohio Aerospace Institute and NASA Lewis Research Center. The main objectives of the research activities at CMOTT are to develop, validate, and implement turbulence and transition models for flows of interest in propulsion systems. Currently, our research covers eddy viscosity one- and two-equation models, Reynolds-stress algebraic equation models, Reynolds-stress transport equation models, nonequilibrium multiple-scale models, bypass transition models, joint scalar probability density function models, and Renormalization Group Theory and Direct Interaction Approximation methods. Some numerical simulations (LES and DNS) have also been carried out to support the development of turbulence modeling. Last year was CMOTT's third year in operation. During this period, in addition to the above mentioned research, CMOTT has also hosted the following programs: an eighteen-hour short course on 'Turbulence--Fundamentals and Computational Modeling (Part I)' given by CMOTT at the NASA Lewis Research Center; a productive summer visitor research program that has generated many encouraging results; collaborative programs with industry customers to help improve their turbulent flow calculations for propulsion system designs; a biweekly CMOTT seminar series with speakers from within and without the NASA Lewis Research Center including foreign speakers. In addition, CMOTT members have been actively involved in the national and international turbulence research activities. The current CMOTT roster and organization are listed in Appendix A. Listed in Appendix B are the abstracts of the biweekly CMOTT seminar. Appendix C lists the papers contributed by CMOTT members.

Liou, William W.↗
DOE OSTI2024

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗
NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗
NASA NTRS1992

Turbulence modeling for hypersonic flight

The objective of the present work is to develop, verify, and incorporate two equation turbulence models which account for the effect of compressibility at high speeds into a three dimensional Reynolds averaged Navier-Stokes code and to provide documented model descriptions and numerical procedures so that they can be implemented into the National Aerospace Plane (NASP) codes. A summary of accomplishments is listed: (1) Four codes have been tested and evaluated against a flat plate boundary layer flow and an external supersonic flow; (2) a code named RANS was chosen because of its speed, accuracy, and versatility; (3) the code was extended from thin boundary layer to full Navier-Stokes; (4) the K-omega two equation turbulence model has been implemented into the base code; (5) a 24 degree laminar compression corner flow has been simulated and compared to other numerical simulations; and (6) work is in progress in writing the numerical method of the base code including the turbulence model.

Bardina, Jorge E.↗
NASA NTRS1984

Development of a turbulence near-wall model and its application to separated and reattached flows

A numerical study is reported of flow and heat transfer in the separated and reattached flows created by an abrupt pipe expansion. In the study attention has been given primarily to the development of turbulence near-wall models based on the wall function. An efficient numerical method is also employed for the computation of high-velocity flows. The computed results are compared with experimental data obtained earlier. Generally, better results are obtained by employing the present near-wall models, and among them the three-layer model is superior to the two-layer one.

Amano, R. S.↗
NASA NTRSDate not supplied

Thermal Modelling of the in-Situ Consolidation of Automated Fiber Placement of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Thermoplastics↗