Search NASA⌕ Search

SEARCH · Search NASA

Results for “Implicit 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Electrification Futures Study: Methodological Approaches for Assessing Long-Term Power System Impacts of End-Use Electrification

By its nature, electrification enhances the coupling between the electric sector and end-use sectors. Assessing the impacts of electrification requires both an examination of the complex interactions between sectors and a broader assessment of multiple parts of the energy system. The Electrification Futures Study (EFS) uses several complementary modeling and analysis tools to analyze the impacts of electrification on the U.S. energy system. In particular, the EFS relies on an overarching scenario analysis approach, but through the use of separate modeling approaches designed to assess various electricity demand- and supply-side futures. The primary model employed to generate the supply-side scenarios is the Regional Energy Deployment System (ReEDS) model, which is a capacity expansion model for the U.S. electricity system through 2050. Traditionally, the model has been primarily exercised in scenario analysis that implicitly assumed limited electrification. Because of this assumption, resulting future load profiles are approximated by historical ones, load growth is driven primarily by population and economic growth only, and changes in direct end-use natural gas consumption do not effect natural gas costs for electricity generation. In this report, we (1) reflect the potential for resource sharing between regions given these changes in demand, (2) represent how changes in natural gas consumption in end-use sectors could impact the economics of natural gas-fired generation, and (3) document a new model representation of demand-side flexibility used for the EFS. These improvements to ReEDS are employed for the EFS supply-wide analysis, which is summarized in a companion EFS report titled Electrification Futures Study: Scenarios of Power System Evolution and Infrastructure Development for the United States (Murphy et al. 2019). The data and methods documented in this report could also be adapted for other models with similar scope and limitations as ReEDS, and these data and methods could be used to assess future electric system scenarios.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1989

Supersonic flow computations by two-equation turbulence modeling

In the present work a solver for the Reynolds averaged compressible Navier-Stokes equations, to compute high speed turbulent flows characterized by interacting shock waves and viscous layers, is presented. A k-epsilon turbulence model that accounts for compressibility effects is developed. The numerical algorithm is based on a finite volume multistage Runge Kutta technique that is explicit for the solution of the mean flow variables, and implicit for the solution of the k-epsilon equations. The model is validated by extensive comparison with experimental results of flows over compression ramps characterized by interacting shock waves/boundary layers.

Grasso, Francesco↗
DOE OSTI2022

The role of surface forces in environment-enhanced cracking of brittle solids

Fracture initiation and propagation in brittle materials is promoted in surface-reactive (sorptive) environments, a phenomenon known as subcritical crack growth (SCG). Laboratory measured crack-propagation velocity vs. stress intensity factor relationships typically exhibit highly nonlinear, multi-stage characteristics that are sensitive to environmental factors such as adsorbate concentration and temperature. For practical purposes, empirical relationships (e.g., a power law) have been used to describe this complex phenomenon. However, how the overall SCG behavior emerges from the underlying fundamental processes near the crack tip, such as the interaction of the crack surfaces separated by only a few nanometers and mass transport within the nano-confined space, is still not well understood. This paper develops a mechanistic, surface-force-based fracture theory (SFFT) which integrates surface force models, fluid transport models, and linear elastic fracture mechanics to quantitatively explain the multi-stage characteristics of SCG in brittle solids. A numerical model is developed based on SFFT and solved through an implicit partitioned scheme for efficiency and modularity. The results are validated by Wiederhorn's data on crack propagation in soda-lime glasses at a wide range of relative humidity levels. We show that, for the first time, the entire range of an SCG curve can be captured by a single physics-based model. The predicted SCG curves reveal that the development of repulsive disjoining pressure behind the crack tip can be responsible for the reduced apparent fracture toughness in a sorptive environment. In conclusion, the shape of the SCG curve, and its changes with respect to the environment, is found to critically depend on the assumed transport models.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Exponential integration algorithms applied to viscoplasticity

Four, linear, exponential, integration algorithms (two implicit, one explicit, and one predictor/corrector) are applied to a viscoplastic model to assess their capabilities. Viscoplasticity comprises a system of coupled, nonlinear, stiff, first order, ordinary differential equations which are a challenge to integrate by any means. Two of the algorithms (the predictor/corrector and one of the implicits) give outstanding results, even for very large time steps.

Freed, Alan D.↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗
NASA NTRS1998

A Conceptual Framework for Predicting Error in Complex Human-Machine Environments

We present a Goals, Operators, Methods, and Selection Rules-Model Human Processor (GOMS-MHP) style model-based approach to the problem of predicting human habit capture errors. Habit captures occur when the model fails to allocate limited cognitive resources to retrieve task-relevant information from memory. Lacking the unretrieved information, decision mechanisms act in accordance with implicit default assumptions, resulting in error when relied upon assumptions prove incorrect. The model helps interface designers identify situations in which such failures are especially likely.

Freed, Michael↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗
NASA NTRS1984

The calculation of transonic rotor noise

It is pointed out that an accurate prediction of the high-speed impulsive rotor noise is very difficult since the noise field depends on many complex factors. Schmitz and Yu (1981) have conducted a study with the aim to include local aerodynamic nonlinearities in a calculation regarding the acoustic shock. In the present investigation, the hover case of the model rotor studied by Schmitz and Yu was recomputed without certain implicit assumptions made in the earlier calculation. The obtained results show that the Schmitz-Yu computations are very much base-support area dependent. Attention is given to the quadrupole integral, and the pressure-time curves for the model rotor in hover.

Aggarwal, H. R.↗
DOE OSTI2020

A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons

Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

Turbulence modeling in supersonic combusting flows

To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.

Chitsomboon, Tawit↗
NASA NTRS1991

Three-dimensional calculations of supersonic reacting flows using an LU scheme

A 3-D numerical program that incorporates comprehensive real gas property models was developed to simulate supersonic reacting flows. The code employs an implicit, finite volume, Lower-Upper (LU), time-marching method to solve the complete Navier-Stokes and species equations in a fully-coupled and efficient manner. A chemistry model with 9 species and 18 reaction steps is adopted in the program to represent the chemical reactions of H2 and air. To demonstrate the capability of the program, flow fields of underexpanded hydrogen jets transversely injected into the supersonic airstream inside the combustors of scramjets are calculated. Results clearly depict the flow characteristics, including the shock structure, the separated flow regions around the injector, and the distribution of the combustion products.

Yu, Sheng-Tao↗
NASA NTRS1992

Implementation of a kappa-epsilon turbulence model to RPLUS3D code

The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.

Chitsomboon, Tawit↗
NASA NTRSDate not supplied

Introduction to Quantum Computing

Quantum computing offers the potential to revolutionize high-performance computing by providing a means to solve certain computational problems asymptotically faster than any classical computer. Quantum computing has advanced recently from merely a theoretical possibility to engineered reality, including commercial entities offering early prototype quantum processors, both special-purpose quantum annealers and general-purpose gate-model processors. The media have been showcasing each new development and implicitly conveying the message that quantum-computing ubiquity is nigh. Here, we will respond to this hype and provide an overview of the exciting but still early state of the field. In this tutorial, we introduce participants to the computational models that give quantum computing its immense computational power. We examine the thought processes that programmers need to map problems to quantum computers. And we discuss hardware and algorithmic challenges that must be overcome before quantum computing becomes a component of every software developer's repertoire.

Quantum computing↗
NASA NTRSDate not supplied

Introduction to Quantum Computing

Quantum computing offers the potential to revolutionize high-performance computing by providing a means to solve certain computational problems asymptotically faster than any classical computer. Quantum computing has advanced recently from merely a theoretical possibility to engineered reality, including commercial entities offering early prototype quantum processors, both special-purpose quantum annealers and general-purpose gate-model processors. The media have been showcasing each new development and implicitly conveying the message that quantum-computing ubiquity is nigh. Here, we will respond to this hype and provide an overview of the exciting but still early state of the field. In this tutorial, we introduce participants to the computational models that give quantum computing its immense computational power. We examine the thought processes that programmers need to map problems to quantum computers. And we discuss hardware and algorithmic challenges that must be overcome before quantum computing becomes a component of every software developer's repertoire. (Update of 2022 slides)

Quantum computing↗
NASA NTRS1989

Three-dimensional calculation of supersonic reacting flows using an LU scheme

A new three-dimensional numerical program incorporated with comprehensive real gas property models was developed to simulate supersonic reacting flows. The code employs an implicit finite volume, Lower-Upper (LU) time-marching method to solve the complete Navier-Stokes and species equations in a fully-coupled and very efficient manner. A chemistry model with nine species and eighteen reaction steps are adopted in the program to represent the chemical reaction of H2 and air. To demonstrate the capability of the program, flow fields of underexpanded hydrogen jets transversely injected into supersonic air stream inside the combustors of scramjets are calculated. Results clearly depict the flow characteristics, including the shock structure, separated flow regions around the injector, and the distribution of the combustion products.

Yu, Sheng-Tao↗
NASA NTRS1994

Computation of Separated and Unsteady Flows with One- and Two-Equation Turbulence Models

The ability of one- and two-equation turbulence models to predict unsteady separated flows over airfoils is evaluated. An implicit, factorized, upwind-biased numerical scheme is used for the integration of the compressible, Reynolds averaged Navier-Stokes equations. The turbulent eddy viscosity is obtained from the computed mean flowfield by integration of the turbulent field equations. The two-equation turbulence models are discretized in space with an upwind-biased, second order accurate total variation diminishing scheme. One and two-equation turbulence models are first tested for a separated airfoil flow at fixed angle of incidence. The same models are then applied to compute the unsteady flowfields about airfoils undergoing oscillatory motion at low subsonic Mach numbers. Experimental cases where the flow has been tripped at the leading edge and where natural transition was allowed to occur naturally are considered. The more recently developed field-equation turbulence models capture the physics of unsteady separated flow significantly better than the standard kappa-epsilon and kappa-omega models. However, certain differences in the hysteresis effects are obtained. For an untripped high-Reynolds-number flow, it was found necessary to take into account the leading edge transitional flow region in order to capture the correct physical mechanism that leads to dynamic stall.

Ekaterinaris, John A.↗
NASA NTRS1991

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.↗