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 487 records · Page 27

Wavenumber selection for small-wavelength Goertler vortices in curved channel flows

The problem of wavenumber selection for fully nonlinear, small-wavelength Goertler vortices in a curved channel flow is considered. These types of Goertler vortices were first considered by Hall & Lakin (1988) for an external boundary layer flow. They proved particularly amenable to asymptotic description, it was possible to consider vortices large enough so that the mean flow correction driven by them is as large as the basic state, and this prompted the authors to consider them in a curved channel flow as an initial application of the phase-equation approach to Goertler vortices. This involves the assumption that the phase variable of these Goertler vortices varies on slow spanwise and time scales, then an analysis of both inside and outside the core region, to which vortex activity is restricted, leads to a system of partial differential equations which can be solved numerically for the wavenumber. The authors consider in particular the effect on the wavenumber of the outer channel wall varying on the same slow spanwise scale as the phase variable.

Dando, Andrew↗

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↗

Nonlinear Waves and Particle Heating in the Auroral Ionosphere

In the two year period, the researchers have developed innovative 2-D simulation codes for modeling the interaction between electron beams and plasma waves and particles in the auroral ionosphere. These new simulations include kinetic wave particles as well as wave-wave interactions near and above approximately 500 km. The principal numerical model consists of partial-differential equations which evolve large amplitude magnetized 2-D Langmuir waves self-consistently with the electron distribution function. Progress was also made in developing 2-D Particle in Cell (PIC) codes.

Goldman, Martin↗

Cellerator: extending a computer algebra system to include biochemical arrows for signal transduction simulations

Cellerator describes single and multi-cellular signal transduction networks (STN) with a compact, optionally palette-driven, arrow-based notation to represent biochemical reactions and transcriptional activation. Multi-compartment systems are represented as graphs with STNs embedded in each node. Interactions include mass-action, enzymatic, allosteric and connectionist models. Reactions are translated into differential equations and can be solved numerically to generate predictive time courses or output as systems of equations that can be read by other programs. Cellerator simulations are fully extensible and portable to any operating system that supports Mathematica, and can be indefinitely nested within larger data structures to produce highly scaleable models.

Non-NASA Center↗

Calculation of the Full Scattering Amplitude without Partial Wave Decomposition: Inclusion of Exchange - 2

The development of a practical method of accurately calculating the full scattering amplitude, without making a partial wave decomposition is continued. The method is developed in the context of electron-hydrogen scattering, and here exchange is dealt with by considering e-H scattering in the static exchange approximation. The Schroedinger equation in this approximation can be simplified to a set of coupled integro-differential equations. The equations are solved numerically for the full scattering wave function. The scattering amplitude can most accurately be calculated from an integral expression for the amplitude; that integral can be formally simplified, and then evaluated using the numerically determined wave function. The results are essentially identical to converged partial wave results.

Shertzer, Janine↗

Application of the Finite Element Method in Atomic and Molecular Physics

The finite element method (FEM) is a numerical algorithm for solving second order differential equations. It has been successfully used to solve many problems in atomic and molecular physics, including bound state and scattering calculations. To illustrate the diversity of the method, we present here details of two applications. First, we calculate the non-adiabatic dipole polarizability of Hi by directly solving the first and second order equations of perturbation theory with FEM. In the second application, we calculate the scattering amplitude for e-H scattering (without partial wave analysis) by reducing the Schrodinger equation to set of integro-differential equations, which are then solved with FEM.

Shertzer, Janine↗

Badhwar-O'Neill 2014 Galactic Cosmic Ray Flux Model Description

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) model is based on GCR measurements from particle detectors. The model has mainly been used by NASA to certify microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BON14 model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model also incorporates an empirical time delay function to account for the lag of the solar activity to reach the boundary of the heliosphere. This technical paper describes the most recent improvements in parameter fits to the BON model (BON14). Using a comprehensive measurement database, it is shown that BON14 is significantly improved over the previous version, BON11.

O'Neill, P. M.↗

NASA Galactic Cosmic Radiation Environment Model: Badhwar - O'Neill (2014)

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) flux model has been used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a comprehensive database of GCR measurements from various particle detectors to determine boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

Golge, S.↗

Thermal history and evolution of the moon.

In this work, theoretical lunar temperature models are computed taking into account different initial conditions to represent possible accretion models and various abundances of heat sources to correspond to different compositions. Differentiation and convection are simulated in the numerical computational scheme. Models of the thermal evolution of the moon that fit the chronology of igneous activity on the lunar surface, the stress history of the lunar lithosphere implied by the presence of mascons, and the surface concentrations of radioactive elements, involve extensive differentiation early in lunar history. This differentiation may be the result of rapid accretion and large-scale melting or of primary chemical layering during accretion. Differences in present-day temperatures for these two possibilities are significant only in the inner 1000 km of the moon and are not resolvable with presently available data.

Toksoz, M. N.↗

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗

Chimpanzee counting and rhesus monkey ordinality judgments

An investigation is conducted to address the questions of whether chimpanzees can count and whether rhesus monkeys can differentiate written numbers. One investigation demonstrates the capacity of a chimpanzee to produce a quantity of responses appropriate to a given Arabic numeral. Rhesus monkeys are shown to have the capability for making fine differentiations between quantities of pellets and Arabic numerals.

Rumbaugh, Duane M.↗

Processing CCD Images to Detect Transits of Earth-Sized Planets: Maximizing Sensitivity While Achieving Reasonable Downlink Requirements

We have performed end-to-end laboratory and numerical simulations to demonstrate the capability of differential photometry under realistic operating conditions to detect transits of Earth-sized planets orbiting solar-like stars. Data acquisition and processing were conducted using the same methods planned for the proposed Kepler Mission. These included performing aperture photometry on large-format CCD images of an artificial star fields obtained without a shutter at a readout rate of 1 megapixel/sec, detecting and removing cosmic rays from individual exposures and making the necessary corrections for nonlinearity and shutterless operation in the absence of darks. We will discuss the image processing tasks performed `on-board' the simulated spacecraft, which yielded raw photometry and ancillary data used to monitor and correct for systematic effects, and the data processing and analysis tasks conducted to obtain lightcurves from the raw data and characterize the detectability of transits. The laboratory results are discussed along with the results of a numerical simulation carried out in parallel with the laboratory simulation. These two simulations demonstrate that a system-level differential photometric precision of 10-5 on five- hour intervals can be achieved under realistic conditions.

Earth-size planets↗

Numerical Stability In Hyperbolic Boundary-Value Problems

Technical memorandum discusses stability of numerical solutions involving semidiscrete approximations to hyperbolic partial differential equations in initial-and-boundary-value problems. Topic of practical significance because hyperbolic partial differential equations arise in mathematical modeling of waves and blasts. Solutions often needed over restricted regions of unbounded spaces. Outer boundaries artificial, introduced only to limit domains of numerical solutions. Conditions at such artificial boundaries cause numerical instabilities that degrade computed solutions.

Warming, Robert F.↗

Finite-analytic numerical solution of heat transfer in two-dimensional cavity flow

Heat transfer in cavity flow is numerically analyzed by a new numerical method called the finite-analytic method. The basic idea of the finite-analytic method is the incorporation of local analytic solutions in the numerical solutions of linear or nonlinear partial differential equations. In the present investigation, the local analytic solutions for temperature, stream function, and vorticity distributions are derived. When the local analytic solution is evaluated at a given nodal point, it gives an algebraic relationship between a nodal value in a subregion and its neighboring nodal points. A system of algebraic equations is solved to provide the numerical solution of the problem. The finite-analytic method is used to solve heat transfer in the cavity flow at high Reynolds number (1000) for Prandtl numbers of 0.1, 1, and 10.

Chen, C.-J.↗

Accuracy of the domain method for the material derivative approach to shape design sensitivities

Numerical accuracy for the boundary and domain methods of the material derivative approach to shape design sensitivities is investigated through the use of mesh refinement. The results show that the domain method is generally more accurate than the boundary method, using the finite element technique. It is also shown that the domain method is equivalent, under certain assumptions, to the implicit differentiation approach not only theoretically but also numerically.

Yang, R. J.↗

A Functional Interpolation Approach to Compute Period Orbits in the Circular Restricted Three-body Problem

In this paper, we develop a method to solve for periodic orbits, i.e., Lyapunov and Halo orbits, using a functional interpolation scheme called the Theory of Func- tional Connections (TFC). Using this technique, a periodic constraint is analyti- cally embedded into the TFC constrained expression. By doing this, the system of differential equations governing the three-body problem is transformed into an unconstrained optimization problem where simple numerical schemes can be used to find a solution, e.g., nonlinear least-squares is used. This allows for a simpler numerical implementation with comparable accuracy and speed to the traditional differential corrector method.

Mortari, Daniele↗