Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained 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 523 records · Page 29

The Herschel DUNES Open Time Key Programme

We will use the unique photometric capabilities provided by Herschel to perform a deep and systematic survey for faint, cold debris disks around nearby stars. Our sensitivity-limited Open Time Key Programme (OTKP) aims at finding and characterizing faint extrasolar analogues to the Edgeworth-Kuiper Belt (EKB) in an unbiased, statistically significant sample of nearby FGK main-sequence stars. Our target set spans a broad range of stellar ages (from 0.1 to 10 Gyr) and is volume-limited (distances < 20 pc). All stars with known extrasolar planets within this distance are included; additionally, some M- and A-type stars will be observed in collaboration with the Herschel DEBRIS OTKP, so that the entire sample covers a decade in stellar mass, from 0.2 to 2 solar masses. We will perform PACS and SPIRE photometric observations covering the wavelength range from 70 to 500 microns. The PACS observations at 100 microns have been designed to detect the stellar photospheres down to the confusion limit with a signal-to-noise ratio > 5. The observations in the other Herschel bands will allow us to characterize, model, and constrain the disks. As a result, it will be possible for us to reach fractional dust luminosities of a few times 10-7, close to the EKB level in the Solar System. This will provide an unprecedented lower limit to the fractional abundance of planetesimal systems and allow us to assess the presence of giant planets, which would play dynamical roles similar to those played by Jupiter and Neptune in the Solar System. The proposed observations will provide new and unique evidence for the presence of mature planetary systems in the solar neighbourhood and, in turn, will address the universality of planet/planetary system formation in disks around young stars.

Danchi, William C.↗

Imaging with HST the Time Evolution of Eta Carinae's Colliding Winds

We report new HST/STIS observations that map the high-ionization forbidden line emission in the inner arcsecond of Eta Car, the first that fully image the extended wind-wind interaction region of the massive colliding wind binary. These observations were obtained after the 2009.0 periastron at orbital phases 0.084, 0.163, and 0.323 of the 5.54-year spectroscopic cycle. We analyze the variations in brightness and morphology of the emission, and find that blue-shifted emission (-400 to -200 km/s is symmetric and elongated along the northeast-southwest axis, while the red-shifted emission (+ 100 to +200 km/s) is asymmetric and extends to the north-northwest. Comparison to synthetic images generated from a 3-D dynamical model strengthens the 3-D orbital orientation found by Madura et al. (2011), with an inclination i approx. 138deg, argument of periapsis omega approx. 270deg, and an orbital axis that is aligned at the same PA on the sky as the symmetry axis of the Homunculus, 312deg. We discuss the potential that these and future mappings have for constraining the stellar parameters of the companion star and the long-term variability of the system.

Gull, Theodore R.↗

Scientific Rationale for Saturn's in Situ Exploration

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases׳ abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere may help constraining its bulk O/H ratio. We compare predictions of Jupiter and Saturn's bulk compositions from different formation scenarios, and highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations, all based on an entry probe that would descend through Saturn's stratosphere and troposphere under parachute down to a minimum of 10 bar of atmospheric pressure. We finally discuss the science payload required on a Saturn probe to match the measurement requirements.

Galileo↗

The Influence of surface binding energy on sputtering in models of the sodium exosphere of Mercury

We have simulated the sodium (Na) exosphere of Mercury to show how the exosphere is affected by the assumed surface binding energy (SBE) of Na in the sputtered component. We constrained ion precipitation onto the surface using distributions for the cusp regions that are consistent with measurements by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging Fast Imaging Plasma Spectrometer instrument. We have simulated sputtering with SBEs of 0.27, 2.6, 4.4, and 7.9 eV, with the lowest value commonly used in exosphere models and the highest from recent molecular dynamics calculations for the Na-bearing feldspar end-member, albite. A gradual change in the exosphere is seen as the yield decreases and the ejecta energy increases with increasing SBE. We describe the corresponding exosphere source functions for ion sputtering (IS), as well as for the previously studied processes of micrometeoroid impact vaporization and photon-stimulated desorption (PSD), along with their release energy distributions and spatial distributions. We have summed the contributions of the various source processes to explain how and when the different sources can be distinguished by observations. The modeled exosphere scale heights range from 72 km for PSD to over 1000 km for IS using a SBE of 7.9 eV. We find that the processes responsible for generating Mercury's Na exosphere are separable by measuring line-of-sight column densities tangent to the planet at various altitudes and positions around the planet. Our initial results are consistent with the Na being sputtered from a high-SBE material such as feldspar, which has been predicted to be abundant on the Mercury's surface.

Rosemary M. Killen↗

The dynamics of a flexible bladed disc on a flexible rotor in a two-rotor system

This paper describes the development of the analysis of the transient dynamic response of a bladed disk on a flexible rotor. The rotating flexible bladed disk is considered as a module in a complete turbine engine structure. The analysis of the flexible bladed disk (FBD) module is developed for the non-equilibrated one-diameter axial mode. The FBD motion is considered as a sum of two standing axial waves constrained to the rotor. The FBD is coupled inertially and gyroscopically to its rotor support, and indirectly through connecting elements, to the adjacent rotor and/or other supporting structures. Incorporated in the basic Turbine Engine Transient Response Analysis program (TETRA), the FBD module is demonstrated with a two-rotor model where the FBD can be excited into resonance by an unbalance in the adjacent rotor and at a frequency equal to the differential rotor speed. The FBD module also allows the analysis of two flexible bladed disks in the same rotor.

Gallardo, V. C.↗

Development of Head Injury Assessment Reference Values Based on NASA Injury Modeling

NASA is developing a new capsule-based, crewed vehicle that will land in the ocean, and the space agency desires to reduce the risk of injury from impact during these landings. Because landing impact occurs for each flight and the crew might need to perform egress tasks, current injury assessment reference values (IARV) were deemed insufficient. Because NASCAR occupant restraint systems are more effective than the systems used to determine the current IARVs and are similar to NASA s proposed restraint system, an analysis of NASCAR impacts was performed to develop new IARVs that may be more relevant to NASA s context of vehicle landing operations. Head IARVs associated with race car impacts were investigated by completing a detailed analysis of all of the 2002-2008 NASCAR impact data. Specific inclusion and exclusion criteria were used to select 4071 impacts from the 4015 recorder files provided (each file could contain multiple impact events). Of the 4071 accepted impacts, 274 were selected for numerical simulation using a custom NASCAR restraint system and Humanetics Hybrid-III 50th percentile numerical dummy model in LS-DYNA. Injury had occurred in 32 of the 274 selected impacts, and 27 of those injuries involved the head. A majority of the head injuries were mild concussions with or without brief loss of consciousness. The 242 non-injury impacts were randomly selected and representative of the range of crash dynamics present in the total set of 4071 impacts. Head dynamics data (head translational acceleration, translational change in velocity, rotational acceleration, rotational velocity, HIC-15, HIC-36, and the Head 3ms clip) were filtered according to SAE J211 specifications and then transformed to a log scale. The probability of head injury was estimated using a separate logistic regression analysis for each log-transformed predictor candidate. Using the log transformation constrains the estimated probability of injury to become negligible as IARVs approach zero. For the parameters head translational acceleration, head translational velocity change, head rotational acceleration, HIC-15, and HIC-36, conservative values (in the lower 95% confidence interval) that gave rise to a 5% risk of any injury occurring were estimated as 40.0 G, 7.9 m/s, 2200 rad/s2, 98.4, and 77.4 respectively. Because NASA is interested in the consequence of any particular injury on the ability of the crew to perform egress tasks, the head injuries that occurred in the NASCAR dataset were classified according to a NASA-developed scale (Classes I - III) for operationally relevant injuries, which classifies injuries on the basis of their operational significance. Additional analysis of the data was performed to determine the probability of each injury class occurring, and this was estimated using an ordered probit model. For head translational acceleration, head translational velocity change, head rotational acceleration, head rotational velocity, HIC-36, and head 3ms clip, conservative values of IARVs that produced a 5% risk of Class II injury were estimated as 50.7 G, 9.5 m/s, 2863 rad/s2, 11.0 rad/s, 30.3, and 46.4 G respectively. The results indicate that head IARVs developed from the NASCAR dataset may be useful to protect crews during landing impact.

Somers, Jeffrey T.↗

Global CO2 Distributions over Land from the Greenhouse Gases Observing Satellite (GOSAT)

January 2009 saw the successful launch of the first space-based mission specifically designed for measuring greenhouse gases, the Japanese Greenhouse gases Observing SATellite (GOSAT). We present global land maps (Level 3 data) of column-averaged CO2 concentrations (X(sub CO2)) derived using observations from the GOSAT ACOS retrieval algorithm, for July through December 2009. The applied geostatistical mapping approach makes it possible to generate maps at high spatial and temporal resolutions that include uncertainty measures and that are derived directly from the Level 2 observations, without invoking an atmospheric transport model or estimates of CO2 uptake and emissions. As such, they are particularly well suited for comparison studies. Results show that the Level 3 maps for July to December 2009 on a lO x 1.250 grid, at six-day resolution capture much of the synoptic scale and regional variability of X(sub CO2), in addition to its overall seasonality. The uncertainty estimates, which reflect local data coverage, X(sub CO2) variability, and retrieval errors, indicate that the Southern latitudes are relatively well-constrained, while the Sahara Desert and the high Northern latitudes are weakly-constrained. A probabilistic comparison to the PCTM/GEOS-5/CASA-GFED model reveals that the most statistically significant discrepancies occur in South America in July and August, and central Asia in September to December. While still preliminary, these results illustrate the usefulness of a high spatiotemporal resolution, data-driven Level 3 data product for direct interpretation and comparison of satellite observations of highly dynamic parameters such as atmospheric CO2.

Hammerling, Dorit M.↗

Comet C2012 S1 (ISON)s Carbon-rich and Micron-size-dominated Coma Dust

Comet C/2012 S1 (ISON) was unique in that it was a dynamically new comet derived from the Nearly Isotropic Oort cloud reservoir of comets with a sun-grazing orbit. We present thermal models for comet ISON (rh approx.1.15 AU, 2013-Oct-25 11:30 UT) that reveal comet ISON's dust was carbon-rich and dominated by a narrow size distribution dominated by approx. micron-sized grains. We constrained the models by our SOFIA FORCAST photometry at 11.1, 19.7 and 31.5 microns and by a silicate feature strength of approx.1.1 and an 8-13microns continuum greybody color temperature of approx. 275-280 K (using Tbb ∝ r−0.5 h and Tbb approx. 260-265 K from Subaru COMICS, 2013-Oct-19 UT)[1,2]. N-band spectra of comet ISON with the BASS instrument on the NASA IRTF (2013-Nov-11-12 UT) show a silicate feature strength of approx. 1.1 and an 11.2microns forsterite peak.[3] Our thermal models yield constraints the dust composition as well as grain size distribution parameters: slope, peak grain size, porosity. Specifically, ISON's dust has a low silicate-to- amorphous carbon ratio (approx. 1:9), and the coma size distribution has a steep slope (N4.5) such that the coma is dominated by micron-sized, moderately porous, carbon-rich dust grains. The N-band continuum color temperature implies submicronto micron-size grains and the steep fall off of the SOFIA far-IR photometry requires the size distribution to have fewer relative numbers of larger and cooler grains compared to smaller and hotter grains. A proxy for the dust production rate is fρ approx.1500 cm, akin to Afρ. ISON has a moderate-to-low dust-to-gas ratio. Comet ISON's dust grain size distribution does not appear similar to the few well-studied long-period Nearly Isotropic Comets (NICs), namely C/1995 O1 (Hale-Bopp) and C/2001 Q4 (NEAT) that had smaller and/or more highly porous grains and larger sizes, or C/2007 N4 (Lulin) and C/2006 P1 (McNaught) that had large and/or compact grains. Radial transport to comet-forming disk distances (≥ 20 AU) is easier for smaller grains (≤1 micron) than for larger grains (approx. 20 microns like Stardust terminal particles). The presence of predominantly micron-sized and smaller grains suggests comet ISON may have formed either earlier in disk evolution whereby larger grains did not have the time to be transported to distances beyond Neptune, or the comet formed so far out in the disk that larger grains did not traverse such large radial distances. The high carbon-content of ISON's refractory dust appears to be complimented by the presence of limitedlifetime organic (CHON-like) grain materials: preliminary analyses of near-IR and high-resolution optical spectra indicate that gas-phase daughter molecules C2, CN, and CH were more abundant than their parent molecules (C2H2, C2H6, measured in the near- IR). Dust composition as well as grain size distribution parameters (slope, peak grain size, and porosity) give clues to comet origins.

NASA/DLR Stratospheric Observatory↗

High-Fidelity Multidisciplinary Design Optimization Methodology with Application to Rotor Blades

A multidisciplinary design optimization procedure has been developed and applied to rotorcraft simulations involving tightly-coupled, high-fidelity computational fluid dynamics and comprehensive analysis. A discretely-consistent, adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows on unstructured, dynamic, overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system for high-fidelity rotorcraft analysis is verified; simulation results exhibit good agreement with established solutions. A constrained gradient-based design optimization for a HART-II rotorcraft configuration is demonstrated. The computational cost for individual components of the multidisciplinary sensitivity analysis is assessed and improved.

Wang, Li↗

The DYNAMO Orbiter Project: High Resolution Mapping of Gravity/Magnetic Fields and In Situ Investigation of Mars Atmospheric Escape

Dynamo is a small Mars orbiter planned to be launched in 2005 or 2007, in the frame of the NASA/CNES Mars exploration program. It is aimed at improving gravity and magnetic field resolution, in order to better understand the magnetic, geologic and thermal history of Mars, and at characterizing current atmospheric escape, which is still poorly constrained. These objectives are achieved by using a low periapsis orbit, similar to the one used by the Mars Global Surveyor spacecraft during its aerobraking phases. The proposed periapsis altitude for Dynamo of 120-130 km, coupled with the global distribution of periapses to be obtained during one Martian year of operation, through about 5000 low passes, will produce a magnetic/gravity field data set with approximately five times the spatial resolution of MGS. Low periapsis provides a unique opportunity to investigate the chemical and dynamical properties of the deep ionosphere, thermosphere, and the interaction between the atmosphere and the solar wind, therefore atmospheric escape, which may have played a crucial role in removing atmosphere, and water, from the planet. There is much room for debate on the importance of current atmosphere escape processes in the evolution of the Martian atmosphere, as early "exotic" processes including hydrodynamic escape and impact erosion are traditionally invoked to explain the apparent sparse inventory of present-day volatiles. Yet, the combination of low surface gravity and the absence of a substantial internally generated magnetic field have undeniable effects on what we observe today. In addition to the current losses in the forms of Jeans and photochemical escape of neutrals, there are solar wind interaction-related erosion mechanisms because the upper atmosphere is directly exposed to the solar wind. The solar wind related loss rates, while now comparable to those of a modest comet, nonetheless occur continuously, with the intriguing possibility of important cumulative and/or enhanced effects over the several billion years of the solar system's life. If the detailed history of the Martian internal field could be traced back, and the current escape processes could be understood well enough to model the expected stronger losses under early Sun conditions, one could go a long way toward constraining this part of the mysterious history of Mars' atmosphere.

Smrekar, S.↗

Robust computer-aided synthesis and optimization of linear multivariable control systems with varying plant dynamics via AUTOCON

AUTOCON is an automated computer-aided design tool for the synthesis and optimization of linear multivariable control systems based upon user-defined control parameter optimization. Violations in stability and performance requirements are computed from constraints on Single Input/Single Output (SISO) open- and closed-loop transfer function frequency responses, and from constraints on the singular-value frequency responses of Multiple Input/Multiple Output (MIMO) transfer functions, for all critical plant variations. Optimum nonlinear programming algorithms are used in the search for local constrained solutions in which violations in stability and performance are caused either to vanish or be minimized for a proper selection of the control parameters. Classical control system stability and performance design can, in this way, be combined with modern multivariable robustness methods to offer general frequency response loop-shaping via a computer-aided design tool. Complete Nichols, Nyquist, Bode, singular-value Bode magnitude and transient response plots are produced, including user-defined boundary responses. AUTOCON is used to synthesize and optimize the lateral/directional flight control system for a typical high-performance aircraft.

Lefkowitz, C. P.↗

Time-controlled descent guidance in uncertain winds

A procedure has been developed for constructing a statistical model of the altitude-dependent mean wind profile from the historical record of wind measurements at particular locations. The model is constructed by fitting a Markov process, with altitude as the stage variable, to the historical wind data. The wind model, together with the aircraft dynamics and the error characteristics of the navigation system, are incorporated in the design of a state estimator, which gives the minimum variance estimate of the aircraft state and the wind vector. The state and wind estimates are used as inputs to a linear feedback law for guiding the aircraft along the nominal trajectory. An example design of a time-constrained (4D RNAV) descent guidance system is presented, showing tracking accuracy, control activity, and probability of arrival time with and without the wind estimator.

Menga, G.↗

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

The Time Evolution of Eta Carinae's Colliding Winds

We report new HST/STIS observations that map the high-ionization forbidden line emission in the inner arc second of Eta Car, the first that fully image the extended wind-wind interaction region of the massive colliding wind binary. These observations were obtained after the 2009.0 periastron at orbital phases 0.084, 0.163, and 0.323 of the 5.54-year spectroscopic cycle. We analyze the variations in brightness and morphology of the emission, and find that blue-shifted emission (-400 to -200 km/s is symmetric and elongated along the northeast-southwest axis, while the red-shifted emission (+ 100 to +200 km/s) is asymmetric and extends to the north-northwest. Comparison to synthetic images generated from a 3-D dynamical model strengthens the 3-D orbital orientation found by Madura et al. (2011), with an inclination i = 138 deg, argument of periapsis w = 270 deg, and an orbital axis that is aligned at the same P A on the sky as the symmetry axis of the Homunculus, 312 deg. We discuss the potential that these and future mappings have for constraining the stellar parameters of the companion star and the long-term variability of the system. Plain-Language Abstract: With HST, we resolved the interacting winds of the binary, Eta Carinae. With a 3-D model, we find the binary orbit axis is aligned to the Homunculus axis. This suggests a connection between the binary and Homunculus ejection mechanism.

Gull, Theodore R.↗

Spatially resolved methane band photometry of Saturn. I - Absolute reflectivity and center-to-limb variations in the 6190-, 7250-, and 8900-A bands

Spatially resolved measurements of Saturn's absolute reflectivity in methane bands at 6190, 7250, and 8900 A and in nearby continuum regions are presented. Images were obtained through narrow-band interference filters with a 500 x 500-pixel charge-coupled device. Band/continuum ratios were measured to high accuracy by referencing to the ring brightness in each image. Several data processing techniques enhanced the quality of the observations. These are the use of the ring symmetry to find center position and orientation, accurate subtraction of ring light, and constrained image deconvolution. Uncertainty in the continuum absolute reflectivity is within 10%. Uncertainties in band/continuum ratios are from one to several percent. The Equatorial Zone was much brighter than any other latitude in the strong 8900 band image. Northern midlatitudes were brighter than southern midlatitudes. The latter observation indicates fewer high-altitude aerosols in the south, a possible result of atmospheric dynamics or seasonal sublimation of NH3 crystals. The data are tabulated and presented in a form suitable for quantitative scattering model analyses.

West, R. A.↗

A parametric LQ approach to multiobjective control system design

The synthesis of a constant parameter output feedback control law of constrained structure is set in a multiple-objective linear quadratic regulator (MOLQR) framework. The use of intuitive objective functions such as model-following ability and closed-loop trajectory sensitivity allows multiple-objective decision-making techniques, such as the surrogate worth tradeoff method, to be applied. For the continuous-time deterministic problem with an infinite time horizon, dynamic compensators as well as static output feedback controllers can be synthesized using a descent Anderson-Moore algorithm modified to impose linear equality constraints on the feedback gains by moving in feasible directions. Results of three different examples are presented, including a unique reformulation of the sensitivity reduction problem.

Kyr, Doug↗

Graphite/epoxy Composite Laminates with Co-cured Interlaminar Damping Layers

Damped composite laminates were fabricated by co-curing viscoelastic damping film with graphite/epoxy prepreg plies. The dynamic response of the damped plates was measured using an impulse response technique and compared with the response of similar undamped laminates. Modal damping was computed from the frequency response data. Micrographs of the damped laminates showed that the damping layers retained their integrity during the fabrication process. The layers significantly increased the damping in the composite laminates. The use of the constrained viscoelastic film as an integral part of composite structures appears to be a feasible approach to passive vibration control. Composite plates manufactured with co-cured damping layers may have commercial applications in cases where light weight, strength, and vibration and noise reduction are important considerations.

Pereira, J. Michael↗

NASP and SDI Spearhead CFD Developments

The National Aerospace Plane (NASP) program's purpose, as stated by the National Space Council, is to "develop and demonstrate hypersonic technologies with the ultimate goal of single stage to orbit." The council has also directed that "performance of the experimental flight vehicle will be constrained to the minimum necessary to meet the highest priority research, as opposed to operational objectives .... The program will be conducted in such a way as to minimize technical and cost uncertainty associated with the experimental vehicle." The purpose of the Strategic Defense Initiative (SDI), as defined by President Bush, is "...protection from limited ballistic missile strikes, whatever their source." Computational fluid dynamics (CFD) plays a vital role in both endeavors.

Mehta, Unmeel B.↗