Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fundamental constants”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Evolution of Structure in the Intergalactic Medium and the Nature of the LY-Alpha Forest

We have performed a detailed statistical study of the evolution of structure in a photoionized intergalactic medium (IGM) using analytical simulations to extend the calculation into the mildly nonlinear density regime found to prevail at z = 3. Our work is based on a simple fundamental conjecture: that the probability distribution function of the density of baryonic diffuse matter in the universe is described by a lognormal (LN) random field. The LN distribution has several attractive features and follows plausibly from the assumption of initial linear Gaussian density and velocity fluctuations at arbitrarily early times. Starting with a suitably normalized power spectrum of primordial fluctuations in a universe dominated by cold dark matter (CDM), we compute the behavior of the baryonic matter, which moves slowly toward minima in the dark matter potential on scales larger than the Jeans length. We have computed two models that succeed in matching observations. One is a nonstandard CDM model with OMEGA = 1, h = 0.5, and GAMMA = 0.3, and the other is a low-density flat model with a cosmological constant (LCDM), with OMEGA = 0.4, OMEGA(sub LAMBDA) = 0.6, and h = 0.65. In both models, the variance of the density distribution function grows with time, reaching unity at about z = 4, where the simulation yields spectra that closely resemble the Ly-alpha forest absorption seen in the spectra of high-z quasars. The calculations also successfully predict the observed properties of the Ly-alpha forest clouds and their evolution from z = 4 down to at least z = 2, assuming a constant intensity for the metagalactic UV background over this redshift range. However, in our model the forest is not due to discrete clouds, but rather to fluctuations in a continuous intergalactic medium. At z = 3; typical clouds with measured neutral hydrogen column densities N(sub H I) = 10(exp 13.3), 10(exp 13.5), and 10(exp 11.5) /sq cm correspond to fluctuations with mean total densities approximately 10, 1, and 0.1 times the universal mean baryon density. Perhaps surprisingly, fluctuations whose amplitudes are less than or equal to the mean density still appear as "clouds" because in our model more than 70% of the volume of the IGM at z = 3 is filled with gas at densities below the mean value.

Bi, Hongguang↗
NASA NTRSDate not supplied

Cal X-1, a Pioneer mission for absolute in-orbit flux calibration for current and future X-ray observatories

X-ray astronomy missions are now making observations capable of answering fundamental physics and cosmology questions. However, those studies are often limited by lack of absolute X-ray flux calibration; the most powerful X-ray observatories reveal discrepancies in the source fluxes at the 10-20% level. These discrepancies are plausibly related to changes in the instrument performance between the ground calibration and the orbit. We propose to overcome this problem by establishing "standard candles" in the X-ray sky. Our Cal X-1 mission, a Pioneer concept, will consist of 2 small satellites flying in tandem. One will have an X-ray telescope and the other, a pair of absolutely-calibrated X-ray sources (E=1.5 keV and 6 keV). The telescope will ping-pong between the source satellite and several well-chosen constant celestial sources, determining the accurate fluxes of those sources and turning them into standard candles for all other X-ray observatories, current and future.

Maxim Markevitch↗
DOE OSTI2024

Jets from shocked metal surfaces with grooves: Missing experiments

Many studies have investigated the mass outflows generated when a planar shock transits an imperfect (“defected”) metal surface, where the defects are symmetric triangular or sinusoidal grooves. Yet a fundamental question remains unanswered: how does the quantity of outflow mass and its maximum velocity vary as a function of the groove cross-sectional aspect ratio? We identify two sets of missing experiments that must be addressed to answer the question. The aspect ratio (groove depth over width) is equivalently represented by θ, the cross-sectional half angle, or by η 0 k, the amplitude multiplied by an effective wavenumber. Low θ (high η 0 k) grooves comprise the first set of missing experiments, which are necessary to determine the validity of theoretical predictions of the nonlinear regime (η 0 k ≥ 1, θ < 57.5°). The second set of missing experiments are those in which the volume of the groove (or equivalently, the axial cross-sectional area) has been held constant as θ or η 0 k are varied. Such experiments are necessary to independently measure the effects of variations in groove volume and groove aspect ratio on the resulting jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Randomized Adiabatic Quantum Linear Solver Algorithm with Optimal Complexity Scaling and Detailed Running Costs

Solving linear systems of equations is a fundamental problem with a wide variety of applications across many fields of science, and there is increasing effort to develop quantum linear solver algorithms. Subaşı et al. [Phys. Rev. Lett. 122, 060504 (2019)] proposed a randomized algorithm inspired by adiabatic quantum computing, based on a sequence of random Hamiltonian simulation steps, with suboptimal scaling in the condition number 𝜅 of the linear system and the target error 𝜖. Here we go beyond these results in several ways. Firstly, using filtering [Lin and Tong, Quantum 4, 361 (2020)] and Poissonization techniques [Cunningham and Roland, ArXiv:2406.03972 (2024)], the algorithm complexity is improved to the optimal scaling 𝑂⁡(𝜅⁢log (1/𝜖))—an exponential improvement in 𝜖, and a shaving of a log 𝜅 scaling factor in 𝜅. Secondly, the algorithm is further modified to achieve constant factor improvements, which are vital as we progress towards hardware implementations on fault-tolerant devices. We introduce a cheaper randomized walk operator method replacing Hamiltonian simulation—which also removes the need for potentially challenging classical precomputations; randomized routines are sampled over optimized random variables; circuit constructions are improved. We obtain a closed formula rigorously upper bounding the expected number of times one needs to apply a block-encoding of the linear system matrix to output a quantum state encoding the solution to the linear system. The upper bound is 837⁢𝜅 at 𝜖 = 10 −10 for Hermitian matrices.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗
DOE OSTI2023

High-resolution vibrational predissociation spectroscopy of I - ·H 2 O by single-mode CW infrared excitation in a 3D cryogenic ion trap

We describe the integration of a tunable, single-mode, continuous wave infrared laser into a cryogenic ion spectroscopy experiment to measure the rovibrational spectrum of the I - ∙H 2 O complex in the OH stretching region. These upper levels lie about 300 cm -1 above the dissociation threshold. The measurements are carried out by loading the ions in a radiofrequency ion trap at 10 Hz and cooling them to 5 K with pulsed He buffer gas. IR photodissociation (PD) of the I - ∙H 2 O complex is monitored by recording the I - product yield by time-of-flight mass spectrometry as a function of laser wavelength. Very narrow (Δν~75 MHz) rotational lines are observed throughout the spectrum, indicating long (ca. 2 ns) lifetimes for the excited metastable rovibrational levels. Rotational analysis of the band arising from the K”=1 to K’=2 transition of the free OH stretching fundamental yields the structure of the complex for the first time. Over 50% of the trapped ion ensemble in the trap can be photodissociated upon excitation of a single rotational line. This enables very high signal-to-noise in the PD spectrum, and is traced to a mechanism in which the ground state rotational levels are rapidly equilibrated by collisions with the buffer gas.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Metal Decorated Multi-Walled Carbon Nanotube/Polyimide Composites with High Dielectric Constants and Low Loss Factors

The measurement of observable electromagnetic phenomena in materials and their derived intrinsic electric material properties are of prime importance in the discovery and development of material systems for electronic and aerospace applications. Nanocomposite materials comprised of metal decorated multi-walled carbon nanotubes (MWCNTs) were prepared by a facile method and characterized. Metal particles such as silver, platinum and palladium with diameters ranging from less than 5 to over 50 nanometers were distributed randomly on the MWCNTs. The metal-containing MWCNTs were then melt mixed into a polymer matrix and the mixture extruded as ribbons. These extruded ribbons exhibited a moderate to high degree of MWCNT alignment as determined by HRSEM. These ribbons were then fabricated into test specimens while maintaining MWCNT alignment and subsequently characterized for electromagnetic properties at 8-12 GHz. The present study is focused on silver decorated MWCNTs dispersed in an Ultem polyimide matrix. The results of the electromagnetic characterization showed that certain sample configurations exhibited a decoupling of the permittivity and loss factor (?? and ??) indicating that these properties could be tailored within certain limits. The decoupling and independent control of these fundamental electrical material parameters offer a new class of materials with potential applications in electronics, microwave engineering and optics.

Ghose, Sayata↗
NASA NTRS2011

Full-Scale Turbofan-Engine Turbine-Transfer Function Determination Using Three Internal Sensors

Existing NASA/Honeywell EVNERT full-scale static engine test data is analyzed by using source-separation techniques in order to determine the turbine transfer of the currently sub-dominant combustor noise. The results are used to assess the combustor-noise prediction capability of the Aircraft Noise Prediction Program (ANOPP). Time-series data from three sensors internal to the Honeywell TECH977 research engine is used in the analysis. The true combustor-noise turbine-transfer function is educed by utilizing a new three-signal approach. The resulting narrowband gain factors are compared with the corresponding constant values obtained from two empirical acoustic-turbine-loss formulas. It is found that a simplified Pratt & Whitney formula agrees better with the experimental results for frequencies of practical importance. The 130 deg downstream-direction far-field 1/3-octave sound-pressure levels (SPL) results of Hultgren & Miles are reexamined using a post-correction of their ANOPP predictions for both the total noise signature and the combustion-noise component. It is found that replacing the standard ANOPP turbine-attenuation function for combustion noise with the simplified Pratt & Whitney formula clearly improves the predictions. It is recommended that the GECOR combustion-noise module in ANOPP be updated to allow for a user-selectable switch between the current transmission-loss model and the simplified Pratt & Whitney formula. The NASA Fundamental Aeronautics Program has the principal objective of overcoming today's national challenges in air transportation. The Subsonic Fixed Wing Project's Reduce-Perceived-Noise Technical Challenge aims to develop concepts and technologies to dramatically reduce the perceived aircraft noise outside of airport boundaries. The reduction of aircraft noise is critical to enabling the anticipated large increase in future air traffic.

Hultgren, Lennart S.↗
NASA NTRS2014

Visualizing Flutter Mechanism as Traveling Wave Through Animation of Simulation Results for the Semi-Span Super-Sonic Transport Wind-Tunnel Model

It has long been recognized that frequency and phasing of structural modes in the presence of airflow play a fundamental role in the occurrence of flutter. Animation of simulation results for the long, slender Semi-Span Super-Sonic Transport (S4T) wind-tunnel model demonstrates that, for the case of mass-ballasted nacelles, the flutter mode can be described as a traveling wave propagating downstream. Such a characterization provides certain insights, such as (1) describing the means by which energy is transferred from the airflow to the structure, (2) identifying airspeed as an upper limit for speed of wave propagation, (3) providing an interpretation for a companion mode that coalesces in frequency with the flutter mode but becomes very well damped, (4) providing an explanation for bursts of response to uniform turbulence, and (5) providing an explanation for loss of low frequency (lead) phase margin with increases in dynamic pressure (at constant Mach number) for feedback systems that use sensors located upstream from active control surfaces. Results from simulation animation, simplified modeling, and wind-tunnel testing are presented for comparison. The simulation animation was generated using double time-integration in Simulink of vertical accelerometer signals distributed over wing and fuselage, along with time histories for actuated control surfaces. Crossing points for a zero-elevation reference plane were tracked along a network of lines connecting the accelerometer locations. Accelerometer signals were used in preference to modal displacement state variables in anticipation that the technique could be used to animate motion of the actual wind-tunnel model using data acquired during testing. Double integration of wind-tunnel accelerometer signals introduced severe drift even with removal of both position and rate biases such that the technique does not currently work. Using wind-tunnel data to drive a Kalman filter based upon fitting coefficients to analytical mode shapes might provide a better means to animate the wind tunnel data.

Christhilf, David M.↗
NASA NTRSDate not supplied

On The Incorporation of Both Function-Driven Requirements and Topology Optimization in The Development of Lunar Launch & Landing Pads

In the design of additively manufactured structures, there can be a constant battle between focusing on function-driven design and topology optimization in an attempt to improve a products ability to meet all requirements. Function-driven optimization is defined as the process of identifying the best design parameters that satisfy project functionality requirements. Function-driven optimization is typically implemented by a designer(s) who, instead of computational limitations, has human bias for the input or inputs determined. On the other hand, topology optimization is defined as a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions, and constraints with the goal of maximizing the performance of the system. In topology optimization, the work is defined by the limited mathematical inputs and computationally available goals/restraints. If not evaluated simultaneously, these two approaches can result in fundamentally different design concepts meant to address the same requirements. NASA’s Artemis program includes the development of a base camp near the lunar South Pole. After site selection is made, one of the first infrastructure elements needed will be a launch and landing pad (LLP). The construction of such a pad will prevent damage to nearby structures from plume-surface interaction (PSI) while also providing a stable platform for the lander. These benefits reduce the required distance between the established landing zones and other settlement infrastructure. In this paper, the authors evaluate the requirements associated with lunar LLP development, application of both functional and topological optimization techniques, methodologies to combine the two techniques, and potential impacts on final design of such a pad.

In-situ↗
NASA NTRS2009

Metal Decorated Multi-Walled Carbon Nanotube/Polyimide Composites with High Dielectric Constants and Low Loss Factors

The measurement of observable electromagnetic phenomena in materials and their derived intrinsic electrical material properties are of prime importance in the discovery and development of material systems for electronic and aerospace applications. Nanocomposite materials comprised of metal decorated multi-walled carbon nanotubes (MWCNTs) were prepared by a facile method and characterized. Metal particles such as silver(Ag), platinum(Pt) and palladium(Pd) with diameters ranging from less than 5 to over 50 nanometers were distributed randomly on the MWCNTs. The present study is focused on silver decorated MWCNTs dispersed in a polyimide matrix. The Ag-containing MWCNTs were melt mixed into Ultem(TradeMark) and the mixture extruded as ribbons. The extruded ribbons exhibited a moderate to high degree of MWCNT alignment as determined by HRSEM. These ribbons were then fabricated into test specimens while maintaining MWCNT alignment and subsequently characterized for electrical and electromagnetic properties at 8-12 GHz. The results of the electromagnetic characterization showed that certain sample configurations exhibited a decoupling of the permittivity (epsilon ) and loss factor (epsilon") indicating that these properties could be tailored within certain limits. The decoupling and independent control of these fundamental electrical material parameters offers a new class of materials with potential applications in electronics, microwave engineering and optics.

Elliott, Holly A.↗
DOE OSTI2020

The acidity of atmospheric particles and clouds

Acidity, defined as pH, is a central component of aqueous chemistry. In the atmosphere, the acidity of condensed phases (aerosol particles, cloud water, and fog droplets) governs the phase partitioning of semivolatile gases such as HNO 3 , NH 3 , HCl, and organic acids and bases as well as chemical reaction rates. It has implications for the atmospheric lifetime of pollutants, deposition, and human health. Despite its fundamental role in atmospheric processes, only recently has this field seen a growth in the number of studies on particle acidity. Even with this growth, many fine-particle pH estimates must be based on thermodynamic model calculations since no operational techniques exist for direct measurements. Current information indicates acidic fine particles are ubiquitous, but observationally constrained pH estimates are limited in spatial and temporal coverage. Clouds and fogs are also generally acidic, but to a lesser degree than particles, and have a range of pH that is quite sensitive to anthropogenic emissions of sulfur and nitrogen oxides, as well as ambient ammonia. Historical measurements indicate that cloud and fog droplet pH has changed in recent decades in response to controls on anthropogenic emissions, while the limited trend data for aerosol particles indicate acidity may be relatively constant due to the semivolatile nature of the key acids and bases and buffering in particles. This paper reviews and synthesizes the current state of knowledge on the acidity of atmospheric condensed phases, specifically particles and cloud droplets. It includes recommendations for estimating acidity and pH, standard nomenclature, a synthesis of current pH estimates based on observations, and new model calculations on the local and global scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Relativistic Hartree–Fock chiral Lagrangians with confinement, nucleon finite size and short-range effects

A relativistic Hartree–Fock Lagrangian including a chiral potential and nucleon polarisation is investigated in hopes of providing a better description of dense nuclear matter. Here, we fully consider the contribution of the exchange Fock term to the energy and the self-energies, and in addition we investigate the nucleon’s compositeness and finite size effects (confinement and form factors) and short range correlations modeled by a Jastrow ansatz. These effects are added step by step, such that their impact on the dense matter properties can be analysed in details. The parameters of the model are adjusted to reproduce fundamental properties related to the QCD theory at low energy, such as the chiral symmetry breaking, nucleon’s quark substructure and Lattice-QCD predictions, as well as two empirical properties at saturation: the binding energy and the density. All other empirical parameters, e.g., symmetry energy and its slope, incompressibility modulus, effective mass, as well as spin–isospin Landau–Midgal parameter are predictions of the models and can be used to evaluate the gain of the different approximation schemes in describing nuclear properties. Bayesian statistics is employed in order to propagate parameter uncertainties into predictions for the nuclear matter properties. We show that the splitting of the effective Landau mass is largely influenced by the value of the ρ T coupling, and we show that the fit to the symmetry energy, which induces an increase of the coupling constant g ρ by about 20–25% compared to the case where it is fixed by the quark model, provides a very good EoS compatible with the present nuclear physics knowledge.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2000

Kinetics of constant gravitropic stimulus responses in Arabidopsis roots using a feedback system

The study of gravitropism is hindered by the fact that as a root responds, the gravitational stimulus changes. Using a feedback system to connect a rotating stage platform to a video digitizer system, we were able to maintain a constant angle of gravistimulation to Arabidopsis roots for long time periods. The rate of curvature approximated the sine rule for angles of stimulation between 20 degrees and 120 degrees. For a given angle of stimulation, the rate of curvature also remained constant, with no observed diminishment of the response. Although previous reports of Arabidopsis root gravitropism suggest latent periods of approximately 30 min, using a smooth mechanical stage to reorient the root, we observed a mean time lag of approximately 10 min. This more rapid onset of curvature can, in part, be explained by reduced mechanical perturbation during the process of gravistimulation. This suggests that mechanical stimulation associated with rapid root re-orientation may confound investigations of early gravitropic events.

NASA Program Fundamental Space Biology↗
DOE OSTI2024

Efficient Anomaly Detection Driven By Different Machine Learning Architectures And Models

The rapid growth and ubiquitous adoption of the internet and cyber-physical systems (CPS) have fundamentally transformed modern communication, work, and human-system interactions. While networks now form the backbone of critical digital ecosystems, enabling seamless data transmission across diverse, interconnected systems, this increased connectivity also expands the attack surface, making real-time detection of network intrusions and anomalies a pressing challenge. Detecting unusual activities within network infrastructure requires advanced data traffic analysis to differentiate between legitimate and malicious interactions. Traditional approaches to network anomaly detectionâ??such as rule-based and signature-based systemsâ??often depend on predefined patterns to identify known anomalies, limiting their effectiveness against emerging, stealthy, or previously unseen threats. These conventional methods suffer from high false alarm rates and fail to adapt to the ever-evolving nature of network traffic, particularly in large-scale, decentralized environments where data volume, velocity, and variety are constantly increasing. This dissertation presents artificial intelligence (AI)-driven approaches to anomaly detection that leverage graphics processing unit (GPU)-enabled high-performance computing (HPC) platforms for processing massive network traffic data and monitoring the components of cyber-physical systems (CPS) for potentially hazardous conditions. The research advances several key contributions: (1) Designing efficient machine learning techniques for CPS condition monitoring and anomaly detection; (2) enabling federated learning (FL) frameworks that enable distributed detection while preserving data privacy and system resilience; (3) exploring graph-based methodologies combining graph neural networks (GNN) and graph machine learning (ML) approaches for the Internet of Things (IoT) and automotive network security, and (4) performing distributed edge computing optimizations that integrate FL with scalable technologies for reduced communication overhead. Through extensive experiments, these methodologies demonstrate that complex anomaly detection and condition monitoring tasks can be achieved while balancing computational efficiency and detection accuracy through fine-grained network information processing. The frameworks developed in this research establish a robust foundation for network anomaly detection, providing scalable, adaptive, and privacy-preserving solutions for safeguarding CPS and IoT networks in an increasingly interconnected digital landscape. The practical implications of these research findings are significant, as they can inform the development of next-generation network security systems and contribute to the protection of critical infrastructure against sophisticated cyber attacks.

Marfo, William↗
NASA NTRS2000

Dynamical Reference Frame: Current Relevance and Future Prospects

Planetary and lunar ephemerides are no longer used for the determination of inertial space. Instead, the new fundamental reference frame, the International Celestial Reference Frame (ICRF), is inherently less susceptible to extraneous, non-inertial rotations than a dynamical reference frame determined by the ephemerides would be. Consequently, the ephemerides are now adjusted onto the ICRF, and they are fit to two modern, accurate observational data types: ranging (radar, lunar laser, spacecraft) and Very Long Baseline Interferometry (VLBI) (of spacecraft near planets). The uncertainties remaining in the inner planet ephemerides are on the order of 1 kilometer, both in relative positions between the bodies and in the orientation of the inner system as a whole. The predictive capabilities of the inner planet ephemerides are limited by the uncertainties in the masses of many asteroids. For this reason, future improvements to the ephemerides must await determinations of many asteroid masses. Until then, it will be necessary to constantly update the ephemerides with a continuous supply of observational data.

VLBI (VERY LONG BASELINE INTERFEROMETRY)↗
NASA NTRS2011

Long Term Missions at the Sun-Earth Libration Point L1: ACE, SOHO, and WIND

Three heliophysics missions -- the Advanced Composition Explorer (ACE), Solar Heliospheric Observatory (SOHO), and the Global Geoscience WIND -- have been orbiting the Sun-Earth interior libration point L1 continuously since 1997, 1996, and 2004, respectively. ACE and WIND (both NASA missions) and SOHO (an ESA-NASA joint mission) are all operated from the NASA Goddard Space Flight Center (GSFC). While ACE and SOHO have been dedicated libration point orbiters since their launches, WIND has had also a remarkable 10-year career flying a deep-space, multiple lunar-flyby trajectory prior to 2004. That era featured 36 targeted lunar flybys with excursions to both L1 and L2 before its final insertion in L1 orbit. A figure depicts the orbits of the three spacecraft, showing projections of the orbits onto the orthographic planes of a solar rotating ecliptic frame of reference. The SOHO orbit is a quasi-periodic halo orbit, where the frequencies of the in-plane and out-of-plane motions are practically equal. Such an orbit is seen to repeat itself with a period of approximately 178 days. For ACE and WIND, the frequencies of the in-plane and out-of-plane motions are unequal, giving rise to the characteristic Lissajous motion. ACE's orbit is of moderately small amplitude, whereas WIND's orbit is a large-amplitude Lissajous of dimensions close to those of the SOHO halo orbit. As motion about the collinear points is inherently unstable, stationkeeping maneuvers are necessary to prevent orbital decay and eventual escape from the L1 region. Though the three spacecraft are dissimilar (SOHO is a 3-axis stabilized Sun pointer, WIND is a spin-stabilized ecliptic pole pointer, and ACE is also spin-stabilized with its spin axis maintained between 4 and 20 degrees of the Sun), the stationkeeping technique for the three is fundamentally the same. The technique consists of correcting the energy of the orbit via a delta-V directed parallel or anti-parallel to the Spacecraft-to-Sun line. SOHO achieves this using thrusters oriented in line with the solar direction. WIND achieves the delta-V via pulsing radial thrusters when aligned with the Sun. ACE uses axial thrusters to apply delta-V with a component that is 94% or more aligned with the ACE-Sun line. Sunward thrust adds energy to the orbit preventing decay back toward Earth. Thrust directed anti-Sunward takes energy out of the L1 orbit, thereby preventing escape from the Earth-Moon system into independent heliocentric orbit. Libration point orbit stationkeeping delta-V costs grow exponentially with time elapsed from the last maneuver performed. The doubling time constant is approximately 16 days. For the sake of fuel conservation, and for limiting the absolute magnitude of propulsion performance errors, stationkeeping maneuvers should be performed before the delta-V grows too large; for our purposes 'too large' is considered to be greater than 0.5 m/sec. In practice, the typical interval between burns for this trio is about three months, and the typical delta-V is much smaller than 0.5 m/sec. Typical annual stationkeeping costs have been around 1.0 m/sec for ACE and WIND, and much less than that for SOHO. All three spacecraft have ample fuel remaining; barring contingencies all three could, in principle, be maintained at L1 for decades to come. This paper will review the L1 orbits and the mission history of ACE, WIND, and SOHO, and describe the stationkeeping techniques and orbit maneuver experience. The Lissajous phase control that was practiced for ACE during the period from 1999 to 2001 will also be briefly discussed. The final section will consider the future of these ongoing missions.

Roberts, Craig E.↗