Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical differentiation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗
NASA NTRSDate not supplied

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗
NASA NTRS2022

Geodetic Investigations of the Mission Concept MAGIC to Reveal Callisto’s Internal Structure

Geodetic and geophysical investigations of the Galilean moon Callisto address fundamental questions regarding the formation and evolution of the Jovian system. Callisto's evolution and internal structure appear to significantly differ from the other Jovian satellites. Similarly-sized Ganymede is a highly evolved ice-rock moon with a differentiated interior, intrinsic magnetic field, and abundant surface evidence of internal activity. In contrast, Callisto's surface is ancient, and Galileo spacecraft data suggest its interior is only incompletely differentiated, despite the presumed presence of a sub-surface ocean. These properties make Callisto uniquely able to constrain the timing and nature of the Jovian system formation. The Magnetics, Altimetry, Gravity, and Imaging of Callisto (MAGIC) mission concept is conceived to fully characterize the properties of this enigmatic moon from its deep interior to the icy shell. Three main instruments are included as a scientific payload. Highly accurate measurements of Callisto's topography, magnetic field, and morphology are obtained by the onboard laser altimeter, magnetometer, and camera, respectively. The telecommunication system supports an additional gravity and radio science investigation. Long- and short-wavelength gravity anomalies afford powerful constraints on internal differentiation and the properties of the hydrosphere (water and ice). Comprehensive numerical simulations and covariance analyses of MAGIC mission scenarios presented in this paper show that the gravitational degree-2 normalized coefficients and the pole obliquity enable the determination of the moment of inertia with an accuracy better than 0.015%. The combination of gravity and altimetry measurements acquired by MAGIC are essential to the characterization of Callisto's interior if – as is likely – the degree-2 gravity includes non-hydrostatic terms. MAGIC's radio science data yield the estimation of Callisto's gravity field with spatial resolutions of <100 km. The combination of gravitational and deformation tides that are retrieved by the radio science and altimetry investigations, respectively, leads to the recovery of the rigid ice shell thickness to within ∼3 km. Together these datasets would resolve ambiguities inherent in Galileo flyby data, revealing Callisto's interior structure as well as the existence and properties of its postulated internal ocean.

Antonio Genova↗
NASA NTRS1995

Chemical and Thermal Analysis

During the past six months we have conducted significant research in several domains in order to clarify and understanding the aging and chemical failure mechanism of thermoplastics (PVDF or Tefzel) for pipes. We organized numerous analytical studies with methods including Fourier Transform Infrared Spectroscopy, Dynamic Mechanical Analysis, Differential Scanning Calorimetry, and Stress Relaxation experiments. In addition we have reanalyzed previous thermogravimetric data concerning the rate of deplasticization of Coflon pipe. We investigated a number of aged samples of both Tefzel and Coflon that were forwarded from MERL. We conducted stress relaxation experiments of Coflon pipe at several temperatures and determined an activation energy. We also examined the dynamic mechanical response PVDF during deplasticization and during methanol plasticization. We performed numerous DSC analyses to research the changing crystalline morphology. We have noted significant changes in crystallinity upon aging for both PVDF and Tefzel. Little variation in elemental composition was noted for many of the aged Coflon and Tefzel samples tested.

Bulluck, J. W.↗
NASA NTRS2009

Stability of Detached Solidification

Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.

Mazuruk, K.↗
NASA NTRS1964

Analytical Solutions for Non-Linear Differential Equations with the Help of a Digital Computer

A technique was developed with the help of a digital computer for analytic (algebraic) solutions of autonomous and nonautonomous equations. Two operational transform techniques have been programmed for the solution of these equations. Only relatively simple nonlinear differential equations have been considered. In the cases considered it has been possible to assimilate the secular terms into the solutions. For cases where f(t) is not a bounded function, a direct series solution is developed which can be shown to be an analytic function. All solutions have been checked against results obtained by numerical integration for given initial conditions and constants. It is evident that certain nonlinear differential equations can be solved with the help of a digital computer.

Cromwell, P. C.↗
NASA NTRS1972

Analytical treatment of gas flows through multilayer insulation, project 1

A theoretical investigation of gas flow inside a multilayer insulation system was made for the case of the broadside pumping process. A set of simultaneous first-order differential equations for the temperature and pressure of the gas molecules through the perforations on the insulation layers. A modified Runge-Kutta method was used for numerical experiment. The numerical stability problem was also investigated. It was shown that when the relaxation time is less than the time period over which the gas properties change appreciably, the set of differential equations can be replaced by a set of algebraic equations for solution. Numerical examples were given and comparisons with experimental data were made.

Lin, J. T.↗
NASA NTRS1996

Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

Mickens, Ronald E.↗
NASA NTRS2006

Constraining Numerical Geodynamo Modeling with Surface Observations

Numerical dynamo solutions have traditionally been generated entirely by a set of self-consistent differential equations that govern the spatial-temporal variation of the magnetic field, velocity field and other fields related to dynamo processes. In particular, those solutions are obtained with parameters very different from those appropriate for the Earth s core. Geophysical application of the numerical results therefore depends on correct understanding of the differences (errors) between the model outputs and the true states (truth) in the outer core. Part of the truth can be observed at the surface in the form of poloidal magnetic field. To understand these differences, or errors, we generate new initial model state (analysis) by assimilating sequentially the model outputs with the surface geomagnetic observations using an optimal interpolation scheme. The time evolution of the core state is then controlled by our MoSST core dynamics model. The final outputs (forecasts) are then compared with the surface observations as a means to test the success of the assimilation. We use the surface geomagnetic data back to year 1900 for our studies, with 5-year forecast and 20-year analysis periods. We intend to use the result; to understand time variation of the errors with the assimilation sequences, and the impact of the assimilation on other unobservable quantities, such as the toroidal field and the fluid velocity in the core.

Kuang, Weijia↗
NASA NTRS1981

Experimental validation of a numerical model predicting the charging characteristics of Teflon and Kapton under electron beam irradiation

The effect differential charging of spacecraft thermal control surfaces is assessed by studying the dynamics of the charging process. A program to experimentally validate a computer model of the charging process was established. Time resolved measurements of the surface potential were obtained for samples of Kapton and Teflon irradiated with a monoenergetic electron beam. Results indicate that the computer model and experimental measurements agree well and that for Teflon, secondary emission is the governing factor. Experimental data indicate that bulk conductivities play a significant role in the charging of Kapton.

Hazelton, R. C.↗
NASA NTRS2008

Simulink Model of the Ares I Upper Stage Main Propulsion System

A numerical model of the Ares I upper stage main propulsion system is formulated based on first principles. Equation's are written as non-linear ordinary differential equations. The GASP fortran code is used to compute thermophysical properties of the working fluids. Complicated algebraic constraints are numerically solved. The model is implemented in Simulink and provides a rudimentary simulation of the time history of important pressures and temperatures during re-pressurization, boost and upper stage firing. The model is validated against an existing reliable code, and typical results are shown.

Burchett, Bradley T.↗
NASA NTRS1981

Advantages of multigrid methods for certifying the accuracy of PDE modeling

Numerical techniques for assessing and certifying the accuracy of the modeling of partial differential equations (PDE) to the user's specifications are analyzed. Examples of the certification process with conventional techniques are summarized for the three dimensional steady state full potential and the two dimensional steady Navier-Stokes equations using fixed grid methods (FG). The advantages of the Full Approximation Storage (FAS) scheme of the multigrid technique of A. Brandt compared with the conventional certification process of modeling PDE are illustrated in one dimension with the transformed potential equation. Inferences are drawn for how MG will improve the certification process of the numerical modeling of two and three dimensional PDE systems. Elements of the error assessment process that are common to FG and MG are analyzed.

Forester, C. K.↗
NASA NTRS1939

Some notes on the numerical solution of shear-lag and mathematically related problems

The analysis of box beams with shear deformation of the flanges can be reduced to the solution of a differential equation. The same equation is met in other problems of stress analysis. No analytical solutions of this equation can be given for practical cases, and numerical methods of evaluation must be used. Available methods are briefly discussed. Two numerical examples show the application of the step-by-step method of integration to shear-lag problems.

Kuhn, Paul↗
NASA NTRS2023

DDFRG2: Double-Differential FRaGmentation Model with Pion Production

A set of Double-Differential FRaGmentation (DDFRG) models for production of fragments from nucleon-nucleus and nucleus-nucleus collisions, relevant to space radiation, is currently being developed. DDFRG1 was a previous model developed for production of protons and light ions. A new model for pion production is developed in the present work and is called DDFRG2, which is based upon thermal production of pions from a central fireball source. The Lorentzinvariant double-differential cross section is calculated in the central fireball frame, and is Lorentz transformed to the laboratory frame, which results in a closed-form analytic formula involving no numerical integration, and which can be run very efficiently in radiation transport codes. The Lorentz-invariant double-differential cross section is then integrated over angle to give the singledifferential energy spectral distribution. Two versions of the DDFRG2 pion production are developed, namely DDFRG2-SIMPLE and DDFRG2-FULL. Both versions produce pion production Lorentz-invariant double-differential cross sections and also pion production single-differential spectral (energy) distributions. The SIMPLE version contains only one adjustable parameter. The FULL version contains three additional adjustable parameters. Both DDFRG2 versions are compared to over 1,000 data points. The SIMPLE version describes the p + 12 C → π − much better than the current HZETRN pion model. For other reactions, the SIMPLE version is of comparable quality to the current HZETRN pion model compared to data. For both the p + 12 C → π − reaction and all other reactions, DDFRG2-FULL describes the data much more accurately than the current HZETRN pion model.

John W. Norbury↗
NASA NTRS2024

DDFRG3: Double-Differential FRaGmentation model for Neutron production

A Double-Differential FRaGmentation (DDFRG) model for neutron production from nucleus - nucleus collisions is developed for space radiation applications. DDFRG1 was a previous model developed for production of protons and light ions, and DDFRG2 was a model developed for pion production. A new model (DDFRG3) for neutron production is developed in the present work, and is based upon thermal production of neutrons from the projectile, the target and central fireball sources. The Lorentz-invariant double-differential cross sections are calculated in the various source frames, and are Lorentz transformed to the laboratory frame, resulting in a closed-form analytic formula involving no numerical integration, and which can be run very efficiently in radiation transport codes. The Lorentz-invariant double-differential cross section is then integrated over angle to give the single-differential energy spectral distribution. The DDFRG3 neutron model compares very well to an extensive experimental data set.

neutrons↗
NASA NTRS1986

Parameter estimation problems for distributed systems using a multigrid method

The problem of estimating spatially varying coefficients of partial differential equations is considered from observation of the solution and of the right hand side of the equation. It is assumed that the observations are distributed in the domain and that enough observations are given. A method of discretization and an efficient multigrid method for solving the resulting discrete systems are described. Numerical results are presented for estimation of coefficients in an elliptic and a parabolic partial differential equation.

Taasan, S.↗