Search NASA⌕ Search

SEARCH · Search NASA

Results for “electromagnetic transient 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.

33 records · Page 2

Finite Element Modeling of Magnetically-Damped Convection during Solidification

A fully 3-D, transient finite element model is developed to represent the magnetic damping effects on complex fluid flow, heat transfer and electromagnetic field distributions in a Sn- 35.5%Pb melt undergoing unidirectional solidification. The model is developed based on our in- house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The numerical model is tested against numerical and experimental results for water as reported in literature. Various numerical simulations are carried out for the melt convection and temperature distribution with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to stabilize melt flow, reduce turbulence and flow levels in the melt and over a certain threshold value a higher magnetic field resulted in a greater reduction in velocity. Also, for the study of melt flow instability, a long enough running time is needed to ensure the final fluid flow recirculation pattern. Moreover, numerical results suggest that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the 0 convection in the melt is actually enhanced.

deGroh, H. C.↗

Swift-XRT Follow-up of Gravitational-wave Triggers in the Second Advanced LIGO/Virgo Observing Run

The Neil Gehrels Swift Observatory carried out prompt searches for gravitational-wave (GW) events detected by the LIGO/Virgo Collaboration (LVC) during the second observing run ("O2"). Swift performed extensive tiling of eight LVC triggers, two of which had very low false-alarm rates (GW170814 and the epochal GW170817), indicating a high confidence of being astrophysical in origin; the latter was the first GW event to have an electromagnetic counterpart detected. In this paper we describe the follow-up performed during O2 and the results of our searches. No GW electromagnetic counterparts were detected; this result is expected, as GW170817 remained the only astrophysical event containing at least one neutron star after LVC's later retraction of some events. A number of X-ray sources were detected, with the majority of identified sources being active galactic nuclei. We discuss the detection rate of transient X-ray sources and their implications in the O2 tiling searches. Finally, we describe the lessons learned during O2 and how these are being used to improve the Swift follow-up of GW events. In particular, we simulate a population of gamma-ray burst afterglows to evaluate our source ranking system's ability to differentiate them from unrelated and uncataloged X-ray sources. We find that ≈60%–70% of afterglows whose jets are oriented toward Earth will be given high rank (i.e., "interesting" designation) by the completion of our second follow-up phase (assuming that their location in the sky was observed), but that this fraction can be increased to nearly 100% by performing a third follow-up observation of sources exhibiting fading behavior.

N. J. Klingler↗

Data-driven Expectations for Electromagnetic Counterpart Searches Based on LIGO/Virgo Public Alerts

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic after glow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo’s third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics ofO3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers

Polina Petrov↗

Mechanisms for the Dissipation of Alfven Waves in Near-Earth Space Plasma

Alfven waves are a major mechanism for the transport of electromagnetic energy from the distant part of the magnetosphere to the near-Earth space. This is especially true for the auroral and polar regions of the Earth. However, the mechanisms for their dissipation have remained illusive. One of the mechanisms is the formation of double layers when the current associated with Alfven waves in the inertial regime interact with density cavities, which either are generated nonlinearly by the waves themselves or are a part of the ambient plasma turbulence. Depending on the strength of the cavities, weak and strong double layers could form. Such double layers are transient; their lifetimes depend on that of the cavities. Thus they impulsively accelerate ions and electrons. Another mechanism is the resonant absorption of broadband Alfven- wave noise by the ions at the ion cyclotron frequencies. But this resonant absorption may not be possible for the very low frequency waves, and it may be more suited for electromagnetic ion cyclotron waves. A third mechanism is the excitation of secondary waves by the drifts of electrons and ions in the Alfven wave fields. It is found that under suitable conditions, the relative drifts between different ion species and/or between electrons and ions are large enough to drive lower hybrid waves, which could cause transverse accelerations of ions and parallel accelerations of electrons. This mechanism is being further studied by means of kinetic simulations using 2.5- and 3-D particle-in-cell codes. The ongoing modeling efforts on space weather require quantitative estimates of energy inputs of various kinds, including the electromagnetic energy. Our studies described here contribute to the methods of determining the estimates of the input from ubiquitous Alfven waves.

Singh, Nagendra↗

Implementation and Testing of the First Prompt Search for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (OW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient OW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize OW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 172009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify OW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency OW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with approx.50% or better probability with a few pointings of wide-field telescopes. Key words. gravitational waves - methods: observational

Camp, Jordan↗

Nonlinear dynamics of attractive magnetic bearings

The nonlinear dynamics of a ferromagnetic shaft suspended by the force of attraction of 1, 2, or 4 independent electromagnets is presented. Each model includes a state variable feedback controller which has been designed using the pole placement method. The constitutive relationships for the magnets are derived analytically from magnetic circuit theory, and the effects of induced eddy currents due to the rotation of the journal are included using Maxwell's field relations. A rotor suspended by four electro-magnets with closed loop feedback is shown to have nine equilibrium points within the bearing clearance space. As the rotor spin speed increases, the system is shown to pass through a Hopf bifurcation (a flutter instability). Using center manifold theory, this bifurcation can be shown to be of the subcritical type, indicating an unstable limit cycle below the critical speed. The bearing is very sensitive to initial conditions, and the equilibrium position is easily upset by transient excitation. The results are confirmed by numerical simulation.

Hebbale, K. V.↗

Implementation and Testing of the First Prompt for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW-event candidates and to reconstruct-maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline s ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with 50% or better probability with a few pointings of wide-field telescopes.

Abadie, J.↗

Analysis and Simulations of Space Radiation Induced Single Event Effects and Transients

Spacecraft electronic are affected by the space radiation environment. Among the different types of radiation effects that can affect spacecraft electronics is the single event transient (SET). The space environment is responsible for many of the single event transients which can upset the performance of the spacecraft avionics hardware. In this chapter we first explore the origins of single event transients, then explore the modeling of a single event transient in digital and analog circuit. The chapter also addresses the concept of crosstalk that could develop among digital circuits in the present of a SET event. The chapter also provides a discussion of SET hardening. We then provide a discussion concerning propagation of a single event transient event at the local, subsystem, and system level in a spacecraft using two different models, one of the models developed by the author, known as the state transition model. The final goal of the chapter is to provide a qualitatively methodology for assessing single event transients and its effects so that spacecraft avionics engineers can develop either hardware or software countermeasures in their designs. SET is not a form of electromagnetic interference (EMI) in its origin, but semantically SET is very similar to EMI because they are both caused a current source not previously accounted for. SET has the same effects as EMI and it can cause interference problems in electronic circuits via multiple coupling mechanisms similar to EMI, and therefore makes such circuits incompatible.

Perez, Reinaldo J.↗

ISS And Space Environment Interactions Without Operating Plasma Contactor

The International Space Station (ISS) will be the largest, highest power spacecraft placed in orbit. Because of this the design of the electrical power system diverged markedly from previous systems. The solar arrays will operate at 160 V and the power distribution voltage will be 120 V. The structure is grounded to the negative side of the solar arrays so under the right circumstances it is possible to drive the ISS potential very negative. A plasma contactor has been added to the ISS to provide control of the ISS structure potential relative to the ambient plasma. The ISS requirement is that the ISS structure not be greater than 40 V positive or negative of local plasma. What are the ramifications of operating large structures with such high voltage power systems? The application of a plasma contactor on ISS controls the potential between the structure and the local plasma, preventing degrading effects. It is conceivable that there can be situations where the plasma contactor might be non-functional. This might be due to lack of power, the need to turn it off during some of the build-up sequences, the loss of functionality for both plasma contactors before a replacement can be installed, similar circumstances. A study was undertaken to understand how important it is to have the contactor functioning and how long it might be off before unacceptable degradation to ISS could occur. The details of interaction effects on spacecraft have not been addressed until driven by design. This was true for ISS. If the structure is allowed to float highly negative impinging ions can sputter exposed conductors which can degrade the primary surface and also generate contamination due to the sputtered material. Arcing has been known to occur on solar arrays that float negative of the ambient plasma. This can also generate electromagnetic interference and voltage transients. Much of the ISS structure and pressure module surfaces exposed to space is anodized aluminum. The anodization thickness is very thin to provide the required solar absorptance and emittance. For conditions where ISS structure can charge negative a large percentage of the array voltage, the dielectric strength of this layer is low, and dielectric breakdown (arcing) can occur. The energy stored capacitively in the structure can be delivered to the arc. The mechanisms by which this energy is delivered and how much of the energy is available hasn't been fully quantified. Questions have been raised regarding the possibility of whether a sustained arc might result due to current collected by the solar arrays from local plasma. It was postulated that even if dielectric breakdown didn't occur, impacts due to micrometeoroids and space debris could penetrate thin layers of dielectric on ISS and initiate an arc due to the coupling provided by the dense local plasma produced by the impact. This was proven in experiments conducted jointly by MSFC and Auburn University. A target chamber with a simulated ionospheric plasma and a biased, anodized aluminum plate and a 1-microfarad capacitor was used. The plate was then impacted by 75-micron particles accelerated to orbital velocity. Arc discharges were sustained for higher voltages but a threshold appears below which no discharge was initiated. Most items without an exposed power system will float electrically near the local plasma potential. This is true of the Space Shuttle, an Astronaut on EVA, and similar items. The structure of ISS might be at a large negative voltage. Therefore, capacitively stored energy can be transferred during docking, installing external boxes and equipment and Astronaut contact with ISS structure. The circumstances of when this can happen and the resulting effects are evaluated in this study. Also, a crewmember on EVA might be in the vicinity of an arc. All safety aspects of such an encounter including charging, molten particles from the arc site and EMI have been evaluated. This paper will report on the total results of this study focussed on the 4A configuration, scheduled to be complete in November, 2000. Interactions such as arcing, debris induced arcs, sustained arcs, sputtering, contamination from sputtering and arcing, docking interactions and Astronaut safety issues will all be addressed.

Carruth, M. R., Jr.↗

Power electronics for low power arcjets

In anticipation of the needs of future light-weight, low-power spacecraft, arcjet power electronics in the 100 to 400 W operating range were developed. Limited spacecraft power and thermal control capacity of these small spacecraft emphasized the need for high efficiency. Power topologies similar to those in the higher 2 kW and 5 to 30 kW power range were implemented, including a four transistor bridge switching circuit, current mode pulse-width modulated control, and an output current averaging inductor with an integral pulse generation winding. Reduction of switching transients was accomplished using a low inductance power distribution network, and no passive snubber circuits were necessary for power switch protection. Phase shift control of the power bridge was accomplished using an improved pulse width modulation to phase shift converter circuit. These features, along with conservative magnetics designs allowed power conversion efficiencies of greater than 92.5 percent to be achieved into resistive loads over the entire operating range of the converter. Electromagnetic compatibility requirements were not considered in this work, and control power for the converter was derived from AC mains. Addition of input filters and control power converters would result in an efficiency of on the order of 90 percent for a flight unit. Due to the developmental nature of arcjet systems at this power level, the exact nature of the thruster/power processor interface was not quantified. Output regulation and current ripple requirements of 1 and 20 percent respectively, as well as starting techniques, were derived from the characteristics of the 2 kW system but an open circuit voltage in excess of 175 V was specified. Arcjet integration tests were performed, resulting in successful starts and stable arcjet operation at power levels as low as 240 W with simulated hydrazine propellants.

Hamley, John A.↗

Foraging with MUSHROOMS: A Mixed-integer Linear Programming Scheduler for Multimessenger Target of Opportunity Searches with the Zwicky Transient Facility

Electromagnetic follow-up of gravitational-wave detections is very resource intensive, taking up hours of limited observation time on dozens of telescopes. Creating more efficient schedules for follow-up will lead to a commensurate increase in counterpart location efficiency without using more telescope time. Widely used in operations research and telescope scheduling, mixed-integer linear programming is a strong candidate to produce these higher-efficiency schedules, as it can make use of powerful commercial solvers that find globally optimal solutions to provided problems. We detail a new target-of-opportunity scheduling algorithm designed with Zwicky Transient Facility in mind that uses mixed-integer linear programming. We compare its performance to gwemopt, the tuned heuristic scheduler used by the Zwicky Transient Facility and other facilities during the third LIGO–Virgo gravitational-wave observing run. This new algorithm uses variable-length observing blocks to enforce cadence requirements and to ensure field observability, along with having a secondary optimization step to minimize slew time. We show that by employing a hybrid method utilizing both this scheduler and gwemopt, the previous scheduler used, in concert, we can achieve an average improvement in detection efficiency of 3%–11% over gwemopt alone for a simulated binary neutron star merger data set consistent with LIGO–Virgo's third observing run, highlighting the potential of mixed-integer target of opportunity schedulers for future multimessenger follow-up surveys.

B Parazin↗

Development of a Loop Antenna Deployment System for TEM-Based Subsurface Mars Water Detection

The subsurface of Mars has the potential to harbor existing deposits of liquid water, which are of great interest both scientifically (in the search for life) and as resources for future astronauts living on the red planet. However, the depth to a potential subsurface aquifer may be kilometers deep, which is well beyond what current surface or orbital approaches can survey while also confirming the unique signature for liquid water. Inspired by critical technology for exploring deep aquifers on Earth, we are developing the Transient H2O Reconnaissance (TH2OR) instrument at NASA’s Jet Propulsion Laboratory (JPL). TH2OR operates on the principle of transient electromagnetics (TEM), which leverages induction and electromotive force to induce a current in a subsurface water body using a loop-shaped antenna that provides both transmit and receive functions when placed on the surface. On Mars, TEM may be even more effective given the relative dryness of the subsurface compared to a more conductive, saline, liquid, water body on Earth. However, to probe deep within the subsurface (below 5 km) a large enclosed antenna is needed – specifically, a 100-m diameter effective loop or a shape with roughly equivalent area that can transmit at low frequencies (kHz-Hz). The deployment of a large-scale structure on the surface is complicated by the fact that the transmit wire must be both lightweight and robust to contact with the surface. Further, TH2OR may be delivered to the surface by a non-mobile vehicle, so it is desired that the deployment can be activated from a static location. This paper provides an overview of the deployment trade study, focusing on our current, favored approach, using a projectile wire launcher. Building on past approaches in the literature, we have developed and fielded an Earth-based, gas-projectile prototype for launching a triangle-shaped antenna onto analogue terrain. Our results compare simulated launch performance to actual field tests conducted under Earth gravity and pressure conditions. We discuss how Earth performance maps to a prospective Mars deployment under reduced gravity and pressure. We also provide lessons learned and next steps towards the development of an integrated TH2OR instrument for finding water on Mars.

Stamenkovic, Vlada↗

Non-linear resonance of fluids in a crystal growth cavity

In the microgravity environment, the effect of gravity on fluid motion is much reduced. Hence, secondary effects such as vibrations, jitters, surface tension, capillary effects, and electromagnetic forces become the dominant mechanism of fluid convection. Numerous studies have been conducted to investigate fluid behavior in microgravity with the ultimate goal of developing processes with minimal influence from convection. Industrial applications such as crystal growth from solidification of melt and protein growth for pharmatheutical application are just a few examples of the vast potential benefit that can be reaped from material processing in space. However, a space laboratory is not immune from all undesirable disturbances and it is imperative that such disturbances be well understood, quantifiable, and controlled. Non-uniform and transient accelerations such as vibrations, jitters, and impulsive accelerations exist as a result of crew activities, space vehicle maneuvering, and the operations of on-board equipment. Measurements conducted on-board a U.S. Spacelab showed the existence of vibrations in the frequency range of 1 to 100 Hz with a dominant mode of 17 Hz and harmonics of 54 Hz. The observed vibration is not limited to any coordinate plane but exists in all directions. Similar situation exists on-board the Russian MIR space station. Due to the large structure of its design, the future International Space Station will have its own characteristic vibration spectrum. It is well known that vibration can exert substantial influence on heat and mass transfer processes, thus hindering any attempts to achieve a diffusion-limited process. Experiments on vibration convection for a liquid-filled enclosure under one-g environment showed the existence of different flow regimes as vibration frequency and intensity changes. Results showed the existence of a resonant frequency, near which the enhancement is the strongest, and the existence of a high frequency asymptote. Numerical simulations of vibration convection have been conducted by Yurkov, Fu and Shieh, and by Wang. These analyses considered a two-dimensional air-filled cell under weightlessness condition and showed results similar to those of the experiments. It is not yet known whether resonance convection can be triggered by jitter alone or whether it requires the interaction of jitter with other convective forces in low gravity. An order of magnitude analysis, however, can be used to show the dependence of the resonance frequency on the fluid Prandtl number. Even though the onset of resonance convection may depend on other factors, results indicates that fluids with low Prandtl numbers are more susceptible to resonance than those with high Prandtl numbers. The current study is aimed at gaining additional insights to this problem using germanium as working fluid. Germanium was chosen for this analysis because of its common usage in solidification process and its relatively low Prandtl number (Pr = 0.02).

Wang, Francis C.↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗