Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

General Relativistic Implicit Monte Carlo Radiation-hydrodynamics

Abstract We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ : an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation–matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1975

Calculation of arbitrary-order diffraction efficiencies of thick gratings with arbitrary grating shape

A method for calculating arbitrary-order diffraction efficiencies of thick, lossless transmission gratings with arbitrary periodic grating shapes has been developed. A Fourier-series representation of the grating is employed, along with a coupled-mode theory of diffraction. For illustration, numerical values of the diffraction efficiencies at the first three Bragg angles are calculated for sinusoidal, square-wave, triangular, and saw-tooth gratings. Numerical results for the same grating shapes with the same parameters are also calculated for comparison, by extending Burckhardt's numerical method for analyzing thick sinusoidal gratings. The comparison shows that the coupled-mode theory provides results with relative computational ease and results that are in agreement with calculations obtained by extending the more-rigorous Burckhardt theory to nonsinusoidal grating shapes and to higher-order Bragg angles.

Su, S. F.↗
NASA NTRS1981

Analysis of starvation effects on hydrodynamic lubrication in nonconforming contacts

The effects of lubricant starvation on minimum film thickness, under conditions of a hydrodynamic point contact, are determined by numerical methods where: (1) starvation is effected by varying the fluid inlet level; (2) the Reynolds boundary conditions are applied at the cavitation boundary; and (3) zero pressure is stipulated at the meniscus or inlet boundary. Seventy-four cases were used to numerically determine a minimum-film-thickness equation, as a function of the ratio of dimensionless load to dimensionless speed for varying degrees of starvation. A film reduction factor was in turn determined as a function of the fluid inlet level, and a starved, fully-flooded boundary was defined along with an expression determining the onset of starvation. It is found that as the degree of starvation increases, the minimum film thickness decreases gradually until the fluid inlet becomes critical. The changes in the inlet pressure buildup due to changing the available lubricant supply are presented in the form of three-dimensional isometric plots.

Brewe, D. E.↗
NASA NTRS1989

The role of topography in geodetic gravity field modelling

Masses associated with the topography, bathymetry, and its isostatic compensation are a dominant source of gravity field variations, especially at shorter wavelengths. On global scales the topographic/isostatic effects are also significant, except for the lowest harmonics. In practice, though, global effects need not be taken into account as such effects are included in the coefficients of the geopotential reference fields. On local scales, the short-wavelength gravity variations due to the topography may, in rugged terrain, be an order of magnitude larger than other effects. In such cases, explicit or implicit terrain reduction procedures are mandatory in order to obtain good prediction results. Such effects may be computed by space-domain integration or by fast Fourier transformation (FFT) methods. Numerical examples are given for areas of the Canadian Rockies. In principle, good knowledge of the topographic densities is required to produce the smoothest residual field. Densities may be determined from sample measurements or by gravimetric means, but both are somewhat troublesome methods in practice. The use of a standard density, e.g., 2.67 g/cu cm, may often yield satisfactory results and may be put within a consistent theoretical framework. The independence of density assumptions is the key point of the classical Molodensky approach to the geodetic boundary value problem. The Molodensky solutions take into account that land gravity field observations are done on a non-level surface. Molodensky's problem may be solved by integral expansions or more effective FFT methods, but the solution should not be intermixed with the use of terrain reductions. The methods are actually complimentary and may both be required in order to obtain the smoothest possible signal, least prone to aliasing and other effects coming from sparse data coverage, typical of rugged topography.

Forsberg, R.↗
DOE OSTI2017

Efficient and robust relaxation procedures for multi-component mixtures including phase transition

We consider a thermodynamic consistent multi-component model in multi-dimensions that is a generalization of the classical two-phase flow model of Baer and Nunziato. The exchange of mass, momentum and energy between the phases is described by additional source terms. Typically these terms are handled by relaxation procedures. Available relaxation procedures suffer from efficiency and robustness resulting in very costly computations that in general only allow for one-dimensional computations. Therefore we focus on the development of new efficient and robust numerical methods for relaxation processes. We derive exact procedures to determine mechanical and thermal equilibrium states. Further we introduce a novel iterative method to treat the mass transfer for a three component mixture. All new procedures can be extended to an arbitrary number of inert ideal gases. We prove existence, uniqueness and physical admissibility of the resulting states and convergence of our new procedures. Efficiency and robustness of the procedures are verified by means of numerical computations in one and two space dimensions. - Highlights: • We develop novel relaxation procedures for a generalized, thermodynamically consistent Baer–Nunziato type model. • Exact procedures for mechanical and thermal relaxation procedures avoid artificial parameters. • Existence, uniqueness and physical admissibility of the equilibrium states are proven for special mixtures. • A novel iterative method for mass transfer is introduced for a three component mixture providing a unique and admissible equilibrium state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Level 2 Milestone: Develop and Incorporate Novel Code Verification of Fundamental Equations in Gemma and Set Up Appropriate Tests

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementations of the electric-, magnetic-, and combinedfield integral equations, are frequently used to solve Maxwell’s equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this report, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches in Gemma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

The anisotropic radiative transfer problem in optically thick, strongly magnetized plasma - A comparison of results

The accuracy of four different methods for solving the coherent radiative transfer problem in a strongly magnetized plasma is compared. Specific attention is given to the case of a semi-infinite homogeneous plasma with a magnetic field perpendicular to the surface. The four approaches evaluated were: diffusion approximation schemes; the generalized H-function method of Silant'ev (1982) the Feautrier method of Nagel (1981) and the integral equation approach of Meszares and Bonazzola (1981). It is found that the coupled diffusion approximation schemes of Nagel (1981) and Kaminker et al. (1982) provided satisfactory results and were more efficient in many cases than direct numerical methods.

Pavlov, G. G.↗
NASA NTRS2016

Development of a High-Order Space-Time Matrix-Free Adjoint Solver

The growth in computational power and algorithm development in the past few decades has granted the science and engineering community the ability to simulate flows over complex geometries, thus making Computational Fluid Dynamics (CFD) tools indispensable in analysis and design. Currently, one of the pacing items limiting the utility of CFD for general problems is the prediction of unsteady turbulent ows.1{3 Reynolds-averaged Navier-Stokes (RANS) methods, which predict a time-invariant mean flowfield, struggle to provide consistent predictions when encountering even mild separation, such as the side-of-body separation at a wing-body junction. NASA's Transformative Tools and Technologies project is developing both numerical methods and physical modeling approaches to improve the prediction of separated flows. A major focus of this e ort is efficient methods for resolving the unsteady fluctuations occurring in these flows to provide valuable engineering data of the time-accurate flow field for buffet analysis, vortex shedding, etc. This approach encompasses unsteady RANS (URANS), large-eddy simulations (LES), and hybrid LES-RANS approaches such as Detached Eddy Simulations (DES). These unsteady approaches are inherently more expensive than traditional engineering RANS approaches, hence every e ort to mitigate this cost must be leveraged. Arguably, the most cost-effective approach to improve the efficiency of unsteady methods is the optimal placement of the spatial and temporal degrees of freedom (DOF) using solution-adaptive methods.

Adjoint↗
DOE OSTI2015

On the numerical treatment of dissipative particle dynamics and related systems

We review and compare numerical methods that simultaneously control temperature while preserving the momentum, a family of particle simulation methods commonly used for the modelling of complex fluids and polymers. The class of methods considered includes dissipative particle dynamics (DPD) as well as extended stochastic-dynamics models incorporating a generalized pairwise thermostat scheme in which stochastic forces are eliminated and the coefficient of dissipation is treated as an additional auxiliary variable subject to a feedback (kinetic energy) control mechanism. In the latter case, we consider the addition of a coupling of the auxiliary variable, as in the Nosé–Hoover–Langevin (NHL) method, with stochastic dynamics to ensure ergodicity, and find that the convergence of ensemble averages is substantially improved. To this end, splitting methods are developed and studied in terms of their thermodynamic accuracy, two-point correlation functions, and convergence. In terms of computational efficiency as measured by the ratio of thermodynamic accuracy to CPU time, we report significant advantages in simulation for the pairwise NHL method compared to popular alternative schemes (up to an 80% improvement), without degradation of convergence rate. The momentum-conserving thermostat technique described here provides a consistent hydrodynamic model in the low-friction regime, but it will also be of use in both equilibrium and nonequilibrium molecular simulation applications owing to its efficiency and simple numerical implementation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Analysis of acoustic radiation in a jet flow environment

The analysis of sound fields from arbitrary source distributions in terms of Legendre and spherical Hankel functions is well known. The purpose of this paper is to extend this classical method of analysis to environments such as jet flows where flow and flow gradients are inherently present. The wave-equation governing the radiation of sound in such an environment is derived. The steady state flow and flow gradients in the axial and transverse directions appear as coefficients in the terms of the wave-equation. A semi-numerical method is used to solve the wave-equation in terms of modified spherical harmonics yielding the phase velocities and the directivities of an infinite set of modes. The directivity of each mode is obtained in terms of modified Legendre functions by numerical integration. Some results of these directivity and phase-velocity calculations are presented for a limited number of frequency and flow parameters. Both convective and shear refraction are shown to be important.

Mungur, P.↗
NASA NTRS1987

Reduction of the radar cross section of arbitrarily shaped cavity structures

The problem of the reduction of the radar cross section (RCS) of open-ended cavities was studied. The issues investigated were reduction through lossy coating materials on the inner cavity wall and reduction through shaping of the cavity. A method was presented to calculate the RCS of any arbitrarily shaped structure in order to study the shaping problem. The limitations of this method were also addressed. The modal attenuation was studied in a multilayered coated waveguide. It was shown that by employing two layers of coating, it was possible to achieve an increase in both the magnitude of attenuation and the frequency band of effectiveness. The numerical method used in finding the roots of the characteristic equation breaks down when the coating thickness is very lossy and large in terms of wavelength. A new method of computing the RCS of an arbitrary cavity was applied to study the effects of longitudinal bending on RCS reduction. The ray and modal descriptions for the fields in a parallel plate waveguide were compared. To extend the range of validity of the Shooting and Bouncing Ray (SBR) method, the simple ray picture must be modified to account for the beam blurring.

Chou, R.↗
NASA NTRS1992

Flowfield of a lifting rotor in hover - A Navier-Stokes simulation

The viscous, three-dimensional flowfield of a lifting helicopter rotor in hover is calculated by using an upwind, implicit, finite-difference numerical method for solving the thin layer Navier-Stokes equations. The induced effects of the wake, including the interaction of tip vortices with successive blades, are calculated as part off the overall flowfield solution without using any ad hoc wake models. Comparison of the numerical results for the subsonic and transonic conditions show good agreement with the experimental data and with the previously published Navier-Stokes calculations using a simple wake model. Some comparisons with Euler calculations are also presented, along with some discussions of the grid refinement studies.

Srinivasan, G. R.↗
NASA NTRS2001

Left ventricular endocardial surface detection based on real-time 3D echocardiographic data

OBJECTIVE: A new computerized semi-automatic method for left ventricular (LV) chamber segmentation is presented. METHODS: The LV is imaged by real-time three-dimensional echocardiography (RT3DE). The surface detection model, based on level set techniques, is applied to RT3DE data for image analysis. The modified level set partial differential equation we use is solved by applying numerical methods for conservation laws. The initial conditions are manually established on some slices of the entire volume. The solution obtained for each slice is a contour line corresponding with the boundary between LV cavity and LV endocardium. RESULTS: The mathematical model has been applied to sequences of frames of human hearts (volume range: 34-109 ml) imaged by 2D and reconstructed off-line and RT3DE data. Volume estimation obtained by this new semi-automatic method shows an excellent correlation with those obtained by manual tracing (r = 0.992). Dynamic change of LV volume during the cardiac cycle is also obtained. CONCLUSION: The volume estimation method is accurate; edge based segmentation, image completion and volume reconstruction can be accomplished. The visualization technique also allows to navigate into the reconstructed volume and to display any section of the volume.

Non-NASA Center↗
DOE OSTI2013

Proceedings of the 2013 International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering - M and C 2013

The Mathematics and Computation Division of the American Nuclear (ANS) and the Idaho Section of the ANS hosted the 2013 International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M and C 2013). This proceedings contains over 250 full papers with topics ranging from reactor physics; radiation transport; materials science; nuclear fuels; core performance and optimization; reactor systems and safety; fluid dynamics; medical applications; analytical and numerical methods; algorithms for advanced architectures; and validation verification, and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

RELAP-7 Beta Release: Summary of Capabilities

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory (INL). Building upon the decades of software development at the INL, we began the development of RELAP-7 in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway. As part of this development, the first lines of RELAP-7 code were committed to the software revision control repository on November 7th, 2011. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. RELAP-7 is built using the INL’s modern scientific software development framework, MOOSE (Multi-physics Object Oriented Simulation Environment). MOOSE provides improved numerical calculations (including higher-order integration in both space and time, yielding converged second-order accuracy). The RELAP-7 code structure is based on multiple physical component models such as pipes, junctions, pumps, etc. Each component can have options for different fluid models such as single- and two-phase flow. This component-based and physics-based software architecture allows RELAP-7 to adopt different physical models for different applications. A relatively new two-phase hydrodynamic model, termed the ''7-Equation model'' for two phasic pressures, velocities, energies, and volumetric fraction, is incorporated into RELAP-7 for liquid-gas (water-steam) flows. This new model allows second-order integration because it is well-posed, which will reduce the numerical error associated with traditional systems analysis codes. In this paper, we provide a RELAP-7 capability list describing analysis features, range of applicability, and reactor components that will be available for the December 15th, 2014 beta release of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1990

Flowfield of a lifting hovering rotor: A Navier-Stokes simulation

The viscous, three-dimensional flowfield of a lifting helicopter rotor in hover is calculated by using an upwind, implicit, finite-difference numerical method for solving the thin layer Navier-Stokes equations. The induced effects of the wake, including the interaction of tip vortices with successive blades, are calculated as a part of the overall flowfield solution without using any ad hoc wake models. Comparison of the numerical results for the subsonic and transonic conditions show good agreement with the experimental data and with the previously published Navier-Stokes calculations using a simple wake model. Some comparisons with Euler calculations are also presented, along with some discussions of the grid refinement studies.

Srinivasan, G. R.↗
NASA NTRS1991

Numerical prediction of the unsteady flowfield around the F-18 aircraft at large incidence

This paper describes a numerical method capable of solving the steady and unsteady viscous flow around complete aircraft configurations at high angles of attack. This method is used to simulate the external flow around the F-18 aircraft, including deflected control surfaces. The current technique employs a generalized overset zonal grid scheme to decompose the computational space around the aircraft. The grid around various components of the aircraft are created numerically using a three-dimensional hyperbolic grid generation procedure. The Reynolds-averaged Navier-Stokes equations are integrated using a time-accurate, implicit procedure. Results for the turbulent flow around the F-18 aircraft at 30 degrees angle of attack show the details of the flowfield structure, including the unsteadiness created by the vortex burst and the resulting fluctuating airloads exerted on the vertical tail. The computed results agree fairly well with flight data for surface pressure, surface flow pattern, vortex burst location, and the dominant frequency for tail load fluctuations.

Rizk, Yehia M.↗
NASA NTRS1995

Optimal guidance law development for an advanced launch system

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Calise, Anthony J.↗