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

HZE beam transport in multilayered materials

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multiple-layered materials. The code is established to operate on the Langley nuclear fragmentation model used in space engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurement. Comparison with 56Fe ion with Pb-Al and Pb-(CH2)x targets shows reasonable agreement.

NASA Discipline Radiation Health↗
NASA NTRS2015

Thermal Testing and Model Correlation of the Magnetospheric Multiscale (MMS) Observatories

International Conference on Envronmental Systems (ICES), Seattle WA NCTS 20964-15. The Magnetospheric Multiscale (MMS) mission is a Solar Terrestrial Probes mission comprising four identically instrumented spacecraft that will use Earths magnetosphere as a laboratory tostudy the microphysics of three fundamental plasma processes: magnetic reconnection, energetic particle acceleration, and turbulence. This paper presents the complete thermal balance (TB) test performed on the first of four observatories to go through thermal vacuum (TV) and the minibalance testing that was performed on the subsequent observatories to provide a comparison of all four. The TV and TB tests were conducted in a thermal vacuum chamber at the Naval Research Laboratory (NRL) in Washington, D.C. with the vacuum level higher than 1.3 x 10-4 Pa (10-6 torr)and the surrounding temperature achieving -180 C. Three TB test cases were performed that included hot operational science, cold operational science and a cold survival case. In addition to the three balance cases a two hour eclipse and a four hour eclipse simulation was performed during the TV test to provide additional transient data points that represent the orbit in eclipse (or Earth's shadow) The goal was to perform testing such that the flight orbital environments could be simulated as closely as possible. A thermal model correlation between the thermal analysis and the test results was completed. Over 400 1-Wire temperature sensors, 200 thermocouples and 125 flight thermistor temperature sensors recorded data during TV and TB testing. These temperatureversus time profiles and their agreements with the analytical results obtained using Thermal Desktop and SINDAFLUINT are discussed. The model correlation for the thermal mathematical model (TMM) is conducted based on the numerical analysis results and the test data. The philosophy of model correlation was to correlate the model to within 3 C of the test data using the standard deviation and mean deviation error calculation. Individual temperature error goal is to be within 5 C and the heater power goal is to be within 5 of test data. The results of the model correlation are discussed and the effect of some material and interface parameters on the temperature profiles are presented.

Thermal↗
NASA NTRS2015

Thermal Testing and Model Correlation of the Magnetospheric Multiscale (MMS) Observatories

The Magnetospheric Multiscale (MMS) mission is a Solar Terrestrial Probes mission comprising four identically instrumented spacecraft that will use Earth's magnetosphere as a laboratory to study the microphysics of three fundamental plasma processes: magnetic reconnection, energetic particle acceleration, and turbulence. This paper presents the complete thermal balance (TB) test performed on the first of four observatories to go through thermal vacuum (TV) and the minibalance testing that was performed on the subsequent observatories to provide a comparison of all four. The TV and TB tests were conducted in a thermal vacuum chamber at the Naval Research Laboratory (NRL) in Washington, D.C. with the vacuum level higher than 1.3 x 10 (sup -4) pascals (10 (sup -6) torr) and the surrounding temperature achieving -180 degrees Centigrade. Three TB test cases were performed that included hot operational science, cold operational science and a cold survival case. In addition to the three balance cases a two hour eclipse and a four hour eclipse simulation was performed during the TV test to provide additional transient data points that represent the orbit in eclipse (or Earth's shadow) The goal was to perform testing such that the flight orbital environments could be simulated as closely as possible. A thermal model correlation between the thermal analysis and the test results was completed. Over 400 1-Wire temperature sensors, 200 thermocouples and 125 flight thermistor temperature sensors recorded data during TV and TB testing. These temperature versus time profiles and their agreements with the analytical results obtained using Thermal Desktop and SINDA/FLUINT are discussed. The model correlation for the thermal mathematical model (TMM) is conducted based on the numerical analysis results and the test data. The philosophy of model correlation was to correlate the model to within 3 degrees Centigrade of the test data using the standard deviation and mean deviation error calculation. Individual temperature error goal is to be within 5 degrees Centigrade and the heater power goal is to be within 5 percent of test data. The results of the model correlation are discussed and the effect of some material and interface parameters on the temperature profiles are presented.

Correlation↗
NASA NTRS2019

Assessment of a One-Dimensional Finite Element Charring Ablation Material Response Model for Phenolic-Impregnated Carbon Ablator

In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.

Thermal Protection System↗
DOE OSTI2023

Direct Numerical Simulation of Flame-Wall Interaction for Low-Carbon Gas Turbine Combustion

Green hydrogen (H2) and ammonia (NH3) are emerging as carbon free alternatives to hydrocarbon fuels. In gas turbines, flames interact with the combustor liner (wall), which affects the pollutant emissions, the burning efficiency, and the thermal load on the liner. We lack understanding of this important flame-wall interaction (FWI) for alternative fuels. FWI occurs at the scale of the flame thickness, and such scales are only resolved by direct numerical simulation. Here, we resolve the FWI for NH3 and H2 flames. Preliminary results of two distinct projects are presented: 2D laminar anchored V-flame, and 3D turbulent swirling flame. For the laminar flame, we show flame quenching at the wall and resulting leakage of NH3 and increased N2O, but decreased NO emissions. For the swirling flame, we show CO emissions, and comparison of flame shapes and location to experimental measurements. Ultimately we show that quenching and pollutant emissions strongly affected by fuel/air ratio, and that an inclined wall can stabilize very lean turbulent CH4/H2 flames, but prone to large CO emissions.

adaptive mesh refinement↗
NASA NTRS2001

The CE/SE Method: a CFD Framework for the Challenges of the New Millennium

The space-time conservation element and solution element (CE/SE) method, which was originated and is continuously being developed at NASA Glenn Research Center, is a high-resolution, genuinely multidimensional and unstructured-mesh compatible numerical method for solving conservation laws. Since its inception in 1991, the CE/SE method has been used to obtain highly accurate numerical solutions for 1D, 2D and 3D flow problems involving shocks, contact discontinuities, acoustic waves, vortices, shock/acoustic waves/vortices interactions, shock/boundary layers interactions and chemical reactions. Without the aid of preconditioning or other special techniques, it has been applied to both steady and unsteady flows with speeds ranging from Mach number = 0.00288 to 10. In addition, the method has unique features that allow for (i) the use of very simple non-reflecting boundary conditions, and (ii) a unified wall boundary treatment for viscous and inviscid flows. The CE/SE method was developed with the conviction that, with a solid foundation in physics, a robust, coherent and accurate numerical framework can be built without involving overly complex mathematics. As a result, the method was constructed using a set of design principles that facilitate simplicity, robustness and accuracy. The most important among them are: (i) enforcing both local and global flux conservation in space and time, with flux evaluation at an interface being an integral part of the solution procedure and requiring no interpolation or extrapolation; (ii) unifying space and time and treating them as a single entity; and (iii) requiring that a numerical scheme be built from a nondissipative core scheme such that the numerical dissipation can be effectively controlled and, as a result, will not overwhelm the physical dissipation. Part I of the workshop will be devoted to a discussion of these principles along with a description of how the ID, 2D and 3D CE/SE schemes are constructed. In Part II, various applications of the CE/SE method, particularly those involving chemical reactions and acoustics, will be presented. The workshop will be concluded with a sketch of the future research directions.

Chang, Sin-Chung↗
NASA NTRS1993

Unified approach for incompressible flows

An unified approach for solving both compressible and incompressible flows was investigated in this study. The difference in CFD code development between incompressible and compressible flows is due to the mathematical characteristics. However, if one can modify the continuity equation for incompressible flows by introducing pseudocompressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of a compressible flow code to solve incompressible flows becomes feasible. Among numerical algorithms developed for compressible flows, the Centered Total Variation Diminishing (CTVD) schemes possess better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that CTVD schemes can equally well solve incompressible flows. In this study, the governing equations for incompressible flows include the continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the boundary conditions including physical and numerical boundary conditions must be properly specified to obtain accurate solution. The CFD code for this research is currently in progress. Flow past a circular cylinder will be used for numerical experiments to determine the accuracy and efficiency of the code before applying this code to more specific applications.

Chang, Tyne-Hsien↗
NASA NTRS1984

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗
NASA NTRS1991

A study of space-rated connectors using a robotic end-effector

The fabrication and testing of a pair of robot fingers designed to grasp two types of Orbital Replacement Unit (ORU) interfaces, the H Handle type and the Micro Square type are discussed. A closed-form solution is given for the force inverse kinematics. A numerical solution using the Newton-Raphson Method for force forward kinematics is given. Mathematical expressions are derived to compute forces/torques applied to the finger. Suggestions are given for improving finger fabrication. The results of numerous experiments conducted to study the characteristics and feasibility of the fingers are given. The first part of the study was devoted to obtaining data on the forces applied by the finger to the interfaces under various translational and rotational misalignments; the second part was devoted to determining the maximum allowed capture angles that would insure successful mating; and the third part was devoted to the processing and interpretation of the forces/torque data.

Nguyen, Charles C.↗
NASA NTRS1992

Numerical simulation of freckle formation in directional solidification of binary alloys

A mathematical model of solidification is presented which simulates the formation of segregation models known as 'freckles' during directional solidification of binary alloys. The growth of the two-phase or dendritic zone is calculated by solving the coupled equations of momentum, energy, and solute transport, as well as maintaining the thermodynamic constraints dictated by the phase diagram of the alloy. Calculations for lead-tin alloys show that the thermosolutal convection in the dendritic zone during solidification can produce heavily localized inhomogeneities in the composition of the final alloy.

Felicelli, Sergio D.↗
NASA NTRS1984

Applications of numerical optimization methods to helicopter design problems: A survey

A survey of applications of mathematical programming methods is used to improve the design of helicopters and their components. Applications of multivariable search techniques in the finite dimensional space are considered. Five categories of helicopter design problems are considered: (1) conceptual and preliminary design, (2) rotor-system design, (3) airframe structures design, (4) control system design, and (5) flight trajectory planning. Key technical progress in numerical optimization methods relevant to rotorcraft applications are summarized.

Miura, H.↗
NASA NTRS1985

Applications of numerical optimization methods to helicopter design problems - A survey

A survey of applications of mathematical programming methods is used to improve the design of helicopters and their components. Applications of multivariable search techniques in the finite dimensional space are considered. Five categories of helicopter design problems are considered: (1) conceptual and preliminary design, (2) rotor-system design, (3) airframe structures design, (4) control system design, and (5) flight trajectory planning. Key technical progress in numerical optimization methods relevant to rotorcraft applications are summarized.

Miura, H.↗
NASA NTRS1984

Applications of numerical optimization methods to helicopter design problems - A survey

A survey of applications of mathematical programming methods is used to improve the design of helicopters and their components. Applications of multivariable search techniques in the finite dimensional space are considered. Five categories of helicopter design problems are considered: (1) conceptual and preliminary design, (2) rotor-system design, (3) airframe structures design, (4) control system design, and (5) flight trajectory planning. Key technical progress in numerical optimization methods relevant to rotorcraft applications are summarized.

Miura, H.↗
NASA NTRS1980

Numerical solution techniques for unsteady transonic aerodynamics problems

Basic concepts of finite difference solution techniques for unsteady transonic flows are presented. The hierarchy of mathematical forumulations that approximate the Navier-Stokes equations are reviewed. The basic concepts involved in constructing numerical algorthms to solve these formulations are given. Semi-implicit and implicit schemes are constructed and analyzed. The discussion focuses primarily on techniques for solving the low frequency transonic small disturbance equation. This is the simplest formulation that contains the essence of inviscid unsteady transonic flow physics. The low frequency formulation is emphasized here because codes based on this theory can be run in minutes of processor time on currently available computers. Furthermore, numerical techniques involved in solving this simple formulation also apply to the more complicated formulations. Extensions to these formulations are briefly described. An indication of the present capability for solving unsteady transonic flows is provided. Important areas of future research for the advancement of computational unsteady transonic aerodynamics are described.

Ballhaus, W. F.↗
NASA NTRS1984

Intersection of parametric surfaces using lookup tables

When primitive structures in the form of parametric surfaces are combined and modified interactively to form complex intersecting surfaces, it becomes important to find the curves of intersection. One must distinguish between finding the shape of the intersection curve, which may only be useful for display purposes, and finding an accurate mathematical representation of the curve, which is important for any meaningful geometric modeling, analysis, design, or manufacturing involving the intersection. The intersection curve between two or more parametric surfaces is important in a variety of computer-aided design and manufacture areas. A few examples are shape design, analysis of groins, design of fillets, and computation of numerically controlled tooling paths. The algorithm presented here provides a mathematical representation of the intersection curve to a specified accuracy. It also provides the database that can simplify operations such as hidden-surface removal, surface rendering, profile identification, and interference or clearance computations.

Hanna, S. L.↗
DOE OSTI2025

Mathematical Foundation for Quantum Computing of Electromagnetic Wave Propagation in Dielectric Media

Can quantum computers effectively simulate the propagation and scattering of electromagnetic waves in a classical plasma? This chapter introduces some of the basic concepts in mathematics and physics essential to answering that question. The numerical simulations of Maxwell equations for wave propagation in dielectrics are constrained by technological limitations of the present-day computers. In contrast, there has been ample fanfare around quantum computers and their potential to far exceed the performance of traditional computers. Whether the enhanced capabilities of a quantum computer can be put to use for simulating topics in classical physics is a source of intrigue and curiosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗