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.

SourceResearch connectionFrom yearThrough year

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

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Antenna feedhorn software upgrade

The HYBRIDHORN computer program was developed to serve as an item of general purpose antenna feedhorn design and analysis software. The formulation contains a small flare angle approximation which is subject to question for designs such as the S- and X-band feedhorn. Additionally, the original formulation did not allow azimuthal variation indexes other than unity. The HYBRIDHORN program was upgraded to correct both of these deficiencies. A large flare angle formulation was found. In the upgrade, all of the major program elements were converted to Univac 1108 compatible structured FORTRAN (SFTRAN) for ease of software maintenance. The small and large angle formulations are described and sample numerical results are presented.

Potter, P. D.↗

NASA NTRS1987

Stochastic estimation of conditional eddies in turbulent channel flow

Several long-standing issues regarding linear estimation and coherent structures are addressed, and some questions that were addressed partially, but never with the benefit of full, three-dimensional information are answered. The objectives were to: determine how well linear estimates approximate the field obtained by true conditional averaging, using events such as those in quadrant analysis; determine the extent to which the three-dimensional linearly estimated fields correspond to coherent structures, and the degree and manner in which they differ; evaluate the type and nature of the structural information gained by employing several different types of events; and learn more about the 3-D structure of important coherent motions that occur in wall turbulence. The results indicate that linear stochastic estimation can be used effectively in the study of numerical data bases consisting of three-dimensional vector fields, both velocity and vorticity. Two-point stochastic estimation yields more structural information and more detail than single-point estimation. The structures observed occurred repeatedly within the flow, but much can not be said about their dominance or the probability of their occurrence without further systematic studies of their frequency.

Adrian, R. J.↗

DOE OSTI2020

Quantitative metrics for evaluating thermonuclear design codes and physics models applied to the National Ignition Campaign

In this study, physics models and design codes for hot dense plasmas undergoing thermonuclear burn are evaluated objectively using statistical metrics that compare the difference between calculations and data relative to the experimental uncertainties. The analysis is applied to the National Ignition Campaign (NIC) because it is relevant, comprehensive, and well documented. The statistics confirm that a key process afflicting NIC performance is mix driven by hydrodynamic instabilities as approximated here using the KL model [G. Dimonte and R. Tipton, Phys. Fluids 18, 085101 (2006)]. New physics models are also presented for instability-driven magnetic fields [B. Srinivasan et al., Phys. Rev. Lett. 108, 165002 (2012)] and the Coulomb logarithm for electron–ion thermal relaxation [G. Dimonte and J. Daligault, Phys. Rev. Lett. 101, 135001 (2008)]. The plasma-generated magnetic fields improve the agreement between code and data in a statistically significant manner by reducing the electron thermal conduction in the hot-spot via the Hall term. The Coulomb logarithm is presented in a numerically practical form that incorporates recent theoretical advances. However, it does not improve the agreement between the code and data because the thermal relaxation is so fast in non-ignited NIC plasmas that the typical 20% error does not change the tight coupling between electrons and ions. Even with these improved physics models, the one-dimensional simulations presented here are not able to describe the high-convergence NIC implosions in a statistically acceptable manner. Of the heroic multi-dimensional simulations of Clark et al. [Phys. Plasmas 23, 056302 (2016)], only those in three-dimensions satisfy the statistical acceptance criterion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1991

Visual analysis of fluid dynamics at NASA's numerical aerodynamic simulation facility

A study aimed at describing and illustrating visualization tools used in Computational Fluid Dynamics (CFD) and indicating how these tools are likely to change by showing a projected resolution of the human computer interface is presented. The following are outlined using a graphically based test format: the revolution of human computer environments for CFD research; comparison of current environments; current environments with the ideal; predictions for the future CFD environments; what can be done to accelerate the improvements. The following comments are given: when acquiring visualization tools, potential rapid changes must be considered; environmental changes over the next ten years due to human computer interface cannot be fathomed; data flow packages such as AVS, apE, Explorer and Data Explorer are easy to learn and use for small problems, excellent for prototyping, but not so efficient for large problems; the approximation techniques used in visualization software must be appropriate for the data; it has become more cost effective to move jobs that fit on workstations and run only memory intensive jobs on the supercomputer; use of three dimensional skills will be maximized when the three dimensional environment is built in from the start.

Watson, Velvin R.↗

NASA NTRS2001

Solidification Dynamics of Metal Drops in a Free Fall Environment

Comparison of experimental observations were made with numerical solutions to a model of the heat transfer and solidification kinetics associated with the cooling of a molten drop during free fall, particularly with regard to the fraction of liquid transformed. Experimentally, silver drops (99.9%, 4-9 mm diameter) were levitated, melted, and released to fall through Marshall Space Flight Center's 105m drop tube in helium - 6% hydrogen and argon atmospheres. By systematically varying the drops initial superheat the extent of solidification prior to impact ranged from complete to none during the approximately 4.6s of free fall time. Analysis reveals the relative importance of the initial parameters affecting the cooling and solidification rates within the drop. A discussion of the conditions under which the actual observations deviate from the assumptions used in the model is presented.

Grugel, R. N.↗

NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

NASA NTRS1993

Optimal design of composite hip implants using NASA technology

Using an adaptation of NASA software, we have investigated the use of numerical optimization techniques for the shape and material optimization of fiber composite hip implants. The original NASA inhouse codes, were originally developed for the optimization of aerospace structures. The adapted code, which was called OPORIM, couples numerical optimization algorithms with finite element analysis and composite laminate theory to perform design optimization using both shape and material design variables. The external and internal geometry of the implant and the surrounding bone is described with quintic spline curves. This geometric representation is then used to create an equivalent 2-D finite element model of the structure. Using laminate theory and the 3-D geometric information, equivalent stiffnesses are generated for each element of the 2-D finite element model, so that the 3-D stiffness of the structure can be approximated. The geometric information to construct the model of the femur was obtained from a CT scan. A variety of test cases were examined, incorporating several implant constructions and design variable sets. Typically the code was able to produce optimized shape and/or material parameters which substantially reduced stress concentrations in the bone adjacent of the implant. The results indicate that this technology can provide meaningful insight into the design of fiber composite hip implants.

Blake, T. A.↗

NASA NTRS2004

Ground Ice at the Phoenix Landing Site: A Preflight Assessment

One of the objectives of the Mars Scout mission, Phoenix, is to characterize the present state of water in the martian environment, in a location where water may play a significant role in the present and past habitability of Mars. Given the generally dry and cold climate of Mars today any substantial amount of water is expected to occur in the form of ground ice (subsurface ice) within the regolith. The Mars Odyssey Gamma Ray Spectrometer has indicated abundant subsurface hydrogen and inferred ground ice at high latitudes. Therefore, the Phoenix mission will be targeted to land in the northern high latitudes (approximately 65 degrees N - 75 degrees N) where ground ice is expected to be abundantly available for analysis. The lander will be capable of excavating, sampling, and analyzing, dry and water-rich/icy soils. The location and depth of excavation necessary to achieve the goals of sampling and analysis of icy material become important parameters to assess. In the present work we ask two key questions: 1) At what depth within the regolith do we expect to find ice? 2) How might this depth vary over the region of potential landing sites? Numerous lines of evidence can be employed to provide an indication of the presence or absence of shallow ground ice at the potential landing sites. For example geomorphology, neutrons, gamma rays, and theory each contribute clues to an overall understanding of the distribution of ice. Orbital observations provide information on a variety of spatial scales, typically 10 s of meters (patterned ground) to 100 s of kilometers (gamma rays). While information on all of these scales are important, of particular interest is how the presence and depth of ground ice might vary on spatial scales comparable to the lander and its work area. While ground ice may be stable (and present) on a regional scale, local-scale slopes and changes in the physical characteristics of soils can result in significant variations in the distribution of ice.

Mellon, M. T.↗

DOE OSTI2021

Contingency Analysis Based on Partitioned and Parallel Holomorphic Embedding

In the steady-state contingency analysis, the traditional Newton-Raphson method suffers from non-convergence issues when solving post-outage power flow problems, which hinders the integrity and accuracy of security assessment. In this paper, we propose a novel robust contingency analysis approach based on holomorphic embedding (HE). Here, the HE-based simulator provides theoretical convergence guarantee, which is desirable because it avoids the influence of numerical issues and provides a credible security assessment conclusion. In addition, based on the multi-area characteristics of real-world power systems, a partitioned HE (PHE) method is proposed with an interfacebased partitioning of HE formulation. The PHE method does not undermine the numerical robustness of HE and significantly reduces the computation burden in large-scale contingency analysis. The PHE method is further enhanced by parallel or distributed computation to become parallel PHE (P2HE). Tests on a 458-bus system, a synthetic 419-bus system and a large-scale 21447-bus system demonstrate the advantages of the proposed methods in robustness and efficiency.

42 ENGINEERING↗

NASA NTRS1981

Time-temperature effect in adhesively bonded joints

The viscoelastic analysis of an adhesively bonded lap joint was reconsidered. The adherends are approximated by essentially Reissner plates and the adhesive is linearly viscoelastic. The hereditary integrals are used to model the adhesive. A linear integral differential equations system for the shear and the tensile stress in the adhesive is applied. The equations have constant coefficients and are solved by using Laplace transforms. It is shown that if the temperature variation in time can be approximated by a piecewise constant function, then the method of Laplace transforms can be used to solve the problem. A numerical example is given for a single lap joint under various loading conditions.

Delale, F.↗

NASA NTRS1988

A near optimal guidance algorithm for aero-assisted orbit transfer

The paper presents a near optimal guidance algorithm for aero-assited orbit plane change, based on minimizing the energy loss during the atmospheric portion of the maneuver. The guidance algorithm makes use of recent results obtained from energy state approximations and singular perturbation analysis of optimal heading change for a hypersonic gliding vehicle. This earlier work ignored the terminal constraint on altitude needed to insure that the vehicle exits that atmosphere. Thus, the resulting guidance algorithm was only appropriate for maneuvering reentry vehicle guidance. In the context of singular perturbation theory, a constraint on final altitude gives rise to a difficult terminal boundary layer problem, which cannot be solved in closed form. This paper will demonstrate the near optimality of a predictive/corrective guidance algorithm for the terminal maneuver. Comparisons are made to numerically optimized trajectories for a range or orbit plane angles.

Calise, Anthony J.↗

NASA NTRS1992

Analysis of orbital perturbations acting on objects in orbits near geosynchronous earth orbit

The paper presents a numerical investigation of orbital evolution for objects started in GEO or in orbits near GEO in order to study potential orbital debris problems in this region. Perturbations simulated include nonspherical terms in the earth's geopotential field, lunar and solar gravity, and solar radiation pressure. Objects simulated include large satellites, for which solar radiation pressure is insignificant, and small particles, for which solar radiation pressure is an important force. Results for large satellites are largely in agreement with previous GEO studies that used classical perturbation techniques. The orbit plane of GEO satellites placed in a stable plane orbit inclined approximately 7.3 deg to the equator experience very little precession, remaining always within 1.2 percent of their initial orientation. Solar radiation pressure generates two major effects on small particles: an orbital eccentricity oscillation anticipated from previous research, and an oscillation in orbital inclination.

Friesen, Larry J.↗

NASA NTRS1996

A wavelet-optimized, very high order adaptive grid and order numerical method

Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.

Jameson, Leland↗

DOE OSTI2021

Long-term Dynamical Stability in the Outer Solar System. I. The Regular and Chaotic Evolution of the 34 Largest Trans-Neptunian Objects

We carried out an extensive analysis of the stability of the outer solar system, making use of the frequency analysis technique over short-term integrations of nearly 100,000 test particles, as well as a statistical analysis of 200 1 Gyr long numerical simulations, which consider the mutual perturbations of the giant planets and the 34 largest trans-Neptunian objects (we have called all 34 objects “dwarf planets,” DPs, even if probably only the largest of them are true DPs). From the frequency analysis, we produced statistical diffusion maps for a wide region of the a–e phase-space plane; we also present the average diffusion time for orbits as a function of perihelion. We later turned our attention to the 34 DPs, making an individualized analysis for each of them and producing a first approximation of their future stability. From the 200 distinct realizations of the orbital evolution of the 34 DPs, we classified the sample into three categories, including 17 stable, 11 unstable, and 6 resonant objects; we also found that, statistically, two objects from the sample will leave the trans-Neptunian region within the next gigayear, most likely being ejected from the solar system, but with a nonnegligible probability of going inside the orbit of Neptune, either to collide with a giant planet or even falling to the inner solar system, where our simulations are no longer able to resolve their continuous evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1990

Structure and evolution of small-amplitude intermediate shock waves

A simplified set of equations is derived that approximates the magnetohydrodynamic (MHD) Navier-Stokes equations for weakly nonlinear disturbances whose speeds are close to the MHD intermediate speed. Its shock structure solutions are then examined. The fast and slow shock solutions are uniquely specified by their Rankine-Hugoniot relations. However, the intermediate shock solutions, which are not unique, are characterized by the integral through the shock of the noncoplanar component of the magnetic field. For situations in which this integral is conserved, the Riemann problem is well defined and predicts the evolution of intermediate shocks. This analysis is substantiated by numerical computations.

Kennel, C. F.↗

NASA NTRS1990

Combined design of structures and controllers for optimal maneuverability

Approaches to the combined design of structures and controllers for achieving optimal maneuverability are presented. A maneuverability index which directly reflects the minimum time required to perform a given set of maneuvers is introduced. By designing the flexible appendages, the maneuver time of the spacecraft is minimized under the constraints of structural properties, and post maneuver spillover is kept within a specified bound. The spillover reduction is achieved by making use of an appropriate reduced order model. The distributed parameter design problem is approached using assumed shape functions, and finite element analysis with dynamic reduction. Solution procedures have been investigated. Approximate design methods have been developed to overcome the computational difficulties. Some new constraints on the modal frequencies of the spacecraft are introduced in the original optimization problem to facilitate the solution process. It is shown that the global optimal design may be obtained by tuning the natural frequencies to satisfy specific constraints. Researchers quantify the difference between a lower bound to the solution for maneuver time associated with the original problem and the estimate obtained from the modified problem, for a specified application requirement. Numerical examples are presented to demonstrate the capability of this approach.

Ling, Jer↗

NASA NTRS1993

Antenna pattern control using impedance surfaces

A hybrid numerical technique is developed for electrically large pyramidal horn antennas radiating in free space. A stepped-waveguide method is used to analyze the interior surfaces of the horn transition. The Electric Field Integral Equation (EFIE) is employed on the outer surfaces of the pyramidal horn including the radiating aperture. Meanwhile, the Magnetic Field Integral Equation (MFIE) is used on the aperture to relate the aperture fields and those in the horn transition. The resultant hybrid field integral equation (HFIE) is solved numerically by the method of moments. This formulation is both accurate and numerically stable so that high-gain microwave pyramidal horns can be analyzed rigorously. Far-field radiation patterns, both computed and measured, are presented for three electrically-large x-band horn antennas. The comparisons demonstrate that this method is accurate enough to predict the fine pattern structure at wide angles and in the back region. Computed far-field patterns and aperture field distribution of two smaller x-band horns are also presented along with a discussion on the validity of the approximate aperture field distributions routinely used in the analysis and design of pyramidal horns.

Balanis, Constantine A.↗

NASA NTRS1996

An MHD Study of the Interaction Between the Solar Wind and the Interstellar Medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results have been compared with observations where possible and with the results from previous analytic and numerical studies. The primary accomplishment of this project has been the development of codes for 2-D models in both spherical and cylindrical coordinates and the application of the codes to the solar wind/interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. The simulations have also been used along with a data analysis study to provide a quantitative estimate of the distance to the termination and bow shocks. Some of the specific topics that have been studied are: (1) gas dynamic models of the solar wind/interstellar medium interaction, (2) termination shock response to large-scale solar wind fluctuations, and (3) distances to the termination shock and heliopause. The main results from each of these studies are summarized. The results were published in three papers which are included as attachments.

Steinolfson, R. S.↗