Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1981

Numerical simulation of flows in curved diffusers with cross-sectional transitioning using a three-dimensional viscous analysis

A three dimensional analysis for fully viscous, subsonic, compressible flow is evaluated. An approximate form of the Navier Stokes equations is solved by an implicit spatial marching technique. Calculations were made for flow in a circular S duct and in the F 16 inlet duct. The computed total pressure contours and secondary flow velocity vectors are presented. Qualitative comparisons with experiment are shown for both ducts. The analysis is used to show how the cross section transitioning in the F 16 inlet suppresses the development of a secondary flow vortex.

Towne, C. E.↗
NASA NTRS1985

Numerical Estimation of the Curvature of Biological Surfaces

Many biological systems may profitably be studied as surface phenomena. A model consisting of isotropic growth of a curved surface from a flat sheet is assumed. With such a model, the Gaussian curvature of the final surface determines whether growth rate of the surface is subharmonic or superharmonic. These properties correspond to notions of convexity and concavity, and thus to local excess growth and local deficiency of growth. In biological models where the major factors controlling surface growth are intrinsic to the surface, researchers thus gained from geometrical study information on the differential growth undergone by the surface. These ideas were applied to an analysis of the folding of the cerebral cortex, a geometrically rather complex surface growth. A numerical surface curvature technique based on an approximation to the Dupin indicatrix of the surface was developed. A metric for comparing curvature estimates is introduced, and considerable numerical testing indicated the reliability of this technique.

Todd, P. H.↗
NASA NTRS1986

Parametric analytical study of instability-related delamination growth

The effect of various parameters on instability-related delamination growth was studied analytially. The configuration studied consisted of a thick composite laminate with a single through-width delamination located near one surface. Both mechanical and thermal loads were considered. All conclusions were based on the assumption that GI and GII govern delamination growth. An approximate superposition stress analysis was developed which gives closed form expressions for GI and GII. The simplicity of the analysis permitted examination of numerous configurations. Both GI and GII were found to be very sensitive to delamination length and location through the thickness. The magnitude of GI was also very sensitive to initial imperfections, which might be the results of an inclusion of finite thickness in the delamination. Critical loads for delamination growth were calculated based on three growth criteria. Large differences in the predictions highlight the need for a verified mixed-mode delamination growth criterion.

Whitcomb, J. D.↗
NASA NTRS1997

Wave Number Selection for Incompressible Parallel Jet Flows Periodic in Space

The temporal instability of a spatially periodic parallel flow of an incompressible inviscid fluid for various jet velocity profiles is studied numerically using Floquet Analysis. The transition matrix at the end of a period is evaluated by direct numerical integration. For verification, a method based on approximating a continuous function by a series of step functions was used. Unstable solutions were found only over a limited range of wave numbers and have a band type structure. The results obtained are analogous to the behavior observed in systems exhibiting complexity at the edge of order and chaos.

Miles, Jeffrey Hilton↗
NASA NTRS1998

Preliminary Work for Modeling the Propellers of an Aircraft as a Noise Source in an Acoustic Boundary Element Analysis

An algorithm for generating appropriate velocity boundary conditions for an acoustic boundary element analysis from the kinematics of an operating propeller is presented. It constitutes the initial phase of Integrating sophisticated rotorcraft models into a conventional boundary element analysis. Currently, the pressure field is computed by a linear approximation. An initial validation of the developed process was performed by comparing numerical results to test data for the external acoustic pressure on the surface of a tilt-rotor aircraft for one flight condition.

Vlahopoulos, Nickolas↗
NASA NTRS2008

Lunar Soil Erosion Physics for Landing Rockets on the Moon

To develop a lunar outpost, we must understand the blowing of soil during launch and landing of the new Altair Lander. For example, the Apollo 12 Lunar Module landed approximately 165 meters from the deactivated Surveyor Ill spacecraft, scouring its surfaces and creating numerous tiny pits. Based on simulations and video analysis from the Apollo missions, blowing lunar soil particles have velocities up to 2000 m/s at low ejection angles relative to the horizon, reach an apogee higher than the orbiting Command and Service Module, and travel nearly the circumference of the Moon [1-3]. The low ejection angle and high velocity are concerns for the lunar outpost.

Clegg, Ryan N.↗
NASA NTRS1990

A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals

The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.

Brush, L. N.↗
NASA NTRS1982

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗
NASA NTRS1984

Applications of finite element and wave envelope element approximations to turbofan engine noise radiation including flight effects

The problem of acoustic radiation from turbofan engine inlets in flow has not lent itself fully to analysis by numerical means because of the large domains and high frequencies involved. The current work has extended the use of finite elements and wave envelope elements, elements which simulate decay and wavelike behaviour in their interpolation functions, from the no-flow case in which they have been proven, to cases incorporating mean flow. By employing an irrotational mean flow assumption, the acoustics problem has been posed in an axisymmetric formulation in terms of acoustic velocity potential, thus minimizing computer solution storage requirements. The results obtained from the numerical procedures agree well with known analytical solutions, static experimental jet engines inflow data, and also with flight test results.

Parrett, A. V.↗
NASA NTRS1981

On the utility of finite element theory for computational fluid dynamics

An implicit finite element numerical solution algorithm is derived for the compressible Navier-Stokes equations expressed in generalized coordinates. The theoretical basis utilizes a Galerkin-Weighted Residuals formulation, and extremization of approximation error within the context of a multipole expansion. A von Neumann analysis for a simplified form indicates the algorithm fourth- to sixth-order phase accurate, with third-order dissipation for the elementary linear element construction. Performance is improved for the algorithm constructed using quadratic interpolation. Numerical experiments for shocked duct flows are employed to optimize the several algorithm parameters. Additional numerical solutions validate algorithm accuracy and utility for aerodynamics applications.

Baker, A. J.↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used, with huge success, to model phenomena arising across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDE models fail to adequately model observed phenomena or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. In this article, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis, and specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference, and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modeling and algorithmic extensions which serve to show the wide applicability of nonlocal modeling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The eXtended virtual element method for elliptic problems with weakly singular solutions

This paper introduces a novel eXtended virtual element method, an extension of the conforming virtual element method. The X-VEM is formulated by incorporating appropriate enrichment functions in the local spaces. The method is designed to handle highly generic enrichment functions, including singularities arising from fractured domains. By achieving consistency on the enrichment space, the method is proven to achieve arbitrary approximation orders even in the presence of singular solutions. The paper includes a complete convergence analysis under general assumptions on mesh regularity, and numerical experiments validating the method’s accuracy on various mesh families, demonstrating optimal convergence rates in the L 2 - and H 1 - norms on fractured or L-shaped domains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

The quasi-rigid rotation of coronal magnetic fields

Spherical harmonic analysis and numerical simulations are used to study the rotational properties of the coronal magnetic field under the assumption that it can be approximated by a current-free extension of the photospheric field. It is found that the rotation rate in the outer corona is determined, principally, by coronal filtering, the global averages of the photospheric rotation rate, and ongoing source eruptions. The present model is able to account for observationally inferred rotational properties. It is suggested that the coronal rotation rate accelerates gradually due to the equatorward migration of sunspots, and that the 27-day equatorial period is approached toward sunspot minimum as the decaying photospheric flux becomes localized near the equator.

Wang, Y.-M.↗
NASA NTRS2014

Mission Design and Analysis for Suborbital Intercept and Fragmentation of an Asteroid with Very Short Warning Time

Small near-Earth objects (NEOs) approximately 50-150 m in size are far more numerous (hundreds of thousands to millions yet to be discovered) than larger NEOs. Small NEOs, which are mostly asteroids rather than comets, are very faint in the night sky due to their small sizes, and are, therefore, difficult to discover far in advance of Earth impact. Furthermore, even small NEOs are capable of creating explosions with energies on the order of tens or hundreds of megatons (Mt). We are, therefore, motivated to prepare to respond effectively to short warning time, small NEO impact scenarios. In this paper we explore the lower bound on actionable warning time by investigating the performance of notional upgraded Intercontinental Ballistic Missiles (ICBMs) to carry Nuclear Explosive Device (NED) payloads to intercept and disrupt a hypothetical incoming NEO at high altitudes (generally at least 2500 km above Earth). We conduct this investigation by developing optimal NEO intercept trajectories for a range of cases and comparing their performances. Our results show that suborbital NEO intercepts using Minuteman III or SM-3 IIA launch vehicles could achieve NEO intercept a few minutes prior to when the NEO would strike Earth. We also find that more powerful versions of the launch vehicles (e.g., total delta V of approximately 9.5-11 km/s) could intercept incoming NEOs several hours prior to when the NEO would strike Earth, if launched at least several days prior to the time of intercept. Finally, we discuss a number of limiting factors and practicalities that affect whether the notional systems we describe could become feasible.

near- Earth object↗
NASA NTRS2006

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries, such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quality of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral waveforms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with the corresponding 3PN and 3.5PN orbital phasing. We find agreement consistent with internal error estimates based on either approach at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger/

Baker, John G.↗