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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A new technique for solving the Parker-type wind equations

Substitution of the novel function Phi for velocity, as one of the dependent variables in Parker-type solar wind equations, removes the critical point, and therefore the numerical difficulties encountered, from the set of coupled differential wind equations. The method has already been successfully used in a study of radiatively-driven mass loss from the surface of X-ray bursting neutron stars. The present technique for solving the equations of time-independent mass loss can be useful in similar applications.

Melia, Fulvio↗
NASA NTRS1990

A parallel iterative solution method for systems of nonlinear hyperbolic equations

An iterative algorithm suitable for the solution of a system of nonlinear hyperbolic partial differentiation equations in multiple dimensions is discussed. Current numerical methods for systems of nonlinear PDEs have limited parallelism due to strong coupling between the equations. This method decouples the PDEs by linearizing the convention coefficient for a space-time domain. This provides large grain parallelism. The linearization also allows the treatment of some terms in the equations as source terms, providing more freedom to choose from a wider variety of numerical methods. Smaller grain parallelism may be exploited within the solves for each equation. Thus, the method has potential for parallelism at several levels.

Scroggs, Jeffrey S.↗
NASA NTRS1991

Heat Transfer to a Thin Solid Combustible in Flame Spreading at Microgravity

The heat transfer rate to a thin solid combustible from an attached diffusion flame, spreading across the surface of the combustible in a quiescent, microgravity environment, was determined from measurements made in the drop tower facility at NASA-Lewis Research Center. With first-order Arrhenius pyrolysis kinetics, the solid-phase mass and energy equations along with the measured spread rate and surface temperature profiles were used to calculate the net heat flux to the surface. Results of the measurements are compared to the numerical solution of the complete set of coupled differential equations that describes the temperature, species, and velocity fields in the gas and solid phases. The theory and experiment agree on the major qualitative features of the heat transfer. Some fundamental differences are attributed to the neglect of radiation in the theoretical model.

Bhattacharjee, S.↗
NASA NTRS2001

Left ventricular endocardial surface detection based on real-time 3D echocardiographic data

OBJECTIVE: A new computerized semi-automatic method for left ventricular (LV) chamber segmentation is presented. METHODS: The LV is imaged by real-time three-dimensional echocardiography (RT3DE). The surface detection model, based on level set techniques, is applied to RT3DE data for image analysis. The modified level set partial differential equation we use is solved by applying numerical methods for conservation laws. The initial conditions are manually established on some slices of the entire volume. The solution obtained for each slice is a contour line corresponding with the boundary between LV cavity and LV endocardium. RESULTS: The mathematical model has been applied to sequences of frames of human hearts (volume range: 34-109 ml) imaged by 2D and reconstructed off-line and RT3DE data. Volume estimation obtained by this new semi-automatic method shows an excellent correlation with those obtained by manual tracing (r = 0.992). Dynamic change of LV volume during the cardiac cycle is also obtained. CONCLUSION: The volume estimation method is accurate; edge based segmentation, image completion and volume reconstruction can be accomplished. The visualization technique also allows to navigate into the reconstructed volume and to display any section of the volume.

Non-NASA Center↗
NASA NTRS2000

Kinetics of DSB rejoining and formation of simple chromosome exchange aberrations

PURPOSE: To investigate the role of kinetics in the processing of DNA double strand breaks (DSB), and the formation of simple chromosome exchange aberrations following X-ray exposures to mammalian cells based on an enzymatic approach. METHODS: Using computer simulations based on a biochemical approach, rate-equations that describe the processing of DSB through the formation of a DNA-enzyme complex were formulated. A second model that allows for competition between two processing pathways was also formulated. The formation of simple exchange aberrations was modelled as misrepair during the recombination of single DSB with undamaged DNA. Non-linear coupled differential equations corresponding to biochemical pathways were solved numerically by fitting to experimental data. RESULTS: When mediated by a DSB repair enzyme complex, the processing of single DSB showed a complex behaviour that gives the appearance of fast and slow components of rejoining. This is due to the time-delay caused by the action time of enzymes in biomolecular reactions. It is shown that the kinetic- and dose-responses of simple chromosome exchange aberrations are well described by a recombination model of DSB interacting with undamaged DNA when aberration formation increases with linear dose-dependence. Competition between two or more recombination processes is shown to lead to the formation of simple exchange aberrations with a dose-dependence similar to that of a linear quadratic model. CONCLUSIONS: Using a minimal number of assumptions, the kinetics and dose response observed experimentally for DSB rejoining and the formation of simple chromosome exchange aberrations are shown to be consistent with kinetic models based on enzymatic reaction approaches. A non-linear dose response for simple exchange aberrations is possible in a model of recombination of DNA containing a DSB with undamaged DNA when two or more pathways compete for DSB repair.

DNA Damage↗
NASA NTRS1993

Issues in the Design of a Multistep Code

This paper summarizes things the author considers important from his past work, and presents some new ideas for use in the implementation of multistep methods. An effort has been made to include material of interest to those of a theoretical bent, as well as to those whose primary interest is in the implementation of methods.

Numerical solution of ordinary differential equati↗
NASA NTRS2018

Atmospheric Chemistry Modeling Using a Regression Forest Model

Atmospheric chemistry is central to many environmental issues such as air pollution, climate change, and stratospheric ozone loss. Chemistry Transport Models (CTM) are a central tool for understanding these issues, whether for research or for forecasting. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O1D to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a CTM. We have investigated a machine learning approach to solving the differential equations instead of solving them numerically. From an annual simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry. This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, centennial scale climate simulations etc. We discuss our approches' speed and accuracy, and highlight some potential future directions for improving this approach.

Keller, Christoph A.↗
NASA NTRS1998

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: The Two Dimensional Time Marching Schemes - 1

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly.

Chang, Sin-Chung↗
NASA NTRS1994

Gradient drift eigenmodes in the equatorial electrojet

The problem of kilometer-scale irregularities in the daytime equatorial electrojet is revisited by means of an eigenmode analysis of the gradient drift instability. Realistic physical parameters are used, including the modeled altitude variations of ion and electron collision frequencies and mobilities. The full fourth-order system of two coupled differential equations (each of second order) for the denisty and electrostatic potential perturbations is solved numerically by a relaxation technique. Under some approximations, the fourth-order system can be shown to reduce to a second-order differential equation for the perturbed potential or density. The latter is solved using a shooting technique and provides initial guesses for numerical solutions to the full problem. It is shown that the linear growth rate peaks for kilometer-scale waves, contrary to the findings of recent initial-value studies. This occurs because the equilibrium velocity shear is much more effective as a damping mechanism for short-wavelength modes than it is for the longer, kilometer-scale modes. These results provide a natural qualitative explanation for the observed dominance of kilometer-scale structures in the daytime electrojet spectrum.

Wang, X.-H.↗
NASA NTRS1976

A differential game solution to the Coplanar tail-chase aerial combat problem

Numerical results obtained in a simplified version of the one on one aerial combat problem are presented. The primary aim of the data is to specify the roles of pursuer and evader as functions of the relative geometry and of the significant physical parameters of the problem. Numerical results are given in a case in which the slower aircraft is more maneuverable than the faster aircraft. A third order dynamic model of the relative motion is described, for which the state variables are relative range, bearing, and heading. The ranges at termination are arbitary in the present version of the problem, so the weapon systems of both aircraft can be visualized as forward firing high velocity weapons, which must be aimed at the tail pipe of the evader. It was found that, for the great majority of the ralative geometries, each aircraft can evade the weapon system of the other.

Merz, A. W.↗
NASA NTRS1982

Recent advances in electromagnetic synthesis and analysis of dual-shaped reflector antennas

Dual-shaped reflectors have been used for many years. Thus, these reflectors have been used as high gain antennas on Voyagers 1 and 2. The objectives of the geometrical optics (GO) dual shaped synthesis are considered. Concerning the synthesis of dual shaped reflectors, it has been shown for circular symmetric reflectors that an exact GO solution can be found to the problem of transforming, by two reflections, any feed pattern into any aperture distribution. This problem involves solving two simultaneous nonlinear ordinary differential equations. The same approach for offset geometry leads to two simultaneous nonlinear partial differential equations. It is shown that these equations could also be integrated numerically, except that in general these equations are not total and therefore, in general, they do not have a 'smooth' solution. It is further shown that the offset partial differentials often very nearly form a total differential in many cases of practical importance.

Galindo-Israel, V.↗
NASA NTRS1973

Analysis of free turbulent shear flows by numerical methods

Studies are described in which the effort was essentially directed to classes of problems where the phenomenologically interpreted effective transport coefficients could be absorbed by, and subsequently extracted from (by comparison with experimental data), appropriate coordinate transformations. The transformed system of differential equations could then be solved without further specifications or assumptions by numerical integration procedures. An attempt was made to delineate different regimes for which specific eddy viscosity models could be formulated. In particular, this would account for the carryover of turbulence from attached boundary layers, the transitory adjustment, and the asymptotic behavior of initially disturbed mixing regions. Such models were subsequently used in seeking solutions for the prescribed two-dimensional test cases, yielding a better insight into overall aspects of the exchange mechanisms.

Korst, H. H.↗
NASA NTRS1994

Post-Gondwana geomorphic evolution of southwestern Africa: Implications for hte controls on landscape development from observations and numerical experiments

The relationship between morphology and surficial geology is used to quantify the denudation that has occurred across southwestern Africa sicne the fragmentation of Gondwana during the Early Mesozoic. Two main points emerge. Signficant denudation, of the order of kilometers, is widespread except in the Kalahari region of the continental interior. The denudation is systematically distributed so that the continental exterior catchment, draining directly to the Cape basin, is denuded to a greater depth than the interior catchment inland of the Great Escarpment. The analysis also implies tha the majority of the denudation occurred before the beginning of the Cenozoic for both teh exerior and interior catchments. Existing models of landscape development are reviewed, and implications of the denudation chronology are incorporated into a revised conceptual model. This revision implies tha thte primary effect of rifting on the subsequent landscape evolution is that it generates two distinct drainage regimes. A marginal upwarp, or rift flank uplift, separates rejuvenated rivers that drain into the subsiding rift from rivers in the continetal interior that are deflected but not rejuvenated. The two catchments evolve independently unless they are integrated by breaching of hte marginal upwarp. If this occurs, the exterior baselevel is communicated to the interior catchment that is denuded accordingly. Denudation rates generally decrease as the margin evolves, and this decrease is reinforced by the exposure of substrate that is resistant to denudation and/or a change to a more arid climate. The observations do not reveal a particular style of smaller-scale landscape evolution, sucha s escarpment retreat, that is responsible for the differential denudation across the region. It is proposed that numerical model experiments, which reflect the observational insights at the large scale, may identify the smaller-scale controls on escarpment development if the model and natural systems are analogous. Four numerical experiments are presented in which the roles of antecedent topography, resistant substrate, climte change, and lowering the baselevel of the interior catchment are investigated for an initially high elevation margin bordered by an escarpment. The model results suggest several styles of landscape evolution that are compatible with the observations. Escarpments may retreat in a regular manner, but they also degrade and are destroyed, only to reform at the drainage divide between exterior and interior catchments.

Gilchrist, Alan R.↗
NASA NTRS1980

An interaction algorithm for prediction of mean and fluctuating velocities in two-dimensional aerodynamic wake flows

A theoretical analysis is presented yielding sets of partial differential equations for determination of turbulent aerodynamic flowfields in the vicinity of an airfoil trailing edge. A four phase interaction algorithm is derived to complete the analysis. Following input, the first computational phase is an elementary viscous corrected two dimensional potential flow solution yielding an estimate of the inviscid-flow induced pressure distribution. Phase C involves solution of the turbulent two dimensional boundary layer equations over the trailing edge, with transition to a two dimensional parabolic Navier-Stokes equation system describing the near-wake merging of the upper and lower surface boundary layers. An iteration provides refinement of the potential flow induced pressure coupling to the viscous flow solutions. The final phase is a complete two dimensional Navier-Stokes analysis of the wake flow in the vicinity of a blunt-bases airfoil. A finite element numerical algorithm is presented which is applicable to solution of all partial differential equation sets of inviscid-viscous aerodynamic interaction algorithm. Numerical results are discussed.

Baker, A. J.↗