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 433 records · Page 24

Numerical simulation of convective response to mesoscale forcing

A two-dimensional mesoscale model was used to initialize a two-dimensional cloud model with both mesoscale thermodynamic and dynamic information utilized in the initialization. The Mesoscale-cloud mdoel linkage was used to examine differential convective response due to mesoscale variations for the April 24, 1982 Atmospheric Variability Experiment-Vertical Atmospheric Sounder (AVE-VAS), Case 4. On this day both data analyses and mesoscale model simulations indicated strong variations in thermodynamic and dynamic structure across the panhandle of Texas where a moderately strong convective line formed. Using the cloud model, the current research has shown a preferred area for strong convection to occur due to concomitant mesoscale convergence and mesoscale destabilization of the atmosphere.

Mcnider, R. T.↗

A new second-order integration algorithm for simulating mechanical dynamic systems

A new integration algorithm which has the simplicity of Euler integration but exhibits second-order accuracy is described. In fixed-step numerical integration of differential equations for mechanical dynamic systems the method represents displacement and acceleration variables at integer step times and velocity variables at half-integer step times. Asymptotic accuracy of the algorithm is twice that of trapezoidal integration and ten times that of second-order Adams-Bashforth integration. The algorithm is also compatible with real-time inputs when used for a real-time simulation. It can be used to produce simulation outputs at double the integration frame rate, i.e., at both half-integer and integer frame times, even though it requires only one evaluation of state-variable derivatives per integration step. The new algorithm is shown to be especially effective in the simulation of lightly-damped structural modes. Both time-domain and frequency-domain accuracy comparisons with traditional integration methods are presented. Stability of the new algorithm is also examined.

Howe, R. M.↗

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗

Simulation of multicorrelated random processes using the FFT algorithm

A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the desired random processes, and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method of simulation over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. The response of a nonlinear string to a distributed noise input is presented as an example.

Wittig, L. E.↗

Development and validation of a hybrid-computer simulator for a transonic cryogenic wind tunnel

A study was undertaken to model the cryogenic wind tunnel process, to validate the model by the use of experimental data from the Langley 0.3 Meter Transonic Cryogenic Tunnel, and to construct an interactive simulator of the cryogenic tunnel using the validated model. Additionally, this model was used for designing closed loop feedback control laws for regulation of temperature and pressure in the 0.3 meter cryogenic tunnel. The global mathematical model of the cryogenic tunnel that were developed consists of coupled, nonlinear differential governing equations based on an energy state concept of the physical cryogenic phenomena. Process equations and comparisons between actual tunnel responses and computer simulation predictions were examined. Also included are the control laws and simulator responses obtained using the feedback schemes for closed loop control of temperature and pressure were also included.

Thibodeaux, J. J.↗

Spaceborne lidar investigations of the atmosphere

Atmospheric investigations with a spaceborne lidar system are discussed. Measurements of aerosols, O3, and H2O with the NASA/LaRC airborne DIAL system are presented as examples of data obtained from space. The NASA/CNES study of an autonomous differential absorption lidar system is described. This system is a precursor to a spaceborne lidar system. Simulations of spaceborne lidar experiments are reviewed, and laser requirements for a spaceborne lidar system are presented.

Browell, E. V.↗

The effects of simulated hypogravity on murine bone marrow cells

Mouse bone marrow cells grown in complete medium at unit gravity were compared with a similar population cultured in conditions that mimic some aspects of microgravity. After the cells adjusted to the conditions that simulated microgravity, they proliferated as fetal or oncogenic populations; their numbers doubled in twelve hour periods. Differentiated subpopulations were depleted from the heterogeneous mixture with time and the undifferentiated hematopoietic stem cells increased in numbers. The cells in the control groups in unit gravity and those in the bioreactors in conditions of microgravity were monitored under a number of parameters. Each were phenotyped as to cell surface antigens using a panel of monoclonal antibodies and flow cytometry. Other parameters compared included: pH, glucose uptake, oxygen consumption and carbon-dioxide production. Nuclear DNA was monitored by flow cytometry. Functional responses were studied by mitogenic stimulation by various lectins. The importance of these findings should have relevance to the space program. Cells should behave predictably in zero gravity; specific populations can be eliminated from diverse populations and other populations isolated. The availability of stem cell populations will enhance both bone marrow and gene transplant programs. Stem cells will permit developmental biologists study the paths of hematopoiesis.

Lawless, Desales↗

Examples of data assimilation in mesoscale models

The keynote address was the problem of physical initialization of mesoscale models. The classic purpose of physical or diabatic initialization is to reduce or eliminate the spin-up error caused by the lack, at the initial time, of the fully developed vertical circulations required to support regions of large rainfall rates. However, even if a model has no spin-up problem, imposition of observed moisture and heating rate information during assimilation can improve quantitative precipitation forecasts, especially early in the forecast. The two key issues in physical initialization are the choice of assimilating technique and sources of hydrologic/hydrometeor data. Another example of data assimilation in mesoscale models was presented in a series of meso-beta scale model experiments with and 11 km version of the MASS model designed to investigate the sensitivity of convective initiation forced by thermally direct circulations resulting from differential surface heating to four dimensional assimilation of surface and radar data. The results of these simulations underscore the need to accurately initialize and simulate grid and sub-grid scale clouds in meso- beta scale models. The status of the application of the CSU-RAMS mesoscale model by the NOAA Forecast Systems Lab for producing real-time forecasts with 10-60 km mesh resolutions over (4000 km)(exp 2) domains for use by the aviation community was reported. Either MAPS or LAPS model data are used to initialize the RAMS model on a 12-h cycle. The use of MAPS (Mesoscale Analysis and Prediction System) model was discussed. Also discussed was the mesobeta-scale data assimilation using a triply-nested nonhydrostatic version of the MM5 model.

Carr, Fred↗

Cryogenic Integration of the 2-14 GHz Eleven Feed in a Wideband Receiver for VLBI2010

The next generation VLBI systems require the design of a wideband receiver covering the 2-14 GHz range, necessitating a wideband feed. Presented here are the 2009 development of a cryogenic 2-14 GHz Eleven feed for reflector radio telescope antennas, including its integration into a cryogenic receiver. The Eleven feed is designed for dual linear polarization and consists of four log-periodic folded dipole arrays. Each pair of arrays is fed by a differential two-wire transmission line connected either to balun or a differential LNA. The present configuration has been measured in many configurations, at various independent labs - corresponding simulations have been done. The results show (across the band) a high polarization efficiency for the feed, with a nearly constant beam width, a reflection coefficient below -10dB, and a constant phase center. Electrical parameters under cryogenic conditions and measured receiver noise temperatures are presented.

Pantaleev, Miroslaw↗

Communication-Aware Orbit Design for Small Spacecraft Swarms around Small Bodies

Exploration of small Solar System bodies has traditionally been performed by single monolithic spacecraft carrying a number of science instruments. However, science instruments typically cannot be operated simultaneously due to the instrument requirements including optimal viewing angle, surface illumination, altitude and ground resolution, power, and data constraints. This observation has motivated interest in multi-spacecraft architectures where a swarm of small spacecraft, each carrying a single science instrument, studies a small body after being deployed by a carrier spacecraft, which then collects data from the vehicles and relays it to Earth. Such architectures hold promise to yield significant improvements in mission efficiency, increases in data quality, and shorter mission duration. A key difficulty in the design of such missions is the selection of orbits for the small spacecraft, which must satisfy not only instrument requirements, but also strict inter-spacecraft communication and on-board storage constraints. To address this, in this paper, we present a novel computationally-efficient optimization algorithm for \emph{communication-aware design} of the orbits of a small spacecraft swarm orbiting a small body. The proposed approach captures constraints including instrument requirements, inter-spacecraft communication bandwidths, and on-board memory usage, and it can accommodate highly irregular gravity field models and surface geometries. We propose an efficient algorithm for optimization of instrument observations and inter-spacecraft communications; we then leverage the differentiable nature of the proposed algorithm to accelerate a gradient-based global search algorithm. Numerical simulations of a six-spacecraft swarm studying 433 Eros show that the proposed approach successfully identifies high-quality orbits, and significantly outperform communication-agnostic optimization techniques, resulting in a 10% increase in scientific returns and a 30% increase in the quality of the collected data.

Rahmani, Amir↗

Estimating the Concentration of Large Raindrops from Polarimetric Radar and Disdrometer Observations

Estimation of rainfall integral parameters, including radar observables, and empirical relations between them are sensitive to the truncation of the drop size distribution (DSD), particularly at the large drop end. The sensitivity of rainfall integral parameters to the maximum drop diameter (D(sub max)) is exacerbated at C‐band since resonance effects are pronounced for large drops in excess of 5 mm diameter (D). Due to sampling limitations, it is often difficult to reliably estimate D(sub max) with disdrometers. The resulting uncertainties in D(sub max0 potentially increase errors in radar retrieval methods, particularly at C‐band, that rely on disdrometer observations for DSD input to radar models. In fact, D(sub max) is typically an assumed DSD parameter in the development of radar retrieval methods. Because of these very uncertainties, it is difficult to independently confirm disdrometer estimates of D(sub max) with polarimetric radar observations. A couple of approaches can be taken to reduce uncertainty in large drop measurement. Longer integration times can be used for the collection of larger disdrometer samples. However, integration periods must be consistent with a radar resolution volume (RRV) and the temporal and spatial scales of the physical processes affecting the DSD therein. Multiple co‐located disdrometers can be combined into a network to increase the sample size within a RRV. However, over a reasonable integration period, a single disdrometer sample volume is many orders of magnitudes less than a RRV so it is not practical to devise a network of disdrometers that has an equivalent volume to a typical RRV. Since knowledge of DSD heterogeneity and large drop occurrence in time and space is lacking, the specific accuracy or even general representativeness of disdrometer based D(sub max) and large drop concentration estimates within a RRV are currently unknown. To address this complex issue, we begin with a simpler question. Is the frequency of occurrence of large rain drops (D > 5 mm) in disdrometer observations, either stand alone or networked, generally representative and consistent with polarimetric radar observations? We first show from simulations that the concentration of large (D > 5 mm) rain drops (N(sub T5)) can be estimated from polarimetric observations of specific differential phase (K(sub dp)) and differential reflectivity (Z(sub dr)), N(sub T5)=F(K(sub dp),Z(sub dr)), or horizontal reflectivity (Z(sub h)) and Z(sub dr), N(sub T5)=(Z(sub h),Z(sub dr)). We assess the error associated with polarimetric retrieval of N(sub T5), including sensitivity to D(sub max) parameterization assumptions and measurement error in the radar simulations. Polarimetric measurements at S‐band and C‐band will then be used to retrieve estimates of N(sub T5) and compared to disdrometer estimates of N(sub T5). After careful consideration of retrieval error, we will check consistency between disdrometer and polarimetric radar estimates of N(sub T5) and the frequency of occurrence of large rain drops in a variety of precipitating regimes using data from NASA's Global Precipitation Measurement (GPM) Ground Validation (GV) program, including field campaigns such as MC3E (Oklahoma) and IFloodS (Iowa) and extended measurements over Huntsville, Alabama and NASA Wallops Flight Facility in coastal Virginia.

Carey, Lawrence D.↗

Improved technique for digital simulation of bending and slosh phenomena

Mathematical model representation of bending and slosh phenomena in the Saturn vehicle results in linear second order differential equations. Improved technique was developed to provide a real-time digital solution of the equations. The technique may also be applied to nonreal time digital simultations, resulting in savings of digital computer time.

Stauffer, N. E.↗

Attitude dynamics simulation subroutines for systems of hinge-connected rigid bodies with nonrigid appendages

Three computer subroutines designed to solve the vector-dyadic differential equations of rotational motion for systems that may be idealized as a collection of hinge-connected rigid bodies assembled in a tree topology, with an optional flexible appendage attached to each body are reported. Deformations of the appendages are mathematically represented by modal coordinates and are assumed small. Within these constraints, the subroutines provide equation solutions for (1) the most general case of unrestricted hinge rotations, with appendage base bodies nominally rotating at a constant speed, (2) the case of unrestricted hinge rotations between rigid bodies, with the restriction that those rigid bodies carrying appendages are nominally nonspinning, and (3) the case of small hinge rotations and nominally nonrotating appendages. Sample problems and their solutions are presented to illustrate the utility of the computer programs.

Fleischer, G. E.↗

Simulation of Dust-Laden Flow in the DLR L2K Facility

An integrated CFD-particle trajectory code called US3D-DUST is used to simulate the dusty flow environments performed in the German Aerospace Center (DLR) L2K arcjet facility. The particle trajectory code uses a Lagrangian approach, where the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The L2K experiments measured mass loss to flat-face shaped samples of the Norcoat-Liège thermal protection system material. The US3D-DUST simulations were able to reproduce particle velocity and stagnation point heating rate data. When applied to the L2K dusty flow conditions, the US3D-DUST predictions of mass loss due to particle impacts and surface recession rates over-predicted the experimental values. It is recommended to perform additional dusty flow tests in the L2Kto provide and improve the surface damage model for Norcoat-Liège.

Thermal protection systems↗

Computer analysis of the negative differential resistance switching phenomenon of double-injection devices

By directly solving the semiconductor differential equations for the double-injection (DI) devices involving two interacting deep levels, the authors studied the negative differential resistance switching characteristic and its relationship with the device dimension, doping level, and dependence on the deep impurity profile. Computer simulation showed that although one can increase the threshold voltage by increasing the device length, the excessive holding voltage that would follow would put this device in a very limited application such as pulse power source. The excessive leakage current in the low conductance state also jeopardizes the attempt to use the device for any practical purpose. Unless there are new materials and deep impurities found that have a great differential hole and electron capture cross sections and a reasonable energy bandgap for low intrinsic carrier concentration, no big improvement in the fate of DI devices is expected in the near future.

Shieh, Tsay-Jiu↗

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

Electrostatic Beneficiation of Lunar Simulant

Electrostatic beneficiation of lunar regolith is a method allowing refinement of specific minerals in the material for processing on the moon. The use of tribocharging the regolith prior to separation was investigated on the lunar simulant MLS-I by passing the dust through static mixers constructed from different materials; aluminum, copper, stainless steel, and polytetrafluoroethylene (PTFE). The amount of charge acquired by the simulant was dependent upon the difference in the work function of the dust and the charging material. XPS and SEM were used to characterize the simulant after it was sieved into five size fractions (> 100 pm, 75-100 pm, 50- 75 pm, 50-25 pm, and < 25 pm), where very little difference in surface composition was observed between the sizes. Samples of the smallest (< 25 pm) and largest (> 100 pm) size fractions were beneficiated through a charge separator using the aluminum (charged the simulant negatively) and PTFE (charged positively) mixers. The mass fractions of the separated simulant revealed that for the larger particle size, significant unipolar charging was observed for both mixers, whereas for the smaller particle sizes, more bipolar charging was observed, probably due to the finer simulant adhering to the inside of the mixers shielding the dust from the charging material. Subsequent XPS analysis of the beneficiated fractions showed the larger particle size fraction having some species differentiation, but very little difference for the smaller.size. Although MLS-1 was made to have similar chemistry to actual lunar dust, its mineralogy is quite different. On-going experiments are using NASA JSC-1 lunar simulant. A vacuum chamber has been constructed, and future experiments are planned in a simulated lunar environment.

Trigwell, Steve↗

Nanoflare vs Footpoint Heating : Observational Signatures

Time lag analysis shows very long time lags between all channel pairs. Impulsive heating cannot address these long time lags. 3D Simulations of footpoint heating shows a similar pattern of time lags (magnitude and distribution) to observations. Time lags and relative peak intensities may be able to differentiate between TNE and impulsive heating solutions. Adding a high temperature channel (like XRT Be-­thin) may improve diagnostics.

corona simulation heating↗