SEARCH · Search NASA
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.
In this paper, we present a mathematical foundation, including a convergence analysis, for cascading architecture neural networks. From this, a mathematical foundation for the casade correlation learning algorithm can also be found. Furthermore, it becomes apparent that the cascade correlation scheme is a special case of an efficient hardware learning algorithm called Cascade Error Projection.
The most sensitive direct neutrino mass searches today are based on measurement of the end point of the 𝛽 spectrum of tritium to infer limits on the mass of the unobserved neutrino. To avoid the smearing associated with the distribution of molecular final states in the T-He molecule, the next generation of these experiments will need to employ atomic (T) rather than molecular (T 2 ) tritium sources, at currents of at least 10 15 atoms per second. Following production, atomic T can be trapped in gravitational and/or magnetic bottles for 𝛽 spectrum experiments, if and only if it can first be cooled to millikelvin temperatures. Accomplishing this cooling presents substantial technological challenges. The Project 8 collaboration is developing a technique based on magnetic evaporative cooling along a beamline (MECB) for the purpose of cooling T to feed a magnetogravitational trap that also serves as a cyclotron radiation emission spectroscope. Initial tests of the approach are planned in a pathfinder apparatus using atomic Li. Here, this paper presents a method for analyzing the dynamics of the MECB technique and applies these calculations to the design of systems for cooling and slowing of atomic Li and T. A scheme is outlined that could provide a current of T at the millikelvin temperatures required for the Project 8 neutrino mass search.
The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored
Task 5 consisted of two parts, wave energy converter (WEC) arrays and current energy converter (CEC) arrays. The two parts of the project were executed independently, and this report covers only the work on WECs. In this part of the task, the team characterized the behavior of and developed control schemes for wave energy converter (WEC) arrays that improved performance (i.e., maximized aggregate power generation and reduced the levelized cost of energy) over baseline, non-coordinated control approaches. Specifically, for WEC arrays: the team developed optimal layouts of arrays that considered the effect of WEC placement within an array on coordinated array control with real-time estimation schemes. The numerical codes developed through this task were made available to users for array design. These tools utilized existing commercial software as well as software developed within the project team.
Measurements in many-body quantum systems can generate non-trivial phenomena, such as preparation of long-range entangled states, dynamical phase transitions, or measurement-altered criticality. Here, we introduce a new measurement scheme that produces an ensemble of mixed states in a subsystem, obtained by measuring a local Hermitian observable on part of its complement. We refer to this as the observable-projected ensemble . Unlike standard projected ensembles-where pure states are generated by projective measurements on the complement-our approach involves projective partial measurements of specific observables. This setup has two main advantages: theoretically, it is amenable to analytical computations, especially within conformal field theories. Experimentally, it requires only a linear number of measurements, rather than an exponential one, to probe the properties of the ensemble. As a first step in exploring the observable-projected ensemble, we investigate its entanglement properties in conformal field theory and perform a detailed analysis of the free compact boson.
In the currently popular orientation-based unified scheme, a radio galaxy appears as a quasar when its principal radio-axis happens to be oriented within a certain cone opening angle around the observer's line of sight. Due to geometrical projection, the observed sizes of quasars should therefore appear smaller than those of radio galaxies. We show that this simple, unambiguous prediction of the unified scheme is not borne out by the actually observed angular sizes of radio galaxies and quasars. Except in the original 3CR sample, based on which the unified scheme was proposed, in other much larger samples no statistically significant difference is apparent in the size distributions of radio galaxies and quasars. The population of low-excitation radio galaxies with apparently no hidden quasars inside, which might explain the observed excess number of radio galaxies at low redshifts, cannot account for the absence of any foreshortening of the sizes of quasars at large redshifts. On the other hand, from infrared and X-ray studies, there is evidence of a hidden quasar within a dusty torus in many radio galaxies, at z > 0.5. It is difficult to reconcile this with the absence of foreshortening of quasar sizes at even these redshifts, and perhaps one has to allow that the major radio axis may not have anything to do with the optical axis of the torus. Otherwise, to resolve the dichotomy of radio galaxies and quasars, a scheme quite different from the present might be required.
We describe the development and application of stochastic optimization models and algorithms to address an issue of critical importance in the strategic allocation of resources; namely, the selection of a portfolio of capital investment projects under the constraints of a limited and uncertain budget. This issue is significant and one that faces decision-makers across all industries. The objective of this strategic decision process is generally self evident - to maximize the value obtained from the portfolio of selected projects (with value usually measured in terms of the portfolio's net present value). However, heretofore, many organizations have developed processes to make these investment decisions using simplistic rule-based rank-ordering schemes. This approach has the significant limitation of not accounting for the (often large) uncertainties in the costs or economic benefits associated with the candidate projects or in the uncertainties in the actual funds available to be expended over the projected period of time. As a result, the simple heuristic approaches that typically are employed in industrial practice generate outcomes that are non-optimal and do not achieve the level of benefits intended. In this paper we describe the results of research performed to utilize stochastic optimization models and algorithms to address this limitation by explicitly incorporating the evaluation of uncertainties in the analysis and decision making process. (authors)
A desktop-scale, computer-controlled display system, initially developed for NASA and now known as the VolumeViewer(TradeMark), generates three-dimensional (3D) images of 3D objects in a display volume. This system differs fundamentally from stereoscopic and holographic display systems: The images generated by this system are truly 3D in that they can be viewed from almost any angle, without the aid of special eyeglasses. It is possible to walk around the system while gazing at its display volume to see a displayed object from a changing perspective, and multiple observers standing at different positions around the display can view the object simultaneously from their individual perspectives, as though the displayed object were a real 3D object. At the time of writing this article, only partial information on the design and principle of operation of the system was available. It is known that the system includes a high-speed, silicon-backplane, ferroelectric-liquid-crystal spatial light modulator (SLM), multiple high-power lasers for projecting images in multiple colors, a rotating helix that serves as a moving screen for displaying voxels [volume cells or volume elements, in analogy to pixels (picture cells or picture elements) in two-dimensional (2D) images], and a host computer. The rotating helix and its motor drive are the only moving parts. Under control by the host computer, a stream of 2D image patterns is generated on the SLM and projected through optics onto the surface of the rotating helix. The system utilizes a parallel pixel/voxel-addressing scheme: All the pixels of the 2D pattern on the SLM are addressed simultaneously by laser beams. This parallel addressing scheme overcomes the difficulty of achieving both high resolution and a high frame rate in a raster scanning or serial addressing scheme. It has been reported that the structure of the system is simple and easy to build, that the optical design and alignment are not difficult, and that the system can be built by use of commercial off-the-shelf products. A prototype of the system displays an image of 1,024 by 768 by 170 (=133,693,440) voxels. In future designs, the resolution could be increased. The maximum number of voxels that can be generated depends upon the spatial resolution of SLM and the speed of rotation of the helix. For example, one could use an available SLM that has 1,024 by 1,024 pixels. Incidentally, this SLM is capable of operation at a switching speed of 300,000 frames per second. Implementation of full-color displays in future versions of the system would be straightforward: One could use three SLMs for red, green, and blue, respectively, and the colors of the voxels could be automatically controlled. An optically simpler alternative would be to use a single red/green/ blue light projector and synchronize the projection of each color with the generation of patterns for that color on a single SLM.
A new adaptive scheme to be used in particle-in-cell codes for carrying out gyroaveraging operations with matrices is presented. This new scheme uses an intermediate velocity grid whose resolution is adapted to the local thermal Larmor radius. The charge density is computed by projecting marker weights in a field-line following manner while preserving the adiabatic magnetic moment μ. These choices permit to improve the accuracy of the gyroaveraging operations performed with matrices even when strong spatial variation of temperature and magnetic field is present. Accuracy of the scheme in different geometries from simple 2D slab geometry to realistic 3D toroidal equilibrium has been studied. As a result, a successful implementation in the gyrokinetic code XGC is presented in the delta-f limit.
The Austin SHINES project and solution is a software management platform, for an electric grid with a high penetration of dispersed photovoltaic (PV) solar generation sites, which maintains the traditional power quality and reliability associated with grid service. This project developed and deployed the platform as a Distributed Energy Resource Management System (DERMS), engaging multiple advanced controls, to evaluate operation and optimization of a fleet of diverse DER assets, installed at several locations among Austin Energy’s customers and distribution system. The project also produced a methodology to create a replicable DERMS template, adaptable to other regions and market structures. Last, Austin SHINES aimed to demonstrate the solution’s methodology would enable the DER grid ecosystem to serve load at a technical cost (System Levelized Cost of Electricity, or System LCOE) of less than the U.S. Department of Energy SHINES program metric of $0.14/kWh, in a defined boundary, while enabling a high penetration of distributed PV. Research was categorized in 6 reports (Final Deliverables = FD) listed below, with titles and descriptions indicating which area of understanding was investigated: FD-1: System Levelized Cost of Electricity (System LCOE) Methodology The creation and use of the System LCOE to Serve Load metric that encompasses the holistic, system-level costs and benefits of all resources, and enables them to be evaluated based on their ability to support an efficient and low-cost integrated grid ecosystem. FD-2: Software Platform Product Description The creation of new DER control methodologies deployable within a utility-grade software platform that enable DER's to maximize their benefit within a grid, that is capable of serving load enabling a high penetration of distributed PV generation. FD-3: Optimal Design Methodology Optimal design methodologies for individual DER installations that enable utilities to determine the optimal combinations and sizing for individual DER sites. FD-4: Austin SHINES Ownership and Operation Models for DER System Performance A comparison of multiple DER aggregation and ownership methodologies including direct utility control, third-party aggregator, and autonomous. FD-5: Economic Modeling & Optimization A comparison of multiple DER technology mixes and configurations within the distribution system, providing insight into an optimal blend of technologies that best enable the distribution system to serve load at the lowest cost at high penetrations of solar. FD-6: Fielded Assets Deployed DER assets within the Austin Energy SHINES circuits. Austin SHINES provided an opening for state-of-the-art technology products to be deployed, providing a rich opportunity for improving how each of the products perform as stand-alone products, and in concert with other complementary products. The Austin SHINES project comprised of two key metrics for System LCOE: SystemLCOE_SHINES<$0.14/kWh Modeled ΔSystemLCOE_SHINES/ΔSystemLCOE_Base≥20% at same solar penetration The System LCOE calculation uses the costs of the utility-owned infrastructure as it exists today, the cost of the DERs that exist in the system today, and the cost of the purchase of energy from ERCOT wholesale markets over the course of the calendar year. All costs are on an annualized basis. The capital and operating costs are derived from the rate case, which produces a yearly cost. The net cost of energy and services imported to the system is integrated over the test year, as is the load served and solar penetration. The first metric was easily achieved by every scenario considered. The goal was set when the Department of Energy’s SHINES Funding Opportunity Announcement was written in 2015 and was a more difficult target at the time. Due mostly to rapidly declining costs for DERs and the significant decrease in the Electric Reliability Council of Texas (ERCOT) energy market prices, which results in lower net cost of energy purchases, the System LCOE is well below this target for all scenarios considered. A fleet of DERs can assume different mixtures, each of which serves the load at a different LCOE. The optimal mixture of DERs serves load at the smallest System LCOE. The second metric (hereinafter %delta metric) asks that the holistic DERMS controls reduce the incremental cost above the baseline of going to a high solar penetration future by at least 20% as compared to the case of a DER deployment with no sophisticated controls (autonomous). Many comparison sets were created throughout this project. Physical technology was installed for informing utility engineering and testing several types of operational control schemes, through the DERMS. The types of operational control which were compared for valuation of the System LCOE Metric were: Holistic control = using the full suite of the DERMS platform to decide and optimize how/why the systems operate depending on weather, market, and reliability signal input. Autonomous control = a local mode at the asset site, wherein a schedule operates the asset, with visibility into performance only No control = the baseline for comparing value against the other two types of control The types of ownership control included: Direct Utility control = the utility dispatches a signal to each asset Third-Party Aggregator = a third party aggregates a fleet of assets and the utility dispatches one signal for all Autonomous = a local mode is set for operation at the asset site, wherein a schedule operates the asset, with visibility into performance only The types of control methodologies deployable within a utility-grade software platform included: Utility Peak Load Reduction = Lower transmission cost obligation Day-Ahead Energy Arbitrage = Realize economic value through price differential Real-Time Price Dispatch = Realize economic value from real-time price spikes Voltage support = Reduce losses and increase solar generation Distribution Congestion Management = Increase local grid reliability Demand Charge Reduction = Lower customer bills and realize system benefit The fielded assets deployed for the project were: Utility Scale Kingsbery Energy Storage System: 1.5 MW / 3 MWh Li-Ion battery storage Mueller Energy Storage System: 1.75 MW / 3.2 MWh Li-Ion battery storage, 7 Energy Storage Units (250 kW each) La Loma Community Solar: 2.6 MW Commercial Scale Aggregated storage installations at 3 sites, with existing solar (300+ kW): One 18 kW / 36 kWh Li-Ion battery storage Two 72 kW / 144 kWh Li-Ion battery storage Residential Scale Aggregated storage installations: -Six stationary battery storage systems (10 kWh each) at homes with existing solar -One Electric Vehicle installed as Vehicle-to-Grid (V2G) Utility-Controlled Solar via Smart Inverters at 12 homes Autonomously-Controlled Smart Inverters at 6 homes Over the course of the project, Austin SHINES undertook installing more than 3 MW of distributed battery energy storage, smart PV inverters, a DER control platform, and other enabling technologies utilizing customer and utility locations and aggregation models. All of these resources were to be integrated and optimized at the utility level. DER assets and control methodologies were designed to achieve a credible pathway to a System LCOE for energy delivered to load of $0.14//kWh or less by 2020, while maximizing distributed solar generation and maintaining acceptable standards of power quality. The project also established a template for other regions to follow, to maximize the adoption of distributed solar PV in support of an economic and efficient grid. In total, the Austin SHINES project added value to the DER subject area in each layer of integration. From utility, to commercial to residential scales, the sheer hierarchy of communication and coordination was a significant accomplishment in addition to learnings from what these communications revealed was unique to each. Economically, the most effective method demonstrated was the criticality of planning phases. Contingencies and multiple projection scenarios helped guide the project to deploy optimal design as close as feasible, in real world conditions. The project and reports will serve public benefit by outlining specific areas of DER strategy and installation where many stakeholders and needs can be addressed with improved efficiency. Overall, communities and utilities should use the results to guide the increasing options available for powering the grid with DER, renewables, and carbon considerate energy.
In this report we explore tradeoffs of top-down and bottom-up methods in their precision and computational burden using NREL's dGen model, an agent-based model of residential and non-residential distributed PV adoption. In particular, we assess the role of agent resolution in instantiating statistically-representative populations in the model--and the resulting variance of model projections at the state, sector and county levels. Running the model for a single state with fixed macroeconomic assumptions, we vary both the number of agents instantiated and the number of unique simulations. Notably, in each run we introduce sources of uncertainty in cumulative DPV adoption results during agent instantiation as we stochastically sample agent attributes from localized building characteristic and solar resource distributions. We find three key trends which guide what can be expected of results. First, precision of the model projections, as measured by the variance of the unique simulations increased as a function of the number of agents sampled, though with diminishing returns to scale. The result is expected given the Law of Large Numbers. Second, the mean or expected value of projections decreased as the number of agents instantiated increased. Empirically, aggregate DPV capacities projected by dGen are 7% lower on average at a 10-agent resolution compared to 1-agent resolutions, and up to 10% lower in the non-residential sector. We attribute this outcome to a system sizing scheme that biases lower cumulative installed capacity as agent resolution increases when roof size and load are not positively correlated. Finally, the variance of model projections differed by sector, where variance was substantially larger in the commercial sector, largely due to differences in parameterization by sector. We explain the larger variance in the non-residential sector as an outcome of the greater variance in building uses and load shapes as compared to the residential sector, and we find that a resolution on the order of 100-agents per county would be necessary to bring non-residential variance in line with residential variance at a 2-agent resolution. We propose that further research explore increasing precision and agent resolution in the non-residential sector. Moreover, given that days of intensive computing may be required to resolve a suitable estimate even at the extent of a single state within acceptable levels of precision, to promote computational efficiency we also suggest more clearly quantifying sources of variability at the county level such that resolution can by dynamically increased where county level results are least precise.
The dynamics of spin projections of the electron shell of an alkali metal on the coordinate axis is considered in the electron paramagnetic resonance scheme with continuous pumping by biharmonic circularly polarized laser radiation. The working region is a cell with alkali vapor metal vapors and a buffer gas at a high concentration at temperature 60°C. It was found that the use of biharmonic pumping causes not only the expected electron-spin precession, but also pulsations of the electron-spin projection on the axis along which the magnetic field is directed. The frequency of these pulsations depends on the nuclear angular momentum of alkali metal atoms. In the case of the transverse electron magnetic resonance, this effect is absent.
The work performed in this project consisted of the derivation, implementation, and testing of a new, computationally advantageous approximation to the 3D Boltzmann transport equation. The solution of the Boltzmann equation is the neutron flux in nuclear reactor cores and shields, but solving this equation is difficult and costly. The new “2D/1D” approximation takes advantage of a special geometric feature of typical 3D reactors to approximate the neutron transport physics in a specific (ax- ial) direction, but not in the other two (radial) directions. The resulting equation is much less expensive to solve computationally, and its solutions are expected to be sufficiently accurate for many practical problems. In this project we formulated the new equation, discretized it using standard methods, developed a stable iteration scheme for solving the equation, implemented the new numerical scheme in the MPACT code, and tested the method on several realistic problems. All the hoped- for features of this new approximation were seen. For large, difficult problems, the resulting 2D/1D solution is highly accurate, and is calculated about 100 times faster than a 3D discrete ordinates simulation.
We extend the Virtual Element Method to a two-dimensional unsteady nonlinear convection-diffusion equation characterized by a fractional-order derivative with respect to the time variable. Our methodology is based on three fundamental technical components: a fractional version of the Grunwald-Letnikov approximation, discrete maximal regularity, and the regularity theory associated with non-linearity. We prove the method's well-posedness, i.e., the approximate solution's existence and uniqueness to the time-fractional convection-diffusion equation with a Lipschitz nonlinear source term. The fully discrete scheme inherently maintains stability and consistency by leveraging the discrete maximal regularity and the energy projection operator. The convergence in the L 2 -norm and H 1 -norm to various mesh configurations is validated by numerical results, underlining the practical effectiveness of the proposed method.
For this project, a version of the GLAS Analysis/Forecast System was developed that includes an objective dealiasing scheme as an integral part of the analysis cycle. With this system the (100 sq km) binned SASS wind data generated by S. Peteherych of AER, Canada corresponding of the period 0000 GMT 7 September 1978 to 1200 GMT 13 September 1978 was objectively dealiased. The dealiased wind fields have been requested and received by JPL, NMC and the British Meteorological Office. The first 3.5 days of objectively dealiased fields were subjectively enhanced on the McIDAS system. Approximately 20% of the wind directions were modified, and of these, about 70% were changed by less than 90 deg. Two SASS forecast impact studies, were performed using the dealiased fields, with the GLAS and the NEPRF (Navy Environmental Prediction Research Facility) analysis/forecast systems.
Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) observations over mesoscale convective systems (MCSs) reveal that there are localized maxima in the rain rate with a scale of about 10 to 20 km that represent thunderstorms (Cbs). Some of these Cbs are developing or intense, while others are decaying or weak. These Cbs constitute only about 20 % of the rain area of a given MCS. Outside of Cbs, the average rain rate is much weaker than that within Cbs. From an analysis of the PR data, we find that the spatial distribution of rain and its character, convective or stratiform, is highly inhomogeneous. This complex nature of rain exists on a scale comparable to that of a Cb. The 85 GHz brightness temperature, T85, observations of the TRMM Microwave Imager (TMI) radiometer taken over an MCS reflect closely the PR rain rate pattern over land. Local maxima in rain rate shown by PR are observed as local minima in T85. Where there are no minima in T85, PR observations indicate there is light rain. However, the TMI brightness temperature measurements (Tbs) have poor ability to discriminate convective rain from stratiform rain. For this reason, a TMI rain retrieval procedure that depends primarily on the magnitude of Tbs performs poorly. In order to retrieve rain rate from TMI data on land one has to include the spatial distribution information deduced from the T85 data in the retrieval method. Then, quantitative estimation of rain rate can be accomplished. A TMI rain retrieval method developed along these lines can yield estimates of rain rate and its frequency distribution which agree closely with that given by PR. We find the current TRMM project TMI (Version 5) rain retrieval algorithm on land could be improved with the retrieval scheme developed here. To support the conceptual frame work of the rain retrieval method developed here, a theoretical analysis of the TMI brightness temperatures in convective and stratiform regions is presented.
To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.
The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbo-pump geometry through numerical simulation will be of significant value toward design. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the Message Passage Interface (MPI) and Multi Level Parallel (MLP) versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology, INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbo-pump simulations, moving boundary capability and efficient time-accurate integration methods are built in the flow solver, To handle the geometric complexity and moving boundary problems, an overset grid scheme is incorporated with the solver so that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two methods is compared by obtaining unsteady solutions for the evolution of twin vortices behind a flat plate. Calculated results are compared with experimental and other numerical results. For an unsteady flow, which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiteration procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy the incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in present computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.