SEARCH · Engineering Papers

Results for “Algebraic models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Computation of incompressible, three-dimensional turbulent boundary layers and comparison with experiment

Incompressible three dimensional, turbulent boundary layer (3DTBL) experiments were simulated numerically by integrating the boundary layer equations together with an algebraic eddy viscosity turbulence model. For the flow treated, the downstream portion, where the crossflow was large, was not predicted with the present computational method; the flow was significantly influenced by elliptic flow field effects. Departures from the boundary layer concept are indicated. Calculations agreed reasonably well with the mean flow development up to separation. In one experiment the normal pressure gradients were found to be neligible in regions with large skewing and allowed testing turbulence models using the boundary layer equations. The simulation of this flow compared favorably with the experimental data throughout the flow field and suggested the applicability of algebraic eddy viscosity models for 3DTBLs.

Mueller, U. R.↗

NASA NTRS2018

kL-based Linear and Nonlinear Two-Equation Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) two-equation turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important differences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The present report documents the formulation of two model levels of turbulence models as implemented in the computational fluid dynamics FUN3D code. The levels are the two-equation linear k-kL and the two-equation algebraic Reynolds stress model (ARSM). Free shear, separated and corner flow cases are documented and compared with experimental, and other turbulence model data. The results show generally very good comparisons with experimental data. The results from this formulation are similar or better than results using the SST two-equation turbulence model. ARSM shows great promise with similar level of computational resources as basic two-equation turbulence models.

Abdol-Hamid, K. S.↗

NASA NTRS1995

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗

NASA NTRS2018

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗

DOE OSTI2025

On infinite tensor networks, complementary recovery and type II factors

We initiate a study of local operator algebras at the boundary of infinite tensor networks, using the mathematical theory of inductive limits. In particular, we consider tensor networks in which each layer acts as a quantum code with complementary recovery, a property that features prominently in the bulk-to-boundary maps intrinsic to holographic quantum error-correcting codes. In this case, we decompose the limiting Hilbert space and the algebras of observables in a way that keeps track of the entanglement in the network. As a specific example, we describe this inductive limit for the holographic Harlow-Pastawski-Preskill-Yoshida code model and relate its algebraic and error-correction features. We find that the local algebras in this model are given by the hyperfinite type II$_\infty$ factor. Next, we discuss other networks that build upon this framework and comment on a connection between type II factors and stabilizer circuits. We conclude with a discussion of multiscale entanglement renormalization ansatz networks in which complementary recovery is broken. We argue that this breaking possibly permits a limiting type III von Neumann algebra, making them more suitable ansätze for approximating subregions of quantum field theories.

holographic dualities↗

DOE OSTI2022

Likelihood Maximization and Moment Matching in Low SNR Gaussian Mixture Models

We derive an asymptotic expansion for the log-likelihood of Gaussian mixture models (GMMs) with equal covariance matrices in the low signal-to-noise regime. The expansion reveals an intimate connection between two types of algorithms for parameter estimation: the method of moments and likelihood optimizing algorithms such as Expectation-Maximization (EM). We show that likelihood optimization in the low SNR regime reduces to a sequence of least squares optimization problems that match the moments of the estimate to the ground truth moments one by one. This connection is a stepping stone towards the analysis of EM and maximum likelihood estimation in a wide range of models. A motivating application for the study of low SNR mixture models is cryo-electron microscopy data, which can be modeled as a GMM with algebraic constraints imposed on the mixture centers. We discuss the application of our expansion to algebraically constrained GMMs, among other example models of interest. © 2022 The Authors. Communications on Pure and Applied Mathematics published by Wiley Periodicals LLC.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Variational quantum and neural quantum states algorithms for the linear complementarity problem

Variational quantum algorithms (VQAs) are promising hybrid quantum-classical methods designed to leverage the computational advantages of quantum computing while mitigating the limitations of current noisy intermediate-scale quantum (NISQ) hardware. Although VQAs have been demonstrated as proofs of concept, their practical utility in solving real-world problems—and whether quantum-inspired classical algorithms can match their performance—remains an open question. We present a novel application of the variational quantum linear solver (VQLS) and its classical neural quantum states-based counterpart, the variational neural linear solver (VNLS), as key components within a minimum map Newton solver for a complementarity-based rigid-body contact model. We demonstrate using the VNLS that our solver accurately simulates the dynamics of rigid spherical bodies during collision events. These results suggest that quantum and quantum-inspired linear algebra algorithms can serve as viable alternatives to standard linear algebra solvers for modelling certain physical systems.

neural quantum states↗

DOE OSTI2024

Axion domain walls, small instantons, and non-invertible symmetry breaking

Non-invertible global symmetry often predicts degeneracy in axion potentials and carries important information about the global form of the gauge group. When these symmetries are spontaneously broken they can lead to the formation of stable axion domain wall networks which support topological degrees of freedom on their worldvolume. Such non-invertible symmetries can be broken by embedding into appropriate larger UV gauge groups where small instanton contributions lift the vacuum degeneracy, and provide a possible solution to the domain wall problem. We explain these ideas in simple illustrative examples and then apply them to the Standard Model, whose gauge algebra and matter content are consistent with several possible global structures. Each possible global structure leads to different selection rules on the axion couplings, and various UV completions of the Standard Model lead to more specific relations. As a proof of principle, we also present an example of a UV embedding of the Standard Model which can solve the axion domain wall problem. The formation and annihilation of the long-lived axion domain walls can lead to observables, such as gravitational wave signals. Observing such signals, in combination with the axion coupling measurements, can provide valuable insight into the global structure of the Standard Model, as well as its UV completion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1981

Numerical models for Halley dust environments

The algebra, parameter values, and results of a reference model for the dust environment of Comet P/Halley are presented. The model employs simplified treatments of the icy conglomerate nucleus and of the hydrodynamic outflow of gas and dust. The model includes a gas production rate derived from observations of Halley's visual magnitude in 1910, gas and dust flux nonuniformly distributed over the nuclear surface, a dust distribution function consistent with large particle inferences from comet antitails, dust density variable with mass, and an approximation to the distribution of particles near the dust envelope. Many parameter values are derived by analogy with the results for other comets. The features are incorporated in a computer program which calculates the concentration, flux, and fluence of dust particles as functions of time, position and particle mass.

Divine, N.↗

NASA NTRS1987

Three-dimensional computations of transverse hydrogen jet combustion in a supersonic airstream

A computational fluid dynamics (CFD) code is being developed to compute the mixing and combustion of hydrogen fuel in the turbulent flow fields of supersonic combustion ramjets (scramjet). The code solves the three-dimensional Reynolds time-averaged complete Navier-Stokes equations including transport equations for a four species, two reaction, global finite rate chemistry model. The code was applied to the case of transverse injection of hydrogen from a sonic circular orifice into a supersonic airstream. The equations were numerically integrated using MacCormack's explicit method, and the algebraic eddy viscosity model of Baldwin-Lomax was used to model the turbulence. In the species transport and energy equations, diffusion coefficients based on Fick's Law and an assumption of unit Lewis number were applied. Computed features of the three-dimensional flow field are depicted by static pressure, static temperature, mass fraction of species, and velocity vectors. For engineering interest, mixing and combustion parameters were examined to assess the effect of injector diameter, injected fuel pressure, fuel-air ratio, and spacing of fuel injectors. The objective of the present paper is to demonstrate the capability of the present three-dimensional spatially elliptic, CFD code for turbulent, reacting flow. Application of the code to specific supersonic combustion configurations is planned.

Uenishi, K.↗

DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗

NASA NTRS1991

Center for Modeling of Turbulence and Transition (CMOTT). Research briefs: 1990

Brief progress reports of the Center for Modeling of Turbulence and Transition (CMOTT) research staff from May 1990 to May 1991 are given. The objectives of the CMOTT are to develop, validate, and implement the models for turbulence and boundary layer transition in the practical engineering flows. The flows of interest are three dimensional, incompressible, and compressible flows with chemistry. The schemes being studied include the two-equation and algebraic Reynolds stress models, the full Reynolds stress (or second moment closure) models, the probability density function models, the Renormalization Group Theory (RNG) and Interaction Approximation (DIA), the Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS).

Povinelli, Louis A.↗

NASA NTRS1984

Ion thruster performance model

A model of ion thruster performance is developed for high flux density, cusped magnetic field thruster designs. This model is formulated in terms of the average energy required to produce an ion in the discharge chamber plasma and the fraction of these ions that are extracted to form the beam. The direct loss of high energy (primary) electrons from the plasma to the anode is shown to have a major effect on thruster performance. The model provides simple algebraic equations enabling one to calculate the beam ion energy cost, the average discharge chamber plasma ion energy cost, the primary electron density, the primary-to-Maxwellian electron density ratio and the Maxwellian electron temperature. Experiments indicate that the model correctly predicts the variation in plasma ion energy cost for changes in propellant gas (Ar, Kr and Xe), grid transparency to neutral atoms, beam extraction area, discharge voltage, and discharge chamber wall temperature. The model and experiments indicate that thruster performance may be described in terms of only four thruster configuration dependent parameters and two operating parameters. The model also suggests that improved performance should be exhibited by thruster designs which extract a large fraction of the ions produced in the discharge chamber, which have good primary electron and neutral atom containment and which operate at high propellant flow rates.

Brophy, J. R.↗

NASA NTRS1991

An algebra of discrete event processes

This report deals with an algebraic framework for modeling and control of discrete event processes. The report consists of two parts. The first part is introductory, and consists of a tutorial survey of the theory of concurrency in the spirit of Hoare's CSP, and an examination of the suitability of such an algebraic framework for dealing with various aspects of discrete event control. To this end a new concurrency operator is introduced and it is shown how the resulting framework can be applied. It is further shown that a suitable theory that deals with the new concurrency operator must be developed. In the second part of the report the formal algebra of discrete event control is developed. At the present time the second part of the report is still an incomplete and occasionally tentative working paper.

Heymann, Michael↗

NASA NTRS2009

Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation

This report gives an overview of physical modeling enhancements to the Wind-US flow solver which were made to improve the capabilities for simulation of hypersonic flows and the reliability of computations to complement hypersonic testing. The improvements include advanced turbulence models, a bypass transition model, a conjugate (or closely coupled to vehicle structure) conduction-convection heat transfer capability, and an upgraded high-speed combustion solver. A Mach 5 shock-wave boundary layer interaction problem is used to investigate the benefits of k- s and k-w based explicit algebraic stress turbulence models relative to linear two-equation models. The bypass transition model is validated using data from experiments for incompressible boundary layers and a Mach 7.9 cone flow. The conjugate heat transfer method is validated for a test case involving reacting H2-O2 rocket exhaust over cooled calorimeter panels. A dual-mode scramjet configuration is investigated using both a simplified 1-step kinetics mechanism and an 8-step mechanism. Additionally, variations in the turbulent Prandtl and Schmidt numbers are considered for this scramjet configuration.

Georgiadis, Nicholas J.↗

NASA NTRS1993

Numerical computation of aerodynamics and heat transfer in a turbine cascade and a turn-around duct using advanced turbulence models

The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.

Lakshminarayana, B.↗

NASA NTRS2019

On the Use of Optimization Techniques for Turbulence Model Calibration

This article discusses the use of numerical optimization procedures to aid in the calibration of turbulence model coefficients. Such methods would increase the rigor and repeatability of the calibration procedure by requiring clearly defined and objective optimization metrics, and could be used to identify unique combinations of coefficient values for specific flow problems. The approach is applied to the re-calibration of an explicit algebraic Reynolds stress model for the incompressible planar mixing layer using the Nelder-Mead simplex algorithm and a micro-genetic algorithm with minimally imposed constraints. Three composite fitness functions, each based upon the error in the mixing layer growth rate and the normal and shear components of the Reynolds stresses, are investigated. The results demonstrate a significant improvement in the target objectives through the adjustment of three pressure-strain coefficients. Adjustments of additional coefficients provide little further benefit. Issues regarding the effectiveness of the fitness functions and the efficiency of the optimization algorithms are also discussed.

Yoder, Dennis A.↗

NASA NTRS2020

Juncture Flow Computations using kL-Based Turbulence Models

The development and implementation of kL-based Reynolds-Average Navier-Stokes (RANS) two-equation turbulence models are reported in this paper. The kL model is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta showed that a reliable transport equation can be formed from the turbulence length scale L, and the turbulence kinetic energy k. Rotta's kL equation is well suited for a term-by-term modeling and displays useful features compared to other scale formulations. One of the important differences is the inclusion of higher-order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The k-kL scheme has been modified to include the Algebraic Reynolds Stress Model (ARSM) and the Quadratic Constitutive Relation (QCR) as nonlinear models. These models have gone through extensive validations using two-dimensional benchmark problems. The present study documents the application of the k-kL schemes to simulate flow around the Juncture Flow Model (JFM). The k-kL prediction results show generally good comparisons with measurements. The results from this formulation are similar to, or better than results using the SA-RC-QCR2000 two-equation turbulence model. The ARSM and the QCR formulations of k-kL show promise with a similar level of computational resources as basic two-equation turbulence models.

K S Abdol-Hamid↗