Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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 271 records · Page 15

Methods for Detection and Correction of Sudden Pixel Sensitivity Drops

PDC 8.0 includes implementation of a new algorithm to detect and correct step discontinuities appearing in roughly one of every twenty stellar light curves during a given quarter. An example of such a discontinuity in an actual light curve is shown in fig. 1. The majority of such discontinuities are believed to result from high-energy particles (either cosmic or solar in origin) striking the photometer and causing permanent local changes (typically -0.5% in summed apertures) in quantum efficiency, though a partial exponential recovery is often observed. Since these features, dubbed sudden pixel sensitivity dropouts (SPSDs), are uncorrelated across targets they cannot be properly accounted for by the current detrending algorithm. PDC de-trending is based on the assumption that features in flux time series are due either to intrinsic stellar phenomena or to systematic errors and that systematics will exhibit measurable correlations across targets. SPSD events violate these assumptions and their successful removal not only rectifies the flux values of affected targets, but demonstrably improves the overall performance of PDC de-trending.

Kolodziejczak, Jeffery↗

Satellite-Based Assessment of Grassland Conversion and Related Fire Disturbance in the Kenai Peninsula, Alaska

Spruce beetle-induced (Dendroctonus rufipennis (Kirby)) mortality on the Kenai Peninsula has been hypothesized by local ecologists to result in the conversion of forest to grassland and subsequent increased fire danger. This hypothesis stands in contrast to empirical studies in the continental US which suggested that beetle mortality has only a negligible effect on fire danger. In response, we conducted a study using Landsat data and modeling techniques to map land cover change in the Kenai Peninsula and to integrate change maps with other geospatial data to predictively map fire danger for the same region. We collected Landsat imagery to map land cover change at roughly five-year intervals following a severe, mid-1990s beetle infestation to the present. Landcover classification was performed at each time step and used to quantify grassland encroachment patterns over time. The maps of land cover change along with digital elevation models (DEMs),temperature, and historical fire data were used to map and assess wildfire danger across the study area. Results indicate the highest wildfire danger tended to occur in herbaceous and black spruce land cover types, suggesting that the relationship between spruce beetle damage and wildfire danger in costal Alaskan forested ecosystems differs from the relationship between the two in the forests of the coterminous United States. These change detection analyses and fire danger predictions provide the Kenai National Wildlife Refuge (KENWR) ecologists and other forest managers a better understanding of the extent and magnitude of grassland conversion and subsequent change in fire danger following the 1990s spruce beetle outbreak.

wildfire↗

Physiological Maturation of Regenerating Hair Cells

The bullfrog saccule, a sensor of gravity and substrate-borne vibration, is a model system for hair cell transduction. Saccular hair cells also increase in number throughout adult life and rapidly recover after hair cell damage, making this organ an ideal system for studying hair cell development, repair, and regeneration. We have used of hair cell and supporting cell immunocytochemical markers to identify damaged hair cells and hair cell precursors in organotypic cultures of the bullfrog saccule. We then used an innovative combination of confocal, electron, and time-lapse microscopy to study the fate of damaged hair cells and the origin of new hair cells after gentamicin ototoxicity in normal and mitotically blocked saccular cultures. These studies have shown that gentamicin ototoxicity produces both lethal and sublethal hair cell damage. They have also shown that hair cell recovery in this organ takes place by both the repair of sublethally damaged hair cells and by the replacement of lost hair cells by mitotic regeneration. In parallel studies, we have used biophysical and molecular biological techniques to study the differentiation and innervation of developing, repairing, and regenerating hair cells. More specifically, we have used RT-PCR to obtain the bullfrog homologues of L-type voltage- gated calcium (L-VGCC) and large-conductance Ca(2+)-activated potassium (BK) channel genes. We have then obtained probes for these genes and, using in situ hybridization, begun to examine their expression in the bullfrog saccule and amphibian papilla. We have also used fluorescent-labeled channel toxins and channel toxin derivatives to determine the time of appearance of L-type voltage-gated calcium (L-VGCC) and Ca(2+)-activated potassium (BK) channels and to study dynamic changes in the number, distribution, and co-localization of these proteins in developing, repairing, and regenerating hair cells. Using time-lapse microscopy, we are also studying the dynamic relationship between ion channel clustering and synaptic formation in hair cells and afferent neurons. In future studies, we will determine when hair cell precursors acquire electrical tuning, and, using whole-cell patch-clamp techniques, identify and characterize their L-VGCC and BK currents. We will also use biophysical techniques to determine the number of L-VGCC and BK channels and the size and gating kinetics of their underlying L-VGCC and BK conductances, correlating these variables with the amplitude and frequency of membrane oscillations produced by intracellular current steps. We expect these studies to determine how hair cells regulate ion channel expression to achieve specific physiological responses.

Baird, Richard A.↗

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

XTCE (XML Telemetric and Command Exchange) Standard Making It Work at NASA. Can It Work For You?

The XML Telemetric and Command Exchange (XTCE) standard is intended as a way to describe telemetry and command databases to be exchanged across centers and space agencies. XTCE usage has the potential to lead to consolidation of the Mission Operations Center (MOC) Monitor and Control displays for mission cross-support, reducing equipment and configuration costs, as well as a decrease in the turnaround time for telemetry and command modifications during all the mission phases. The adoption of XTCE will reduce software maintenance costs by reducing the variation between our existing mission dictionaries. The main objective of this poster is to show how powerful XTCE is in terms of interoperability across centers and missions. We will provide results for a use case where two centers can use their local tools to process and display the same mission telemetry in their MOC independently of one another. In our use case we have first quantified the ability for XTCE to capture the telemetry definitions of the mission by use of our suite of support tools (Conversion, Validation, and Compliance measurement). The next step was to show processing and monitoring of the same telemetry in two mission centers. Once the database was converted to XTCE using our tool, the XTCE file became our primary database and was shared among the various tool chains through their XTCE importers and ultimately configured to ingest the telemetry stream and display or capture the telemetered information in similar ways.Summary results include the ability to take a real mission database and real mission telemetry and display them on various tools from two centers, as well as using commercially free COTS.

CCSDS↗

Numerical Computation of a Continuous-thrust State Transition Matrix Incorporating Accurate Hardware and Ephemeris Models

A significant body of work exists showing that providing a nonlinear programming (NLP) solver with expressions for the problem constraint gradient substantially increases the speed of program execution and can also improve the robustness of convergence, especially for local optimizers. Calculation of these derivatives is often accomplished through the computation of spacecraft's state transition matrix (STM). If the two-body gravitational model is employed as is often done in the context of preliminary design, closed form expressions for these derivatives may be provided. If a high fidelity dynamics model, that might include perturbing forces such as the gravitational effect from multiple third bodies and solar radiation pressure is used then these STM's must be computed numerically. We present a method for the power hardward model and a full ephemeris model. An adaptive-step embedded eight order Dormand-Prince numerical integrator is discussed and a method for the computation of the time of flight derivatives in this framework is presented. The use of these numerically calculated derivatieves offer a substantial improvement over finite differencing in the context of a global optimizer. Specifically the inclusion of these STM's into the low thrust missiondesign tool chain in use at NASA Goddard Spaceflight Center allows for an increased preliminary mission design cadence.

Low Thrust↗

CFD Assessment of Forward Booster Separation Motor Ignition Overpressure on ET XT 718 Ice/Frost Ramp

Computational fluid dynamics assessment of the forward booster separation motor ignition over-pressure was performed on the space shuttle external tank X(sub T) 718 ice/frost ramp using the flow solver OVERFLOW. The main objective of this study was the investigation of the over-pressure during solid rocket booster separation and its affect on the local pressure and air-load environments. Delta pressure and plume impingement were investigated as a possible contributing factor to the cause of the debris loss on shuttle missions STS-125 and STS-127. A simplified computational model of the Space Shuttle Launch Vehicle was developed consisting of just the external tank and the solid rocket boosters with separation motor nozzles and plumes. The simplified model was validated by comparison to full fidelity computational model of the Space Shuttle without the separation motors. Quasi steady-state plume solutions were used to calibrate the thrust of the separation motors. Time-accurate simulations of the firing of the booster-separation motors were performed. Parametric studies of the time-step size and the number of sub-iterations were used to find the best converged solution. The computed solutions were compared to previous OVERFLOW steady-state runs of the separation motors with reaction control system jets and to ground test data. The results indicated that delta pressure from the overpressure was small and within design limits, and thus was unlikely to have contributed to the foam losses.

Tejnil, Edward↗

Accurate Assessment of Land-Atmosphere Coupling in Climate Models Requires High Frequency Data Output

Land-atmosphere (L-A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centred climate anomalies. Local L-A coupling (LoCo) metrics capture relevant L-A processes, highlighting the impact of soil and vegetation states on surface flux partitioning, and the impact of surface fluxes on boundary layer (BL) growth, development, and entrainment of air above the BL. A primary goal of the Climate Process Team on Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L-A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land-atmosphere interactions span time scales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability of behavioural regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyse, and better understand L-A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behaviour. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant time scales, (ii) enable demonstration of inter-model spread in L-A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

Kirsten L. Findell↗

Central difference TVD and TVB schemes for time dependent and steady state problems

We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.

Jorgenson, P.↗

Central difference TVD and TVB schemes for time dependent and steady state problems

We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.

Jorgenson, P.↗

The Two Micron All Sky Survey

The Two Micron All Sky Survey (2MASS) will provide a uniform survey of the entire sky at three near-infrared wavebands: J(lambda(sub eff) = 1.25 micrometers), H(lambda(sub eff) = 1.65 micrometers), and K(sub s)(lambda(sub eff) = 2.16 micrometers). A major goal of the survey is to probe large scale structures in the Milky Way and in the Local Universe, exploiting the relatively high transparency of the interstellar medium in the near-infrared, and the high near-infrared luminosities of evolved low- and intermediate-mass stars. A sensitive overview of the near-infrared sky is also an essential next step to maximize the gains achievable with infrared array technology. Our assessment of the astrophysical questions that might be addressed with these new arrays is currently limited by the very bright flux limit of the only preceding large scale near-infrared sky survey, the Two Micron Sky Survey carried out at Caltech in the late 1960's. Near-infrared instruments based on the new array technology have already obtained spectra of objects 1 million times fainter than the limit of the TMSS! This paper summarizes the essential parameters of the 2MASS project and the rationale behind those choices, and gives an overview of results obtained with a prototype camera that has been in operation since May 1992. We conclude with a list of expected data products and a statement of the data release policy.

Kleinmann, S. G.↗

Phase holograms in PMMA with proximity effect correction

Complex computer generated phase holograms (CGPH's) have been fabricated in PMMA by partial e-beam exposure and subsequent partial development. The CGPH was encoded as a sequence of phase delay pixels and written by the JEOL JBX-5D2 E-beam lithography system, a different dose being assigned to each value of phase delay. Following carefully controlled partial development, the pattern appeared rendered in relief in the PMMA, which then acts as the phase-delay medium. The exposure dose was in the range 20-200 micro-C/sq cm, and very aggressive development in pure acetone led to low contrast. This enabled etch depth control to better than plus or minus lambda(sub vis)/60. That result was obtained by exposing isolated 50 micron square patches and measuring resist removal over the central area where the proximity effect dose was uniform and related only to the local exposure. For complex CGPH's with pixel size of the order of the e-beam proximity effect radius, the patterns must be corrected for the extra exposure caused by electrons scattered back up out of the substrate. This has been accomplished by deconvolving the two-dimensional dose deposition function with the desired dose pattern. The deposition function, which plays much the same role as an instrument response function, was carefully measured under the exact conditions used to expose the samples. The devices fabricated were designed with 16 equal phase steps per retardation cycle, were up to 1 cm square, and consisted of up to 100 million 0.3-2.0 micron square pixels. Data files were up to 500 MB long and exposure times ranged to tens of hours. A Fresnel phase lens was fabricated that had diffraction limited optical performance with better than 85 percent efficiency.

Maker, Paul D.↗

LPV Modeling of a Flexible Wing Aircraft Using Modal Alignment and Adaptive Gridding Methods

One of the earliest approaches in gain-scheduling control is the gridding based approach, in which a set of local linear time-invariant models are obtained at various gridded points corresponding to the varying parameters within the flight envelop. In order to ensure smooth and effective Linear Parameter-Varying control, aligning all the flexible modes within each local model and maintaining small number of representative local models over the gridded parameter space are crucial. In addition, since the flexible structural models tend to have large dimensions, a tractable model reduction process is necessary. In this paper, the notion of s-shifted H2- and H Infinity-norm are introduced and used as a metric to measure the model mismatch. A new modal alignment algorithm is developed which utilizes the defined metric for aligning all the local models over the entire gridded parameter space. Furthermore, an Adaptive Grid Step Size Determination algorithm is developed to minimize the number of local models required to represent the gridded parameter space. For model reduction, we propose to utilize the concept of Composite Modal Cost Analysis, through which the collective contribution of each flexible mode is computed and ranked. Therefore, a reduced-order model is constructed by retaining only those modes with significant contribution. The NASA Generic Transport Model operating at various flight speeds is studied for verification purpose, and the analysis and simulation results demonstrate the effectiveness of the proposed modeling approach.

LPV Modeling↗

Effect of Space Radiation on TES Detectors Performance

The Athena mission and its X-IFU instrument (X-ray Integral Field Unit) will be positioned at the Sun-Earth Lagrangian point L1, where it will be subject to solar wind (low flux plasma of 95% protons and 5% alpha particles) with energy below 0.1 MeV, and to galactic cosmic rays and solar flares (energetic protons, alpha particles, and electrons) with energies up to hundreds of MeV for protons and the GeV for heavier ions. Some of these particles will go through the satellite and hit the focal plane assembly and hence the detectors. These detectors will be TES (Transition-edge sensor) microcalorimeters, flown for the first time in such an environment. In order to ensure the performance of this type of detectors throughout the duration of such mission, it is critical to study the impact of the radiation on their behavior. Indeed, although a lot of reference material exist for semiconductor detectors such as CCDs, little is currently known about the impact of radiation on TES detectors. These energetic events could cause local heating or damage to the detectors and affect their performance. In this work, we describe how we designed a test campaign to assess the impact of L1 radiation on TES detectors for Athena/X-IFU-like missions and present the results of the tests. Analyses includes assessing changes in the pulse shapes and energy resolution of the detectors measured at 55 mK after several radiation dose steps performed at 4 K.

Sophie Beaumont↗

Improved Atmospheric Soundings and Error Estimates from Analysis of AIRS/AMSU Data

The AIRS Science Team Version 5.0 retrieval algorithm became operational at the Goddard DAAC in July 2007 generating near real-time products from analysis of AIRS/AMSU sounding data. This algorithm contains many significant theoretical advances over the AIRS Science Team Version 4.0 retrieval algorithm used previously. Three very significant developments of Version 5 are: 1) the development and implementation of an improved Radiative Transfer Algorithm (RTA) which allows for accurate treatment of non-Local Thermodynamic Equilibrium (non-LTE) effects on shortwave sounding channels; 2) the development of methodology to obtain very accurate case by case product error estimates which are in turn used for quality control; and 3) development of an accurate AIRS only cloud clearing and retrieval system. These theoretical improvements taken together enabled a new methodology to be developed which further improves soundings in partially cloudy conditions, without the need for microwave observations in the cloud clearing step as has been done previously. In this methodology, longwave C02 channel observations in the spectral region 700 cm-' to 750 cm-' are used exclusively for cloud clearing purposes, while shortwave C02 channels in the spectral region 2195 cm-' to 2395 cm-' are used for temperature sounding purposes. The new methodology for improved error estimates and their use in quality control is described briefly and results are shown indicative of their accuracy. Results are also shown of forecast impact experiments assimilating AIRS Version 5.0 retrieval products in the Goddard GEOS 5 Data Assimilation System using different quality control thresholds.

Susskind, Joel↗

Ground Based Studies of Thermocapillary Flows in Levitated Drops

Ground-based experiments together with analytical studies are presently being conducted for levitated drops. Both acoustic and electrostatic techniques are being employed to achieve levitation of drops in a gaseous environment. The scientific effort is principally on the thermal and the fluid phenomena associated with the local heating of levitated drops, both at 1-g and at low-g. In particular, the thermocapillary flow associated with local spot heating is being studied. Fairly stable acoustic levitation of drops has been achieved with some exceptions when random rotational motion of the drop persists. The flow visualization has been carried out by light scattering from smoke particles for the exterior flow and fluorescent tracer particles in the drop. The results indicate a lack of axial symmetry in the internal flow even though the apparatus and the heating are symmetric. The theoretical studies for the past year have included fundamental analyses of acoustically levitated spherical drops. The flow associated with a particle near the velocity antinode is being investigated by the singular perturbation technique. As a first step towards understanding the effect of the particle displacement from the antinode, the flow field about the node has been calculated for the first time. The effect of the acoustic field on the interior of a liquid drop has also been investigated. The results predict that the internal flow field is very weak.

Sadhal, Satwindar Singh↗

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Dejong, Frederik J.↗

Non-Linear Analysis of the NASA Super Pressure Balloons: Whole Flight Simulations

Tensys have a long-established background in the shape generation and load analysis of architectural stressed membrane structures. Founded upon their inTENS finite element analysis suite, these activities have broadened to encompass 'lighter than air' structures such as aerostats, airships, hybrid air vehicles and stratospheric balloons. Tensys have acted as consultants to the NASA Super Pressure Balloon (SPB) Program since 2004. Previous papers have focussed upon the application of inTENS to the overall structural and stability analysis of pumpkin type balloons as used in the SPB Program. Particular emphasis has been placed upon the ability to study both stress and stability at all stages of a flight. As the program has developed, increasing modelling fidelity has been introduced as design refinement has moved emphasis from the overall shape to the performance of individual details. Examples include the introduction of a contact capability into inTENS to better represent the separate cap layer of film. Localised investigations have considered the consequence of local geometric anomalies introduced during fabrication and the effects of debonding of the PBO tendon within its sleeve. Analysis to date has used a material model for the polyethylene shell film developed by Dr Rand of Winzen Engineering, supported by a program of fine resolution material tests at the Balloon Research and Development Laboratory (BRDL) at the GSFC Wallops Island Facility. Based upon work first reported by Schapery, this model has been incorporated into so-called 'snapshot' analyses within inTENS. For a given elapsed time point and temperature distribution, non-linear material properties are iteratively updated for the current stress state in each individual element until changes in those material properties are insignificant. This process has to be iterative as the film stresses and material properties are interdependent. Attempts to derive a time-stepping incremental viscoelastic capability using the Schapery Rand model encountered problems when dealing with the varying temperatures associated with balloon deployment, pressurisation and diurnal behaviour. An alternative approach has been developed by Pellegrino et al at Caltech, again with support from the NASA Balloon Program Office. This is a large strain non-linear viscoelastic model that includes film out-of-plane mechanical and thermal effects. Working from a unified base, this model can be utilised for time-stepping analyses in both a modulus or compliance mode, or a combination of both. This new capability is intended to enhance the current NASA balloon design process. This paper presents the implementation of the Caltech model into a specialist finite element suite, inTENS and its application to whole flight simulations.

Wakefield, David↗