Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed 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 271 records · Page 15

Propulsion Electric Grid Simulator (PEGS) for Future Turboelectric Distributed Propulsion Aircraft

NASA Glenn Research Center, in collaboration with the aerospace industry and academia, has begun the development of technology for a future hybrid-wing body electric airplane with a turboelectric distributed propulsion (TeDP) system. It is essential to design a subscale system to emulate the TeDP power grid, which would enable rapid analysis and demonstration of the proof-of-concept of the TeDP electrical system. This paper describes how small electrical machines with their controllers can emulate all the components in a TeDP power train. The whole system model in Matlab/Simulink was first developed and tested in simulation, and the simulation results showed that system dynamic characteristics could be implemented by using the closed-loop control of the electric motor drive systems. Then we designed a subscale experimental system to emulate the entire power system from the turbine engine to the propulsive fans. Firstly, we built a system to emulate a gas turbine engine driving a generator, consisting of two permanent magnet (PM) motors with brushless motor drives, coupled by a shaft. We programmed the first motor and its drive to mimic the speed-torque characteristic of the gas turbine engine, while the second motor and drive act as a generator and produce a torque load on the first motor. Secondly, we built another system of two PM motors and drives to emulate a motor driving a propulsive fan. We programmed the first motor and drive to emulate a wound-rotor synchronous motor. The propulsive fan was emulated by implementing fan maps and flight conditions into the fourth motor and drive, which produce a torque load on the driving motor. The stator of each PM motor is designed to travel axially to change the coupling between rotor and stator. This feature allows the PM motor to more closely emulate a wound-rotor synchronous machine. These techniques can convert the plain motor system into a unique TeDP power grid emulator that enables real-time simulation performance using hardware-in-the-loop (HIL).

electric power grid emulator↗

Background and Pickup Ion Velocity Distribution Dynamics in Titan's Plasma Environment: 3D Hybrid Simulation and Comparison with CAPS T9 Observations

In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.

Lipatov, A. S.↗

A Comparison of Filter-based Approaches for Model-based Prognostics

Model-based prognostics approaches use domain knowledge about a system and its failure modes through the use of physics-based models. Model-based prognosis is generally divided into two sequential problems: a joint state-parameter estimation problem, in which, using the model, the health of a system or component is determined based on the observations; and a prediction problem, in which, using the model, the stateparameter distribution is simulated forward in time to compute end of life and remaining useful life. The first problem is typically solved through the use of a state observer, or filter. The choice of filter depends on the assumptions that may be made about the system, and on the desired algorithm performance. In this paper, we review three separate filters for the solution to the first problem: the Daum filter, an exact nonlinear filter; the unscented Kalman filter, which approximates nonlinearities through the use of a deterministic sampling method known as the unscented transform; and the particle filter, which approximates the state distribution using a finite set of discrete, weighted samples, called particles. Using a centrifugal pump as a case study, we conduct a number of simulation-based experiments investigating the performance of the different algorithms as applied to prognostics.

prognostics↗

Protecting Venus from Asteroids, Comets, and Meteors

It is well accepted that the dense, thick atmosphere of Venus prevents most small cosmic bodies from reaching the surface and forming craters. We have examined this atmospheric intervention in detail, incorporating the lessons learned from the extensive modeling of impactor deceleration and flattening motivated by the SL-9 impacts with Jupiter. We employ a "pancake" model, which best matches detailed code simulations of atmospheric energy deposition, and Schmidt-Holsapple crater scaling modified for complex (flattened) craters. We adopt the distributions of Venus-crossing asteroids and comets determined by E.M. Shoemaker and co-workers, as well as generalizations of these distributions. Our nominal simulation of the venusian crater record is shown below, calibrated to the total number of venusian craters (940). As nearly all craters on Venus are well-preserved and relatively uniformly distributed, such simulations constrain the age of the surface. The fit is reasonable, with a nominal crater retention age of approx. 700 Ma. The fit at the large-crater end is improved if the number of large asteroids is increased, which Shoemaker argues is in fact more representative of the long-term (over several 100 Ma) average, and if Halley-family comets are included. The ages we obtain under a variety of modeling choices that produce good fits (including using Shoemaker's preferred crater scaling) are approx. 700-900 Ma, substantially greater than the most widely cited age estimate in the literature (-300 Ma). The key difference is that we find very large depletions in the production of 20-30-km craters (see figure) compared with previous estimates, the size range at which atmospheric effects are often calibrated or assumed nearly negligible. As venusian global resurfacing recedes deeper into history, the likelihood that Venus is resting between bouts of activity diminishes. Venus, like Mars, may instead be dying or dead.

McKinnon, William B.↗

Dynamically Evolving Sectors for Convective Weather Impact

A new strategy for altering existing sector boundaries in response to blocking convective weather is presented. This method seeks to improve the reduced capacity of sectors directly affected by weather by moving boundaries in a direction that offers the greatest capacity improvement. The boundary deformations are shared by neighboring sectors within the region in a manner that preserves their shapes and sizes as much as possible. This reduces the controller workload involved with learning new sector designs. The algorithm that produces the altered sectors is based on a force-deflection mesh model that needs only nominal traffic patterns and the shape of the blocking weather for input. It does not require weather-affected traffic patterns that would have to be predicted by simulation. When compared to an existing optimal sector design method, the sectors produced by the new algorithm are more similar to the original sector shapes, resulting in sectors that may be more suitable for operational use because the change is not as drastic. Also, preliminary results show that this method produces sectors that can equitably distribute the workload of rerouted weather-affected traffic throughout the region where inclement weather is present. This is demonstrated by sector aircraft count distributions of simulated traffic in weather-affected regions.

Drew, Michael C.↗

Intercomparison of Large-Eddy Simulations of Arctic Mixed-Phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds.

clouds↗

Radon-222 as a test of convective transport in a general circulation model

A three-dimensional tracer model based on the Goddard Institude of Space Studies GCM is used to simulate the distribution of Rn-222 over North America to test the ability of the model to describe the transport of pollutants in the boundary layer and the exchange of mass between the boundary layer and the free troposphere. The model results are compared with surface observations from five sites in the U.S., showing that Rn-222 concentrations are primarily regulated by dry convection. The simulations show satisfactory agreement with observations although the model underpredicts observations at night and the simulated Rn-222 concentrations over the northeastern U.S. are too high in the spring and too low in the fall.

Jacob, Daniel J.↗

Spectroscopic Observation of the Stardust Re-Entry in the Near UV with SLIT: Deduction of Surface Temperatures and Plasma Radiation

Thermal radiation of the heat-shield and the emission of the post-shock layer around the Stardust capsule, during its re-entry, were detected by a NASA-led observation campaign aboard NASA's DC-8 airborne observatory involving teams from several nations. The German SLIT experiment used a conventional spectrometer, in a Czerny-Turner configuration (300 mm focal length and a 600 lines/mm grating), fed by fiber optics, to cover a wavelength range from 324 nm to 456 nm with a pixel resolution of 0.08 nm. The reentering spacecraft was tracked m uansuinaglly a camera with a view angle of 20 degrees, and light from the capsule was collected using a small mirror telescope with a view angle of only 0.45 degrees. Data were gathered with a measurement frequency of 5 Hz in a 30-second time interval around the point of maximum heating until the capsule left the field of view. The emission of CN (as a major ablation product), N2(+) and different atoms were monitored successfully during that time. Due to the nature of the experimental set up, spatial resolution of the radiation field was not possible. Therefore, all measured values represent an integration of radiation from the visible part of the glowing heat shield, and from the plasma in the post-shock region. Further, due to challenges in tracking not every spectrum gathered contained data. The measured spectra can be split up into two parts: (i) continuum spectra which represent a superposition of the heat shield radiation and the continuum radiation of potential dust particles in the plasma, and (ii) line spectra from the plasma in the shock layer. Planck temperatures (interpreted as the surface temperatures of the Stardust heat shield) were determined assuming either a constant surface temperature, or a temperature distribution deduced from numerical simulation. The constant surface temperatures are in good agreement with numerical simulations, but the peak values at the stagnation point are significantly lower than those in the numerical simulation if a temperature distribution over the surface is assumed. Emission bands of CN and N2(+) were tracked along the visible trajectory and compared to a spectral simulation with satisfying agreement. Values for the integrated radiation of the transitions of interest for these species were extracted from this comparison.

Winter, Michael W.↗

The coupling of fluids, dynamics, and controls on advanced architecture computers

This grant provided for the demonstration of coupled controls, body dynamics, and fluids computations in a workstation cluster environment; and an investigation of the impact of peer-peer communication on flow solver performance and robustness. The findings of these investigations were documented in the conference articles.The attached publication, 'Towards Distributed Fluids/Controls Simulations', documents the solution and scaling of the coupled Navier-Stokes, Euler rigid-body dynamics, and state feedback control equations for a two-dimensional canard-wing. The poor scaling shown was due to serialized grid connectivity computation and Ethernet bandwidth limits. The scaling of a peer-to-peer communication flow code on an IBM SP-2 was also shown. The scaling of the code on the switched fabric-linked nodes was good, with a 2.4 percent loss due to communication of intergrid boundary point information. The code performance on 30 worker nodes was 1.7 (mu)s/point/iteration, or a factor of three over a Cray C-90 head. The attached paper, 'Nonlinear Fluid Computations in a Distributed Environment', documents the effect of several computational rate enhancing methods on convergence. For the cases shown, the highest throughput was achieved using boundary updates at each step, with the manager process performing communication tasks only. Constrained domain decomposition of the implicit fluid equations did not degrade the convergence rate or final solution. The scaling of a coupled body/fluid dynamics problem on an Ethernet-linked cluster was also shown.

Atwood, Christopher↗

Dispersion and Aging of Volcanic Aerosols after the La Soufriere Eruption in April 2021

Volcanic aerosols change the atmospheric composition and thereby affect weather and climate. Aerosol dynamic processes such as nucleation, condensation, and coagulation modify the shape, size, and mass of aerosol particles, which influence their atmospheric lifetime and radiative properties. Nevertheless, most models omit these processes for ash particles. In this work, we explore the ash aerosol aging and sulfate production during the first 4 days following the 2021 La Soufrière (St. Vincent) eruption with the ICON-ART model (ICOsahedral Nonhydrostatic model with Aerosol and Reactive Trace gases). Online coupling of ICON-ART with a one-dimensional volcanic plume model calculates volcanic emission, which makes it possible to resolve the different eruption phases of the noncontinuous La Soufrière eruption. We compared our simulated aerosol distribution and composition with observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument, the Multiangle Imaging SpectroRadiometer (MISR) Research Aerosol (RA) Algorithm, and the Barbados Cloud Observatory (BCO). We show that online coupling is essential to adequately model the emissions and plume development close to the volcano. The modeled aerosol aging is in very good agreement with observations from MISR near the emission source and with CALIOP at larger distances. Furthermore, particle aging occurs faster in the troposphere than in the stratosphere due to the availability of water vapor and OH, but a layer of coated ash appears at the plume top due to faster oxidation of SO2 and lofting by aerosol-radiation interaction. This paper gives the first direct comparison of aerosol aging in volcanic eruption plumes between simulations and observations

volcanic aerosols↗

Ion Velocity Distributions in Dipolarization Events: Beams in the Vicinity of the Plasma Sheet Boundary

Using combined MHD/test particle simulations, we further explore characteristic ion velocity distributions in relation to magnetotail reconnection and dipolarization events, focusing on distributions at and near the plasma sheet boundary layer (PSBL). Simulated distributions right at the boundary are characterized by a single earthward beam, as discussed earlier. However, farther inside, the distributions consist of multiple beams parallel and antiparallel to the magnetic field, remarkably similar to recent Magnetospheric Multiscale observations. The simulations provide insight into the mechanisms: the lowest earthward beam results from direct acceleration at an earthward propagating dipolarization front (DF), with a return beam at somewhat higher energy. A higher-energy earthward beam results from dual acceleration, first near the reconnection site and then at the DF, again with a corresponding return beam resulting from mirroring closer to Earth. Multiple acceleration at the X line or the propagating DF with intermediate bounces may produce even higher-energy beams. Particles contributing to the lower energy beams are found to originate from the PSBL with thermal source energies, increasing with increasing beam energy. In contrast, the highest-energy beams consist mostly of particles that have entered the acceleration region via cross-tail drift with source energies in the suprathermal range.

Birn, J.↗

Modeling Impacts of Dust Mineralogy on Fast Climate Response

Mineralogical composition drives dust impacts on Earth's climate systems. However, most climate models still use homogeneous dust, without accounting for the temporal and spatial variation in mineralogy. To quantify the radiative impact of resolving dust mineralogy on Earth's climate, we implement and simulate the distribution of dust minerals (i.e., illite, kaolinite, smectite, hematite, calcite, feldspar, quartz, and gypsum) from Claquin et al. (1999) (C1999) and activate their interaction with radiation in the GFDL AM4.0 model. Resolving mineralogy reduces dust absorption compared to the homogeneous dust used in the standard GFDL AM4.0 model that assumes a globally uniform hematite volume content of 2.7 % (HD27). The reduction in dust absorption results in improved agreement with observation-based single-scattering albedo (SSA), radiative fluxes from CERES (the Clouds and the Earth's Radiant Energy System), and land surface temperature from the CRU (Climatic Research Unit) compared to the baseline HD27 model version. It also results in distinct radiative impacts on Earth's climate over North Africa. Over the 19-year (from 2001 to 2019) modeled period during JJA (June–July–August), the reduction in dust absorption in AM4.0 leads to a reduction of over 50 % in net downward radiation across the Sahara and approximately 20 % over the Sahel at the top of the atmosphere (TOA) compared to the baseline HD27 model version. The reduced dust absorption weakens the atmospheric warming effect of dust aerosols and leads to an alteration in land surface temperature, resulting in a decrease of 0.66 K over the Sahara and an increase of 0.7 K over the Sahel. The less warming in the atmosphere suppresses ascent and weakens the monsoon inflow from the Gulf of Guinea. This brings less moisture to the Sahel, which combined with decreased ascent induces a reduction of precipitation. To isolate the effect of reduced absorption compared to resolving spatial and temporal mineralogy, we carry out a simulation where the hematite volume content of homogeneous dust is reduced from 2.7 % to 0.9 % (HD09). The dust absorption (e.g., single-scattering albedo) of HD09 is comparable to that of the mineralogically speciated model on a global mean scale, albeit with a lower spatial variation that arises solely from particle size. Comparison of the two models indicates that the spatial inhomogeneity in dust absorption resulting from resolving mineralogy does not have significant impacts on Earth's radiation and climate, provided there is a similar level of dust absorption on a global mean scale before and after resolving dust mineralogy. However, uncertainties related to emission and distribution of minerals may blur the advantages of resolving minerals to study their impact on radiation, cloud properties, ocean biogeochemistry, air quality, and photochemistry. On the other hand, lumping together clay minerals (i.e., illite, kaolinite, and smectite), but excluding externally mixed hematite and gypsum, appears to provide both computational efficiency and relative accuracy. Nevertheless, for specific research, it may be necessary to fully resolve mineralogy to achieve accuracy.

Soil dust↗

Evolving the SpaceFOM: Lessons Learned and Future Development

For over fifty years, simulation has been a cornerstone technology for space missions. In many cases, one single simulator, developed by one team, is not enough to meet the simulation objectives. Astronaut crew training, docking, vehicle testing and complex mission planning all require several simulations to be connected and to interoperate, in order to achieve the goal of the simulation. To facilitate the development of interoperable simulations, an open standard, the SISO SpaceFOM was developed in collaboration between government, industry and academia. The standard was released in 2020. The initial release of SpaceFOM focuses on the core interoperability topics for space simulation: (i) specifying spatial data using several, well defined reference frames, (ii) managing time (including real-time, hard real-time, and as-fast-as-possible execution), (iii) composition of systems with subsystems such as space vehicles, and (iv) execution control, such as initializing, running, pausing, and terminating simulations. A number of requirements for the next version of SpaceFOM have been identified and proposed solutions are already under development. One key area is the object classes used for data exchange, where a richer set of classes have been proposed, based on practical use cases. For time management, the use of variable time resolutions and variable time steps is also proposed. For execution control, save and restore of simulations is considered as well as improved robustness, to handle failing components. A number of clarifications are also under consideration. This paper provides both a forward-looking view of the technical developments underway and a retrospective on the standard's evolution. The authors emphasize that prioritizing core challenges and delivering a practical first version were crucial steps. While not every idea was implemented in the initial release, the timely delivery of a functional standard enabled real-world applications, the growth of a user community, and the collection of valuable insights to shape future versions. The authors look forward to working with the next version of SpaceFOM together with the growing user community.

FOMs↗

Models of Small-Scale Patchiness

Patchiness is perhaps the most salient characteristic of plankton populations in the ocean. The scale of this heterogeneity spans many orders of magnitude in its spatial extent, ranging from planetary down to microscale. It has been argued that patchiness plays a fundamental role in the functioning of marine ecosystems, insofar as the mean conditions may not reflect the environment to which organisms are adapted. Understanding the nature of this patchiness is thus one of the major challenges of oceanographic ecology. The patchiness problem is fundamentally one of physical-biological-chemical interactions. This interconnection arises from three basic sources: (1) ocean currents continually redistribute dissolved and suspended constituents by advection; (2) space-time fluctuations in the flows themselves impact biological and chemical processes, and (3) organisms are capable of directed motion through the water. This tripartite linkage poses a difficult challenge to understanding oceanic ecosystems: differentiation between the three sources of variability requires accurate assessment of property distributions in space and time, in addition to detailed knowledge of organismal repertoires and the processes by which ambient conditions control the rates of biological and chemical reactions. Various methods of observing the ocean tend to lie parallel to the axes of the space/time domain in which these physical-biological-chemical interactions take place. Given that a purely observational approach to the patchiness problem is not tractable with finite resources, the coupling of models with observations offers an alternative which provides a context for synthesis of sparse data with articulations of fundamental principles assumed to govern functionality of the system. In a sense, models can be used to fill the gaps in the space/time domain, yielding a framework for exploring the controls on spatially and temporally intermittent processes. The following discussion highlights only a few of the multitude of models which have yielded insight into the dynamics of plankton patchiness. In addition, this particular collection of examples is intended to furnish some exposure to the diversity of modeling approaches which can be brought to bear on the problem. These approaches range from abstract theoretical models intended to elucidate specific processes, to complex numerical formulations which can be used to actually simulate observed distributions in detail.

McGillicuddy, D. J.↗

Space simulator Patent

Space simulator with uniform test region radiation distribution, adapted to simulate Venus solar radiations

Thiele, C.↗

Determination of the geometrical factor and of the spatial distributions of particles for a cosmic ray telescope using a Monte Carlo approach

A program is presented which uses a Monte Carlo technique to compute the geometrical factor of a detector, and the spatial and angular distributions of simulated particles at designated layers within the detector. As the program is written, each layer is assumed to be rectangular, parallel to every other layer with its edges parallel to the edges of every other layer, and centered about a common axis. This program is written in Fortran 4 for the IBM 360 computers and is currently being run under the Fortran 4 (h) compiler. Same features such as the subroutine ERRSET, RANDU, and REMTIM and the type declarations REAL *8 and LOGICAL *1 may need to be changed for use with a different computer system. The use of input data cards, and a sample data set are presented. Lists of the symbol definition, card images of the FORTRAN object deck and subroutines are included along with a flow diagram.

Hemmer, M.↗

An analytical study of the effects of jets located more than one jet diameter above a wing at subsonic speeds

A procedure has been developed to calculate the effects of blowing two jets over a swept tapered wing at low subsonic speeds. The algorithm used is based on a vortex lattice representation of the wing lifting surface and a line sink-source distribution to simulate the effects of the jet exhaust on the wing lift and drag. The method is limited to those cases where the jet exhaust does not intersect or wash the wing. The predictions of this relatively simple procedure are in remarkably good agreement with experimentally measured interference lift and interference induced drag.

Putnam, L. E.↗