Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modelling and simulations”

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 199 records · Page 11

In-Flight Verification and Validation of Colloid Microthruster Performance

Colloid Micronewton Thrusters (CMNTs) use an electrospray to provide precision spacecraft position and pointing control. They were demonstrated in space for the first time as part of NASA’s Space Technology 7 (ST7) payload hosted by the European Space Agency’s (ESA’s) LISA Pathfinder (LPF) technology demonstration mission in January, 2016. CMNTs were the actuator in the disturbance reduction system (DRS) that provided drag-free operation of the LPF spacecraft, which will be necessary for future gravity wave observatories such as ESA’s Laser Interferometer Space Antenna (LISA) mission, currently in Phase A and scheduled for launch in 2034. The CMNT technology met performance requirements operating at 5-30 µN of thrust with 0.1 µN resolution and ≤0.1 µN/ÖHz thrust noise to deliver the required nanometer-level precision spacecraft control measured by the gravitational reference sensor (GRS) in the ESA LISA Technology Package (LTP). The performance of seven of the eight CMNTs in flight was consistent with ground test results, and as a system, all eight thrusters met mission-level performance requirements. The colloid microthruster performance model of thrust and thrust noise as a function of operational parameters (i.e. beam current, voltage, temperature, etc.) was validated in flight over a wide range of conditions. A model and simulation of the thruster control algorithm was developed and validated with flight data to predict thrust noise. This capability is important for future missions because it relates directly to the acceleration noise on the test masses, which provide the gravity wave measurements. The CMNT thruster model data and validation with LISA Pathfinder/ST7-DRS flight experiments are presented in this paper.

"O'Donnell, James"↗

State of the Art Assessment of Simulation in Advanced Materials Development

Advances in both the underlying theory and in the practical implementation of molecular modeling techniques have increased their value in the advanced materials development process. The objective is to accelerate the maturation of emerging materials by tightly integrating modeling with the other critical processes: synthesis, processing, and characterization. The aims of this report are to summarize the state of the art of existing modeling tools and to highlight a number of areas in which additional development is required. In an effort to maintain focus and limit length, this survey is restricted to classical simulation techniques including molecular dynamics and Monte Carlo simulations.

Wise, Kristopher E.↗

Megafires in a Warming World: What Wildfire Risk Factors Led to California’s Largest Recorded Wildfire

Massive wildfires and extreme fire behavior are becoming more frequent across the westernUnited States, creating a need to better understand how megafire behavior will evolve in our warmingworld. Here, the fire spread model Prometheus is used to simulate the initial explosive growth ofthe 2020 August Complex, which occurred in northern California (CA) mixed conifer forests. Hightemperatures, low relative humidity, and daytime southerly winds were all highly correlated withextreme rates of modeled spread. Fine fuels reached very dry levels, which accelerated simulationgrowth and heightened fire heat release (HR). Model sensitivity tests indicate that fire growth andHR are most sensitive to aridity and fuel moisture content. Despite the impressive early observedgrowth of the fire, shifting the simulation ignition to a very dry September 2020 heatwave predicted a>50% increase in growth and HR, as well as increased nighttime fire activity. Detailed model analysesof how extreme fire behavior develops can help fire personnel prepare for problematic ignitions.

Kevin Varga↗

Experiments in sensing transient rotational acceleration cues on a flight simulator

Results are presented for two transient motion sensing experiments which were motivated by the identification of an anomalous roll cue (a 'jerk' attributed to an acceleration spike) in a prior investigation of realistic fighter motion simulation. The experimental results suggest the consideration of several issues for motion washout and challenge current sensory system modeling efforts. Although no sensory modeling effort is made it is argued that such models must incorporate the ability to handle transient inputs of short duration (some of which are less than the accepted latency times for sensing), and must represent separate channels for rotational acceleration and velocity sensing.

Parrish, R. V.↗

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗

Modelling the solution growth of TGS crystals in low gravity

The experimental growth of triglycine sulfate (TGS) crystals from aqueous solution is modeled here in two dimensions using the PHOENICS finite volume code. Simulations are carried out for steady, impulsive, and periodic accelerations in order to determine tolerable acceleration levels. Scaling arguments are used to estimate the times required for thermal and solutal variations from the initial equilibrium state to be diffusively transported throughout the system, and to obtain order of magnitude information on the relative magnitudes of diffusive and convective transport. The computed concentration fields reflect the features of the concentration distributions found experimentally during experiments conducted aboard Spacelab 3 in 1985.

Nadarajah, Arunan↗

Three-Dimensional MHD Modeling of The Solar Corona and Solar Wind: Comparison with The Wang-Sheeley Model

We present simulation results from a tilted-dipole steady-state MHD model of the solar corona and solar wind and compare the output from our model with the Wang-Sheeley model which relates the divergence rate of magnetic flux tubes near the Sun (inferred from solar magnetograms) to the solar wind speed observed near Earth and at Ulysses. The boundary conditions in our model specified at the coronal base and our simulation region extends out to 10 AU. We assumed that a flux of Alfven waves with amplitude of 35 km per second emanates from the Sun and provides additional heating and acceleration for the coronal outflow in the open field regions. The waves are treated in the WKB approximation. The incorporation of wave acceleration allows us to reproduce the fast wind measurements obtained by Ulysses, while preserving reasonable agreement with plasma densities typically found at the coronal base. We find that our simulation results agree well with Wang and Sheeley's empirical model.

Usmanov, A. V.↗

The Brewer–Dobson Circulation in CMIP6

The Brewer–Dobson circulation (BDC) is a key feature of the stratosphere that models need to accurately represent in order to simulate surface climate variability and change adequately. For the first time, the Climate Model Intercomparison Project includes in its phase 6 (CMIP6) a set of diagnostics that allow for careful evaluation of the BDC. Here, the BDC is evaluated against observations and reanalyses using historical simulations. CMIP6 results confirm the well-known inconsistency in the sign of BDC trends between observations and models in the middle and upper stratosphere. Nevertheless, the large uncertainty in the observational trend estimates opens the door to compatibility. In particular, when accounting for the limited sampling of the observations, model and observational trend error bars overlap in 40 % of the simulations with available output. The increasing CO2 simulations feature an acceleration of the BDC but reveal a large spread in the middle-to-upper stratospheric trends, possibly related to the parameterized gravity wave forcing. The very close connection between the shallow branch of the residual circulation and surface temperature is highlighted, which is absent in the deep branch. The trends in mean age of air are shown to be more robust throughout the stratosphere than those in the residual circulation.

Brewer–Dobson circulation↗

Development and Calibration of a System-Integrated Rotorcraft Finite Element Model for Impact Scenarios

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber (DEA) under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. The presence of this energy absorbing device reduced the peak impact acceleration levels by a factor of three. Accelerations and kinematic data collected from the crash tests were compared to a system-integrated finite element model of the test article developed in parallel with the test program. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests were conducted to evaluate the impact performances of various components and subsystems, including new crush tubes and the DEA blocks. Parameters defined for the system-integrated finite element model were determined from these tests. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the full-scale crash test without the DEA. This combination of heuristic and quantitative methods identified modeling deficiencies, evaluated parameter importance, and proposed required model changes. The multidimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and copilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. One lesson learned was that this approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and pretest predictions. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, potentially reducing overall development costs.

Annett, Martin S.↗

Spacecraft Thermal and Optical Modeling Impacts on Estimation of the GRAIL Lunar Gravity Field

We summarize work performed involving thermo-optical modeling of the two Gravity Recovery And Interior Laboratory (GRAIL) spacecraft. We derived several reconciled spacecraft thermo-optical models having varying detail. We used the simplest in calculating SRP acceleration, and used the most detailed to calculate acceleration due to thermal re-radiation. For the latter, we used both the output of pre-launch finite-element-based thermal simulations and downlinked temperature sensor telemetry. The estimation process to recover the lunar gravity field utilizes both a nominal thermal re-radiation accleration history and an apriori error model derived from that plus an off-nominal history, which bounds parameter uncertainties as informed by sensitivity studies.

thermal modeling↗

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

Computational Study of Flow Establishment in a Ram Accelerator

The temporal evolution of the combustion process established during projectile transition from the launch tube into the ram accelerator section containing an explosive hydrogen-oxygen-argon gas mixture is studied. The Navier-Stokes equations for chemically reacting flow are solved in a fully coupled manner, using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order total variation diminishing scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of two ram accelerator configurations are then presented. In particular, the temporal developments of shock-induced combustion and thrust forces are followed. Positive thrust is established in both cases; however, in one of the ram accelerator configurations studied, combustion in the boundary layer enhances its separation, ultimately resulting in unstart.

Yungster, S.↗

Flight Simulator and Training Human Factors Validation

Loss of control has been identified as the leading cause of aircraft accidents in recent years. Efforts have been made to better equip pilots to deal with these types of events, commonly referred to as upsets. A major challenge in these endeavors has been recreating the motion environments found in flight as the majority of upsets take place well beyond the normal operating envelope of large aircraft. The Environmental Tectonics Corporation has developed a simulator motion base, called GYROLAB, that is capable of recreating the sustained accelerations, or G-forces, and motions of flight. A two part research study was accomplished that coupled NASA's Generic Transport Model with a GYROLAB device. The goal of the study was to characterize physiological effects of the upset environment and to demonstrate that a sustained motion based simulator can be an effective means for upset recovery training. Two groups of 25 Air Transport Pilots participated in the study. The results showed reliable signs of pilot arousal at specific stages of similar upsets. Further validation also demonstrated that sustained motion technology was successful in improving pilot performance during recovery following an extensive training program using GYROLAB technology.

Glaser, Scott T.↗

A numerical solution of the Navier-Stokes equations for supercritical fluid thermodynamic analysis

An explicit numerical solution of the compressible Navier-Stokes equations is applied to the thermodynamic analysis of supercritical oxygen in the Apollo cryogenic storage system. The wave character is retained in the conservation equations which are written in the basic fluid variables for a two-dimensional Cartesian coordinate system. Control-volume cells are employed to simplify imposition of boundary conditions and to ensure strict observance of local and global conservation principles. Non-linear real-gas thermodynamic properties responsible for the pressure collapse phenomonon in supercritical fluids are represented by tabular and empirical functions relating pressure and temperature to density and internal energy. Wall boundary conditions are adjusted at one cell face to emit a prescribed mass flowrate. Scaling principles are invoked to achieve acceptable computer execution times for very low Mach number convection problems. Detailed simulations of thermal stratification and fluid mixing occurring under low acceleration in the Apollo 12 supercritical oxygen tank are presented which model the pressure decay associated with de-stratification induced by an ordinary vehicle maneuver and heater cycle operation.

Heinmiller, P. J.↗

Vibration of a liquid with a free surface in a spinning spherical tank

The problem of forced fluid vibrations in a partially filled spinning spherical tank is solved numerically by using the finite element method. The governing equations include Coriolis acceleration and spatially homogeneous vorticity. An exponential instability is detected in the present simulation for fill ratios below 0.5 and centrifugal acceleration to thrust ratios less than 1.7. This fictitious instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.

El-Raheb, M.↗

Deictic primitives for general purpose navigation

A visually-based deictic primative used as an elementary command set for general purpose navigation was investigated. It was shown that a simple 'follow your eyes' scenario is sufficient for tracking a moving target. Limitations of velocity, acceleration, and modeling of the response of the mechanical systems were enforced. Realistic paths of the robots were produced during the simulation. Scientists could remotely command a planetary rover to go to a particular rock formation that may be interesting. Similarly an expert at plant maintenance could obtain diagnostic information remotely by using deictic primitives on a mobile are used in the deictic primitives, we could imagine that the exact same control software could be used for all of these applications.

Crismann, Jill D.↗

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗