SEARCH · Search NASA

Results for “Methods: numerical”

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 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Numerical Simulation of Unsteady Aerodynamic Models

This report documents the results of the numerical simulations of unsteady aerodynamic models. The results focus on numerical accuracy and efficiency, and the robustness of the numerical methods. The aerodynamic models includes the classical Wagner and Kussner functions and the Leishman-Beddoes dynamic stall model. The simulations includes the numerical approximations of the Duhamel's integrals using both indicial (step) and impulse responses, the numerical integrations of the state-space models, and the exact solutions. The report also presents the conversion among different model representations.

Nguyen, Khanh Q.↗

NASA NTRS2001

Development of a Hybrid RANS/LES Method for Turbulent Mixing Layers

Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Georgiadis, Nicholas J.↗

DOE OSTI2024

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗

DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗

NASA NTRS1987

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.↗

NASA NTRS2009

The Capillary Flow Experiments Aboard the International Space Station: Increments 9-15

This report provides a summary of the experimental, analytical, and numerical results of the Capillary Flow Experiment (CFE) performed aboard the International Space Station (ISS). The experiments were conducted in space beginning with Increment 9 through Increment 16, beginning August 2004 and ending December 2007. Both primary and extra science experiments were conducted during 19 operations performed by 7 astronauts including: M. Fincke, W. McArthur, J. Williams, S. Williams, M. Lopez-Alegria, C. Anderson, and P. Whitson. CFE consists of 6 approximately 1 to 2 kg handheld experiment units designed to investigate a selection of capillary phenomena of fundamental and applied importance, such as large length scale contact line dynamics (CFE-Contact Line), critical wetting in discontinuous structures (CFE-Vane Gap), and capillary flows and passive phase separations in complex containers (CFE-Interior Corner Flow). Highly quantitative video from the simply performed flight experiments provide data helpful in benchmarking numerical methods, confirming theoretical models, and guiding new model development. In an extensive executive summary, a brief history of the experiment is reviewed before introducing the science investigated. A selection of experimental results and comparisons with both analytic and numerical predictions is given. The subsequent chapters provide additional details of the experimental and analytical methods developed and employed. These include current presentations of the state of the data reduction which we anticipate will continue throughout the year and culminate in several more publications. An extensive appendix is used to provide support material such as an experiment history, dissemination items to date (CFE publication, etc.), detailed design drawings, and crew procedures. Despite the simple nature of the experiments and procedures, many of the experimental results may be practically employed to enhance the design of spacecraft engineering systems involving capillary interface dynamics.

Jenson, Ryan M.↗

NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗

NASA NTRS1991

Accretion and evolution of solar system bodies

We use a combination of analytical and numerical methods to study dynamical processes involved in the formation of planets and smaller bodies in the solar system. Our goal was to identify and understand critical processes and to link them in a numerical model of planetesimal accretion. We study effects of these processes by applying them in the context of the standard model of solar system formation, which involves accretion of the terrestrial planets and cores of the giant planet from small planetesimals. The principal focus of our research effort is the numerical simulation of accretion of a swarm of planetesimals into bodies of planetary size. Our computer code uses a Monte Carlo method to determine collisional interactions within the swarm. These interactions are not determined simply by a relative velocity, but rather by explicit distributions of keplerian orbital elements. The planetesimal swarm is divided into a number of zones in semimajor axis, which are allowed to interact. The present version of our code has the capability of following detailed distributions of size, eccentricity, and inclination in each zone.

Weidenschilling, S. J.↗

DOE OSTI2020

i- flow: High-dimensional integration and sampling with normalizing flows

In many fields of science, high-dimensional integration is required. Numerical methods have been developed to evaluate these complex integrals. We introduce the code i-flow, a python package that performs high-dimensional numerical integration utilizing normalizing flows. Normalizing flows are machine-learned, bijective mappings between two distributions. i-flow can also be used to sample random points according to complicated distributions in high dimensions. We compare i-flow to other algorithms for high-dimensional numerical integration and show that i-flow outperforms them for high dimensional correlated integrals. The i-flow code is publicly available on gitlab at https://gitlab.com/i-flow/i-flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Calculating fluctuations and self-correlations numerically for causal charge diffusion in relativistic heavy-ion collisions

Here, we study the propagation and diffusion of electric charge fluctuations in the Bjorken hydrodynamic model with both white and Catteneo noise using purely numerical methods. We show that a lattice of noise fluctuations is required to fully calculate the two-point correlators of charge. We introduce a numerical procedure to solve the stochastic differential equations that arise from the charge conservation equation on the lattice event by event. We explicitly identify the self-correlation term in the case of Catteneo noise and provide a physical interpretation. We provide a numerical recipe to remove this contribution from the full two-point correlators. Finally, we calculate the balance functions for charged hadrons. By limiting the speed of signal propagation, we observe the expected narrowing of the balance functions after removing the self-correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1975

New methods for the numerical integration of ordinary differential equations and their application to the equations of motion of spacecraft

A number of hybrid methods for solving Cauchy problems are described on the basis of an evaluation of advantages of single and multiple-point numerical integration methods. The selection criterion is the principle of minimizing computer time. The methods discussed include the Nordsieck method, the Bulirsch-Stoer extrapolation method, and the method of recursive Taylor-Steffensen power series.

Banyukevich, A.↗

NASA NTRS1987

Upwind and symmetric shock-capturing schemes

The development of numerical methods for hyperbolic conservation laws has been a rapidly growing area for the last ten years. Many of the fundamental concepts and state-of-the-art developments can only be found in meeting proceedings or internal reports. This review paper attempts to give an overview and a unified formulation of a class of shock-capturing methods. Special emphasis is on the construction of the basic nonlinear scalar second-order schemes and the methods of extending these nonlinear scalar schemes to nonlinear systems via the extact Riemann solver, approximate Riemann solvers, and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows is discussed. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas dynamics problems.

Yee, H. C.↗

NASA NTRS1985

Implementation of design sensitivity analysis with existing finite element codes

A numerical method is presented to implement structural design sensitivity analysis theory, using the versatility and convenience of existing finite element structural analysis programs. Design variables such as thickness and cross-sectional areas of components of individual members and built-up structures are considered. Structural performance functionals considered include displacement and stress. The method is also applicable for eigenvalue problem design sensitivity analysis. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. Thus, design sensitivity analysis software does not have to be embedded in an existing finite element code. Feasability of the method is shown through analysis of several problems, including a built-up structure. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of finite difference perturbations.

Choi, Kyung K.↗

NASA NTRS1990

Collisionless drift-tearing modes in the magnetopause

The linear stability properties of collisionless drift-tearing modes are analyzed in a modified Harris equilibrium model of the magnetopause. Particular attention is paid to the relevance of the parametric behavior of growth rates to the 'magnetic percolation' theory of flux transfer event formation (Galeev et al., 1986). Numerical methods are used to solve the drift-tearing eigenmode equations and the results are compared with those previously obtained by analytical methods. The analytical results are found to correctly model important parametric dependencies but to typically overestimate the rate of growth. The eigenmode equations are numerically difficult, and an integration scheme utilizing Ricatti transforms is developed to affect their solution.

Gladd, N. T.↗

NASA NTRS2000

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Convection-Diffusion Problems

In the space-time conservation element and solution element (CE/SE) method, the independent marching variables used comprise not only the mesh value of the physical dependent variables but also, in contrast to it typical numerical method, the Mesh values of the spatial derivatives of the physical variables The use of the extra marching variables results from the need to construct the two-level explicit and nondissipative schemes which are at the core of the CE/SE development. It also results from the need to minimize the stencil while maintaining accuracy. In this paper using the 1D(sub (alpha)-mu) scheme as an example, the effect of this added complication on consistency, accuracy and operation count is assessed. As part of this effort, an equivalent yet more efficient form of the alpha-mu scheme in which the independent marching variables are the local fluxes tied to each mesh point is introduced. Also, the intriguing relations that exist among the alpha-mu. Leapfrog, and DuFort-Frankel schemes are further explored. In addition, the redundance of the Leapfrog, DUFort-Frankel, and Lax scheme and the remedy for this redundance are discussed. This paper is concluded with the construction and evaluation of a CE/SE solver for the inviscid Burger equation.

Chang, Sin-Chung↗

NASA NTRS2002

Improved Multi-Axial, Temperature and Time Dependent (MATT) Failure Model

An extensive effort has recently been completed by the Space Shuttle's Reusable Solid Rocket Motor (RSRM) nozzle program to completely characterize the effects of multi-axial loading, temperature and time on the failure characteristics of three filled epoxy adhesives (TIGA 321, EA913NA, EA946). As part of this effort, a single general failure criterion was developed that accounted for these effects simultaneously. This model was named the Multi- Axial, Temperature, and Time Dependent or MATT failure criterion. Due to the intricate nature of the failure criterion, some parameters were required to be calculated using complex equations or numerical methods. This paper documents some simple but accurate modifications to the failure criterion to allow for calculations of failure conditions without complex equations or numerical techniques.

Richardson, D. E.↗

NASA NTRS1995

Linear Wave Motion from Concentrated Harmonic Sources in Blasius Flow

The motion of individual linear instability waves in shear flows is well described by existing theoretical and numerical methods. However, naturally occuring sources produce coherent wave motions with broadband spanwise wavenumber and frequency spectra, and the different spectral components interact both linearly and nonlinearily. This paper describes a series of calculations for the parameters of three different experiments using locally parallel linear stability theory (LST), the Parabolized Stability Equations (PSE), and Direct Numerical Simulation (DNS). The calculations illustrate the strengths and weaknesses of the different methods, the extent to which the methods agree or disagree, and, finally the extent to which agreement with the measurements can be attained, given that the experiments also have their own difficulties.

linear↗