Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric fluid dynamics”

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 325 records · Page 18

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Integrating an Industrial Source and Commercial Algae Farm with Innovative CO 2 Transfer Membrane and Improved Strain Technologies

This report describes the overall findings of the research project “Integrating an Industrial Source and Commercial Algae Farm with Innovative CO 2 Transfer Membrane and Improved Strain Technologies”. The motivation and goal for this project were to increase the carbon utilization efficiency (CUE) and areal productivity for algal cultivations, thereby reducing CO 2 costs to cultivation operations and improving economics. This was achieved through a combination of enhanced delivery of inorganic carbon and improved strains of Nannochloropsis oceanica capable of higher rates of bicarbonate uptake and metabolism. The project succeeded in these goals. First, a bubble-free, membrane-based technology was developed for delivering CO 2 to cultivations, which increased the CUE from the 15-20% that is standard in the industry to more than 65%. In addition, N. oceanica was modified to express a bicarbonate transporter protein, BicA, which enabled the cells to grow more quickly. In addition, advances were made in the protein engineering of carbonic anhydrase (CA) for enhanced stability and catalytic performance, the computational fluid dynamics modeling of algal cultivation systems, and in the life-cycle assessment and technoeconomic analysis of algal production. Together, these outcomes contribute to advancing algal cultivation as an economically viable platform for production of fuels, materials, and other chemical products. The project results aid in addressing the dual challenges of reducing atmospheric CO 2 levels and achieving green energy solutions.

09 BIOMASS FUELS↗

Assessment of the Reconstructed Aerodynamics of the Mars Science Laboratory Entry Vehicle

On August 5, 2012, the Mars Science Laboratory entry vehicle successfully entered Mars atmosphere, flying a guided entry until parachute deploy. The Curiosity rover landed safely in Gale crater upon completion of the Entry Descent and Landing sequence. This paper compares the aerodynamics of the entry capsule extracted from onboard flight data, including Inertial Measurement Unit (IMU) accelerometer and rate gyro information, and heatshield surface pressure measurements. From the onboard data, static force and moment data has been extracted. This data is compared to preflight predictions. The information collected by MSL represents the most complete set of information collected during Mars entry to date. It allows the separation of aerodynamic performance from atmospheric conditions. The comparisons show the MSL aerodynamic characteristics have been identified and resolved to an accuracy better than the aerodynamic database uncertainties used in preflight simulations. A number of small anomalies have been identified and are discussed. This data will help revise aerodynamic databases for future missions and will guide computational fluid dynamics (CFD) development to improved prediction codes.

Schoenenberger, Mark↗

Quantifying Turbine-Level Risk to Golden Eagles Using a High-Fidelity Updraft Model and a Stochastic Behavioral Model

To minimize the effects of wind farms on Golden Eagle (Aquila chrysaetos) populations while enabling sustainable development of renewable energy resources, it is important to understand how eagles interact with atmospheric flows, terrain features, and anthropogenic structures. Models that predict migratory flight paths provide one tool that helps us grasp how the location of wind farms may influence interactions and impacts on migrating Golden Eagles. The current state-of-the-art in predicting migratory flight paths uses a deterministic fluid-flow analogy to predict eagle trajectory using only an orographic updraft potential computed from topographical features. This model does not take into account variables, such as thermal updrafts and time varying atmospheric conditions that are known to influence migratory behavior. In this work, we improve on the model with the objective of developing tools that advance our understanding of how atmospheric flows and terrain features affect migratory eagle behavior and their interactions with wind farms. Specifically, we 1) incorporate both orographic and thermal updraft information in simulating eagle flight paths; 2) incorporate stochasticity into eagle travel patterns to better capture the influence of exogenous factors on, and the inherent stochasticity of eagle behavior; 3) consider spatio-temporal atmospheric data at wind-farm-scale when computing updraft potential; and 4) account for how atmospheric conditions and the direction of migration change seasonally and how these changes affect eagle migratory flight behavior. We tested the model using a 50km by 50km region with 50 m resolution in the western United States. We simulated 900 independent, probabilistic eagle tracks during southerly and northerly migration, assuming eagles solely rely on orographic updrafts. The preliminary results indicate that the inclusion of finer resolution atmospheric data allows for the inclusion of realistic conditions that an eagle experiences. The stochasticity in eagle tracks provides a platform to include uncertainty in eagle decision making and help produce robust eagle presence maps. We will deploy updraft and downdraft velocities computed using a high-fidelity, wind farm scale, computational fluid dynamics solver under development at National Renewable Energy Laboratory. This work is a first step in the development of a predictive and generalizable eagle behavior model at the wind farm scale that does not rely on empirical data collection. Although the current model is intended for migratory eagles, we will extend and refine this model to inform the development of additional behavioral modes, including resident eagle behavior. This modeling approach improves our ability to understand eagle use of the landscape at a fine scale, and it is our hope that this work will ultimately help advance strategies that minimize the impact of wind development on Golden Eagle populations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗

The Effect of Pressure on NO x Entitlement and Reaction Timescales in a Premixed Axial Jet-In-Crossflow

This paper investigates the pressure dependency of a lean premixed jet injected into a lean vitiated crossflow with an experimentally verified detailed chemistry computational fluid dynamics (CFD) model and 53 species considered. Experimental data were taken in an axially staged combustor with an optically accessible test section, allowing the use of particle image velocimetry (PIV) and CH* chemiluminescence techniques as well as point measurement of species concentration, temperature, and pressure. The experimental data cases at one, three, and five atmospheres were selected to describe the flame stabilization dependency on pressure and gain the required knowledge for an extrapolation to engine condition. Simulated exit nitrogen oxide levels were validated with experimental emission data, and a global emission trend for the NO reduction at elevated pressure and constant turbine inlet temperature level was defined. The nitrogen oxide benefit at elevated operating pressure was justified with the significantly smaller flame surface area: the analysis of the simulated spanwise and top-view profiles showed a relatively short receded core flame with nitrogen oxide production in the center at high pressure relative to a longer and larger shear layer flame at atmospheric condition that produced NO toward the inner and outer side of the flame. Furthermore, decomposition of the Damköhler number revealed the strong influence of the reaction timescales with higher reaction rates at elevated pressure, along with a moderate influence of the turbulent timescales, showing higher turbulence intensity in the lee-side recirculation zone at lower pressure.

03 NATURAL GAS↗

Computation of wake/exhaust mixing downstream of advanced transport aircraft

The mixing of engine exhaust with the vortical wake of high speed aircraft operating in the stratosphere can play an important role in the formation of chemical products that deplete atmospheric ozone. An accurate analysis of this type of interaction is therefore necessary as a part of the assessment of the impact of proposed High Speed Civil Transport (HSCT) designs on atmospheric chemistry. This paper describes modifications to the parabolic Navier-Stokes flow field analysis in the UNIWAKE unified aircraft wake model to accommodate the computation of wake/exhaust mixing and the simulation of reacting flow. The present implementation uses a passive chemistry model in which the reacting species are convected and diffused by the fluid dynamic solution but in which the evolution of the species does not affect the flow field. The resulting analysis, UNIWAKE/PCHEM (Passive CHEMistry) has been applied to the analysis of wake/exhaust flows downstream of representative HSCT configurations. The major elements of the flow field model are described, as are the results of sample calculations illustrating the behavior of the thermal exhaust plume and the production of species important to the modeling of condensation in the wake. Appropriate steps for further development of the UNIWAKE/PCHEM model are also outlined.

Quackenbush, Todd R.↗

Aerodynamics of the Mars Microprobe Entry Vehicles

The selection of the unique aeroshell shape for the Mars Microprobes is discussed. A description of its aerodynamics in hypersonic rarefied, hypersonic continuum, supersonic and transonic flow regimes is then presented. This description is based on Direct Simulation Monte Carlo analyses in the rarefied-flow regime, thermochemical nonequilibrium Computational Fluid Dynamics in the hypersonic regime, existing wind tunnel data in the supersonic and transonic regime, additional computational work in the transonic regime, and finally, ballistic range data. The aeroshell is shown to possess the correct combination of aerodynamic stability and drag to convert the probe's initial tumbling attitude and high velocity at atmospheric-interface into the desired surface-impact orientation and velocity.

Mitcheltree, R. A.↗

Space Launch Vehicle Transient Particle Redistribution Modeling and Implications for Optically Sensitive Payloads

The performance of contamination sensitive components—such as optical components—can be degraded by particulate matter depositing on the surfaces. Particles can accumulate during manufacturing, handling and operation. For a space-based system, particles can shed from the fairing and redistribute onto sensitive surfaces during launch. An engineering modeling approach has been developed for modeling particle migration during launch. The approach involves particle detachment from the fairing, particle transport through the venting atmosphere inside the fairing, and attachment to the receiving surface. Particle size and amounts on the fairing surface can be modeled using distributions from standards, such as IEST-STD-CC1246E, as well as from empirical data obtained from tape lifts. Surface interactions are modeled using theoretical as well as empirical data. Commercial computational fluid dynamics codes are used to calculate the gas flow in the fairing during depressurization during launch. This approach not only provides insight into particle redistribution during launch but also can be used to establish fairing cleanliness requirements.

Brieda, Lubos↗

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

Atmospheric dynamics of the outer planets

Despite major differences in the solar and internal energy inputs, the atmospheres of the four Jovian planets all exhibit latitudinal banding and high-speed jet streams. Neptune and Saturn are the windiest planets, Jupiter is the most active, and Uranus is a tipped-over version of the others. Large oval storm systems exhibit complicated time-dependent behavior that can be simulated in numerical models and laboratory experiments. The largest storm system, the Great Red Spot of Jupiter, has survived for more than 300 years in a chaotic shear zone where smaller structures appear and dissipate every few days. Future space missions will add to the understanding of small-scale processes, chemical composition, and vertical structure. Theoretical hypotheses about the interiors provide input for fluid dynamical models that reproduce many observed features of the winds, temperatures, and cloud patterns. In one set of models the winds are confined to the thin layer where clouds form. In other models, the winds extend deep into the planetary fluid interiors. Hypotheses will be tested further as observations and theories become more exact and detailed comparisons are made.

Ingersoll, Andrew P.↗

Ascent Aborts 2 Flight, Reconstruction of Structural Dynamics (Couple Loads Transients)

- This work describes the elements of AA2 post-flight data analysis, couple loads analysis, and comparisons of the flight data with the results of analytical methods transiently. - The measured transient responses ranged from rigid body frequencies (zero Hz) to 150 Hz. - The analytical methods consisted of defining and determining all forcing functions, i.e., buffet, aerodynamics, thrusts from the abort motor and the carrier lunch vehicle motor, and the initial conditions at pre-separation. - Mathematical models representing the mechanics of Orion abort is relatively complex because 1) two motors (carrier’s and LAS’s) are operating at the same time at pre-separation, and 2) the separation occurs in atmosphere. - Other challenges consist of determining the aerodynamics with plume-on condition at every millisecond of the flight ranging from transonic to supersonic (i.e., Mach 1.1 to Mach 1.6), and defining the mechanics of separation of LAS from the launch vehicle.

Launch Vehicle↗

Risk Analysis Associated with Loss of Toxic Gases During Orion Landing and Recovery Operations

Mission, landing and recovery operations for the Orion crew module involve reentry into the Earth's atmosphere and the deployment of three Nomex parachutes to slow the descent before landing along the west coast of the United States. Orion may have residual fuel (hydrazine, N2H4) or coolant (ammonia, NH3) on board which are both highly toxic to crew in the event of exposure. These risks were evaluated using a first principles analysis approach through fluid dynamics modeling. Plume calculations were first performed with the ANSYS Fluent computational fluid dynamics code. Data were then extracted at locations relevant to crew safety such as the snorkel fan inlet and the egress hatch. Mixing calculations were performed to quantify exposure concentrations within the crew bay before and during egress and departure. Finally, results included herein were used to inform the Orion post-landing Concept of Operations (ConOps) so that strategies could be formulated to maintain crew safety in the event of the loss of fuel or coolant.

Orion↗

Numerical Evaluation of Entry System Trajectory Control via Active Porosity Control of Transpiration Cooled Thermal Protection System

As space exploration expands, atmospheric descent requirements will tighten to improve landing accuracy to areas of interest. One way to meet higher landing accuracy requirements is by increasing vehicle maneuverability. This study explores the concept of manipulating a transpiration cooling scheme for entry capsules that could provide both active aerodynamic control as well as thermal protection for external surfaces. Numerical methods of fluid dynamics and heat transfer are presented to study the effectiveness of the concept in both thermal protection and entry vehicle attitude control. This effort works through development of an entry vehicle model at aeromaneuvering conditions with designated boundary activated for effects of transpiration cooling. Preliminary results show calculated values of heat flux and aerodynamic forces as a function of varying transpiration cooling flow rates. This contributes evidence of asymmetric cooling utilization as an additional mechanism for vehicle attitude control.

Caroline Anderson↗

The Aqua-Planet Experiment (APE): CONTROL SST Simulation

Climate simulations by 16 atmospheric general circulation models (AGCMs) are compared on an aqua-planet, a water-covered Earth with prescribed sea surface temperature varying only in latitude. The idealised configuration is designed to expose differences in the circulation simulated by different models. Basic features of the aqua-planet climate are characterised by comparison with Earth. The models display a wide range of behaviour. The balanced component of the tropospheric mean flow, and mid-latitude eddy covariances subject to budget constraints, vary relatively little among the models. In contrast, differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical uncertainty in modelled lower stratospheric temperatures persists in APE.Aspects of the circulation generated more directly by interactions between the resolved fluid dynamics and parameterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter-tropical convergence zone (ITCZ) at the equator, with large variations in mean precipitation. The equatorial wave spectrum shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propagation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive than the equatorial Rossby modes, dominates in a few models or occurs within an eastward propagating envelope in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies.The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top of atmosphere global net flux, dominated by differences in shortwave reflection by clouds. A number of newly developed models, not optimised for Earth climate, contribute to this. Possible reasons for differences in the optimised models are discussed.The aqua-planet configuration is intended as one component of an experimental hierarchy used to evaluate AGCMs. This comparison does suggest that the range of model behaviour could be better understood and reduced in conjunction with Earth climate simulations. Controlled experimentation is required to explore individual model behavior and investigate convergence of the aqua-planet climate with increasing resolution.

AGCM↗

A simulation of saltation dynamics with application to terrestrial and Martian model surfaces.

The feasibility of predicting the saltation threshold from a dynamic model of the saltation process is demonstrated. The model is based on mechanical equations that describe the trajectories of particles (sand grains) near the ground-atmosphere interface. The equations of particle motion are formulated in terms of forces and moments. The forces employed are based on theoretical and experimental results obtained for continuum flow. Surface structuring analytics and simulated fluid turbulence further modify the particle behavior. The model is used to predict the threshold velocities on Mars.

Andres, R. M.↗

Experimental investigation of unsteady flow dynamics in therminol-66 using edge tone-induced oscillations

Recent research has demonstrated that flow oscillations induced by a wedge positioned opposite the jet outlet in edge tones can be effectively utilized for flow measurement purposes. This highlights the necessity for developing advanced instrumentation capable of precise flow measurements in experimental setups exposed to radiation and in challenging fluids such as molten salts. Here, this study contributes to the existing body of knowledge by conducting experiments using Therminol-66, an effective surrogate for FLiBe under appropriate temperature conditions for nuclear applications, across a broad range of Reynolds numbers from 17,700 to 555,000. Experiments were conducted with Therminol-66 at approximately 368.8 K under atmospheric pressure, with varying flow rates. A reliable correlation was established between the generated frequencies, corresponding flow rates, and device geometry parameters. The data were subjected to non-dimensional analysis to reveal the interrelationships among flow rate, fluid properties, device geometry, and their impact on the acoustic tones produced by the edge tone device. These findings contribute to the development of precise flow measurement instrumentation for use in challenging environments.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Saturn's Magnetosphere Interaction with Titan for T9 Encounter: 3D Hybrid Modeling and Comparison with CAPS Observations

Global dynamics of ionized and neutral gases in the environment of Titan plays an important role in the interaction of Saturn s magnetosphere with Titan. Several hybrid simulations of this problem have already been done (Brecht et al., 2000; Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Modolo and Chanteur, 2008). Observational data from CAPS for the T9 encounter (Sittler et al., 2009) indicates an absence of O(+) heavy ions in the upstream that change the models of interaction which were discussed in current publications (Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Ma et al., 2007; Szego et al., 2007). Further analysis of the CAPS data shows very low density or even an absence of H(+) ions in upstream. In this paper we discuss two models of the interaction of Saturn s magnetosphere with Titan: (A) high density of H(+) ions in the upstream flow (0.1/cu cm), and (B) low density of H(+) ions in the upstream flow (0.02/cu cm). The hybrid model employs a fluid description for electrons and neutrals, whereas a particle approach is used for ions. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. The model atmosphere includes exospheric H(+), H(2+), N(2+)and CH(4+) pickup ion production as well as an immobile background ionosphere and a shell distribution for active ionospheric ions (M(sub i)=28 amu). The hybrid model allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures. The results of the ion dynamics in Titan s environment are compared with Cassini T9 encounter data (CAPS).

Lipatov, A. S.↗