Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Radiation and scattering from cylindrically conformal printed antennas

Microstrip patch antennas offer considerable advantages in terms of weight, aerodynamic drag, cost, flexibility, and observables over more conventional protruding antennas. These flat patch antennas were first proposed over thirty years ago by Deschamps in the United States and Gutton and Baisinot in France. Such antennas have been analyzed and developed for planar as well as curved platforms. However, the methods used in these designs employ gross approximations, suffer from extreme computational burden, or require expensive physical experiments. The goal of this thesis is to develop accurate and efficient numerical modeling techniques which represent actual antenna structures mounted on curved surfaces with a high degree of fidelity. In this thesis, the finite element method is extended to cavity-backed conformal antenna arrays embedded in a circular, metallic, infinite cylinder. Both the boundary integral and absorbing boundary mesh closure conditions will be used for terminating the mesh. These two approaches will be contrasted and used to study the scattering and radiation behavior of several useful antenna configurations. An important feature of this study will be to examine the effect of curvature and cavity size on the scattering and radiation properties of wraparound conformal antenna arrays.

Kempel, Leo C.↗
NASA NTRS1992

Reliability analysis of composite structures

A probabilistic static stress analysis methodology has been developed to estimate the reliability of a composite structure. Closed form stress analysis methods are the primary analytical tools used in this methodology. These structural mechanics methods are used to identify independent variables whose variations significantly affect the performance of the structure. Once these variables are identified, scatter in their values is evaluated and statistically characterized. The scatter in applied loads and the structural parameters are then fitted to appropriate probabilistic distribution functions. Numerical integration techniques are applied to compute the structural reliability. The predicted reliability accounts for scatter due to variability in material strength, applied load, fabrication and assembly processes. The influence of structural geometry and mode of failure are also considerations in the evaluation. Example problems are given to illustrate various levels of analytical complexity.

Kan, Han-Pin↗
NASA NTRS1994

The stability of thick, self-gravitating disks in protostellar systems

Using numerical hydrodynamics techniques, we perform a nonlinear stability analysis of accretion disk systems that contain thick, self-gravitating disks. The systems are initially represented by a point mass M(sub c) at the center and a geometrically thick, axisymmetric disk of mass M(sub d) that supports uniform specific angular momentum and obeys an n = 3/2, polytropic equation of state. The equilibrium disk structure is uniquely defined upon the specification of two key dimensionless system parameters: M(sub d)/M(sub c) and T/absolute value of W (the ratio of rotational kinetic energy if the disk to the gravitational potential energy of the system). The focus of this work is on the identification of systems within this two-dimensional parameter space that are marginally unstable toward the development of nonaxisymmetric distortions. The geometric form and relative pattern speed of the disk's distortion as well as the likelihood of disk fragmentation as a result of such instabilities is examined, particularly in the context of protostellar systems. The value of T/absolute value of W at which thick disks first become dynamically unstable to nonaxisymmetric distortions is found to vary significantly with the mass ratio of the accertion disk system. Nonaxisymmetric eigenmodes with four distinctly different characters are identified in systems with mass ratios in the range 0.2 less than or equal to (M(sub d)/M(sub c)) less than or equal to 5.

Woodward, John W.↗
NASA NTRS1995

Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis

An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.

Burgreen, Gregory W.↗
NASA NTRS1995

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗
NASA NTRS1995

Rapid Mars transits with exhaust-modulated plasma propulsion

The operational characteristics of the Exhaust-Modulated Plasma Rocket are described. Four basic human and robotic mission scenarios to Mars are analyzed using numerical optimization techniques at variable specific impulse and constant power. The device is well suited for 'split-sprint' missions, allowing fast, one-way low-payload human transits of 90 to 104 days, as well as slower, 180-day, high-payload robotic precursor flights. Abort capabilities, essential for human missions, are also explored.

Chang-Diaz, Franklin R.↗
NASA NTRS1995

Low-energy ion acceleration at quasi-perpendicular shocks: Transverse diffusion

The problem of ion injection and acceleration at quasi perpendicular shocks has been the subject of some debate over the past two decades. It is widely known that these shocks efficiently accelerate particles that are well in the high-energy tail of the distribution. However, the issue of injection, or the acceleration of low-energy ions, has yet to reach a consensus. The fundamental issue is whether there is enough diffusion normal to the magnetic field for the particles to remain near the shock. Since transverse diffusion is a physical process that is not well understood in space plasmas, this is an important, and difficult issue to address. In this report, we will investigate the ion injection problem by performing test particle orbit integrations using synthesized turbulent fields. These fields are fully three-dimensional so that transverse diffusion is possible (cross-field diffusion is not possible in geometries where the electromagnetic fields are less than three dimensional). The synthesized fields are produced by superimposing a three-dimensional wave field on a background field. For completeness, we will compare the results from this model with the more well-established theories, such as the diffusive approximation and scatter-free shock drift acceleration. We will also compare these results with other numerical simulation techniques such as the well known hybrid simulation, and other test-particle calculations in which the shock fields are specified to have less than three dimensions. We will also discuss some recent relevant observations and how these compare with our results.

Giacalone, J.↗
NASA NTRS2019

Numerical analysis of unsteady flow in a converging-diverging nozzle

This study investigates the physics of jet noise by the numerical solution of the three dimensional unsteady Navier-Stokes equations for a low-supersonic jet (Mach numbers on the order of 1.4) from a converging-diverging (C-D) nozzle and mixing with a subsonic surrounding flow. The numerical method used for the solution of these equations is the second order accurate MacCormack explicit predictor-corrector scheme. The numerical results are compared to the experimental results of Yamamoto for the same C-D nozzle geometry and similar operating conditions. The goal of the study is to use numerical simulation techniques to resolve accurately physical processes which contribute to jet noise, and to assess whether performing these simulations in a three dimensional geometry is worth the additional cost in computing resources as compared to a comparable simulation on a two dimensional planar or axisymmetric geometry. Scientific visualization is used to observe the spreading characteristics of the jet, as well as to examine the structure of the shocks within the jet core. The time average velocity and turbulence intensity level at various locations are computed and compared with experimental results. A two point space-time cross correlation is used to compute the mean convection speed through the shear layer, and these predictions are compared with results from theory.

James N Scott↗
NASA NTRS1997

Specific Volumes of the Zr(41.2)Ti(13.8)Cu(12.5)Ni(10.0)Be(22.5) Alloy in the Liquid, Glass, and Crystalline States

The specific volumes of the Zr(41.2)Ti(3.8)Cu(2.5)Ni(10.0)Be(22.5) alloy as a function of temperature, T, are determined by employing an image digitizing technique and numerical calculation methods applied to the electrostatically levitated spherical alloy. The linear fitting of the volumes of the alloy in the liquid, V(sub l), glass, V(sub g) and crystalline V(sub c), states in the temperature ranges shown in parentheses are V(sub l)(T) = 0.1583 + 8.877 x 10(exp -6) T(cu cm/g) (700-1300 K);V(sub g)(T) = 0.1603 + 5.528 x 10(exp -6) T (400-550 K);V(sub c)(T) = 0.1583 + 6.21 x 10(exp -6)T(400-850 K). The average volume thermal expansion coefficients within the temperature ranges are determined to be 5.32, 3.39. and 3.83 x 10(exp -5) (1/K) for the liquid, glass, and crystalline states, respectively.

Ohsaka, K.↗
NASA NTRS2000

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗
NASA NTRS2000

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗
NASA NTRS2000

Actinic Flux Calculations: A Model Sensitivity Study

calculate direct and diffuse surface irradiance and actinic flux (downwelling (2p) and total (4p)) for the reference model. Sensitivity analysis has shown that the accuracy of the radiative transfer flux calculations for a unit ETS (i.e. atmospheric transmittance) together with a numerical interpolation technique for the constituents' vertical profiles is better than 1% for SZA less than 70(sub o) and wavelengths longer than 310 nm. The differences increase for shorter wavelengths and larger SZA, due to the differences in pseudo-spherical correction techniques and vertical discretetization among the codes. Our sensitivity study includes variation of ozone cross-sections, ETS spectra and the effects of wavelength shifts between vacuum and air scales. We also investigate the effects of aerosols on the spectral flux components in the UV and visible spectral regions. The "aerosol correction factors" (ACFs) were calculated at discrete wavelengths and different SZAs for each flux component (direct, diffuse, reflected) and prescribed IPMMI aerosol parameters. Finally, the sensitivity study was extended to calculation of selected photolysis rates coefficients.

Krotkov, Nickolay A.↗
NASA NTRS2000

Parallel Implementation of a High Order Implicit Collocation Method for the Heat Equation

We combine a high order compact finite difference approximation and collocation techniques to numerically solve the two dimensional heat equation. The resulting method is implicit arid can be parallelized with a strategy that allows parallelization across both time and space. We compare the parallel implementation of the new method with a classical implicit method, namely the Crank-Nicolson method, where the parallelization is done across space only. Numerical experiments are carried out on the SGI Origin 2000.

Kouatchou, Jules↗
NASA NTRS1999

Fast Quantum Algorithm for Predicting Descriptive Statistics of Stochastic Processes

Stochastic processes are used as a modeling tool in several sub-fields of physics, biology, and finance. Analytic understanding of the long term behavior of such processes is only tractable for very simple types of stochastic processes such as Markovian processes. However, in real world applications more complex stochastic processes often arise. In physics, the complicating factor might be nonlinearities; in biology it might be memory effects; and in finance is might be the non-random intentional behavior of participants in a market. In the absence of analytic insight, one is forced to understand these more complex stochastic processes via numerical simulation techniques. In this paper we present a quantum algorithm for performing such simulations. In particular, we show how a quantum algorithm can predict arbitrary descriptive statistics (moments) of N-step stochastic processes in just O(square root of N) time. That is, the quantum complexity is the square root of the classical complexity for performing such simulations. This is a significant speedup in comparison to the current state of the art.

Williams Colin P.↗
NASA NTRS2001

Geomagnetic Cutoff Rigidity Computer Program: Theory, Software Description and Example

The access of charged particles to the earth from space through the geomagnetic field has been of interest since the discovery of the cosmic radiation. The early cosmic ray measurements found that cosmic ray intensity was ordered by the magnetic latitude and the concept of cutoff rigidity was developed. The pioneering work of Stoermer resulted in the theory of particle motion in the geomagnetic field, but the fundamental mathematical equations developed have 'no solution in closed form'. This difficulty has forced researchers to use the 'brute force' technique of numerical integration of individual trajectories to ascertain the behavior of trajectory families or groups. This requires that many of the trajectories must be traced in order to determine what energy (or rigidity) a charged particle must have to penetrate the magnetic field and arrive at a specified position. It turned out the cutoff rigidity was not a simple quantity but had many unanticipated complexities that required many hundreds if not thousands of individual trajectory calculations to solve. The accurate calculation of particle trajectories in the earth's magnetic field is a fundamental problem that limited the efficient utilization of cosmic ray measurements during the early years of cosmic ray research. As the power of computers has improved over the decades, the numerical integration procedure has grown more tractable, and magnetic field models of increasing accuracy and complexity have been utilized. This report is documentation of a general FORTRAN computer program to trace the trajectory of a charged particle of a specified rigidity from a specified position and direction through a model of the geomagnetic field.

Smart, D. F.↗
NASA NTRS1997

Some Experiences with High Order Discretization Schemes for Compressible Fluid Flow

It is widely recognized that the structure and complexity of compressible fluid flow places high demands on numerical discretization techniques for the fluid flow equations. Fluid flows arising in external aerodynamics often contain both flow field discontinuities and fluid boundary-layers. Both must be accurately resolved to provide useful information to aerodynamic design and analysis engineers. These accuracy requirements motivated the present author to examine a class of finite-volume techniques on arbitrary triangulated domains based on linear or quadratic reconstruction of integral-averaged data followed by upwind flux function evaluation and small time evolution. More recently, we have considered some new upwind techniques which yield compact discretizations while maintaining higher order accuracy. In the mini-symposium talk we will discuss both of these techniques as well as demonstrate the relative merits of each method by computing a number of aerodynamic flows containing shock waves and boundary-layers.

Barth, Timothy J.↗
NASA NTRS1998

Three Program Architecture for Design Optimization

In this presentation, I would like to review historical perspective on the program architecture used to build design optimization capabilities based on mathematical programming and other numerical search techniques. It is rather straightforward to classify the program architecture in three categories as shown above. However, the relative importance of each of the three approaches has not been static, instead dynamically changing as the capabilities of available computational resource increases. For example, we considered that the direct coupling architecture would never be used for practical problems, but availability of such computer systems as multi-processor. In this presentation, I would like to review the roles of three architecture from historical as well as current and future perspective. There may also be some possibility for emergence of hybrid architecture. I hope to provide some seeds for active discussion where we are heading to in the very dynamic environment for high speed computing and communication.

Miura, Hirokazu↗
NASA NTRS2002

The Impact of Urbanization on the Precipitation Component of the Water Cycle: A New Perspective

It is estimated that by the year 2025, 60% of the world s population will live in cities (UNFP, 1999). As cities continue to grow, urban sprawl (e.g., the expansion of urban surfaces outward into rural surroundings) creates unique problems related to land use, transportation, agriculture, housing, pollution, and development. Urban expansion also has measurable impacts on environmental processes. Urban areas modify boundary layer processes through the creation of an urban heat island (UHI). The literature indicates that the signature of the urban heat island effect may be resolvable in rainfall patterns over and downwind of metropolitan areas. However, a recent U.S. Weather Research Program panel concluded that more observational and modeling research is needed in this area (Dabberdt et al. 2000). NASA and other agencies initiated programs such as the Atlanta Land-use Analysis: Temperature and Air Quality Project (ATLANTA) (Quattrochi et al. 1998) which aimed to identify and understand how urban heat islands impact the environment. However, a comprehensive assessment of the role of urban-induced rainfall in the global water and energy cycle (GWEC) and cycling of freshwater was not a primary focus of these efforts. NASA's Earth Science Enterprise (ESE) seeks to develop a scientific understanding of the Earth system and its response to natural or human-induced changes to enable improved prediction capability for climate, weather, and natural hazards (NASA, 2000). Within this mission, the ESE has three basic thrusts: science research to increase Earth system knowledge; an applications program to transfer science knowledge to practical use in society; and a technology program to enable new, better, and cheaper capabilities for observing the earth. Within this framework, a research program is underway to further address the co-relationship between land cover use and change (e.g. urban development) and its impact on key components of the GWEC (e.g., precipitation). This presentation discusses the feasibility of using the TRMM or GPM satellite to identify precipitation anomalies likely caused by urbanization (Shepherd et al. 2002). Recent results from analyses of TRMM data around several major U.S. cities (e.g. Dallas, Atlanta, Houston) will be discussed. The presentation also summarizes a NASA-funded research effort to investigate the phenomenon of urban-induced precipitation anomalies using TRMM (future GPM) satellite-based remote sensing, an intensive ground observation/validation effort near Atlanta, and coupled atmosphere-land numerical modeling techniques.

Shephard, J. Marshal↗