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 163 records · Page 9

Numerical calculations of high-altitude differential charging: Preliminary results

A two dimensional simulation program was constructed in order to obtain theoretical predictions of floating potential distributions on geostationary spacecraft. The geometry was infinite-cylindrical with angle dependence. Effects of finite spacecraft length on sheath potential profiles can be included in an approximate way. The program can treat either steady-state conditions or slowly time-varying situations, involving external time scales much larger than particle transit times. Approximate, locally dependent expressions were used to provide space charge, density profiles, but numerical orbit-following is used to calculate surface currents. Ambient velocity distributions were assumed to be isotropic, beam-like, or some superposition of these.

Laframboise, J. G.↗

Optimum spectral bands for rock discrimination

Using stepwise discriminant analysis on spectral reflectance and spectral emissivity data collected by a Multispectral Scanner and Data System, mounted in an NC-130B aircraft and flown at an altitude of approximately 3 km, spectral bands were ranked as to their usefulness in separating specific rock types and rock alteration products in seven geologically diverse Utah sites. The optimum band for rock discrimination included the 1.18 to 1.30 micron interval, and the optimum combination of bands comprised the 1.18 to 1.30, 4.50 to 4.75, 0.46 to 0.50, 1.52 to 1.73, and 2.10 to 2.36 micron intervals. It is concluded that the spectral interval combination was more successful in differentiating geologic materials than either simulated Multispectral Scanner bands or simulated Thematic Mapper bands.

Siegrist, A. W.↗

Numerical observations of dynamic behavior in two-dimensional compressible convection

Using numerical simulations of the anelastic partial differential equations, convective flow in a two-dimensional polytropic atmosphere is investigated. Fixing the aspect ratio at unity and choosing the viscosity as the order parameter, the numerical solutions have indicated a variety of dynamic behavior including (1) oscillatory motion with a period of order the eddy turnover time, (2) quasi-periodic motion with a long period much longer than an eddy turnover time, and (3) at least two branches of solutions in certain regimes.

Ginet, G. P.↗

Effects of the space environment on space-based radar phased-array antenna; status and preliminary observations (LDEF Experiment A0133)

The overall objective of the experiment was to evaluate the effect of the space environment on components considered for a Space-Based Radar (SBR) Phased-Array Antenna. Of primary interest was a study of the degradation of the polyimide film Kapton (DuPont trademark), the material considered for use in the antenna plane. The most striking result of the experiment was the overall good condition of the Kapton antenna planes and Kapton tensile specimens, despite nearly six years of exposure to the space environment. This was largely attributable to the orientation of the Kapton (parallel and flush on the space end) and the stability of the Long Duration Exposure Facility (LDEF) in orbit. However, weathering of exposed Kapton surfaces was not insignificant. Results on elongation and mechanical properties of the plain and the fiberglass-reinforced Kapton are presented. Reduction in strain to failure of flight-exposed Kapton is attributed to surface defects of these specimens. Physical property testing of the materials to date reveals no significant difference between flight-exposed and control material. The second objective was to investigate the interaction between high-voltage electrodes and typical spacecraft contaminants in simulation of discharge triggering across differentially charged dielectric surfaces (spacecraft charging conditions). Electronic data acquisition and memory systems appeared to operate correctly, but very few discharges were recorded. Induced radioactivity, contamination, impacts, and orientation features of atomic oxygen erosion were observed.

J B Whiteside↗

Low-Power SOI CMOS Transceiver

The work aims at developing a low-power Silicon on Insulator Complementary Metal Oxide Semiconductor (SOI CMOS) Transceiver for deep-space communications. RF Receiver must accomplish the following tasks: (a) Select the desired radio channel and reject other radio signals, (b) Amplify the desired radio signal and translate them back to baseband, and (c) Detect and decode the information with Low BER. In order to minimize cost and achieve high level of integration, receiver architecture should use least number of external filters and passive components. It should also consume least amount of power to minimize battery cost, size, and weight. One of the most stringent requirements for deep-space communication is the low-power operation. Our study identified that two candidate architectures listed in the following meet these requirements: (1) Low-IF receiver, (2) Sub-sampling receiver. The low-IF receiver uses minimum number of external components. Compared to Zero-IF (Direct conversion) architecture, it has less severe offset and flicker noise problems. The Sub-sampling receiver amplifies the RF signal and samples it using track-and-hold Subsampling mixer. These architectures provide low-power solution for the short- range communications missions on Mars. Accomplishments to date include: (1) System-level design and simulation of a Double-Differential PSK receiver, (2) Implementation of Honeywell SOI CMOS process design kit (PDK) in Cadence design tools, (3) Design of test circuits to investigate relationships between layout techniques, geometry, and low-frequency noise in SOI CMOS, (4) Model development and verification of on-chip spiral inductors in SOI CMOS process, (5) Design/implementation of low-power low-noise amplifier (LNA) and mixer for low-IF receiver, and (6) Design/implementation of high-gain LNA for sub-sampling receiver. Our initial results show that substantial improvement in power consumption is achieved using SOI CMOS as compared to standard CMOS process. Potential advantages of SOI CMOS for deep-space communication electronics include: (1) Radiation hardness, (2) Low-power operation, and (3) System-on-Chip (SOC) solutions.

Fujikawa, Gene↗

Aeroelastic Analysis of Helicopter Rotor Blades Incorporating Anisotropic Piezoelectric Twist Actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consists of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for dynamics simulation using numerical integration. The twist actuation responses for three conceptual fullscale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗

Template Matching Approach to Signal Prediction

A new approach to signal prediction and prognostic assessment of spacecraft health resolves an inherent difficulty in fusing sensor data with simulated data. This technique builds upon previous work that demonstrated the importance of physics-based transient models to accurate prediction of signal dynamics and system performance. While models can greatly improve predictive accuracy, they are difficult to apply in general because of variations in model type, accuracy, or intended purpose. However, virtually any flight project will have at least some modeling capability at its disposal, whether a full-blown simulation, partial physics models, dynamic look-up tables, a brassboard analogue system, or simple hand-driven calculation by a team of experts. Many models can be used to develop a predict, or an estimate of the next day s or next cycle s behavior, which is typically used for planning purposes. The fidelity of a predict varies from one project to another, depending on the complexity of the simulation (i.e. linearized or full differential equations) and the level of detail in anticipated system operation, but typically any predict cannot be adapted to changing conditions or adjusted spacecraft command execution. Applying a predict blindly, without adapting the predict to current conditions, produces mixed results at best, primarily due to mismatches between assumed execution of spacecraft activities and actual times of execution. This results in the predict becoming useless during periods of complicated behavior, exactly when the predict would be most valuable. Each spacecraft operation tends to show up as a transient in the data, and if the transients are misaligned, using the predict can actually harm forecasting performance. To address this problem, the approach here expresses the predict in terms of a baseline function superposed with one or more transient functions. These transients serve as signal templates, which can be relocated in time and space against the signal background. One then has the ability to reconstruct a signal regardless of the precise timing of the transients. During operation, one applies the actual start times of spacecraft activities as they occur, and produces a reconstructed, accurate predict in real-time. This general approach is valid under two important conditions. First, the transients themselves must be time-invariant. Second, the transients must be reasonably consistent with respect to different operating points. Both of these assumptions are generally valid, but in the case of a complicated system with numerous types of overlapping transients, this approach may not be effective. Fortunately, there tend to be few transients in spacecraft telemetry of sensor quantities or low-level health and status information because these signals rarely reflect multiple different types of operation. Furthermore, if the predict is at least reasonably close to actual operation, the shift either in time, or in the operating point when the transient occurs - is likely to be small. The proposed approach considers three ways to recognize a transient. The first and most reliable is to use a different signal that identifies operating mode - often the transient will be correlated to a change in operating mode, and this change can usually be detected positively from discrete signals in spacecraft telemetry. If there is no useful mode signal, the second option is to identify the transient template by hand, and detect the onset of the transient using a curve-fitting approach. Finally, if one elects not to choose by hand, there is an option for automatic selection. This technique has been applied to sensor data from several JPL missions and industrial applications, demonstrating an improvement in accurate prediction of future behavior and early detection of system problems. This technique is applicable to practically any time-varying, quantitative sensor measurement.

Mackey, Ryan↗

An Observational Approach to Low Lunar Frozen Orbit Design

Low lunar frozen orbits can be used to enable long-duration missions at the Moon while also minimizing the fuel usage necessary for orbit maintenance. There is no analytical expression for determining these frozen orbits, so a high-fidelity simulation approach is required. A differential correction targeting approach was used to find initial conditions that zeroed out argument of periapsis and eccentricity growth. This approach was used during the Lunar Reconnaissance Orbiter extended mission period for four different frozen orbit instances. Observations of the evolution of key orbit elements (e.g. periapsis altitude and the eccentricity vector) led to an understanding of the orbit dynamics and helped to formulate a simple observational approach to finding frozen or quasi-frozen low lunar orbits that looks at minimizing the spread in the periselene altitude evolution.

Moon Frozen Orbit LRO↗

An Observational Approach to Low Lunar Frozen Orbit Design

Low lunar frozen orbits can be used to enable long-duration missions at the Moon while also minimizing the fuel usage necessary for orbit maintenance. There is no analytical expression for determining these frozen orbits, so a high-fidelity simulation approach is required. A differential correction targeting approach was used to find initial conditions that zeroed out argument of periapsis and eccentricity growth. This approach was used during the Lunar Reconnaissance Orbiter extended mission period for four different frozen orbit instances. Observations of the evolution of key orbit elements (e.g. periapsis altitude and the eccentricity vector) led to an understanding of the orbit dynamics and helped to formulate a simple observational approach to finding frozen or quasi-frozen low lunar orbits that looks at minimizing the spread in the periselene altitude evolution.

Moon Frozen Orbit LRO↗

Geothermal reservoir simulation

The prediction of long-term geothermal reservoir performance and the environmental impact of exploiting this resource are two important problems associated with the utilization of geothermal energy for power production. Our research effort addresses these problems through numerical simulation. Computer codes based on the solution of partial-differential equations using finite-element techniques are being prepared to simulate multiphase energy transport, energy transport in fractured porous reservoirs, well bore phenomena, and subsidence.

Mercer, J. W., Jr.↗

Trellis coding with Continuous Phase Modulation (CPM) for satellite-based land-mobile communications

This volume of the final report summarizes the results of our studies on the satellite-based mobile communications project. It includes: a detailed analysis, design, and simulations of trellis coded, full/partial response CPM signals with/without interleaving over various Rician fading channels; analysis and simulation of computational cutoff rates for coherent, noncoherent, and differential detection of CPM signals; optimization of the complete transmission system; analysis and simulation of power spectrum of the CPM signals; design and development of a class of Doppler frequency shift estimators; design and development of a symbol timing recovery circuit; and breadboard implementation of the transmission system. Studies prove the suitability of the CPM system for mobile communications.

Source record↗

Influence of Rayleigh-Doppler broadening on the selection of H2O dial system parameters

Computer simulations have enabled the performance of a H2O Differential Absorption Lidar (DIAL) system to be studied by spectrally analyzing the forward propagating and backscattered laser energy. The simulations were done for a high altitude (21 km) DIAL system operating in a nadir-viewing mode. The influence of Rayleigh Doppler broadening on DIAL measurement accuracies were evaluated and show that the Rayleigh broadening influence, which can be corrected to first order in regions free of large aerosol gradients, reduces the sensitivity of DIAL H2O measurement errors in the upper tropospheric region. The ability to correct the Rayleigh broadening and the selection of H2O DIAL parameters when all the systematic effects are combined, were discussed.

Ismail, S.↗

Wave equations and computational models for sonic boom propagation through a turbulent atmosphere

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by atmospheric turbulence. A nonlinear partial differential equation that can be used to simulate the effects of smaller-scale atmospheric turbulence on sonic boom waveforms is described. The equation is first order in the time derivative and involves an extension of geometrical acoustics to include diffraction phenomena. Various terms in the equation are explained in physical terms. Such terms include those representing convection at the wave speed, diffraction, molecular relaxation, classical dissipation, and nonlinear steepening. The atmospheric turbulence enters through an effective sound speed, which varies with all three spatial coordinates, and which is the sum of the local sound speed and the component of the turbulent flow velocity projected along a central ray that connects the aircraft trajectory with the listener.

Pierce, Allan D.↗

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗

Simulating Bone Loss in Microgravity Using Mathematical Formulations of Bone Remodeling

Most mathematical models of bone remodeling are used to simulate a specific bone disease, by disrupting the steady state or balance in the normal remodeling process, and to simulate a therapeutic strategy. In this work, the ability of a mathematical model of bone remodeling to simulate bone loss as a function of time under the conditions of microgravity is investigated. The model is formed by combining a previously developed set of biochemical, cellular dynamics, and mechanical stimulus equations in the literature with two newly proposed equations; one governing the rate of change of the area of cortical bone tissue in a cross section of a cylindrical section of bone and one governing the rate of change of calcium in the bone fluid. The mechanical stimulus comes from a simple model of stress due to a compressive force on a cylindrical section of bone which can be reduced to zero to mimic the effects of skeletal unloading in microgravity. The complete set of equations formed is a system of first order ordinary differential equations. The results of selected simulations are displayed and discussed. Limitations and deficiencies of the model are also discussed as well as suggestions for further research.

Pennline, James A.↗