Engineering Papers⌕ Search

Engineering topics

Wu, S. T.

Publications and source records attributed to Wu, S. T..

At least 73 records · Page 4

Study of basic physical processes in liquid rocket engines

Inconsistencies between analytical results and measurements for liquid rocket thrust chamber performance, which escape suitable explanations, have motivated the examination of the basic phys ical modeling formulations as to their unlimited application. The publication of Prof. D. Straub's book, 'Thermofluid-dynamics of Optimized Rocket Propulsions,' further stimulated the interest of understanding the gas dynamic relationships in chemically reacting mixtures. A review of other concepts proposed by Falk-Ruppel (Gibbsian Thermodynamics), Straub (Alternative Theory, AT), Prigogine (Non-Equilibrium Thermodynamics), Boltzmann (Kinetic Theory), and Truesdell (Rational Mechanism) has been made to obtain a better understanding of the Navier-Stokes equation, which is now used extensively for chemically reacting flow treatment in combustion chambers. In addition to the study of the different concepts, two workshops were conducted to clarify some of the issues. The first workshop centered on Falk-Ruppel's new 'dynamics' concept, while the second one concentrated on Straub's AT. In this report brief summaries of the reviewed philosophies are presented and compared with the classical Navier-Stokes formulation in a tabular arrangement. Also the highlights of both workshops are addressed.

Wu, S. T.↗

Three-dimensional numerical simulation of interplanetary magnetic field changes at 1 AU as a consequence of simulated solar flares

In order to study the relationship between different forms of activity and transient variations of the north-south component, B(sub z), of the Interplanetary Magnetic Field (IMF), at 1 AU, a three dimensional numerical simulation code was employed to study several aspects of this problem. One form of solar activity, the flare, is simulated by using a pressure pulse at different locations near the solar surface and observing the simulated IMF evolution of B(sub theta)(-B(sub z)) at 1 AU. It was found, for a given pressure pulse, that the orientation of the corresponding transient variation of B(sub z) has a strong relationship with the location of the pressure pulse and the initial condition of IMF.

Wu, S. T.↗

3-dimensional current collection model

A three-dimensional, time dependent current collection model of a satellite has been developed for the TSS-1 system. The system has been simulated particularly for the Research of Plasma Electrodynamics (ROPE) experiment. The Maxwellian distributed particles with the geomagnetic field effects are applied in this numerical simulation. The preliminary results indicate that a ring current is observed surrounding the satellite in the equatorial plane. This ring current is found between the plasma sheath and the satellite surface and is oscillating with a time scale of approximately 1 microsec. This is equivalent to the electron plasma frequency. An hour glass shape of electron distribution was observed when the viewing direction is perpendicular to the equatorial plane. This result is consistent with previous findings from Linson (1969) and Antoniades et al. (1990). Electrons that are absorbed by the satellite are limited from the background ionosphere as indicated by Parker and Murphy (1967).

Hwang, Kai-Shen↗

Development of generalized pressure velocity coupling scheme for the analysis of compressible and incompressible combusting flows

The objective of this investigation has been to develop an algorithm (or algorithms) for the improvement of the accuracy and efficiency of the computer fluid dynamics (CFD) models to study the fundamental physics of combustion chamber flows, which are necessary ultimately for the design of propulsion systems such as SSME and STME. During this three year study (May 19, 1978 - May 18, 1992), a unique algorithm was developed for all speed flows. This newly developed algorithm basically consists of two pressure-based algorithms (i.e. PISOC and MFICE). This PISOC is a non-iterative scheme and the FICE is an iterative scheme where PISOC has the characteristic advantages on low and high speed flows and the modified FICE has shown its efficiency and accuracy to compute the flows in the transonic region. A new algorithm is born from a combination of these two algorithms. This newly developed algorithm has general application in both time-accurate and steady state flows, and also was tested extensively for various flow conditions, such as turbulent flows, chemically reacting flows, and multiphase flows.

Chen, C. P.↗

Forecasting the arrival of fast coronal mass ejecta at Earth by the detection of 2-20 keV neutral atoms

Studies have shown that Earth passages of fast coronal mass ejections (CMEs) trigger geomagnetic storms. Early identification of fast earth-directed CME can help provide storm warnings, but detection of such by coronagraphs is extremely difficult. We suggest that energetic hydrogen atoms (EHA) between 2 and 10 keV produced during the transit phase of an Earth-directed CME by recombination between protons and electrons in the CME can travel ahead of the CME and act as harbingers of a magnetic storm. This forecasting scheme should work if enough EHA are produced, because while CMEs decelerate continuously after their ejection, the EHA fluxes produced in the initial phase of fast CMEs propagate at their initial high speeds. Model simulations support this proposed mechanism.

Hsieh, K. C.↗

A two-dimensional MHD global coronal model - Steady-state streamers

A 2D, time-dependent, numerical, MHD model for the simulation of coronal streamers from the solar surface to 15 solar is presented. Three examples are given; for dipole, quadrupole and hexapole (Legendre polynomials P1, P2, and P3) initial field topologies. The computed properties are density, temperature, velocity, and magnetic field. The calculation is set up as an initial-boundary value problem wherein a relaxation in time produces the steady state solution. In addition to the properties of the solutions, their accuracy is discussed. Besides solutions for dipole, quadrupole, and hexapole geometries, the model use of realistic values for the density and Alfven speed while still meeting the requirement that the flow speed be super-Alfvenic at the outer boundary by extending the outer boundary to 15 solar radii.

Wang, A.-H.↗

Shear-induced instability and arch filament eruption - A magnetohydrodynamic (MHD) numerical simulation

A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.

Wu, S. T.↗

Modified ICED-ALE method for astrogeophysical plasma flows

The Implicit-Continuous-Eulerian-Difference Mesh-Arbitrary-Lagrangian-Eulerian (ICED-ALE) algorithm of Brackbill and Pracht (1973) is modified for the study of astrophysical plasma flows in which dynamical effects are important. In the present study the general-energy-conservation law is directly applied to the iteration process, with the total (kinetic-, specific-internal-, and magnetic-) energy density being obtained implicitly at the end of the process. An example is computed in which the convergence speed of the latter method is substantially faster than that of the former. The initializing calculation, or explicit phase, in which the electric current density, magnetic diffusion of field, energy augmentation, and zero-order approximation of the flow velocity are given, is described. Consideration is given to the iteration process, or implicit phase, from which the exact Lagrangian solution for energy density, velocity, and a magnetic field is obtained.

Wu, S. T.↗

A study of physical mechanisms for filament eruption and coronal mass ejection via numerical simulation

It is well-known from both observation and theory that photospheric shear motion has played a key role in physical processes of the energy build-up and release for the solar flare. In order to further our understanding of the occurrence of solar flares it is necessary to investigate the triggering mechanism. One popular scenario for the onset of eruptive solar flares is that in response to photospheric shear motions the magnetic field evolves slowly through a series of magnetohydrodynamic-equilibria until a threshold is reached where magnetohydrodynamic (MHD) non-equilibrium sets in. Thus, a magnetic eruption occurs, causing the solar flare. To substantiate this claim we have employed our newly developed three-dimensional time-dependent MHD code with gravity to simulate the evolution of the coronal field. We use plasma beta = 0.1 to closely approximate the condition in the actual corona. Some preliminary results are presented.

Wu, S. T.↗

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.↗

Data evaluation, analysis, and scientific study

Extensive work was performed in data analysis and modeling of solar active phenomena. The work consisted in the study of UV data from the Ultraviolet Spectrometer and Polarimeter (UVSP) instrument on board the Solar Maximum Mission satellite. These data were studied in conjunction with X-rays from the Hard X-ray Imaging Spectrometer (HXIS) instrument, and with H-alpha and magnetographic data from ground-based observatories. The processes we studied are the active phenomena which result from the interaction of the solar magnetic fields with the plasma in the outer regions of the solar atmosphere. These processes include some very dynamic processes such as the prominence eruptions and the 'microflares'. Our research aimed at characterizing the following: the observed phenomena, the possible physical models, and the relevance to the chromospheric and coronal heating.

Wu, S. T.↗

A signal preprocessing algorithm applicable to both DTF and SSME spectral diagnostics

An account is given of research conducted by NASA-Stennis' Diagnostic Test Facility (DTF), which is designed to obtain spectra approximating those of the SSME while allowing full control of dopant and thruster conditions. A simple preprocessing technique which emphasizes the essential spectral components of plume spectra that are indicative of dopant levels has been used and found to remain sensitive to wavelength shift even for high-resolution spectra. Mathematical analyses of the problem are presented which are based on the Taylor's series linearization procedure.

Nail, Bert↗

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗

The role of condensation and heat conduction in the formation of prominences - An MHD simulation

The effects of condensation and thermal conduction on the formation of Kippenhahn-Schlueter (K-S) type prominences in quiet regions (QP) due to symmetric mass injection are investigated. To implement this investigation a self-consistent, two-dimensional, nonplanar, time-dependent MHD simulation model is developed. In the model, various values of the injection velocity, density, and magnetic field strength are used to determine the most favorable conditions for the QP formation. Based on these simulation results, it is found that the formation of a K-S type field configuration should be considered as a dynamic process which needs both condensation amd mass injection to supply enough mass to maintain such a configuration to complete the quiescent prominence formation process.

Wu, S. T.↗

Current collection model on tether satellite system

Results are reported from numerical simulations of current collection by a tethered satellite moving in the earth magnetic field. The focus is on the instruments of the Research on Orbital Plasma Electrodynamics (ROPE) experiment planned for the TSS-1 mission (scheduled launch 1991). The derivation of the governing equations is outlined, and the results of simulations with and without the geomagnetic field are presented in graphs. It is predicted that beam-beam interaction will occur when no magnetic effect is present, especially when the instrument bias is lower than the local potential; with the magnetic effect, there should be complex multiple-stream interactions, magnetic-field shielding covering 30 percent or more of the satellite surface when the B field is greater than 0.35 G, and nearly complete collection of the electrons by the instrument arm.

Hwang, K. S.↗

A comparison between progressive extension method (PEM) and iterative method (IM) for magnetic field extrapolations in the solar atmosphere

This paper presents a comparison between two numerical methods for the extrapolation of nonlinear force-free magnetic fields, viz the Iterative Method (IM) and the Progressive Extension Method (PEM). The advantages and disadvantages of these two methods are summarized, and the accuracy and numerical instability are discussed. On the basis of this investigation, it is claimed that the two methods do resemble each other qualitatively.

Wu, S. T.↗

Signal processing of spectroscopic data for engine diagnostic applications

A description is given of research based on a series of spectra taken from the Diagnostic Test Facility (DTF) plume at NASA-Stennis Space Center. The DTF was designed to obtain spectra scaling approximately to that of the SSME, yet allow full control of dopant and thruster conditions. DTF spectra were used in this research to develop a technique for quantitative identification of an alloy. A mathematically based identification technique was tried, with reasonable sucess for one-element alloys, and only fair success for a two-element alloy. This approach was found to have the simplicity required for DSP implementation, thus allowing real-time diagnosis as a possibility. Spreadsheets were organized to handle the proliferation of development equations required, which number about 50,000. The identification work reinforced the presence of nonlinearity to dopant level.

Nail, Bert↗

Numerical simulations of solar disturbances and their interplanetary consequences

Time-dependent MHD numerical simulations are used to study responses of the solar atmosphere and interplanetary medium to simulated solar disturbances. A number of 2D and 3D examples of coronal mass ejection (CME) simulations and some current controversies concerning the basic processes of CME initiation are discussed. Footpoint shearing motion is tested to determine whether it can provide a reasonable mechanism for CME development from arch filament configurations. Possible interplanetary consequences to CME-like disturbances are demonstrated by using 3D simulations to determine the dynamic response of the solar wind to a plasmoid injection from an eruptive filament or prominence. The possibility that a plasmoid may be generated in the interplanetary medium by a solar-generated shock that propagates through a heliospheric current sheet is discussed. Application of the 3D model for the interpretation of interplanetary scintillation observations is addressed.

Dryer, M.↗