SEARCH · Engineering Papers

Results for “phase field method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Cloud properties inferred from 8-12 micron data

A trispectral combination of observations at 8-, 11-, and 12-micron bands is suggested for detecting cloud and cloud properties in the infrared. Atmospheric ice and water vapor absorption peak in opposite halves of the window region so that positive 8-minus-11-micron brightness temperature differences indicate cloud, while near-zero or negative differences indicate clear regions. The absorption coefficient for water increases more between 11 and 12 microns than between 8 and 11 microns, while for ice, the reverse is true. Cloud phases is determined by a scatter diagram of 8-minus-11-micron versus 11-minus-12-micron brightness temperature differences; ice cloud shows a slope greater than 1 and water cloud less than 1. The trispectral brightness temperature method was tested upon high-resolution interferometer data resulting in clear-cloud and cloud-phase delineation. Simulations using differing 8-micron bandwidths revealed no significant degradation of cloud property detection. Thus, the 8-micron bandwidth for future satellites can be selected based on the requirements of other applications, such as surface characterization studies. Application of the technique to current polar-orbiting High-Resolution Infrared Sounder (HIRS)-Advanced Very High Resolution Radiometer (AVHRR) datasets is constrained by the nonuniformity of the cloud scenes sensed within the large HIRS field of view. Analysis of MAS (MODIS Airborne Simulator) high-spatial resolution (500 m) data with all three 8-, 11-, and 12-micron bands revealed sharp delineation of differing cloud and background scenes, from which a simple automated threshold technique was developed. Cloud phase, clear-sky, and qualitative differences in cloud emissivity and cloud height were identified on a case study segment from 24 November 1991, consistent with the scene. More rigorous techniques would allow further cloud parameter clarification. The opportunities for global cloud delineation with the Moderate-Resolution Imaging Spectrometer (MODIS) appear excellent. The spectral selection, the spatial resolution, and the global coverage are all well suited for significant advances.

Strabala, Kathleen I.↗

NASA NTRS2007

The Spatiotemporal Structure of 20th Century Climate Variations in Observations and Reanalyses: Pacific Pan-Decadal Variability - Part 2

The dominant interannual El Nino-Southern Oscillation phenomenon (ENSO) and the short length of climate observation records make it difficult to study long-term climate variations in the spatiotemporal domain. Based on the fact that the ENS0 signal spreads to remote regions and induces delayed climate variation through atmospheric teleconnections, we develop an ENSO-removal method through which the ENS0 signal can be approximately removed at the grid box level from the spatiotemporal field of a climate parameter. After this signal is removed, long-term climate variations, namely, the global warming trend (GW) and the Pacific pan-decadal variability (PDV), are isolated at middle and low latitudes in the climate parameter fields from observed and reanalyses datasets. In this study, we show that one of several PDV interdecadal regime shifts occurred during the 1990s. This significant change in the Pacific basin is comparable but opposite in phase to the 1976 climate regime shift, which results persisting warming in the central-eastern Pacific, and cooling in the North and South Pacific. The 1990s PDV regime shift is consistent with observed changes in ocean biosphere and ocean circulation. A comprehensive picture of PDV as manifested in the troposphere and at the surface is described. In general, the PDV spatial patterns in different parameter fields share some similarities with the patterns associated with ENSO, but important differences exist. First, the PDV atmospheric circulation pattern is shifted westward by about 20deg and its zonal extent is limited to approx.60deg compared to approx.110deg for ENS0 pattern. The westward shift of the PDV wave train produces a different, more west-east oriented, North American teleconnection pattern. The lack of a strong PDV surface temperature (ST) signal in the western equatorial Pacific and the relatively strong ST signal in the subtropical regions are consistent with an atmospheric overturning circulation response that differs from the one associated with ENSO.

Chen, Junye↗

NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗

NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

processing↗

NASA NTRS1995

Boundary layer transition studies

A small-scale wind tunnel previously used for turbulent boundary layer experiments was modified for two sets of boundary layer transition studies. The first study concerns a laminar separation/turbulent reattachment. The pressure gradient and unit Reynolds number are the same as the fully turbulent flow of Spalart and Watmuff. Without the trip wire, a laminar layer asymptotes to a Falkner & Skan similarity solution in the FPG. Application of the APG causes the layer to separate and a highly turbulent and approximately 2D mean flow reattachment occurs downstream. In an effort to gain some physical insight into the flow processes a small impulsive disturbance was introduced at the C(sub p) minimum. The facility is totally automated and phase-averaged data are measured on a point-by-point basis using unprecedently large grids. The evolution of the disturbance has been tracked all the way into the reattachment region and beyond into the fully turbulent boundary layer. At first, the amplitude decays exponentially with streamwise distance in the APG region, where the layer remains attached, i.e. the layer is viscously stable. After separation, the rate of decay slows, and a point of minimum amplitude is reached where the contours of the wave packet exhibit dispersive characteristics. From this point, exponential growth of the amplitude of the disturbance is observed in the detached shear layer, i.e. the dominant instability mechanism is inviscid. A group of large-scale 3D vortex loops emerges in the vicinity of the reattachment. Remarkably, the second loop retains its identify far downstream in the turbulent boundary layer. The results provide a level of detail usually associated with CFD. Substantial modifications were made to the facility for the second study concerning disturbances generated by Suction Holes for laminar flow Control (LFC). The test section incorporates suction through interchangeable porous test surfaces. Detailed studies have been made using isolated holes in the impervious test plate that used to establish the Blasius base flow. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance, for conditions corresponding to strong suction and without suction. The technique was enhanced by using up to nine multiple probes to reduce the experimental run-time. In both cases, 3D contour surfaces in the vicinity of the hole show highly 3D TS waves which fan out in the spanwise direction forming bow-shaped waves downstream. The case without suction has proved useful for evaluating calculation methods. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the TS waves in the far field are similar to the case without suction. Downstream, the contour surfaces of the TS waves develop spanwise irregularities which eventually form into clumps. The spanwise clumping is evidence of a secondary instability that could be associated with suction vortices. Designers of porous surfaces use Goldsmith's Criterion to minimize cross-stream interactions. It is shown that partial TS wave cancellation is possible, depending on the hole spacing, disturbance frequency and free-stream velocity. New high-performance Constant Temperature Hot-Wire Anemometers were designed and built, based on a linear system theory analysis that can be extended to arbitrary order. The motivation was to achieve the highest possible frequency reponse while ensuring overall system stability. The performance is equal to or superior to commercially available instruments at about 10% of the cost. Details, such as fabrication drawings and a parts list, have been published to enable the instrument to be construced by others.

Watmuff, Jonathan H.↗

NASA NTRS2010

Characterizing the Space Debris Environment with a Variety of SSA Sensors

Damaging space debris spans a wide range of sizes and altitudes. Therefore no single method or sensor can fully characterize the space debris environment. Space debris researchers use a variety of radars and optical telescopes to characterize the space debris environment in terms of number, altitude, and inclination distributions. Some sensors, such as phased array radars, are designed to search a large volume of the sky and can be instrumental in detecting new breakups and cataloging and precise tracking of relatively large debris. For smaller debris sizes more sensitivity is needed which can be provided, in part, by large antenna gains. Larger antenna gains, however, produce smaller fields of view. Statistical measurements of the debris environment with less precise orbital parameters result. At higher altitudes, optical telescopes become the more sensitive instrument and present their own measurement difficulties. Space Situational Awareness, or SSA, is concerned with more than the number and orbits of satellites. SSA also seeks to understand such parameters as the function, shape, and composition of operational satellites. Similarly, debris researchers are seeking to characterize similar parameters for space debris to improve our knowledge of the risks debris poses to operational satellites as well as determine sources of debris for future mitigation. This paper will discuss different sensor and sensor types and the role that each plays in fully characterizing the space debris environment.

Stansbery, Eugene G.↗

NASA NTRS2018

High-Efficiency Flexible Multilevel Photon Sieves by Single-Step Laser-Based Fabrication and Optical Analysis

Over the past several decades, the need for high-resolution, high-efficiency, lightweight, high contrast focusing optics has continued to increase due to their applications in fields such as astronomy, spectroscopy, free-space optical communications, defense, and remote sensing. In recent years, photon sieve planar diffractive optics have been developed on flexible, lightweight polyimide substrates. However, transmission efficiencies have continuously been very low (~1-4%), thus impeding the widespread use of photon sieves in practical applications. Here, we present a flexible, lightweight, N = 4 level phase photon sieve with over 25% transmission efficiency; nearly triple that of any other photon sieve reported thus far. In this study, the photon sieve was fabricated via a novel pulsed laser ablation method using an ultraviolet source, and the total time to fabricate a ~3 sq cm sample was tens of seconds. Theoretical analysis of the photon sieve was carried out via the FDTD method, and was in very good agreement with experimental results. We also calculated via FDTD modeling the behavior of higher step (N = 8, 16, 32) photon sieves for further enhanced efficiencies, and show a fundamental limit on photon sieve efficiency of 70% in the limit of increasing step number. This multilevel photon sieve represents a new step in high-resolution diffractive optics, moving towards efficiencies suitable for widespread applications.

Julian, Matthew N.↗

NASA NTRS2006

Jet Noise Diagnostics Supporting Statistical Noise Prediction Methods

The primary focus of my presentation is the development of the jet noise prediction code JeNo with most examples coming from the experimental work that drove the theoretical development and validation. JeNo is a statistical jet noise prediction code, based upon the Lilley acoustic analogy. Our approach uses time-average 2-D or 3-D mean and turbulent statistics of the flow as input. The output is source distributions and spectral directivity. NASA has been investing in development of statistical jet noise prediction tools because these seem to fit the middle ground that allows enough flexibility and fidelity for jet noise source diagnostics while having reasonable computational requirements. These tools rely on Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) solutions as input for computing far-field spectral directivity using an acoustic analogy. There are many ways acoustic analogies can be created, each with a series of assumptions and models, many often taken unknowingly. And the resulting prediction can be easily reverse-engineered by altering the models contained within. However, only an approach which is mathematically sound, with assumptions validated and modeled quantities checked against direct measurement will give consistently correct answers. Many quantities are modeled in acoustic analogies precisely because they have been impossible to measure or calculate, making this requirement a difficult task. The NASA team has spent considerable effort identifying all the assumptions and models used to take the Navier-Stokes equations to the point of a statistical calculation via an acoustic analogy very similar to that proposed by Lilley. Assumptions have been identified and experiments have been developed to test these assumptions. In some cases this has resulted in assumptions being changed. Beginning with the CFD used as input to the acoustic analogy, models for turbulence closure used in RANS CFD codes have been explored and compared against measurements of mean and rms velocity statistics over a range of jet speeds and temperatures. Models for flow parameters used in the acoustic analogy, most notably the space-time correlations of velocity, have been compared against direct measurements, and modified to better fit the observed data. These measurements have been extremely challenging for hot, high speed jets, and represent a sizeable investment in instrumentation development. As an intermediate check that the analysis is predicting the physics intended, phased arrays have been employed to measure source distributions for a wide range of jet cases. And finally, careful far-field spectral directivity measurements have been taken for final validation of the prediction code. Examples of each of these experimental efforts will be presented. The main result of these efforts is a noise prediction code, named JeNo, which is in middevelopment. JeNo is able to consistently predict spectral directivity, including aft angle directivity, for subsonic cold jets of most geometries. Current development on JeNo is focused on extending its capability to hot jets, requiring inclusion of a previously neglected second source associated with thermal fluctuations. A secondary result of the intensive experimentation is the archiving of various flow statistics applicable to other acoustic analogies and to development of time-resolved prediction methods. These will be of lasting value as we look ahead at future challenges to the aeroacoustic experimentalist.

Bridges, James E.↗

NASA NTRS1982

Initial stages of cavitation damage and erosion on copper and brass tested in a rotating disk device

In view of the differences in flow and experimental conditions, there has been a continuing debate as to whether or not the ultrasonic method of producing cavitation damage is similar to the damage occurring in cavitating flow systems, namely, venturi and rotating disk devices. In this paper, the progress of cavitation damage during incubation periods on polycrystalline copper and brass tested in a rotating disk device is presented. The results indicate several similarities and differences in the damage mechanism encountered in a rotating disk device (which simulates field rotary devices) and a magnetostriction apparatus. The macroscopic erosion appears similar to that in the vibratory device except for nonuniform erosion and apparent plastic flow during the initial damage phase.

Rao, P. V.↗

NASA NTRS1984

Analysis of the surface load and radiated sound of a vibrating airfoil with application to the experiment of Brooks

A method is developed for calculating the surface load and radiated sound from a vibrating surface in a compressible viscous fluid. The method is applied to a thin two-dimensional elliptic cross-section. For large values of the viscous diffusion parameter, the surface load tends to an elliptic distribution in agreement with the results of inviscid theory when edge pressure continuity is enforced. For thin surfaces, the surface load is insensitive to variations in the thickness ratio. A three-dimensional spectral technique is developed to calculate the inviscid surface load and radiated sound from a thin vibrating airfoil. The inviscid theory predicts the correct form of the far field sound pressure and its phase. The actual levels are somewhat sensitive to the choice of theoretical spanwise surface pressure mode but are in better agreement with the experiment than the surface pressure. The comparison of theoretical and experimental surface pressure indicates that the viscous theory, used to validate the inviscid theory, is either inadequate or there is a source of experimental error.

Yates, J. E.↗

NASA NTRS1992

A comparative evaluation of in-plane shear test methods for laminated graphite-epoxy composites

The objectives were to evaluate popular shear test methods for various forms of graphite-epoxy composite materials and to determine the shear response of graphite-epoxy composites with various forms of fiber architecture. Numerical and full-field experimental stress analyses were performed on four shear test configurations for unidirectional and bidirectional graphite-epoxy laminates to assess the uniformity and purity of the shear stress (strain) fields produced in the specimen test section and to determine the material in-plane shear modulus and shear response. The test methods were the 10 deg off-axis, the +/- 45 deg tension, the Iosipescu V-notch, and a compact U-notch specimen. Specimens were prepared from AS4/3501-6 graphite-epoxy panels, instrumented with conventional strain gage rosettes and with a cross-line moire grating, and loaded in a convenient testing machine. The shear responses obtained for each test method and the two methods of specimen instrumentation were compared. In a second phase of the program the shear responses obtained from Iosipescu V-notch beam specimens were determined for woven fabric geometries of different weave and fiber architectures. Again the responses of specimens obtained from strain gage rosettes and moire interferometry were compared. Additional experiments were performed on a bidirectional cruciform specimen which was also instrumented with strain gages and a moire grating.

Morton, John↗

NASA NTRS2011

Method for Implementing Optical Phase Adjustment

A method has been developed to mechanically implement the optical phase shift by adjusting the polarization of the pump and probe beams in an NMOR (nonlinear magneto-optical rotation) magnetometer as the proper phase shift is necessary to induce self-oscillation. This innovation consists of mounting the pump beam on a ring that surrounds the atomic vapor sample. The propagation of the probe beam is perpendicular to that of the pump beam. The probe beam can be considered as defining the axis of a cylinder, while the pump beam is directed radially. The magnetic field to be measured defines a third vector, but it is also taken to lie along the cylinder axis. Both the pump and probe beams are polarized such that their electric field vectors are substantially perpendicular to the magnet field. By rotation of the ring supporting the pump beam, its direction can be varied relative to the plane defined by the probe electric field and the magnetic field to be measured.

Hovde, David C.↗

NASA NTRS1990

The effects of forcing on a single stream shear layer and its parent boundary layer

Forcing and its effect on fluid flows has become an accepted tool in the study and control of flow systems. It has been used both as a diagnostic tool, to explore the development and interaction of coherent structures, and as a method of controlling the behavior of the flow. A number of forcing methods have been used in order to provide a perturbation to the flow; among these are the use of an oscillating trailing edge, acoustically driven slots, external acoustic forcing, and mechanical piston methods. The effect of a planar mechanical piston forcing on a single stream shear layer is presented; it can be noted that this is one of the lesser studied free shear layers. The single stream shear layer can be characterized by its primary flow velocity scale and the thickness of the separating boundary layer. The velocity scale is constant over the length of the flow field; theta (x) can be used as a width scale to characterize the unforced shear layer. In the case of the forced shear layer the velocity field is a function of phase time and definition of a width measure becomes somewhat problematic.

Haw, Richard C.↗

NASA NTRS2000

Reliability Modeling of Microelectromechanical Systems Using Neural Networks

Microelectromechanical systems (MEMS) are a broad and rapidly expanding field that is currently receiving a great deal of attention because of the potential to significantly improve the ability to sense, analyze, and control a variety of processes, such as heating and ventilation systems, automobiles, medicine, aeronautical flight, military surveillance, weather forecasting, and space exploration. MEMS are very small and are a blend of electrical and mechanical components, with electrical and mechanical systems on one chip. This research establishes reliability estimation and prediction for MEMS devices at the conceptual design phase using neural networks. At the conceptual design phase, before devices are built and tested, traditional methods of quantifying reliability are inadequate because the device is not in existence and cannot be tested to establish the reliability distributions. A novel approach using neural networks is created to predict the overall reliability of a MEMS device based on its components and each component's attributes. The methodology begins with collecting attribute data (fabrication process, physical specifications, operating environment, property characteristics, packaging, etc.) and reliability data for many types of microengines. The data are partitioned into training data (the majority) and validation data (the remainder). A neural network is applied to the training data (both attribute and reliability); the attributes become the system inputs and reliability data (cycles to failure), the system output. After the neural network is trained with sufficient data. the validation data are used to verify the neural networks provided accurate reliability estimates. Now, the reliability of a new proposed MEMS device can be estimated by using the appropriate trained neural networks developed in this work.

Perera. J. Sebastian↗

NASA NTRS2005

System concepts for transmit arrays of parabolic antennas for deep space uplinks

Phased arrays of parabolic antennas are a potentially lower-cost way to provide uplink transmission to distant spacecraft, compared to the 34-m and 70-m antennas now used by the NASA Deep Space Network. A large transmit array could provide very high EIRP when needed for spacecraft emergencies, such as the equivalent of 1 MW radiated from a 70-m antenna. Cost-effectiveness is realized by dividing the array into smaller arrays to provide routine support to many spacecraft simultaneously. The antennas might be as small as 12-m in diameter, with as many as 100 antennas covering an area of 0.5 km to 1 km in extent. Such arrays present significant technical challenges in phase alignment, which must be maintained at close to 1 mm. The concept requires a very stable system with accurately known antenna phase center locations. The system is first calibrated by transmitting from all antennas, and observing the signals at a target located in the far fields of the individual antennas. The antennas are then pointed to the operational targets, with the signal phases and time delays set to reinforce in the target directions. This requires accurate knowledge of the target directions and calculation of the required phases. The system must be phase-stable for all directions and over the time between calibrations, which is desired to be at least one day. In this paper, a system concept is presented, the major error sources are identified, a rough error budget is established, and key elements of the system are discussed. A calibration method is recommended which uses satellites as radar targets. The performance goal is to achieve a combining loss of less than 0.2 dB in good weather, and of less than 1 dB in all but extremely bad weather.

Effective Isotropic Radiated Power (EIRP)↗

NASA NTRS2004

Fast-Acquisition/Weak-Signal-Tracking GPS Receiver for HEO

A report discusses the technical background and design of the Navigator Global Positioning System (GPS) receiver -- . a radiation-hardened receiver intended for use aboard spacecraft. Navigator is capable of weak signal acquisition and tracking as well as much faster acquisition of strong or weak signals with no a priori knowledge or external aiding. Weak-signal acquisition and tracking enables GPS use in high Earth orbits (HEO), and fast acquisition allows for the receiver to remain without power until needed in any orbit. Signal acquisition and signal tracking are, respectively, the processes of finding and demodulating a signal. Acquisition is the more computationally difficult process. Previous GPS receivers employ the method of sequentially searching the two-dimensional signal parameter space (code phase and Doppler). Navigator exploits properties of the Fourier transform in a massively parallel search for the GPS signal. This method results in far faster acquisition times [in the lab, 12 GPS satellites have been acquired with no a priori knowledge in a Low-Earth-Orbit (LEO) scenario in less than one second]. Modeling has shown that Navigator will be capable of acquiring signals down to 25 dB-Hz, appropriate for HEO missions. Navigator is built using the radiation-hardened ColdFire microprocessor and housing the most computationally intense functions in dedicated field-programmable gate arrays. The high performance of the algorithm and of the receiver as a whole are made possible by optimizing computational efficiency and carefully weighing tradeoffs among the sampling rate, data format, and data-path bit width.

Wintemitz, Luke↗

NASA NTRS1996

LCTV Holographic Imaging

Astronauts are required to interface with complex systems that require sophisticated displays to communicate effectively. Lightweight, head-mounted real-time displays that present holographic images for comfortable viewing may be the ideal solution. We describe an implementation of a liquid crystal television (LCTV) as a spatial light modulator (SLM) for the display of holograms. The implementation required the solution of a complex set of problems. These include field calculations, determination of the LCTV-SLM complex transmittance characteristics and a precise knowledge of the signal mapping between the LCTV and frame grabbing board that controls it. Realizing the hologram is further complicated by the coupling that occurs between the phase and amplitude in the LCTV transmittance. A single drive signal (a gray level signal from a framegrabber) determines both amplitude and phase. Since they are not independently controllable (as is true in the ideal SLM) one must deal with the problem of optimizing (in some sense) the hologram based on this constraint. Solutions for the above problems have been found. An algorithm has been for field calculations that uses an efficient outer product formulation. Juday's MEDOF 7 (Minimum Euclidean Distance Optimal Filter) algorithm used for originally for filter calculations has been successfully adapted to handle metrics appropriate for holography. This has solved the problem of optimizing the hologram to the constraints imposed by coupling. Two laboratory methods have been developed for determining an accurate mapping of framegrabber pixels to LCTV pixels. A friendly software system has been developed that integrates the hologram calculation and realization process using a simple set of instructions. The computer code and all the laboratory measurement techniques determining SLM parameters have been proven with the production of a high quality test image.

Knopp, Jerome↗

NASA NTRS1996

Parallelization of Rocket Engine System Software (Press)

The main goal is to assess parallelization requirements for the Rocket Engine Numeric Simulator (RENS) project which, aside from gathering information on liquid-propelled rocket engines and setting forth requirements, involve a large FORTRAN based package at NASA Lewis Research Center and TDK software developed by SUBR/UWF. The ultimate aim is to develop, test, integrate, and suitably deploy a family of software packages on various aspects and facets of rocket engines using liquid-propellants. At present, all project efforts by the funding agency, NASA Lewis Research Center, and the HBCU participants are disseminated over the internet using world wide web home pages. Considering obviously expensive methods of actual field trails, the benefits of software simulators are potentially enormous. When realized, these benefits will be analogous to those provided by numerous CAD/CAM packages and flight-training simulators. According to the overall task assignments, Hampton University's role is to collect all available software, place them in a common format, assess and evaluate, define interfaces, and provide integration. Most importantly, the HU's mission is to see to it that the real-time performance is assured. This involves source code translations, porting, and distribution. The porting will be done in two phases: First, place all software on Cray XMP platform using FORTRAN. After testing and evaluation on the Cray X-MP, the code will be translated to C + + and ported to the parallel nCUBE platform. At present, we are evaluating another option of distributed processing over local area networks using Sun NFS, Ethernet, TCP/IP. Considering the heterogeneous nature of the present software (e.g., first started as an expert system using LISP machines) which now involve FORTRAN code, the effort is expected to be quite challenging.

Cezzar, Ruknet↗