Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical mathematics”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A comparison of analytical and numerical models for steadily driven magnetic reconnection

The effect of boundary conditions in both analytical and numerical solutions of steadily driven reconnection on the reconnection of antiparallel magnetic fields is considered, with special attention given to the mathematical problem of properly specifying boundary conditions for the MHD equation. A unified formulation developed by Priest and Forbes (1986) for steady state reconnection, which includes the Petschek solution and the Sonnerup (1970) solution as special cases, is used to reinterpret the previous numerical experiments of steadily driven reconnection. It is shown that many contradictory features of these experiments were caused by the use of boundary conditions which are different from those required by Petschek's (1964) theory.

Forbes, T. G.↗
NASA NTRS1988

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗
NASA NTRS1989

An investigation of a mathematical model of an optically pumped Ti(3+):Al2O3 laser system

During the last several years, solid state lasers were developed that have the potential for meeting rigorous performance requirements for space-based remote sensing of the atmosphere. In order to design a stable and efficient laser and to understand the effect on laser output of changes in the physical and design parameters, an understanding of the development of the dynamical processes of the laser is necessary. Typically, the dynamical processes in a laser system are investigated via rate equations describing the evolution of the occupancy in the electronic levels and of the photon density in the laser cavity. There are two approaches to this type of study. Most often, for the sake of simplicity, the spatial variations of the dynamic variables in the laser system are disregarded and the mathematical model consists of a system of first order nonlinear ordinary differential equations (ODE). The second approach is to take into account both spatial and temporal variations in the dynamic variables in the laser cavity. The resulting model consists of a first order semilinear system of partial differential equations (PDE). The model which was studied was studied was generic in the sense that it was a four-level laser system, but the parameters used in the numerical study were specific to Titanium-doped sapphire. For simplicity, a constant, spatially uniform pumping scheme was considered. In addition, a simplification of the model was made so that it treats a single lasing wavelength with a narrow bandwidth. The purpose was to investigate both versions of the mathematical model and to determine whether the numerical solutions are similar both qualitatively and quantitatively. The systems of ordinary differential equations were solved numerically using a Runge-Kutta-Fehlberg algorithm which was very efficient for typical values of the physical parameters. A numerical scheme, based on the Modified Euler method, for computing solutions to the system of partial differential equations was developed and implemented. The PDE model was solved numerically at the expense of greatly increased computer time.

Roberts, Lila F.↗
NASA NTRS1984

Mathematical Modelling of High-speed Ribbbon Systems: a Case Study of Edge-defined Film-fed Growth

Finite element numerical analysis was used to solve the coupled problem of heat transfer and capillarity to describe low and high speed silicon sheet growth in meniscus defined systems. Heat transfer models which neglect the details of convective heat flow in the melt are used to establish operating limits for an EFG system in terms of the growth rate, die temperature and the static head acting on the meniscus. It is shown that convective heat transfer in the melt becomes important only at high growth rates or for materials with low thermal conductivities.

Ettouney, H. M.↗
NASA NTRS1929

Mathematical and experimental investigation of heat control and power increase in air-cooled aircraft engines

In order to understand the numerical relations between the air velocity, temperature of the cylinder walls, heat dissipation, cylinder dimensions and type of construction an experimental plant was installed in the Siemens and Halske laboratory. The experimental cylinder was exposed to the air stream of a wind tunnel. The compression chamber was heated by an electrically heated oil bath kept constantly in motion by a stirrer. The wall temperatures were measured by thermocouples. The air stream was produced a seven-watt blower. The air flowed through a current rectifier (honeycomb), diffuser, air chamber with quieting sieves and a nozzle.

ENGINES - AIR COOLED - POWER INCREASES↗
NASA NTRS1997

Efficient Parallel Kernel Solvers for Computational Fluid Dynamics Applications

Distributed-memory parallel computers dominate today's parallel computing arena. These machines, such as Intel Paragon, IBM SP2, and Cray Origin2OO, have successfully delivered high performance computing power for solving some of the so-called "grand-challenge" problems. Despite initial success, parallel machines have not been widely accepted in production engineering environments due to the complexity of parallel programming. On a parallel computing system, a task has to be partitioned and distributed appropriately among processors to reduce communication cost and to attain load balance. More importantly, even with careful partitioning and mapping, the performance of an algorithm may still be unsatisfactory, since conventional sequential algorithms may be serial in nature and may not be implemented efficiently on parallel machines. In many cases, new algorithms have to be introduced to increase parallel performance. In order to achieve optimal performance, in addition to partitioning and mapping, a careful performance study should be conducted for a given application to find a good algorithm-machine combination. This process, however, is usually painful and elusive. The goal of this project is to design and develop efficient parallel algorithms for highly accurate Computational Fluid Dynamics (CFD) simulations and other engineering applications. The work plan is 1) developing highly accurate parallel numerical algorithms, 2) conduct preliminary testing to verify the effectiveness and potential of these algorithms, 3) incorporate newly developed algorithms into actual simulation packages. The work plan has well achieved. Two highly accurate, efficient Poisson solvers have been developed and tested based on two different approaches: (1) Adopting a mathematical geometry which has a better capacity to describe the fluid, (2) Using compact scheme to gain high order accuracy in numerical discretization. The previously developed Parallel Diagonal Dominant (PDD) algorithm and Reduced Parallel Diagonal Dominant (RPDD) algorithm have been carefully studied on different parallel platforms for different applications, and a NASA simulation code developed by Man M. Rai and his colleagues has been parallelized and implemented based on data dependency analysis. These achievements are addressed in detail in the paper.

Sun, Xian-He↗
NASA NTRS2017

Mathematical Modeling of Electrodynamics Near the Surface of Earth and Planetary Water Worlds

An interesting feature of planetary bodies with hydrospheres is the presence of an electrically conducting shell near the global surface. This conducting shell may typically lie between relatively insulating rock, ice, or atmosphere, creating a strong constraint on the flow of large-scale electric currents. All or parts of the shell may be in fluid motion relative to main components of the rotating planetary magnetic field (as well as the magnetic fields due to external bodies), creating motionally-induced electric currents that would not otherwise be present. As such, one may expect distinguishing features in the types of electrodynamic processes that occur, as well as an opportunity for imposing specialized mathematical methods that efficiently address this class of application. The purpose of this paper is to present and discuss such specialized methods. Specifically, thin-shell approximations for both the electrodynamics and fluid dynamics are combined to derive simplified mathematical formulations describing the behavior of these electric currents as well as their associated electric and magnetic fields. These simplified formulae allow analytical solutions featuring distinct aspects of the thin-shell electrodynamics in idealized cases. A highly efficient numerical method is also presented that is useful for calculations under inhomogeneous parameter distributions. Finally, the advantages as well as limitations in using this mathematical approach are evaluated. This evaluation is presented primarily for the generic case of bodies with water worlds or other thin spherical conducting shells. More specific discussion is given for the case of Earth, but also Europa and other satellites with suspected oceans.

Tyler, Robert H.↗
NASA NTRS1988

Positioning Rotors In Turbine Flowmeters

Lengths of wakes roughly proportional to thickness of vanes. Mathematical model simplifies analysis of effects of flow-straightening vanes in turbine flowmeter. Yields numerical solution of differential equations of flow for quick examination of efforts of thicknesses of vanes and rate of flow on extent of wake behind vanes. From examination, minimum distance at which flowmeter rotor placed behind vanes determined.

Lynch, Edward D.↗
NASA NTRS2002

Minimizing Segregation During the Controlled Directional Solidification of Dendritic Alloys Publication: Metallurgical and Materials Transactions

Gravity-driven thermosolutal convection that arises during controlled directional solidification (DS) of dendritic alloys promotes detrimental macro-segregation (e.g. freckles and steepling) in products such as turbine blades. Considerable time and effort has been spent to experimentally and theoretically investigate this phenomena; although our knowledge has advanced to the point where convection can be modeled and accurately compared to experimental results, little has been done to minimize its onset and deleterious effects. The experimental work demonstrates that segregation can be. minimized and microstructural uniformity promoted when a slow axial rotation is applied to the sample crucible during controlled directional solidification processing. Numerical modeling utilizing continuation and bifurcation methods have been employed to develop accurate physical and mathematical models with the intent of identifying and optimizing processing parameters.

Grugel, R. N.↗
NASA NTRS2002

A Lightning Channel Retrieval Algorithm for the North Alabama Lightning Mapping Array (LMA)

A new multi-station VHF time-of-arrival (TOA) antenna network is, at the time of this writing, coming on-line in Northern Alabama. The network, called the Lightning Mapping Array (LMA), employs GPS timing and detects VHF radiation from discrete segments (effectively point emitters) that comprise the channel of lightning strokes within cloud and ground flashes. The network will support on-going ground validation activities of the low Earth orbiting Lightning Imaging Sensor (LIS) satellite developed at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama. It will also provide for many interesting and detailed studies of the distribution and evolution of thunderstorms and lightning in the Tennessee Valley, and will offer many interesting comparisons with other meteorological/geophysical wets associated with lightning and thunderstorms. In order to take full advantage of these benefits, it is essential that the LMA channel mapping accuracy (in both space and time) be fully characterized and optimized. In this study, a new revised channel mapping retrieval algorithm is introduced. The algorithm is an extension of earlier work provided in Koshak and Solakiewicz (1996) in the analysis of the NASA Kennedy Space Center (KSC) Lightning Detection and Ranging (LDAR) system. As in the 1996 study, direct algebraic solutions are obtained by inverting a simple linear system of equations, thereby making computer searches through a multi-dimensional parameter domain of a Chi-Squared function unnecessary. However, the new algorithm is developed completely in spherical Earth-centered coordinates (longitude, latitude, altitude), rather than in the (x, y, z) cartesian coordinates employed in the 1996 study. Hence, no mathematical transformations from (x, y, z) into spherical coordinates are required (such transformations involve more numerical error propagation, more computer program coding, and slightly more CPU computing time). The new algorithm also has a more realistic definition of source altitude that accounts for Earth oblateness (this can become important for sources that are hundreds of kilometers away from the network). In addition, the new algorithm is being applied to analyze computer simulated LMA datasets in order to obtain detailed location/time retrieval error maps for sources in and around the LMA network. These maps will provide a more comprehensive analysis of retrieval errors for LMA than the 1996 study did of LDAR retrieval errors. Finally, we note that the new algorithm can be applied to LDAR, and essentially any other multi-station TWA network that depends on direct line-of-site antenna excitation.

Koshak, William↗
NASA NTRS2013

Mixed Integer Programming and Heuristic Scheduling for Space Communication

Optimal planning and scheduling for a communication network was created where the nodes within the network are communicating at the highest possible rates while meeting the mission requirements and operational constraints. The planning and scheduling problem was formulated in the framework of Mixed Integer Programming (MIP) to introduce a special penalty function to convert the MIP problem into a continuous optimization problem, and to solve the constrained optimization problem using heuristic optimization. The communication network consists of space and ground assets with the link dynamics between any two assets varying with respect to time, distance, and telecom configurations. One asset could be communicating with another at very high data rates at one time, and at other times, communication is impossible, as the asset could be inaccessible from the network due to planetary occultation. Based on the network's geometric dynamics and link capabilities, the start time, end time, and link configuration of each view period are selected to maximize the communication efficiency within the network. Mathematical formulations for the constrained mixed integer optimization problem were derived, and efficient analytical and numerical techniques were developed to find the optimal solution. By setting up the problem using MIP, the search space for the optimization problem is reduced significantly, thereby speeding up the solution process. The ratio of the dimension of the traditional method over the proposed formulation is approximately an order N (single) to 2*N (arraying), where N is the number of receiving antennas of a node. By introducing a special penalty function, the MIP problem with non-differentiable cost function and nonlinear constraints can be converted into a continuous variable problem, whose solution is possible.

Lee, Charles H.↗
NASA NTRSDate not supplied

Clifford Neural Operators on Atmospheric Data Influenced Partial Differential Equations

Mathematical representations of the atmosphere are key to forecasting and research tasks across Earth science. Numerically solving the underlying partial differential equations(PDEs) of the atmosphere, however, can be difficult and computationally expensive with numerous trade-offs between computing efficiency and accuracy. Utilizing neural net-works to learn approximations of the PDE solutions from the data can help us model complex phenomena more efficiently than traditional numerical schemes. Here, we have applied Clifford algebra-based neural operators for predicting atmospheric variables. Clifford Fourier neural operators are used with two different backbone architectures, ResNet and UNet, on custom data of U10, V10 and surface pressure as well as U500, V500 and Z500. Clifford Fourier neural operators, coupled with ResNet and UNet architectures, are applied to a key reanalysis dataset. Model performance is initially strong, but we observe increasing errors, resulting in the model becoming highly unstable.

Sujit Roy↗
NASA NTRS2002

Annual Research Briefs, 2002: Center for Turbulence Research

Turbulent combustion remains the largest component of the CTR's core program. This program and several related activities at CTR are supported by NASA's Ultra Efficient Engine Technology Program. It is also intimately connected with the Department of Energy's ASCI program at Stanford which develops the technology for numerical simulation of realistic aircraft engines using state of the art massively parallel computers. In combustion modeling the attention has been directed to the modeling of higher levels of complexity such as spray dynamics, radiation and soot formation. Major aircraft engine manufacturers have shown considerable interest in this program; in particular, a significant active collaboration exists between CTR and the Pratt & Whitney Corporation. CTR's combustion program is essentially based on the large-eddy simulation technique, LES, which is actively being pursued at CTR for this and many other applications. Important accomplishments in LES included simulations with three-dimensional filters which result in grid independent calculations (that is why we call it "true" LES), and the development of the methodology for integration of LES and Reynolds Averaged computations. Optimization techniques are being studied and used for the important problem of wall boundary conditions for LES as well as for optimal shape design for aeroacoustic and aerodynamic performance gains.

ULTRA EFFICIENT ENGINE TECHNOLOGY↗
NASA NTRS1975

The finite element method in low speed aerodynamics

The finite element procedure is shown to be of significant impact in design of the 'computational wind tunnel' for low speed aerodynamics. The uniformity of the mathematical differential equation description, for viscous and/or inviscid, multi-dimensional subsonic flows about practical aerodynamic system configurations, is utilized to establish the general form of the finite element algorithm. Numerical results for inviscid flow analysis, as well as viscous boundary layer, parabolic, and full Navier Stokes flow descriptions verify the capabilities and overall versatility of the fundamental algorithm for aerodynamics. The proven mathematical basis, coupled with the distinct user-orientation features of the computer program embodiment, indicate near-term evolution of a highly useful analytical design tool to support computational configuration studies in low speed aerodynamics.

Baker, A. J.↗
NASA NTRS1994

Numerical Analysis of an Impinging Jet Reactor for the CVD and Gas-Phase Nucleation of Titania

We model a cold-wall atmospheric pressure impinging jet reactor to study the CVD and gas-phase nucleation of TiO2 from a titanium tetra-iso-propoxide (TTIP)/oxygen dilute source gas mixture in nitrogen. The mathematical model uses the computational code FIDAP and complements our recent asymptotic theory for high activation energy gas-phase reactions in thin chemically reacting sublayers. The numerical predictions highlight deviations from ideality in various regions inside the experimental reactor. Model predictions of deposition rates and the onset of gas-phase nucleation compare favorably with experiments. Although variable property effects on deposition rates are not significant (approximately 11 percent at 1000 K), the reduction rates due to Soret transport is substantial (approximately 75 percent at 1000 K).

Gokoglu, Suleyman A.↗