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 613 records · Page 34

The diurnal variability of marine stratocumulus clouds

Observations of the diurnal variation of marine stratocumulus clouds obtained during July 1987 as part of the First ISCCP Regional Experiment are described. The observations show a clear diurnal signature in the structure of the cloud-topped marine boundary layer, with a large variation in the height of the cloud base, with lower heights observed at night. The observational results are numerically simulated using a 1D second-moment turbulence-closure model of Koracin and Rogers (1990) for a hydrostatically incompressible fluid. The model results indicate that the differential heating of the cloud layer promotes mixing in the upper part of the cloud and stabilizes the lower part of the layer.

Rogers, David P.↗

Investigation of a Nonlinear Control System

A discontinuous variation of coefficients of the differential equation describing the linear control system before nonlinear elements are added is studied in detail. The nonlinear feedback is applied to a second-order system. Simulation techniques are used to study performance of the nonlinear control system and to compare it with the linear system for a wide variety of inputs. A detailed quantitative study of the influence of relay delays and of a transport delay is presented.

Flugge-Lotz, I↗

Testing the Equivalence Principle in an Einstein Elevator: Detector Dynamics and Gravity Perturbations

We discuss specific, recent advances in the analysis of an experiment to test the Equivalence Principle (EP) in free fall. A differential accelerometer detector with two proof masses of different materials free falls inside an evacuated capsule previously released from a stratospheric balloon. The detector spins slowly about its horizontal axis during the fall. An EP violation signal (if present) will manifest itself at the rotational frequency of the detector. The detector operates in a quiet environment as it slowly moves with respect to the co-moving capsule. There are, however, gravitational and dynamical noise contributions that need to be evaluated in order to define key requirements for this experiment. Specifically, higher-order mass moments of the capsule contribute errors to the differential acceleration output with components at the spin frequency which need to be minimized. The dynamics of the free falling detector (in its present design) has been simulated in order to estimate the tolerable errors at release which, in turn, define the release mechanism requirements. Moreover, the study of the higher-order mass moments for a worst-case position of the detector package relative to the cryostat has led to the definition of requirements on the shape and size of the proof masses.

Hubbard, Dorthy↗

Numerical Study of the Generation of Linear Energy Transfer Spectra for Space Radiation Applications

In analyzing charged particle spectra in space due to galactic cosmic rays (GCR) and solar particle events (SPE), the conversion of particle energy spectra into linear energy transfer (LET) distributions is a convenient guide in assessing biologically significant components of these spectra. The mapping of LET to energy is triple valued and can be defined only on open energy subintervals where the derivative of LET with respect to energy is not zero. Presented here is a well-defined numerical procedure which allows for the generation of LET spectra on the open energy subintervals that are integrable in spite of their singular nature. The efficiency and accuracy of the numerical procedures is demonstrated by providing examples of computed differential and integral LET spectra and their equilibrium components for historically large SPEs and 1977 solar minimum GCR environments. Due to the biological significance of tissue, all simulations are done with tissue as the target material.

Badavi, Francis F.↗

Design of synthetic collagens that assemble into supramolecular banded fibers as a functional biomaterial testbed

Collagens are the most abundant proteins of the extracellular matrix, and the hierarchical folding and supramolecular assembly of collagens into banded fibers is essential for mediating cell-matrix interactions and tissue mechanics. Collagen extracted from animal tissues is a valuable commodity, but suffers from safety and purity issues, limiting its biomaterials applications. Synthetic collagen biomaterials could address these issues, but their construction requires molecular-level control of folding and supramolecular assembly into ordered banded fibers, comparable to those of natural collagens. Here, we show an innovative class of banded fiber-forming synthetic collagens that recapitulate the morphology and some biological properties of natural collagens. The synthetic collagens comprise a functional-driver module that is flanked by adhesive modules that effectively promote their supramolecular assembly. Multiscale simulations support a plausible molecular-level mechanism of supramolecular assembly, allowing precise design of banded fiber morphology. We also experimentally demonstrate that synthetic fibers stimulate osteoblast differentiation at levels comparable to natural collagen. This work thus deepens understanding of collagen biology and disease by providing a ready source of safe, functional biomaterials that bridge the current gap between the simplicity of peptide biophysical models and the complexity of in vivo animal systems.

Jinyuan Hu↗

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗

Particle/fluid simulations of an eruptive flare: Identifying the field-aligned currents responsible for the hard x-rays

While magnetohydrodynamics (MHD) can provide a reasonable description of the overall magnetic reconnection that is believed to drive flares, additional, and often separate processes have to be envoked to in order to explain the electron acceleration that is responsible for many of the observed flare emissions. A new model that incorporates the dynamic coronal current sheets, the reconnection site, and possible electron acceleration processes is developed through the use of two-dimensional particle and modified two-fluid simulations. The specific example of an eruptive flare driven by the coalescence of flux tubes supported by prescribed photospheric current elements is evaluated. It is shown that the electrons and ions have differential trajectories through the coronal current sheet which leads to the development of additonal plasma currents that flow around the surface of the current sheet. These surface currents are explicitly neglected in MHD but they are vital to the flare dynamics because they divert current from the coronal current sheet into the chromosphere, producing an effective resistivity that aids the development of fast reconnection. Because the surface currents are in the plane of the magnetic field, electrons in them experience strong acceleration and can account for the observed hard X-ray emissions. Model predictions are compared with observed time profiles of hard X-ray emissions and Doppler shifts seen in soft X-ray line emissions and are able to account for such features as (1) the asymmetry in the rise and decay time of the hard X-rays, (2) the apparent delay between the largest Doppler shifts and the hard X-ray peak, and (3) the relatively low intensity of the blue-shifted component. The use of particle and fluid simulations is important because it provides different, but complementary treatments of the electron acceleration, the global magnetic morphology, and the flare current system.

Winglee, R. M.↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

Retrieving DSD Moments from GPM-DPR: A Simulation Study Based on the Full DSD Spectra Characterized by the Generalized Gamma Model

Recently, a method to retrieve rain drop size distribution (DSD) moments from X-band dual polarization radar measurements has been developed using copolar reflectivity, differential reflectivity and specific attenuation. Two reference moments are retrieved first, followed by reconstructing the DSDs and calculating other moments using a generalized gamma model to represent the underlying shape, h(x), corresponding to the pair of chosen reference moments. Here we explore a similar approach but the retrieval method in this study uses dual-frequency radar measurements. In the case of GPM-DPR, the two frequencies are 13.8and 35 GHz, and their products include attenuation corrected reflecivities at both frequencies (Z and Z) and the specific attenuation (k, k, if available). Our approach is to use these products to determine two reference moments, namely M3 andM6 representing the third and the sixth moments respectively. Then, as with the polarimetric radar retrievals, we use the most probable h(x) to reconstruct the full DSD spectra, from which other moments are calculated. The best two DPR products for estimating M3 and M6 appear to be A and Z. Figure 1 shows the retrievals versus the ‘true’ moments. 2930 three-minute DSDs were used for the (T-matrix) scattering calculations at Ku and Ka bands but only cases with A > 0.5 dB/km were chosen. For the retrieved moments, the integration was performed only up to 6 mm drop diameter. The [1:1] line is included in Fig. 1. Even with just two DPR products as input the results seem promising, although the lower order moments show somewhat more scatter, especially the zeroth moment, M0. We will quantify the retrieval errors, and additionally examine the stability of h(x). For the latter, data from (i) Greeley, Colorado, (ii) Huntsville, Alabama, and (iii)Wallops, Virginia, will be used. Finally, a GPM overpass case over Huntsville, during a widespread rain event on 11 April 2016, will be considered as an initial test case.

precipitation↗

Quantitative EEG patterns of differential in-flight workload

Four test pilots were instrumented for in-flight EEG recordings using a custom portable recording system. Each flew six, two minute tracking tasks in the Calspan NT-33 experimental trainer at Edwards AFB. With the canopy blacked out, pilots used a HUD display to chase a simulated aircraft through a random flight course. Three configurations of flight controls altered the flight characteristics to achieve low, moderate, and high workload, as determined by normative Cooper-Harper ratings. The test protocol was administered by a command pilot in the back seat. Corresponding EEG and tracking data were compared off-line. Tracking performance was measured as deviation from the target aircraft and combined with control difficulty to achieve an estimate of 'cognitive workload'. Trended patterns of parietal EEG activity at 8-12 Hz were sorted according to this classification. In all cases, high workload produced a significantly greater suppression of 8-12 Hz activity than low workload. Further, a clear differentiation of EEG trend patterns was obtained in 80 percent of the cases. High workload produced a sustained suppression of 8-12 Hz activity, while moderate workload resulted in an initial suppression followed by a gradual increment. Low workload was associated with a modulated pattern lacking any periods of marked or sustained suppression. These findings suggest that quantitative analysis of appropriate EEG measures may provide an objective and reliable in-flight index of cognitive effort that could facilitate workload assessment.

Sterman, M. B.↗

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗

The test facility requirements for the thermal vacuum thermal balance test of the Cosmic Background Explorer Observatory

The Cosmic Background Explorer Observatory (COBE) underwent a thermal vacuum thermal balance test in the Space Environment Simulator (SES). This was the largest and most complex test ever conducted at this facility. The 4 x 4 m (13 x 13 ft) spacecraft weighed approx. 2223 kg (4900 lbs) for the test. The test set up included simulator panels for the inboard solar array panels, simulator panels for the flight cowlings, Sun and Earth Sensor stimuli, Thermal Radio Frequency Shield heater stimuli and a cryopanel for thermal control in the Attitude Control System Shunt Dissipator area. The fixturing also included a unique 4.3 m (14 ft) diameter Gaseous Helium Cryopanel which provided a 20 K environment for the calibration of one of the spacecraft's instruments, the Differential Microwave Radiometer. This cryogenic panel caused extra contamination concerns and a special method was developed and written into the test procedure to prevent the high buildup of condensibles on the panel which could have led to backstreaming of the thermal vacuum chamber. The test was completed with a high quality simulated space environment provided to the spacecraft. The test requirements, test set up, and special fixturing are described.

Milam, Laura J.↗

Test facility requirements for the thermal vacuum thermal balance test of the Cosmic Background Explorer Observatory

The Cosmic Background Explorer Observatory (COBE) underwant a thermal vacuum thermal balance test in the Space Environment Simulator (SES). This was the largest and most complex test ever conducted at this facility. The 4 x 4 m (13 x 13 ft) spacecraft weighed approx. 2223 kg (4900 lbs) for the test. The test set up included simulator panels for the inboard solar array panels, simulator panels for the flight cowlings, Sun and Earth Sensor stimuli, Thermal Radio Frequency Shield heater stimuli and a cryopanel for thermal control in the Attitude Control System Shunt Dissipator area. The fixturing also included a unique 4.3 m (14 ft) diameter Gaseous Helium Cryopanel which provided a 20 K environment for the calibration of one of the spacecraft's instruments, the Differential Microwave Radiometer. This cryogenic panel caused extra contamination concerns and a special method was developed and written into the test procedure to prevent the high buildup of condensibles on the panel which could have led to backstreaming of the thermal vacuum chamber. The test was completed with a high quality simulated space environment provided to the spacecraft. The test requirements, test set up, and special fixturing are described.

Milam, Laura J.↗

Numerical determination of the transmissibility characteristics of a squeeze film damped forced vibration system

Numerical solutions of the governing equations of motion of a liquid squeeze film damped forced vibration system were carried out to examine the feasibility of using a liquid squeeze film to cushion and protect large structures, such as buildings, located in areas of high seismic activity. The mathematical model used was that for a single degree of freedom squeeze film damped spring mass system. The input disturbance was simulated by curve fitting actual seismic data with an eleventh order Lagranging polynomial technique. Only the normal component of the seismic input was considered. The nonlinear, nonhomogeneous governing differential equation of motion was solved numerically to determine the transmissibility over a wide range of physical parameters using a fourth-order Runge-Kutta technique. It is determined that a liquid squeeze film used as a damping agent in a spring-mass system can significantly reduce the response amplitude for a seismic input disturbance.

Sutton, M. A.↗

Numerical study of baroclinic instability

The effect of a power law gravity field on baroclinic instability is examined with emphasis on the case of inverse fifth power gravity, since this is the power law produced when terrestrial gravity is simulated in spherical geometry by electrostatic means. Growth rates of unstable normal modes were obtained as a function of parameters of the problem by solving a second order differential equation numerically. Results are compared with those from an earlier study where gravity was a constant. The conclusion is that, over the range of parameter space explored here, there is no significant change in the character of theoretical regime diagrams if the vertically averaged gravity is used as a parameter.

Geisler, J. E.↗

Spectral simulation of unsteady compressible flow past a circular cylinder

An unsteady compressible viscous wake flow past a circular cylinder was successfully simulated using spectral methods. A new approach in using the Chebyshev collocation method for periodic problems is introduced. It was further proved that the eigenvalues associated with the differentiation matrix are purely imaginary, reflecting the periodicity of the problem. It was been shown that the solution of a model problem has exponential growth in time if improper boundary conditions are used. A characteristic boundary condition, which is based on the characteristics of the Euler equations of gas dynamics, was derived for the spectral code. The primary vortex shedding frequency computed agrees well with the results in the literature for Mach = 0.4, Re = 80. No secondary frequency is observed in the power spectrum analysis of the pressure data.

Don, Wai-Sun↗

A variational approach to understanding white dwarf evolution

A variational approach is used to map out the effects that uncertainties in the theoretical model parameters have upon the derived ages near the observed cutoff in the white dwarf luminosity function. Two representative sequences are assessed, simulating a white dwarf with a 50/50 C/O mixture in the core and an outer helium layer and a white dwarf with a C/O convective overshooting profile. The differential effects that the variation of selected model input quantities has on the ages are reported, and it is concluded that internal theoretical uncertainties are small and getting smaller with time, and the results underscore the power of using the observed white dwarf luminosity function for studying the history of star formation in the Galaxy.

Wood, M. A.↗

Bioelectric Control of a 757 Class High Fidelity Aircraft Simulation

This paper presents results of a recent experiment in fine grain Electromyographic (EMG) signal recognition, We demonstrate bioelectric flight control of 757 class simulation aircraft landing at San Francisco International Airport. The physical instrumentality of a pilot control stick is not used. A pilot closes a fist in empty air and performs control movements which are captured by a dry electrode array on the arm, analyzed and routed through a flight director permitting full pilot outer loop control of the simulation. A Vision Dome immersive display is used to create a VR world for the aircraft body mechanics and flight changes to pilot movements. Inner loop surfaces and differential aircraft thrust is controlled using a hybrid neural network architecture that combines a damage adaptive controller (Jorgensen 1998, Totah 1998) with a propulsion only based control system (Bull & Kaneshige 1997). Thus the 757 aircraft is not only being flown bioelectrically at the pilot level but also demonstrates damage adaptive neural network control permitting adaptation to severe changes in the physical flight characteristics of the aircraft at the inner loop level. To compensate for accident scenarios, the aircraft uses remaining control surface authority and differential thrust from the engines. To the best of our knowledge this is the first time real time bioelectric fine-grained control, differential thrust based control, and neural network damage adaptive control have been integrated into a single flight demonstration. The paper describes the EMG pattern recognition system and the bioelectric pattern recognition methodology.

Jorgensen, Charles↗