Search NASA⌕ Search

SEARCH · Search NASA

Results for “differential drag”

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 145 records · Page 8

Supersonic Aerodynamic Characteristics of a Low-Drag Aircraft Configuration having an Arrow Wing of Aspect Ratio 1.86 and a Body of Fineness Ratio 20

A free-flight rocket-propelled-model investigation was conducted at Mach numbers of 1.2 to 1.9 to determine the longitudinal and lateral aero-dynamic characteristics of a low-drag aircraft configuration. The model consisted of an aspect-ratio -1.86 arrow wing with 67.5 deg. leading-edge sweep and NACA 65A004 airfoil section and a triangular vertical tail with 60 deg. sweep and NACA 65A003 section in combination with a body of fineness ratio 20. Aerodynamic data in pitch, yaw, and roll were obtained from transient motions induced by small pulse rockets firing at intervals in the pitch and yaw directions. From the results of this brief aerodynamic investigation, it is observed that very slender body shapes can provide increased volumetric capacity with little or no increase in zero-lift drag and that body fineness ratios of the order of 20 should be considered in the design of long-range supersonic aircraft. The zero-lift drag and the drag-due-to-lift parameter of the test configuration varied linearly with Mach number. The maximum lift-drag ratio was 7.0 at a Mach number of 1.25 and decreased slightly to a value of 6.6 at a Mach number of 1.81. The optimum lift coefficient, normal-force-curve slope, lateral-force-curve slope, static stability in pitch and yaw, time to damp to one-half amplitude in pitch and yaw, the sum of the rotary damping derivatives in pitch and also in yaw, and the static rolling derivatives all decreased with an increase in Mach number. Values of certain rolling derivatives were obtained by application of the least-squares method to the differential equation of rolling motion. A comparison of the experimental and calculated total rolling-moment-coefficient variation during transient oscillations of the model indicated good agreement when the damping-in-roll contribution was included with the static rolling-moment terms.

Gillespie, Warren, Jr.↗

Permeability of Two Parachute Fabrics - Measurements, Modeling, and Application

Two parachute fabrics, described by Parachute Industry Specifications PIA-C-7020D Type I and PIA-C-44378D Type I, were tested to obtain their permeabilities in air (i.e., flow-through volume of air per area per time) over the range of differential pressures from 0.146 psf (7 Pa) to 25 psf (1197 Pa). Both fabrics met their specification permeabilities at the standard differential pressure of 0.5 inch of water (2.60 psf, 124 Pa). The permeability results were transformed into an effective porosity for use in calculations related to parachutes. Models were created that related the effective porosity to the unit Reynolds number for each of the fabrics. As an application example, these models were used to calculate the total porosities for two geometrically-equivalent subscale Disk-Gap-Band (DGB) parachutes fabricated from each of the two fabrics, and tested at the same operating conditions in a wind tunnel. Using the calculated total porosities and the results of the wind tunnel tests, the drag coefficient of a geometrically-equivalent full-scale DGB operating on Mars was estimated.

Cruz, Juan R.↗

Permeability of Two Parachute Fabrics: Measurements, Modeling, and Application

Two parachute fabrics, described by Parachute Industry Specifications PIA-C-7020D Type I and PIA-C-44378D Type I, were tested to obtain their permeabilities in air (i.e., flow-through volume of air per area per time) over the range of differential pressures from 0.146 psf (7 Pa) to 25 psf (1197 Pa). Both fabrics met their specification permeabilities at the standard differential pressure of 0.5 inch of water (2.60 psf, 124 Pa). The permeability results were transformed into an effective porosity for use in calculations related to parachutes. Models were created that related the effective porosity to the unit Reynolds number for each of the fabrics. As an application example, these models were used to calculate the total porosities for two geometrically-equivalent subscale Disk-Gap-Band (DGB) parachutes fabricated from each of the two fabrics, and tested at the same operating conditions in a wind tunnel. Using the calculated total porosities and the results of the wind tunnel tests, the drag coefficient of a geometrically-equivalent full-scale DGB operating on Mars was estimated.

Cruz, Juan R.↗

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗

Applying Physics-Informed Neural Networks to Solve Navier–Stokes Equations for Laminar Flow around a Particle

In recent years, Physics-Informed Neural Networks (PINNs) have drawn great interest among researchers as a tool to solve computational physics problems. Unlike conventional neural networks, which are black-box models that “blindly” establish a correlation between input and output variables using a large quantity of labeled data, PINNs directly embed physical laws (primarily partial differential equations) within the loss function of neural networks. By minimizing the loss function, this approach allows the output variables to automatically satisfy physical equations without the need for labeled data. The Navier–Stokes equation is one of the most classic governing equations in thermal fluid engineering. This study constructs a PINN to solve the Navier–Stokes equations for a 2D incompressible laminar flow problem. Flows passing around a 2D circular particle are chosen as the benchmark case, and an elliptical particle is also examined to enrich the research. The velocity and pressure fields are predicted by the PINNs, and the results are compared with those derived from Computational Fluid Dynamics (CFD). Additionally, the particle drag force coefficient is calculated to quantify the discrepancy in the results of the PINNs as compared to CFD outcomes. The drag coefficient maintained an error within 10% across all test scenarios.

Hu, Beichao (ORCID:0009000163215151)↗

Dust particles from comets and asteroids collected at the Earth's orbit: Parent-daughter relationships

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not well known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of Asteroids in the main belt. It is believed that one may combine lab analysis of the physics and chemistry of captured particles with orbital data in order to identify comet and asteroid parent bodies. It is possible to use the collected orbits of the dust to connect with its source in two ways. One is to consider the long time orbit evolution of the dust under Poynting-Robertson drag. The other is to look at the prompt orbit change of dust from comets onto trajectories that intersect the earth's orbit. In order to characterize the orbits of dust particles evolved over a long period of time, a study of its orbital evolution was undertaken. Various parameters associated with these dust orbits as they cross the Earth's orbit were considered in order to see if one may discriminate between particles evolved from comets and asteroids. The method was to calculate by a numerical procedure the orbits of dust particles after they left their parent bodies. It appears that as the particles pass the Earth's orbit, asteroidal grains and cometary grains can be differentiated on the basis of their measured orbital eccentricities even after much planetary perturbation. Broad parent daughter associations can be made on this basis from measurement of their trajectories intercepted in earth orbit.

Jackson, A. A.↗

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental iterative approach for differentiating advanced flow codes is successfully demonstrated on a two-dimensional inviscid model problem. The method employs the reverse-mode capability of the automatic differentiation software tool ADIFOR 3.0 and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straightforward, black-box reverse-mode applicaiton of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-rder aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoinct) procedures; then, a very efficient noniterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hesian matrices) of lift, wave drag, and pitching-moment coefficients are calculated with respect to geometric shape, angle of attack, and freestream Mach number.

Taylor, Arthur C., III↗

Potential Vorticity Mixing in a Tangled Magnetic Field

A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is introduced. The analysis is in the context of β -plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfvénizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field–resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfvénization—i.e., before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose here a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗

Mass property control of a spin stabilized spacecraft with restrictive mission and weight constraints

In the primary experiment of the spin stabilized San Marco D/L spacecraft, the drag effects on a light spherical shell coupled to a relatively massive center body will be measured. To achieve the precise mass property control, the centroids of both the shell and the center body must coincide with each other and with the center of pressure of the shell. Precise spin balancing is needed for launch and orbital stability, and the deployable antennas need accurate alignment. Corrective measures, developed after the preliminary mass measurements showed flaws in the mass property control, are described in detail. Inertia booms and a yo-yo despin system were developed as add-on units, spin balance measurements about all three geometric axes were used to define weight minimized correction within the outer shell envelope, and boom tip mass differentials were optimized for the most favorable inertia ratio margin achievable within the mission weight constraints. The weight versus orbital lifetime trade-offs were also considered.

Lang, W. E.↗

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental-iterative approach for differentiating advanced flow codes is successfully demonstrated on a 2D inviscid model problem. The method employs the reverse-mode capability of the automatic- differentiation software tool ADIFOR 3.0, and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straight-forward, black-box reverse- mode application of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-order aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoint) procedures; then, a very efficient non-iterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hessian matrices) of lift, wave-drag, and pitching-moment coefficients are calculated with respect to geometric- shape, angle-of-attack, and freestream Mach number

Taylor, Arthur C., III↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗

Machine Learning based Aircraft Performance Model Estimation for Trajectory Prediction

The accurate prediction of aircraft trajectory by ground-based decision support tools is a critical component of air traffic management in the US National Airspace System (NAS). Accurate predictions of where the aircraft will be in the future or when they will arrive at specific locations (e.g., fixes) is a key enabler for sequencing and efficient arrival management of flights. Traditional physics based aircraft trajectory prediction relies on a simplified point-mass total energy model whose parameters are referred to as Aircraft Performance Model (APM) parameters. Even though the performance coefficients and weight of an aircraft are a vital part of the aircraft performance model’s predictions and accuracy, these coefficients are proprietary in nature and therefore, unavailable to decision-support tools. Current approaches freeze some coefficients to default base of aircraft data (BADA) values and optimize others. However, the APM parameters are highly coupled by the flight dynamics and prioritizing one parameter over others leads to bias and skewed predictions. To alleviate this problem, we provide a combined optimization framework to predict all the critical (thrust, drag and weight) APM parameters. This paper is focused on training Machine Learning (ML) models that map historical flights to optimized APM parameters that provide the best fit (in terms of prediction error). Our dataset obtained from NASA’s Sherlock data warehouse is comprised of thousands of historical flights and includes weather and track data collected from 2019. Using different subsets of relevant features (e.g., aircraft type), we trained several ML models to estimate the aircraft’s take off weight, drag polar coefficients (both parasitic and lift induced), and thrust settings (multiplier applied to the maximum engine thrust). The chosen flights are from three of the most common aircraft types (B738, B737, and A320) arriving at four airports (LAX, DEN, MSP, and DFW). Our ML approach is comprised of two different solutions: 1- using a subset of features that are known prior to the flight departure and do not change during flight (such as engine type, current temperature at departure & destination airports, aircraft type) and 2 - using a subset of temporal features of the flight trajectory (such as cruise altitude, Mach, airspeed, and rate of climb) in addition to the pre-departure features from the first solution. The labels or target variables are the APM parameters that were obtained by an optimized ordinary differential equations (ODE) fitting process (applied to individual flights). The ODE-fitting is very time intensive and is therefore performed offline. Thus, training an ML model to learn the relationship between the flight features and ODE-generated labels enables faster estimation of the APM parameters and is therefore amenable to real-time prediction. Various ML models including linear regression, random forest, XGBoost, and neural network were trained, and the results are compared. After model validation and hyperparameter-tuning, we observed that the Random Forest model outperformed the other three models by the overall mean square error (MSE) of 2% for the first solution and 1.5% for the second solution. Finally, the ML-derived parameters are compared against default BADA APM parameters using NASA’s Autonomy Development toolkit (ADK) simulation software. The simulation results for one of each aircraft type is shown and discussed.

Aida Sharif Rohani↗

Differential correction capability of the GTDS using TDRSS data

A differential correction (DC) capability was implemented in the Goddard Trajectory Determination System (GTDS) to process satellite tracking data acquired via the Tracking and Data Relay Satellite System (TRDRSS). Configuration of the TDRSS is reviewed, observation modeling is presented, and major features of the capability are discussed. The following types of TDRSS data can be processed by GTDS: two way relay range and Doppler measurements, hybrid relay range and Doppler measurements, one way relay Doppler measurements, and differenced one way relay Doppler measurements. These data may be combined with conventional ground based direct tracking data. By using Bayesian weighted least squares techniques, the software allows the simultaneous determination of the trajectories of up to four different satellites - one user satellite and three relay satellites. In addition to satellite trajectories, the following parameters can be optionally solved: for drag coefficient, reflectivity of a satellite for solar radiation pressure, transponder delay, station position, and biases.

Liu, S. Y.↗

Procedures for minimizing the effects of high solar activity on satellite tracking and ephemeris generation

We are currently experiencing a period of high solar radiation combined with wide short-term fluctuations in the radiation. The short-term fluctuations, especially when combined with highly energetic solar flares, can adversely affect the mission of U.S. Space Command's Space Surveillance Center (SSC) which catalogs and tracks the satellites in orbit around the Earth. Rapidly increasing levels of solar electromagnetic and/or particle radiation (solar wind) causes atmospheric warming, which, in turn, causes the upper-most portions of the atmosphere to expand outward, into the regime of low altitude satellites. The increased drag on satellites from this expansion can cause large, unmodeled, in-track displacements, thus undermining the SSC's ability to track and predict satellite position. On 13 March 1989, high solar radiation levels, combined with a high-energy solar flare, caused an exceptional amount of short-term atmospheric warming. The SSC temporarily lost track of over 1300 low altitude satellites--nearly half of the low altitude satellite population. Observational data on satellites that became lost during the days following the 13 March 'solar event' was analyzed and compared with the satellites' last element set prior to the event (referred to as a geomagnetic storm because of the large increase in magnetic flux in the upper atmosphere). The analysis led to a set of procedures for reducing the impact of future geomagnetic storms. These procedures adjust selected software limit parameters in the differential correction of element sets and in the observation association process and must be manually initiated at the onset of a geomagnetic storm. Sensor tasking procedures must be adjusted to ensure that a minimum of four observations per day are received for low altitude satellites. These procedures have been implemented and, thus far, appear to be successful in minimizing the effect of subsequent geomagnetic storms on satellite tracking and ephemeris computation.

Bredvik, Gordon D.↗

The Enigmatic Longevity of Granular Materials on Mars: The Case for Geologically Episodic Dune Formation

Martian sand dunes are concentrated in vast sand seas in the circumpolar belt of the planet's northern hemisphere, but they are also pervasive over the whole planet. Their occurrence is to be expected on a super-arid planetary surface subjected to boundary layer drag from a continually active atmosphere. Whilst their occurrence is to be expected, their survival is enigmatic. But the enigma only arises if the martian system is considered similar to Earth's --where sand is moved highly frequently, more or less on a seasonal basis. Experimentally it is readily demonstrated that active sand will soon wear down to small grains and eventually diminish to below the critical sand size required to sustain dune formation. According to conventional wisdom, sand moves at higher speeds on Mars than on Earth, and if it were to move as frequently as it does on Earth, then the dune-forming sand population should have long since disappeared, given the great longevity of the martian aeolian system (Sagan coined the term "kamikaze" grains to express this disappearance). No supply of sand could keep pace with this depletion, especially in light of the fact that Mars does not have very active weathering, nor significant crustal differentiation. On Earth, plate tectonics, magmatic activity, and general crustal differentiation over geological time have produced great concentrations of quartz crystals in the continental crustal masses. Not only are these quartz grains chemically and mechanically resilient, they are about the right size for being transported by either wind or water. Add to this, the geologically recent contribution of glacial grinding, and it is easy to see why there are dune field on Earth. So what are the martian dunes composed of, and how does the material survive the eons of attrition? In addition to experimental demonstrations of sand comminution in laboratory aeolian simulations, the problem can be approached from first principles. Sagan showed that by simple considerations of material strength versus mechanical work applied to the material, comminution to sub-sand size would be inevitable. Another semi-analytical approach might be taken by considering that the archetypal aeolian sand surface texture is an irregularly pitted ("frosted") surface composed of chipping hollows approximately 10 microns in diameter, 5 microns deep. Their volume = about 250 cubic microns, or about 1/25000 of the volume of a 100 micron diameter dune grain. Because a saltating grain always strikes another grain, then two surfaces are impacted. Thus each grain undergoes two impacts for every one saltation leap, when the impact statistics are considered for a closed dune system (it can be calculated that a grain can never undergo <1 impact, and never >2 per saltation leap). Hence, if we conservatively assume that there is damage to a grain each time it bounces, but with the minimum damage of only 2 microscopic craters per impact, then approximately 12,500 impacts are required to completely eliminate the grain. Of course, it would require only a fraction of this amount to reduce the grain to below sand size. A grain will make only several tens of saltation leaps on the stoss side of a dune before becoming buried on the lee slope. The dune then has to move its full length before the grain is exhumed again for abrasion. Even with this hiatus in transport, it is easy to see that terrestrial dunes need resupplying with sand in order to survive. In recent theoretical work it has been shown that martian aeolian transport may be initiated with high-speed grains, but this converts to a lower energy dynamic transport equilibrium in which a reptation population dominates grain transport (on Earth, at least half of the flux is by reptation and creep). On Mars, therefore, average grain speeds may be lower than those on Earth, or at least comparable. This would permit greater longevity for martian sands, but it would not go far enough to solve the survival problem. It may, however, explain why martian dunes are about the same size as terrestrial dunes. If martian saltation leaps were significantly longer than on Earth (as usually assumed), then a dune's lee slope would have to be correspondingly longer in order to trap the sand; this would scale up the whole dune structure. But with shorter trajectories in a reptation population, larger dunes would be unnecessary. Additional information is contained in the original.

Marshall, J.↗

Flow field Reconstruction for Inhomogeneous Turbulence using Data and Physics Driven Models

A methodology combining Large Eddy Simulation (LES) trained data and a physics driven wave packet model to obtain a reduced order reconstruction for broadband, three-dimensional, temporally stationary but spatially inhomogeneous, incompressible turbulence. Wake turbulence generated by an axisymmetric dragging disk with a turbulent co-flow serves as the benchmark test case. We begin by studying the proper-orthogonal decomposition of the turbulent fluctuations taken from a high-resolution LES to first identify whether the fields demonstrate a low-rank character. It is argued that the presence of the turbulent co-flow results in a largely broadband character lacking any tonal properties. This is especially true for Strouhal numbers greater than 1 and only a small fraction of energy is contained in the leading order Kelvin-Helmholtz modes. As such reconstructions and reduced order modeling purely relying on data from LES does not appear to be a lucrative solution - contrary to problems with strongly tonal character. To supplement the missing energy from a low order truncated mode expansion, we utilize a physics based super-resolution (enrichment) algorithm that relies on spatio-temporally localized Gabor wave packets whose time evolution is described using a set of ordinary differential equations. The reconstructed flow has single- and two-point correlations that are consistent with the reference high resolution simulation data.

SLS↗