Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable simulation”

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 505 records · Page 28

The rollup of a vortex layer near a wall

The behavior of an inviscid vortex layer of non-zero thickness near a wall is studied, both through direct numerical simulation of the two-dimensional vorticity equation at high Reynolds numbers, and using an approximate ordinary nonlinear integro-differential equation which is satisfied in the limit of a thin layer under long-wavelength perturbations. For appropriate initial conditions the layer rolls up and breaks into compact vortices which move along the wall at constant speed. Because of the effect of the wall, they correspond to equilibrium counter-rotating vortex dipoles. This breakup can be related to the disintegration of the initial conditions of the approximate nonlinear dispersive equation into solitary waves. The study is motivated by the formation of longitudinal vortices from vortex sheets in the wall region of a turbulent channel.

Jimenez, Javier↗

Some aspects of high-order numerical solutions of the linear convection equation with forced boundary conditions

A six-stage low-storage Runge-Kutta time-marching method is presented and shown to be an efficient method for use with high-accuracy spatial difference operators for wave propagation problems. The accuracy of the method for inhomogeneous ordinary differential equations is demonstrated through numerical solutions of the linear convection equation with forced boundary conditions. Numerical experiments are presented simulating a sine wave and a Gaussian pulse propagating into and through the domain. For practical levels of mesh refinement corresponding to roughly ten points per wavelength, the six-stage Runge-Kutta method is more accurate than the popular fourth-order Runge-Kutta method. Further numerical experiments are presented which show that the numerical boundary scheme at an inflow boundary can be a significant source of error when high-accuracy spatial discretizations are used.

Zingg, D. W.↗

Effects of non-unity Lewis numbers in diffusion flames

The purpose of this work is to carry out direct numerical simulations of diffusion controlled combustion with non-unity Lewis numbers for the reactants and products, thus accounting for the differential diffusion effects of the temperature and concentration fields. We use a formulation based on combining the conservation equations in a way to eliminate the reaction terms similar to the method used by Burke and Schumann (1928) for unity Lewis numbers. We present calculations for an axisymmetric fuel jet and for a planar, time evolving mixing layer, leaving out the effects of thermal expansion and variations of the transport coefficients due to the heat release. Our results show that the front of the flame shifts toward the fuel or oxygen sides owing to the effect of the differential diffusion and that the location of maximum temperature may not coincide with the flame. The dependence of the distribution of the reaction products on their Lewis number has been investigated.

Linan, A.↗

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This talk describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗

Reliability of Complex Nonlinear Numerical Simulations

This work describes some of the procedure to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations.

Yee, H. C.↗

The Photochemical Reflectance Index from Directional Cornfield Reflectances: Observations and Simulations

The two-layer Markov chain Analytical Canopy Reflectance Model (ACRM) was linked with in situ hyperspectral leaf optical properties to simulate the Photochemical Reflectance Index (PRI) for a corn crop canopy at three different growth stages. This is an extended study after a successful demonstration of PRI simulations for a cornfield previously conducted at an early vegetative growth stage. Consistent with previous in situ studies, sunlit leaves exhibited lower PRI values than shaded leaves. Since sunlit (shaded) foliage dominates the canopy in the reflectance hotspot (coldspot), the canopy PRI derived from field hyperspectral observations displayed sensitivity to both view zenith angle and relative azimuth angle at all growth stages. Consequently, sunlit and shaded canopy sectors were most differentiated when viewed along the azimuth matching the solar principal plane. These directional PRI responses associated with sunlit/shaded foliage were successfully reproduced by the ACRM. As before, the simulated PRI values from the current study were closer to in situ values when both sunlit and shaded leaves were utilized as model input data in a two-layer mode, instead of a one-layer mode with sunlit leaves only. Model performance as judged by correlation between in situ and simulated values was strongest for the mature corn crop (r = 0.87, RMSE = 0.0048), followed by the early vegetative stage (r = 0.78; RMSE = 0.0051) and the early senescent stage (r = 0.65; RMSE = 0.0104). Since the benefit of including shaded leaves in the scheme varied across different growth stages, a further analysis was conducted to investigate how variable fractions of sunlit/shaded leaves affect the canopy PRI values expected for a cornfield, with implications for 20 remote sensing monitoring options. Simulations of the sunlit to shaded canopy ratio near 50/50 +/- 10 (e.g., 60/40) matching field observations at all growth stages were examined. Our results suggest in the importance of the sunlit/shaded fraction and canopy structure in understanding and interpreting PRI.

Cheng, Yen-Ben↗

Utilization of O4 Slant Column Density to Derive Aerosol Layer Height from a Space-Borne UV-Visible Hyperspectral Sensor: Sensitivity and Case Study

The sensitivities of oxygen-dimer (O4) slant column densities (SCDs) to changes in aerosol layer height are investigated using the simulated radiances by a radiative transfer model, the linearized pseudo-spherical vector discrete ordinate radiative transfer (VLIDORT), and the differential optical absorption spectroscopy (DOAS) technique. The sensitivities of the O4 index (O4I), which is defined as dividing O4 SCD by 10(sup 40) molecules (sup 2) per centimeters(sup -5), to aerosol types and optical properties are also evaluated and compared. Among the O4 absorption bands at 340, 360, 380, and 477 nanometers, the O4 absorption band at 477 nanometers is found to be the most suitable to retrieve the aerosol effective height. However, the O4I at 477 nanometers is significantly influenced not only by the aerosol layer effective height but also by aerosol vertical profiles, optical properties including single scattering albedo (SSA), aerosol optical depth (AOD), particle size, and surface albedo. Overall, the error of the retrieved aerosol effective height is estimated to be 1276, 846, and 739 meters for dust, non-absorbing, and absorbing aerosol, respectively, assuming knowledge on the aerosol vertical distribution shape. Using radiance data from the Ozone Monitoring Instrument (OMI), a new algorithm is developed to derive the aerosol effective height over East Asia after the determination of the aerosol type and AOD from the MODerate resolution Imaging Spectroradiometer (MODIS). About 80 percent of retrieved aerosol effective heights are within the error range of 1 kilometer compared to those obtained from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) measurements on thick aerosol layer cases.

Air Pollution↗

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We review the singular value decomposition (SVD) framework and use it for quantifying and discerning vertical information in greenhouse gas retrievals from column integrated absorption measurements. While the commonly used traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We review the SVD method, explicitly recognize the use of an uninformative prior and show it to incur no bias from the choice of the prior. We also make the connection between the SVD method and the pseudo-inverse, which makes it more intuitive and easy to understand. We illustrate the use of the SVD method on an integrated path differential absorption CO2 lidar measurement model and verify our derivations and bias-free properties versus optimal estimation using numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. Hence, the SVD method is particularly useful for situations in which knowledge of the prior mean and prior covariance of the true state (e.g., greenhouse gas profiles) is inadequate.

Ramanathan, Anand K.↗

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗

Low-dimensional description of the dynamics in separated flow past thick airfoils

Results are presented for the numerical simulation of unsteady viscous incompressible flow past thick airfoils. Specifically, flow past a NACA 4424 at an angle of attack of 2.5 deg and Reynolds numbers in the range of 1700-4000 has been simulated using the spectral element method. At these conditions the flow is separatedd and an unsteady wake is formed. Application of the method of empirical eigenfunction reveals the structure of the most energetic components of the flow. These are found to occur in pairs that, through phase exchange, are responsible for the vortex shedding. A set of ordinary differential equations is obtained for the amplitudes of these eigenfunctions by a Galerkin projection of the Navier-Stokes equations. The solutions of the model system are compared with the full simulation. The work is of relevance to the transition process and observed routes to chaos in airfoil wakes.

Deane, Anil E.↗

NASA Experimental Program to Stimulate Competitive Research: South Carolina

The use of an appropriate relationship model is critical for reliable prediction of future urban growth. Identification of proper variables and mathematic functions and determination of the weights or coefficients are the key tasks for building such a model. Although the conventional logistic regression model is appropriate for handing land use problems, it appears insufficient to address the issue of interdependency of the predictor variables. This study used an alternative approach to simulation and modeling urban growth using artificial neural networks. It developed an operational neural network model trained using a robust backpropagation method. The model was applied in the Myrtle Beach region of South Carolina, and tested with both global datasets and areal datasets to examine the strength of both regional models and areal models. The results indicate that the neural network model not only has many theoretic advantages over other conventional mathematic models in representing the complex urban systems, but also is practically superior to the logistic model in its capability to predict urban growth with better - accuracy and less variation. The neural network model is particularly effective in terms of successfully identifying urban patterns in the rural areas where the logistic model often falls short. It was also found from the area-based tests that there are significant intra-regional differentiations in urban growth with different rules and rates. This suggests that the global modeling approach, or one model for the entire region, may not be adequate for simulation of a urban growth at the regional scale. Future research should develop methods for identification and subdivision of these areas and use a set of area-based models to address the issues of multi-centered, intra- regionally differentiated urban growth.

Sutton, Michael A.↗

Design of an Auto-zeroed, Differential, Organic Thin-film Field-effect Transistor Amplifier for Sensor Applications

Organic strain gauge and other sensors require high-gain, precision dc amplification to process their low-level output signals. Ideally, amplifiers would be fabricated using organic thin-film field-effect transistors (OTFT's) adjacent to the sensors. However, OTFT amplifiers exhibit low gain and high input-referred dc offsets that must be effectively managed. This paper presents a four-stage, cascaded differential OTFT amplifier utilizing switched capacitor auto-zeroing. Each stage provides a nominal voltage gain of four through a differential pair driving low-impedance active loads, which provide common-mode output voltage control. p-type pentacence OTFT's are used for the amplifier devices and auto-zero switches. Simulations indicate the amplifier provides a nominal voltage gain of 280 V/V and effectively amplifies a 1-mV dc signal in the presence of 500-mV amplifier input-referred dc offset voltages. Future work could include the addition of digital gain calibration and offset correction of residual offsets associated with charge injection imbalance in the differential circuits.

large-area electronics↗

T-Cell Shenanigans: The Impact of MHC Pathway & Lipid Metabolism Genes on T-Cell Differentiation in the Thymus

Astronaut health and proper immune function are key to the success and sustainability of long term missions in space. While previous studies have observed adaptive immune dysfunction such as diminished WBC counts and cytokine production during spaceflight, the underlying mechanisms behind why immune function worsens are poorly understood. Using transcriptomic data analyzed from mouse thymus tissues available from the GLDS-421 experiment on the Rodent Research-9 mission, we used the GeneLab standardized RNA-Seq pipeline and identified four genes that were largely dysregulated (p<0.05): lipid metabolic genes Hpgd and Pgr and cell cycle regulator genes Cenpe and Kif11. We analyzed how these genes regulate the MHC-TCR interaction, and used this understanding to propose a novel mechanism for alterations to T-cell differentiation. In our experimental methodology, we propose to use 36 transgenic CD4-CD8 mice and 12 wild-type mice from the Jackson Lab to create a microgravity-simulating hindlimb unloading model and test our hypothesis. By silencing the aforementioned genes using custom siRNA based primers, we will use flow cytometry and RNA-Seq to quantify Treg counts in comparison to naive CD4+ cells, Alamar Blue assay for observing the metabolic activity in T-cells, and Enzyme-Linked Immunosorbent Assay (ELISA) to observe the general immune response in terms of T-cell differentiation. We hope this methodology will improve our knowledge of MHC-TCR interaction, contribute to a better understanding of the mechanisms of autoimmune diseases such as lupus, and improve the success rate of newer cancer treatments such as CAR-T cell therapy.

GL4HS↗

A two-dimensional air-to-air combat game - Toward an air-combat advisory system

Air-to-air combat is modeled as a discrete differential game, and by constraining the game to searching for the best guidance laws from the sets of those considered for each opponent, feedback and outcome charts are obtained which can be used to turn one of the automatic opponents into an intelligent opponent against a human pilot. A one-on-one two-dimensional fully automatic, or manned versus automatic, air-to-air combat game has been designed which includes both attack and evasion alternatives for both aircraft. Guidance law selection occurs by flooding the initial-condition space with four simulated fights for each initial condition, depicting the various attack/evasion strategies for the two opponents, and recording the outcomes. For each initial condition, the minimax method from differential games is employed to determine the best choice from the available strategies.

Neuman, Frank↗

A program for identification of linear systems

A program has been written for the identification of parameters in certain linear systems. These systems appear in biomedical problems, particularly in compartmental models of pharmacokinetics. The method presented here assumes that some of the state variables are regularly modified by jump conditions. This simulates administration of drugs following some prescribed drug regime. Parameters are identified by a least-square fit of the linear differential system to a set of experimental observations. The method is especially suited when the interval of observation of the system is very long.

Buell, J.↗

Numerical calculation of the three dimensional transonic flow over a yawed wing.

Results are presented of calculations of the three dimensional steady transonic flow over a finite yawed wing. The full potential flow equation is solved in a transformed coordinate system which permits the boundary conditions to be satisfied exactly. The correct differential properties are enforced by rotating the difference scheme to conform with the flow direction, and fast convergence is assured by simulating a time dependent equation designed to settle quickly to a steady state. Computed lift drag ratios are consistent with the results of wind tunnel tests of a yawed wing conducted by R. T. Jones (1972).-

Jameson, A.↗