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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A numerical study of the effects of a large sandbar upon sea breeze development

Two-dimensional numerical simulations of sea breeze development over a large sandbar on the North Sea coast of Germany are reported. The numerical model used in these experiments contains a detailed treatment of soil moisture, which allows evaluation of the effects of differential surface characteristics on the airflow pattern. Results of the simulations indicate that the contrast between the moist sandbar and adjacent dry land, the tidal inundation of the sandbar, and the westward penetration of the Baltic sea breeze play important roles in the development of mesoscale airflow patterns in the sandbar region.

Kessler, R. C.↗
NASA NTRS1986

Concurrent processing in nonlinear structural stability

A concurrent processing algorithm is developed for materially nonlinear stability analysis of imperfect columns with biaxial partial rotational end restraints. The algorithm for solving the governing nonlinear ordinary differential equations is implemented on a multiprocessor computer called the 'Finite Element Machine', developed at the NASA Langley Research Center. Numerical results are obtained on up to nine concurrent processors. A substantial computational gain is achieved in using the parallel processing approach.

Darbhamulla, S. P.↗
NASA NTRS1987

The origin of rigidly rotating magnetic field patterns on the sun

Using analytical calculations and numerical simulations, it is shown that a meridional component of magnetic-flux transport will offset the shearing effect of differential rotation and give rise to rigidly rotating patterns of large-scale magnetic field. The nonaxisymmetric field attains a striped polarity pattern which rotates rigidly like a barber pole while its individual small-scale flux elements rotate at the differential rate of the latitudes they are crossing. On the sun, the meridional transport is provided by supergranular diffusion possibly assisted by a small poleward flow. New sources of flux retard this process and exclude the rigid rotation from the sunspot belts until well into the declining phase of the sunspot cycle. This mechanism accounts for a number of heretofore unexplained phenomena including the tendency for coronal holes to rotate rigidly during the declining phase of the sunspot cycle.

Sheeley, N. R., Jr.↗
NASA NTRS1990

A numerical simulation of the microburst phenomenon - Comparison of a discrete and a parameterized microphysical representation

Two versions of a microburst model that differ in how they represent microphysical processes are described and comparison of the numerical solutions is made. This study focuses on how the versions differ in representing the evaporation of raindrops and the differential sedimentation of drops that develops because of their fallspeeds. The two versions are then compared to demonstrate how predicted microburst characteristics are influenced by the choice of microphysical scheme.

Coen, Janice L.↗
NASA NTRS1993

Accurate computation and continuation of homoclinic and heteroclinic orbits for singular perturbation problems

In earlier papers, Doedel and the authors have developed a numerical method and derived error estimates for the computation of branches of heteroclinic orbits for a system of autonomous ordinary differential equations in R(exp n). The idea of the method is to reduce a boundary value problem on the real line to a boundary value problem on a finite interval by using a local (linear or higher order) approximation of the stable and unstable manifolds. A practical limitation for the computation of homoclinic and heteroclinic orbits has been the difficulty in obtaining starting orbits. Typically these were obtained from a closed form solution or via a homotopy from a known solution. Here we consider extensions of our algorithm which allow us to obtain starting orbits on the continuation branch in a more systematic way as well as make the continuation algorithm more flexible. In applications, we use the continuation software package AUTO in combination with some initial value software. The examples considered include computation of homoclinic orbits in a singular perturbation problem and in a turbulent fluid boundary layer in the wall region problem.

Vaughan, William W.↗
NASA NTRS1993

An analytic study of a two-phase laminar airfoil in simulated heavy rain

A mathematical model for a two-phase flow laminar airfoil in simulated heavy rain has been established. The set of non-linear partial differential equations has been converted into a set of finite difference equations; appropriate initial and boundary conditions are provided. The numerical results are compared with the experimental measurements. They show good agreement in quality.

Hsu, Yu-Kao↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS1996

Chemical and Thermal Analysis

Work has included significant research in several areas aimed at further clarification of the aging and chemical failure mechanism of thermoplastics (PVDF or Tefzel) for pipes. Among the areas investigated were the crystallinity changes associated with both the Coflon and Tefzel after various simulated environmental exposures using X-Ray diffraction analysis. We have found that significant changes in polymer crystallinity levels occur as a function of the exposures. These crystallinity changes may have important consequences on the fracture, fatigue, tensile, and chemical resistance of the materials. We have also noted changes in the molecular weight distribution of the Coflon material using a dual detector Gel Permeation Analysis. Again these changes may result in variation in the mechanical and chemical properties in the material. We conducted numerous analytical studies with methods including X-Ray Diffraction, Gel Permeation Chromatography, Fourier Transform Infrared Spectroscopy, Thermogravimetric Analysis, and Differential Scanning Calorimetry. We investigated a number of aged samples of both Tefzel and Coflon that were forwarded from MERL. Pressurized tests were performed in a modified Fluid G, which we will call G2. In this case the ethylene diamine concentration was increased to 3 percent in methanol. Coflon pipe sections and powdered Coflon were exposed in pressure cells at 1700 psi at three separate test temperatures, 70 C, 110 C, and 130 C. The primary purpose of the pressure tests in Fluid G2 was to further elucidate the aging mechanism of PVDF degradation.

Bulluck, J. W.↗
NASA NTRS2001

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to develop and implement new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. We have just completed two (2) years of Phase I of this research. This annual report provides a brief and up-to-date summary of our activities during the period: September 1, 2000 through August 31, 2001. In the work within the past year, a methodology termed "velocity-scalar filtered density function" (VSFDF) is developed and implemented for large eddy simulation (LES) of turbulent flows. In this methodology the effects of the unresolved subgrid scales (SGS) are taken into account by considering the joint probability density function (PDF) of all of the components of the velocity and scalar vectors. An exact transport equation is derived for the VSFDF in which the effects of the unresolved SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source terms appear in closed form. The remaining unclosed terms in this equation are modeled. A system of stochastic differential equations (SDEs) which yields statistically equivalent results to the modeled VSFDF transport equation is constructed. These SDEs are solved numerically by a Lagrangian Monte Carlo procedure. The consistency of the proposed SDEs and the convergence of the Monte Carlo solution are assessed by comparison with results obtained by an Eulerian LES procedure in which the corresponding transport equations for the first two SGS moments are solved. The unclosed SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source in the Eulerian LES are replaced by corresponding terms from VSFDF equation. The consistency of the results is then analyzed for a case of two dimensional mixing layer.

Drozda, T. G.↗
NASA NTRS2013

Semi-Empirical Validation of the Cross-Band Relative Absorption Technique for the Measurement of Molecular Mixing Ratios

Studies were performed to carry out semi-empirical validation of a new measurement approach we propose for molecular mixing ratios determination. The approach is based on relative measurements in bands of O2 and other molecules and as such may be best described as cross band relative absorption (CoBRA). . The current validation studies rely upon well verified and established theoretical and experimental databases, satellite data assimilations and modeling codes such as HITRAN, line-by-line radiative transfer model (LBLRTM), and the modern-era retrospective analysis for research and applications (MERRA). The approach holds promise for atmospheric mixing ratio measurements of CO2 and a variety of other molecules currently under investigation for several future satellite lidar missions. One of the advantages of the method is a significant reduction of the temperature sensitivity uncertainties which is illustrated with application to the ASCENDS mission for the measurement of CO2 mixing ratios (XCO2). Additional advantages of the method include the possibility to closely match cross-band weighting function combinations which is harder to achieve using conventional differential absorption techniques and the potential for additional corrections for water vapor and other interferences without using the data from numerical weather prediction (NWP) models.

Pliutau, Denis↗
DOE OSTI2023

Leverage Score Sampling for Parametric PDEs (Final Technical Report)

This final technical report summarizes the accomplishments of work performed under DOE Office of Science Award DE-SC0022266, which is titled “Leverage Score Sampling for Parametric PDEs”. The goal of the project was to extend methods from Randomized Numerical Linear Algebra (RandNLA) to tackle central computational challenges in model order reduction and uncertainty quantification (UQ) for parametric partial differential equations (PDEs). In particular, we sought to use importance sampling methods originally developed for RandNLA to develop sample efficient active learning algorithms for approximating high-dimensional scalar functions, e.g. by polynomials, Gaussian process models, and simple neural networks. Such methods can be immediately applied to developing surrogate models or to approximating quantity of interest (QoI) surfaces. In the context of PDEs, each sample used for learning equates to the solution of the differential equation for a particular set of parameters, so sample efficiency translates to improved computational efficiency for a variety of downstream tasks.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2017

Microbiological Horticultural Internship Final Abstract

GMO dwarf plum (Prunus domestica) is being evaluated as a candidate food crop for long duration space flight missions. A project was undertaken to develop a protocol for transferring selected genetic lines of GMO plum (previously maintained in pots and propagated by cuttings at NASA's Kennedy Space Center in Florida) into in vitro tissue culture. In vitro culture may reduce the space, materials, and labor required to maintain the current lines of GMO plum and better preserve them for future study. Fresh plant material from three selected GMO plum lines (NASA-5, NASA-10, and NASA-11) and a non-modified control line (Control-5) were processed aseptically into in vitro culture on four separate occasions. The impact of multiple treatments on the successful growth of GMO plum tissue in vitro were tested: Parent explant tissue type (leaf petioles, stem nodes containing buds and internodes without buds), tissue sterilization method [soaking in 10 bleach only (5 min for petioles or 10 min for nodesinternodes), or soaking in 70 EtOH (30 sec) followed by 10 bleach (5 min for petioles and 10 min for nodesinternodes)], and media type [three Murashige and Skoog-based medias (SGM, SRM, and SRM+2,4-D) and one recipe containing woody plant media (WPM)]. 22.2 of the plates containing tissue sterilized with bleach alone developed microbial contamination after two weeks, while only 11.8 of plates containing tissue sterilized sequentially with EtOH and bleach developed contamination. Node bud tissue from all four genetic lines of plum produced leafy plantlets on SGM and SRM media after 4-6 weeks. The most numerous and well-developed plantlets were present on SGM. Upon reaching suitable size, plantlets were transferred to larger media containers for further growth. Some node bud growth occurred on SRM+2,4-D and WPM 2.5 weeks after plating, however as of yet no pieces on SRM+2,4-D have adequate development for transferring. Tissue pieces from NASA-5 plated on WPM are developing leaves and will be ready for transferring soon. Petioles and internode tissue lacking bud meristem failed to produce any plantlets on any plates, however petioles developed large masses of undifferentiated callus tissue on SRM+2,4-D media. These callused pieces were then transferred to SRM+TDZ media, which resulted in even larger callus growth but no differentiation. All four selected plum lines were successfully transitioned into in vitro culture. Nodes from NASA-5 and NASA-10 lines produced the most numerous and well-developed leafy plantlets in vitro, while those from NASA-11 and Control-5 were generally smaller, slower growing and less numerous. The best method overall was to use young stem node tissue with buds, surface sterilize the pieces sequentially with 70 EtOH and 10 bleach, and then plate them onto SGM media. Future areas of study will include introducing additional genetic lines of GMO plum into in vitro culture, attempting to induce shoot growth in petiole callus tissue, testing methods (such as cold storage) that extend the time interval between transferring explants into new media, and testing viability of plantlets transferred from in vitro culture back to traditional pot culture.

bioregenerative life support↗
DOE OSTI2026

A flexible and differentiable coil proxy for stellarator equilibrium optimization

Balancing plasma performance and coil cost is a significant challenge when designing a stellarator power plant. Most current stellarator designs are produced through two-stage optimization: stage-1 for the equilibrium and stage-2 for a coil design that reproduces its magnetic configuration. Because few proxies connect both stages, two-stage optimization can produce plasmas that have high-quality physical properties but overly complex coils. In recent years, single-stage optimization has increasingly been used to optimize the plasma and coils simultaneously in order to improve the plasma–coil balance. However, all existing single-stage tools are specialized for filament coils, cannot model coil systems containing permanent magnets or dipole arrays, and continue to be challenged by numerical problems. The quasi-single-stage optimization finds a middle ground by integrating a coil optimization subproblem into stage-1 optimization. We present a flexible, differentiable coil complexity proxy based on the newly developed QUADCOIL coil optimization code. QUADCOIL is fast and can target realistic coil metrics and constraints that are unavailable to codes with comparable speed. We demonstrate the effectiveness and flexibility of the QUADCOIL proxy by presenting two QSS optimization studies. The first study produces a permanent magnet solution for the MUSE stellarator with 34% fewer magnets than previous solutions. The second study produces a coil solution for the ARIES-CS stellarator with 30% reductions in both peak and root-mean-square force.

Fu, Lanke [Princeton Plasma Physics Laboratory (PP↗
NASA NTRS1974

On the passage of radiation through inhomogeneous, moving media. XI - Nonlinear effects on ray paths in the geometrical optics approximation

Study illustrating how the presence of a high-intensity pulse of radiation can distort its own passage through a plane differentially shearing medium. It is demonstrated that the distortion is a sensitive function of the precise, and detailed, variation of the medium's refractive index by considering a couple of simple examples which are worked out numerically. In view of the high-intensity pulses observed from pulsars (approximately 10 to the 30th ergs per pulse), it is believed that the present calculations are of more than academic interest in helping unravel the fundamental properties of pulse production in, and propagating through, differentially sheared media - such as pulsars' magnetospheres within the so-called speed-of-light circle.

Lee, M. A.↗
NASA NTRS1981

On numerical nonlinear analysis of highly flexible spinning cantilevers

The general nonlinear discretized equations of motion of spinning elastic solids and structures are derived as a set of nonlinear ordinary differential equations for the case when the strain-displacement and velocity-displacement relations are nonlinear up to the second order. It is shown that the cost of generation of such equations is proportional to the fourth power of the number of degrees of freedom. A computer program is written to automatically generate the equations for the case of spinning cantilevers with initial imperfections. The types and the number of the coordinate functions used in the trial solution are parameters of the program.

Utku, S.↗
DOE OSTI2022

Running Primal-Dual Gradient Method for Time-Varying Nonconvex Problems

This paper focuses on a time-varying constrained nonconvex optimization problem, and considers the synthesis and analysis of online regularized primal-dual gradient methods to track a Karush-Kuhn-Tucker (KKT) trajectory. The proposed regularized primal-dual gradient method is implemented in a running fashion, in the sense that the underlying optimization problem changes during the execution of the algorithms. In order to study its performance, we first derive its continuous-time limit as a system of differential inclusions. We then study sufficient conditions for tracking a KKT trajectory, and also derive asymptotic bounds for the tracking error (as a function of the time-variability of a KKT trajectory). Further, we provide a set of sufficient conditions for the KKT trajectories not to bifurcate or merge, and also investigate the optimal choice of the parameters of the algorithm. Illustrative numerical results for a time-varying nonconvex problem are provided.

differential inclusion↗
DOE OSTI2024

Weak baselines and reporting biases lead to overoptimism in machine learning for fluid-related partial differential equations

One of the most promising applications of machine learning in computational physics is to accelerate the solution of partial differential equations (PDEs). The key objective of machine-learning-based PDE solvers is to output a sufficiently accurate solution faster than standard numerical methods, which are used as a baseline comparison. Here, we first perform a systematic review of the ML-for-PDE-solving literature. Out of all of the articles that report using ML to solve a fluid-related PDE and claim to outperform a standard numerical method, we determine that 79% (60/76) make a comparison with a weak baseline. Second, we find evidence that reporting biases are widespread, especially outcome reporting and publication biases. We conclude that ML-for-PDE-solving research is overoptimistic: weak baselines lead to overly positive results, while reporting biases lead to under-reporting of negative results. To a large extent, these issues seem to be caused by factors similar to those of past reproducibility crises: researcher degrees of freedom and a bias towards positive results. We call for bottom-up cultural changes to minimize biased reporting as well as top-down structural reforms to reduce perverse incentives for doing so.

97 MATHEMATICS AND COMPUTING↗