Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization 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.

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

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Quantitative Measurements of Nitric Oxide Concentration in High-Pressure, Swirl-Stabilized Spray Flames

Lean direct-injection (LDI) spray flames offer the possibility of reducing NO(sub x) emissions from gas turbines by rapid mixing of the liquid fuel and air so as to drive the flame structure toward partially-premixed conditions. We consider the technical approaches required to utilize laser-induced fluorescence methods for quantitatively measuring NO concentrations in high-pressure LDI spray flames. In the progression from atmospheric to high-pressure measurements, the LIF method requires a shift from the saturated to the linear regime of fluorescence measurements. As such, we discuss quantitative, spatially resolved laser-saturated fluorescence (LSF), linear laser-induced fluorescence (LIF), and planar laser-induced fluorescence (PLIF) measurements of NO concentration in LDI spray flames. Spatially-resolved LIF measurements of NO concentration (ppm) are reported for preheated, LDI spray flames at pressures of two to five atmospheres. The spray is produced by a hollow-cone, pressure-atomized nozzle supplied with liquid heptane. NO is excited via the Q(sub 2)(26.5) transition of the gamma(0,0) band. Detection is performed in a two nanometer region centered on the gamma(0,1) band. A complete scheme is developed by which quantitative NO concentrations in high-pressure LDI spray flames can be measured by applying linear LIF. NO is doped into the reactants and convected through the flame with no apparent destruction, thus allowing a NO fluorescence calibration to be taken inside the flame environment. The in-situ calibration scheme is validated by comparisons to a reference flame. Quantitative NO profiles are presented and analyzed so as to better understand the operation of lean-direct injectors for gas turbine combustors. Moreover, parametric studies are provided for variations in pressure, air-preheat temperature, and equivalence ratio. Similar parametric studies are performed for lean, premixed-prevaporized flames to permit comparisons to those for LDI flames. Finally, PLIF is expanded to high pressure in an effort to quantify the detected fluorescence image for LDI flames. Success is achieved by correcting the PLIF calibration via a single-point LIF measurement. This procedure removes the influence of any preferential background that occurs in the PLIF detection window. In general, both the LIF and PLIF measurements verify that the LDI strategy could be used to reduce NO(sub x) emissions in future gas turbine combustors.

Cooper, Clayton S.↗
DOE OSTI2019

Direct and indirect measurements of the magnetic and magnetocaloric properties of Ni 0.895 Cr 0.105 MnGe 1.05 melt-spun ribbons in high magnetic fields

We report the magnetic and magnetocaloric properties of rapidly solidified Ni 0.895 Cr 0.105 MnGe 1.05 melt-spun ribbons studied by both direct (adiabatic temperature change) and indirect (isothermal magnetic entropy change) methods in intermediate and high magnetic fields up to 10 T. The maximum values of the adiabatic temperature changes (ΔT ad ) and magnetic entropy changes ( Δ S M ) were found to be ~2.6 K (µ o H = 10 T) and 4.4 J kg -1 K -1 (µ o ΔH = 5 T), respectively, near the Curie temperature (T C ). The ΔT ad curves and magnetization isotherms were found to be completely reversible, which indicates the high degree of reversibility of the MCEs in this system. A large temperature span (of about 61 K) and a non-saturating behavior of ΔT ad were observed at magnetic fields up to 10 T. The adiabatic temperature change was found to be a linear function of (µ o H) 2/3 near T C , in accordance with the Landau theory of phase transitions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗
NASA NTRS1959

Supersonic and Moment-of-Area Rules Combined for Rapid Zero-Lift Wave-Drag Calculations

The concepts of the supersonic area rule and the moment-of-area rule are combined to develop a new method for calculating zero-lift wave drag which is amenable to the use of ordinary desk calculators. The total zero-lift wave drag of a configuration is calculated by the new method as the sum of the wave drag of each component alone plus the interference between components. In calculating the separate contributions each component or pair of components is analyzed over the smallest allowable length in order to improve the convergence of the series expression for the wave drag. The accuracy of the present method is evaluated by comparing the total zero-lift wave-drag solutions for several simplified configurations obtained by the present method with solutions given by slender-body and linearized theory. The accuracy and computational time required by the present method are also evaluated relative to the supersonic area rule and the moment-of-area rule. The results of the evaluation indicate that total zero-lift wave-drag solutions for simplified configurations can be obtained by the present method which differ from solutions given by slender-body and linearized theory by less than 6 percent. This accuracy for simplified configurations was obtained from only nine terms of the series expression for the wave drag as a result of calculating the total zero-lift wave drag by parts. For the same number of terms these results represent an accuracy greater than that for solutions obtained by either of the two methods upon which the present method is based, except in a few isolated cases. For the excepted cases, solutions by the present method and the supersonic area rule are identical. Solutions by the present method are obtained in one fifth the computing time required by the supersonic area rule. This difference in computing time of course would be substantially reduced if the complete procedures for both methods were programmed on electronic computing machines.

Levy, Lionel L., Jr.↗
DOE OSTI2022

The effects of 8-week complex training on lower-limb strength and power of Chinese elite female modern pentathlon athletes

Complex training (CT) is a combination training method that alternates between performing high-load resistance training (RT) and plyometric training within one single session. The study aimed to examine the effects of CT on lower-limb strength and power of elite female modern pentathlon athletes under the new modern pentathlon format and competition rules. Ten female participants (age: 23.55 ± 2.22 years, weight: 60.59 ± 3.87 kg, height: 169.44 ± 4.57 cm, and training experience: 6.90 ± 2.08 years) of the national modern pentathlon team completed 8 weeks of RT as followed by 8 weeks of CT, with 2 weeks of break. Then, the participants conducted 8 weeks of CT, which included RT combined with plyometric training (e.g., drop jump and continuous jump). All stages of training were designed by the linear strength training period theories, requiring participants to train twice for the first 4 weeks and three times for the second 4 weeks. The one-repetition maximum (1RM) of squat, isometric mid-thigh pull (IMTP), counter-movement jump (CMJ), squat jump (SJ), pre-stretch augmentation percentage (PSAP), and reaction strength index (RSI) were assessed before and after both RT and CT training. One-way repeated-measure ANOVA models revealed that the 1RM of squat was significantly improved ( p < 0.001) after RT as compared to pre-RT. No significant improvement in IMTP ( p = 0.055), CMJ ( p = 0.194), SJ ( p = 0.692), PSAP ( p = 0.087), and RSI ( p = 0.238) was not observed. After CT, 1RM of squat ( p < 0.001), IMTP ( p < 0.035), CMJ ( p < 0.001), SJ ( p < 0.008), RSI ( p < 0.006) were significant improved as compared to pre-RT, post-RT and pre-CT, while significant improvements in PSAP were observed as compared to pre-RT ( p = 0.003) and pre-CT ( p = 0.027), but not to post-RT ( p = 0.156). This pilot study showed the promise of CT following RT to improve lower-limb strength and power in elite female modern pentathlon athletes. The findings are worthwhile to be confirmed in future studies with larger sample size and randomized design.

Qiao, Zining↗
NASA NTRS1990

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 1: Theory document (version 3.0)

An outline of the derivation of the differential equation governing linear subsonic and supersonic potential flow is given. The use of Green's Theorem to obtain an integral equation over the boundary surface is discussed. The engineering techniques incorporated in the Panel Aerodynamics (PAN AIR) program (a discretization method which solves the integral equation for arbitrary first order boundary conditions) are then discussed in detail. Items discussed include the construction of the compressibility transformation, splining techniques, imposition of the boundary conditions, influence coefficient computation (including the concept of the finite part of an integral), computation of pressure coefficients, and computation of forces and moments. Principal revisions to version 3.0 are the following: (1) appendices H and K more fully describe the Aerodynamic Influence Coefficient (AIC) construction; (2) appendix L now provides a complete description of the AIC solution process; (3) appendix P is new and discusses the theory for the new FDP module (which calculates streamlines and offbody points); and (4) numerous small corrections and revisions reflecting the MAG module rewrite.

Epton, Michael A.↗
DOE OSTI2022

Initial Development of Variable Design Lifetimes and Creep-Fatigue Evaluations for ASME Section III, Division 5, Class B Code Rules

This report summarizes the initial development of new ASME Section III, Division 5, Class B rules to address the gaps identified for new high temperature reactor designs. The overall objective is to introduce the design-by-analysis approach, the variable design lifetimes for primary load design, and the strain limits and creep-fatigue damage evaluations for Class B components. It is also desired that the new evaluation procedures do not require stress classification and stress linearization. Further, the data requirements to support the development of design parameters for the new Class B rules would not be overly burdensome, and less extensive than those supporting Class A rules. We have identified the Elastic, Perfectly Plastic (EPP) method as the design by-analysis approach for primary load design, and developed allowable stress criteria for variable design lifetimes, up to 500,000 hours, that only require data typical for the allowable stresses in Section II, Part D, Tables 1A and 1B. We have also adapted the elastic analysis approach from HBB-T-1400 and Section III, Mandatory Appendix XIII, XIII-2400 and XIII-3520(e) for the new Class B creep-fatigue evaluation procedure. A new elastic follow-up-based Isochronous Stress Strain Curve (ISSC) stress relaxation procedure has been developed to provide adequate conservatism commensurate with Class B constructions. No stress classification and stress linearization are required in the new procedure. Future work needed to complete the development is also summarized.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Subcontract B635525 (Final Report)

During my visits to LLNL July 8–12 and July 22–26, 2019, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $\textit{hp}$ adaptive finite element calculations, and is implemented in the $\textit{PLTMG}$ software package, version 13. This preconditioner typically requires 3–20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $\textit{p}$ as $\textit{p} → ∞$, but remains strong for $\textit{p}$ ≤ 9, which is the highest polynomial degree allowed in $\textit{PLTMG}$, due to limitations of the numerical quadrature rules implemented in the soft- ware package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

3D spherical functional expansion tallies in Serpent 2 Monte Carlo code

This work extends the application of functional expansion tallies to 3D spherical geometries. The 3D Zernike polynomials are set as an orthonormal polynomials basis for the functional reconstruction. The study describes the construction of the complete set of polynomials, a natural expansion of the spherical harmonics polynomials where 3D Zernike moments can be evaluated as a linear combination of the geometrical moments. The 3D Zernike polynomials formulation and the computational approach implemented in Serpent 2 are presented and tested through the Godiva model from the ICSBEP criticality benchmark test cases. The implementation results are in agreement with a reference solution described in a fine-resolution mesh, enhancing also the performance and memory demand. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2014

Missing Data and Multiple Imputation: An Unbiased Approach

The default method of dealing with missing data in statistical analyses is to only use the complete observations (complete case analysis), which can lead to unexpected bias when data do not meet the assumption of missing completely at random (MCAR). For the assumption of MCAR to be met, missingness cannot be related to either the observed or unobserved variables. A less stringent assumption, missing at random (MAR), requires that missingness not be associated with the value of the missing variable itself, but can be associated with the other observed variables. When data are truly MAR as opposed to MCAR, the default complete case analysis method can lead to biased results. There are statistical options available to adjust for data that are MAR, including multiple imputation (MI) which is consistent and efficient at estimating effects. Multiple imputation uses informing variables to determine statistical distributions for each piece of missing data. Then multiple datasets are created by randomly drawing on the distributions for each piece of missing data. Since MI is efficient, only a limited number, usually less than 20, of imputed datasets are required to get stable estimates. Each imputed dataset is analyzed using standard statistical techniques, and then results are combined to get overall estimates of effect. A simulation study will be demonstrated to show the results of using the default complete case analysis, and MI in a linear regression of MCAR and MAR simulated data. Further, MI was successfully applied to the association study of CO2 levels and headaches when initial analysis showed there may be an underlying association between missing CO2 levels and reported headaches. Through MI, we were able to show that there is a strong association between average CO2 levels and the risk of headaches. Each unit increase in CO2 (mmHg) resulted in a doubling in the odds of reported headaches.

Foy, M.↗
NASA NTRS2014

Variable Acceleration Force Calibration System (VACS)

Conventionally, force balances have been calibrated manually, using a complex system of free hanging precision weights, bell cranks, and/or other mechanical components. Conventional methods may provide sufficient accuracy in some instances, but are often quite complex and labor-intensive, requiring three to four man-weeks to complete each full calibration. To ensure accuracy, gravity-based loading is typically utilized. However, this often causes difficulty when applying loads in three simultaneous, orthogonal axes. A complex system of levers, cranks, and cables must be used, introducing increased sources of systematic error, and significantly increasing the time and labor intensity required to complete the calibration. One aspect of the VACS is a method wherein the mass utilized for calibration is held constant, and the acceleration is changed to thereby generate relatively large forces with relatively small test masses. Multiple forces can be applied to a force balance without changing the test mass, and dynamic forces can be applied by rotation or oscillating acceleration. If rotational motion is utilized, a mass is rigidly attached to a force balance, and the mass is exposed to a rotational field. A large force can be applied by utilizing a large rotational velocity. A centrifuge or rotating table can be used to create the rotational field, and fixtures can be utilized to position the force balance. The acceleration may also be linear. For example, a table that moves linearly and accelerates in a sinusoidal manner may also be utilized. The test mass does not have to move in a path that is parallel to the ground, and no re-leveling is therefore required. Balance deflection corrections may be applied passively by monitoring the orientation of the force balance with a three-axis accelerometer package. Deflections are measured during each test run, and adjustments with respect to the true applied load can be made during the post-processing stage. This paper will present the development and testing of the VASC concept.

Rhew, Ray D.↗
NASA NTRS2013

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗
NASA NTRS1989

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

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

Lefkowitz, C. P.↗
NASA NTRS2013

Covariance Analysis Tool (G-CAT) for Computing Ascent, Descent, and Landing Errors

G-CAT is a covariance analysis tool that enables fast and accurate computation of error ellipses for descent, landing, ascent, and rendezvous scenarios, and quantifies knowledge error contributions needed for error budgeting purposes. Because GCAT supports hardware/system trade studies in spacecraft and mission design, it is useful in both early and late mission/ proposal phases where Monte Carlo simulation capability is not mature, Monte Carlo simulation takes too long to run, and/or there is a need to perform multiple parametric system design trades that would require an unwieldy number of Monte Carlo runs. G-CAT is formulated as a variable-order square-root linearized Kalman filter (LKF), typically using over 120 filter states. An important property of G-CAT is that it is based on a 6-DOF (degrees of freedom) formulation that completely captures the combined effects of both attitude and translation errors on the propagated trajectories. This ensures its accuracy for guidance, navigation, and control (GN&C) analysis. G-CAT provides the desired fast turnaround analysis needed for error budgeting in support of mission concept formulations, design trade studies, and proposal development efforts. The main usefulness of a covariance analysis tool such as G-CAT is its ability to calculate the performance envelope directly from a single run. This is in sharp contrast to running thousands of simulations to obtain similar information using Monte Carlo methods. It does this by propagating the "statistics" of the overall design, rather than simulating individual trajectories. G-CAT supports applications to lunar, planetary, and small body missions. It characterizes onboard knowledge propagation errors associated with inertial measurement unit (IMU) errors (gyro and accelerometer), gravity errors/dispersions (spherical harmonics, masscons), and radar errors (multiple altimeter beams, multiple Doppler velocimeter beams). G-CAT is a standalone MATLAB- based tool intended to run on any engineer's desktop computer.

Boussalis, Dhemetrios↗
NASA NTRS1983

Effect of boundary layers on solid walls in three-dimensional subsonic wind tunnels

A linear method is developed which accounts for the effects of boundary layers on solid walls in subsonic three-dimensional wind tunnels. The streamwise gradient of the displacement thickness for a solid-wall boundary layer is expressed in terms of the von Karman momentum integral. The growth of the boundary layer due to the wall shearing stress is small compared to the variation caused by the model-induced pressure gradient. The viscous boundary condition can be expressed in terms of the edge velocity gradient and the gradient of the inviscid velocity potential function at the wall. Utilizing this analysis on the solid walls of several three-dimensional wind tunnel configurations shows that the most pronounced wall boundary-layer effect is on solid blockage for completely closed wind tunnels. For solid-wall tunnel configurations, the streamline curvature interference factor is reduced by a significant amount, while the lift interference factor at the model station does not depend on the boundary-layer parameter. For combination wall configurations, the slot effect of the horizontal walls dominates the viscous effect of the solid sidewalls.

Adcock, J. B.↗
DOE OSTI2019

Earthquake Arrival Association with Backprojection and Graph Theory

The association of seismic-wave arrivals with causative earthquakes becomes progressively more challenging as arrival detection methods become more sensitive, and particularly when earthquake rates are high. For instance, seismic waves arriving across a monitoring network from several sources may overlap in time, false arrivals may be detected, and some arrivals may be of unknown phase (e.g., P or S waves). We propose an automated method to associate arrivals with earthquake sources and obtain source locations applicable to such situations. To do so, we use a pattern detection metric based on the principle of back-projection to reveal candidate sources followed by graph-theory-based clustering and an integer linear optimization routine to associate arrivals with the minimum number of sources necessary to explain the data. We find that this method solves for all sources and phase assignments simultaneously, rather than in a sequential greedy procedure as is common in other association routines. We demonstrate our method on both synthetic and real data from the Integrated Plate Boundary Observatory Chile seismic network of northern Chile. For the synthetic tests, we report results for cases with varying complexity, including rates of 500 earthquakes/day and 500 false arrivals/station/day, for which we measure true positive detection accuracy of >95%. For the real data, we develop a new catalog between 1 January 2010 and 31 December 2017 containing 817,548 earthquakes, with detection rates on average 279 earthquakes/day and a magnitude-of-completion of M~1.8. A subset of detections are identified as sources related to quarry and industrial site activity, and we also detect thousands of foreshocks and aftershocks of the 1 April 2014 Mw 8.2 Iquique earthquake. During the highest rate of aftershock activity, >600 earthquakes/day are detected in the vicinity of the Iquique earthquake rupture zone.

58 GEOSCIENCES↗
NASA NTRS1961

Methods for Determining Blood Flow Through Intact Vessels of Experimental Animals Under Conditions of Gravitational Stress and in Extraterrestrial Space Capsules

Work during the initial six months of this program has proceeded along three lines: (1) technical development of blood flowmeter instrumentation: (2) surgical considerations related to special problems of chronic implantation of flowmeter probes: (3) study of the relationship between levels of blood flow and organ activity. Principal engineering effort has been focused on the fundamental phenomena operative in the flowmeter probe and on probe design. Two probes have been constructed and a third is nearing completion. A previously unrecognized source of zero drift has been discovered and corrected. has been designed, assembled, tested and successfully employed in work on the several probes. A significant innovation which has been introduced is the use of a transformer coupled input for suppression of unwanted common mode signals. New low noise transistor amplifiers with characteristics suitable for space vehicle packaging are being developed. Linearity of the present apparatus has been found to extend at least into the near turbulent flow region. Attention has been directed to the two vascular beds of physiologic interest that present particular problems with regard to chronic implantation of flowmeter probes, the myocardium and the liver. Surgical techniques have been perfected that will permit myocardial and liver blood flow measurement. Two groups of animals have successfully undergone surgery and await the implantation of probes. Studies of the relationship between organ blood flow and levels of physiologic activity have been initiated. The first organ system studied was the heart, For the first time autoregulatory mechanisms have been demonstrated within the myocardial vascular bed which rapidly and precisely maintain coronary blood flow at a level i n keeping with the magnitude of cardiac work, Following a summary statement of the status of the program, the work performed is described, A financial review of expenditures incurred and a proposal for future funding are appended,

Shaw, R. F.↗
NASA NTRS1998

On Wave Processes in the Solar Atmosphere

This grant was awarded by NASA/MSFC to The University of Alabama in Huntsville (UAH) to investigate the physical processes responsible for heating and wind acceleration in the solar atmosphere, and to construct theoretical, self-consistent and time-dependent solar wind models based on the momentum deposition by finite amplitude and nonlinear Alfven waves. In summary, there are three main goals of the proposed research: (1) Calculate the wave energy spectra and wave energy fluxes carried by magnetic non- magnetic waves. (2) Find out which mechanism dominates in supplying the wave energy to different parts of the solar atmosphere. (3) Use the results obtained in (1) and (2) to construct theoretical, self-consistent and time- dependent models of the solar wind. We have completed the first goal by calculating the amount of non-radiative energy generated in the solar convection zone as acoustic waves and as magnetic tube waves. To calculate the amount of wave energy carried by acoustic waves, we have used the Lighthill-Stein theory for sound generation modified by Musielak, Rosner, Stein & Ulmschneider (1994). The acoustic wave energy fluxes for stars located in different regions of the Hertzsprung-Russell (H-R) diagram have also been computed. The wave energy fluxes carried by longitudinal and transverse waves along magnetic flux tubes have been calculated by using both analytical and numerical methods. Our analytical approach is based a theory developed by Musielak, Rosner & Ulmschnelder and Musielak, Rosner, Gall & Ulmschneider, which allows computing the wave energy fluxes for linear tube waves. A numerical approach has been developed by Huang, Musielak & Ulmschneider and Ulmschneider & Musielak to compute the energy fluxes for nonlinear tube waves. Both methods have been used to calculate the wave energy fluxes for stars located in different regions of the HR diagram (Musielak, Rosner & Ulmschneider 1998; Ulmschneider, Musielak & Fawzy 1998). Having obtained the wave energy fluxes for acoustic and magnetic tube waves, we have investigated the behavior of these waves in the solar and stellar atmospheres. The results of our extensive studies have been published in many papers and presented at numerous scientific meetings. In these studies we have investigated different aspects of propagation of acoustic and magnetic waves, the efficiency of energy transfer along magnetic structures in the solar atmosphere, and behavior of Alfven waves in stgeady and expanding solar and stellar atmospheres. Recently, we have used some of these results to construct first purely theoretical, two component and time-dependent models of solar and stellar chromospheres. Finally, to address the third goal, we have constructed first fully theoretical, self-consistent and time dependent wind models based on the momentum deposition by non-linear Alfven waves. The full set of single-fluid MHD equations with the background flow has been solved by using a modified version of the ZEUS MHD code. The constructed wind models are radially symmetric with the magnetic field decreasing radially and the initial outflow is described by the standard Parker wind solution. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; the models thus naturally account for the backreaction of the wind on the waves as well as for the nonlinear interaction between different types of MHD waves. The models have been used to explain the origin of fast speed streams in solar coronal holes. The obtained results clearly demonstrate that the momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast stream components of the solar wind. The range of wave amplitudes required to obtain the desired results seems to be in good agreement with recent observations.

Musielak, Z. E.↗