SEARCH · Search NASA

Results for “Variational formulation”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1987

The effect of along-front temperature variation in a two-dimensional primitive equation model of a surface frontogenesis

The effects of along-front temperature variations on the structure and evolution of surface fronts are investigated using a two-dimensional frontogenesis model. The formulation and characteristics of the confluence-horizontal shear model are described. The structures of three frontal zones forced by confluence and characterized by zero, negative, and positive along-front potential temperature variations are analyzed; particular attention is given to the ageostrophic vertical circulation and their forcings. The temporal evolution of the model frontal zones is examined in terms of the low-level maxima of vorticity, cross-front potential temperature gradient, and convergence. The kinematic and dynamic processes influencing the frontogenesis are discussed. It is noted that the two-dimensional model is applicable to the study of the dynamical properties of frontogenesis.

Keyser, Daniel↗

NASA NTRS1992

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗

NASA NTRS2007

Projection Method for Flows with Large Density Variations

Numerical models of solidification including a mushy zone are notoriously inefficient; most of them are based on formulations that require the coupled solution of the velocity components in the momentum equation greatly restricting the range of applicability of the models. There are only two models known to the authors that have used a projection or fractional step formulation, but none of these were used to model problems of any significant size. A third model was only applied to a partial mushy zone with no all-fluid region. Our initial attempts at modeling directional solidification in the presence of a developing mushy zone using a projection formulation encountered very serious difficulties once solidification starts. These difficulties were traced to the inability of the method to deal with large local density differences in the vicinity of the fluid-mush interface. As a result, a modified formulation of the projection method has been developed, that maintains the coupling between the body force and the pressure gradient and is presented in this work. The new formulation is shown to be robust and efficient, and can be applied to problems involving very large meshes. This is illustrated in this work through its application to simulations involving Pb-Sb and Pb-Sn alloys.

Heinrich, Juan C.↗

NASA NTRS1959

APPLICATION OF VARIATIONAL METHODS TO RADIATION HEAT-TRANSFER CALCULATIONS

A variational method is presented for solving a class of integral equations which arise in radiation heat-transfer problems. First, to demonstrate the formulation of radiation problems in terms of integral equations, consideration is given to a system consisting of two nonblack, finite, parallel plates. After a general description of the variational method, its use is illustrated by application to the parallel-plate system. Comparisons are made which show very good agreement with exact solutions.

E. M. Sparrow↗

NASA NTRS1974

Some examples of exact and approximate solutions in small particle scattering - A progress report

The formulation of basic equations from which the scattering of radiation by a localized variation in a medium is discussed. These equations are developed in both the differential and the integral form. Primary interest is in the scattering of electromagnetic waves for which the solution of the vector wave equation with appropriate boundary conditions must be considered. Scalar scattering by an infinite homogeneous isotropic circular cylinder, and scattering of electromagnetic waves by infinite circular cylinders are treated, and the case of the finite circular cylinder is considered. A procedure is given for obtaining angular scattering distributions from spheroids.

Greenberg, J. M.↗

NASA NTRS1981

Effect of choked flow on terminal characteristics of MPD thrusters

A one-dimensional steady real gas MPD accelerator model based on data for various electrode lengths and anode orifice radii is presented. The model relates the electric field to the boundary conditions in a form that agrees with experiment over parametric variations in current and mass flow. The electric field is sensitive to the formulations and assumptions bound into the analysis, and is useful in the modeling of outflow properties. It is shown that the location of the sonic point in the arc chamber shifts upstream for high current operation, and the relation between terminal voltage and current is demonstrated. Experiments using an anode suggested by the theory yield a thrust efficiency 1.5 times that of the unmodified thruster.

King, D. Q.↗

NASA NTRS1983

Interpretation of vector magnetograph data including magneto-optic effects. I - Azimuth angle of the transverse field

In this paper, the presence of Faraday rotation in measurements of the orientation of a sunspot's transverse magnetic field is investigated. Using observations obtained with the Marshall Space Flight Center's (MSFC) vector magnetograph, the derived vector magnetic field of a simple, symmetric sunspot is used to calculate the degree of Faraday rotation in the azimuth of the transverse field as a function of wavelength from analytical expressions for the Stokes parameters. These results are then compared with the observed rotation of the field's azimuth which is derived from observations at different wavelengths within the Fe I 5250 A spectral line. From these comparisons, it is found: the observed rotation of the azimuth is simulated to a reasonable degree by the theoretical formulations if the line-formation parameter is varied over the sunspot; these variations are substantiated by the line-intensity data; for the MSFC system, Faraday rotation can be neglected for field strengths less than 1800 G and field inclinations greater than 45 deg; to minimize the effects of Faraday rotation in sunspot umbrae, MSFC magnetograph measurements must be made in the far wings of the Zeeman-sensitive spectral line.

West, E. A.↗

NASA NTRS1987

Grace DAKASEP alkaline battery separator

The Grace DAKASEP separator was originally developed as a wicking layer for nickel-zinc alkaline batteries. The DAKASEP is a filled non-woven separator which is flexible and heat sealable. Through modification of formulation and processing variables, products with a variety of properties can be produced. Variations of DAKASEP were tested in Ni-H2, Ni-Zn, Ni-Cd, and primary alkaline batteries with good results. The properties of DAKASEP which are optimized for Hg-Zn primary batteries are shown in tabular form. This separator has high tensile strength, 12 micron average pore size, relatively low porosity at 46-48 percent, and consequently moderately high resistivity. Versions were produced with greater than 70 percent porosity and resistivities in 33 wt percent KOH as low as 3 ohm cm. Performance data for Hg-Zn E-1 size cells containing DAKASEP with the properties shown in tabular form, are more reproducible than data obtained with a competitive polypropylene non-woven separator. In addition, utilization of active material is in general considerably improved.

Giovannoni, R. T.↗

NASA NTRS1989

Numerical optimization in Hilbert space using inexact function and gradient evaluations

Trust region algorithms provide a robust iterative technique for solving non-convex unstrained optimization problems, but in many instances it is prohibitively expensive to compute high accuracy function and gradient values for the method. Of particular interest are inverse and parameter estimation problems, since function and gradient evaluations involve numerically solving large systems of differential equations. A global convergence theory is presented for trust region algorithms in which neither function nor gradient values are known exactly. The theory is formulated in a Hilbert space setting so that it can be applied to variational problems as well as the finite dimensional problems normally seen in trust region literature. The conditions concerning allowable error are remarkably relaxed: relative errors in the gradient error condition is automatically satisfied if the error is orthogonal to the gradient approximation. A technique for estimating gradient error and improving the approximation is also presented.

Carter, Richard G.↗

NASA NTRS1992

Vibration attenuation of the NASA Langley evolutionary structure experiment using H(sub infinity) and structured singular value (micron) robust multivariable control techniques

The use is studied of active control to attenuate structural vibrations of the NASA Langley Phase Zero Evolutionary Structure due to external disturbance excitations. H sub infinity and structured singular value (mu) based control techniques are used to analyze and synthesize control laws for the NASA Langley Controls Structures Interaction (CSI) Evolutionary Model (CEM). The CEM structure experiment provides an excellent test bed to address control design issues for large space structures. Specifically, control design for structures with numerous lightly damped, coupled flexible modes, collocated and noncollocated sensors and actuators and stringent performance specifications. The performance objectives are to attenuate the vibration of the structure due to external disturbances, and minimize the actuator control force. The control design problem formulation for the CEM Structure uses a mathematical model developed with finite element techniques. A reduced order state space model for the control design is formulated from the finite element model. It is noted that there are significant variations between the design model and the experimentally derived transfer function data.

Balas, Gary J.↗

NASA NTRS1993

Robust integrated flight/propulsion control design for a STOVL aircraft using H-infinity control design techniques

Results are presented from an application of H-infinity control design methodology to a centralized integrated flight/propulsion control (IFPC) system design for a supersonic STOVL fighter aircraft in transition flight. The emphasis is on formulating the H-infinity optimal control synthesis problem such that the critical requirements for the flight and propulsion systems are adequately reflected within the linear, centralized control problem formulation and the resulting controller provides robustness to modeling uncertainties and model parameter variations with flight condition. Detailed evaluation results are presented for a reduced order controller obtained from the improved H-infinity control design showing that the control design meets the specified nominal performance objective as well as provides stability robustness for variations in plant system dynamics with changes in aircraft trim speed within the transition flight envelope.

Garg, Sanjay↗

NASA NTRS1993

A numerical study of the unsteady processes associated with the type IV shock interaction

A time-dependent ADI formulation of the two dimensional Navier-Stokes equations coupled with the total variation diminishing (TVD) scheme is used to examine the type IV shock interaction time accurately. The experimentally measured and computationally calculated unsteadiness of the type IV interaction is shown to be related to the formation and shedding of shear layers in the shock layer, and the frequency of oscillation of the jet is a strong function of the shedding rate of these shear layers. For one particular geometry, an oscillation frequency of 1.4 kHz was calculated. The effects of shock impingement angle and shock strength on the location and magnitude of the peak pressure are also identified. In one case it is shown that a 1-deg change in the impinging shock angle can alter the location of the associated jet by almost 11 deg with a change in the peak surface pressure of 32 percent.

Lind, Charles A.↗

NASA NTRS1998

A Monte Carlo Uncertainty Analysis of Ozone Trend Predictions in a Two Dimensional Model

We use Monte Carlo analysis to estimate the uncertainty in predictions of total O3 trends between 1979 and 1995 made by the Goddard Space Flight Center (GSFC) two-dimensional (2D) model of stratospheric photochemistry and dynamics. The uncertainty is caused by gas-phase chemical reaction rates, photolysis coefficients, and heterogeneous reaction parameters which are model inputs. The uncertainty represents a lower bound to the total model uncertainty assuming the input parameter uncertainties are characterized correctly. Each of the Monte Carlo runs was initialized in 1970 and integrated for 26 model years through the end of 1995. This was repeated 419 times using input parameter sets generated by Latin Hypercube Sampling. The standard deviation (a) of the Monte Carlo ensemble of total 03 trend predictions is used to quantify the model uncertainty. The 34% difference between the model trend in globally and annually averaged total O3 using nominal inputs and atmospheric trends calculated from Nimbus 7 and Meteor 3 total ozone mapping spectrometer (TOMS) version 7 data is less than the 46% calculated 1 (sigma), model uncertainty, so there is no significant difference between the modeled and observed trends. In the northern hemisphere midlatitude spring the modeled and observed total 03 trends differ by more than 1(sigma) but less than 2(sigma), which we refer to as marginal significance. We perform a multiple linear regression analysis of the runs which suggests that only a few of the model reactions contribute significantly to the variance in the model predictions. The lack of significance in these comparisons suggests that they are of questionable use as guides for continuing model development. Large model/measurement differences which are many multiples of the input parameter uncertainty are seen in the meridional gradients of the trend and the peak-to-peak variations in the trends over an annual cycle. These discrepancies unambiguously indicate model formulation problems and provide a measure of model performance which can be used in attempts to improve such models.

Considine, D. B.↗

NASA NTRS1987

An optimal output feedback gain variation scheme for the control of plants exhibiting gross parameter changes

A concept for optimally designing output feedback controllers for plants whose dynamics exhibit gross changes over their operating regimes was developed. This was to formulate the design problem in such a way that the implemented feedback gains vary as the output of a dynamical system whose independent variable is a scalar parameterization of the plant operating point. The results of this effort include derivation of necessary conditions for optimality for the general problem formulation, and for several simplified cases. The question of existence of a solution to the design problem was also examined, and it was shown that the class of gain variation schemes developed are capable of achieving gain variation histories which are arbitrarily close to the unconstrained gain solution for each point in the plant operating range. The theory was implemented in a feedback design algorithm, which was exercised in a numerical example. The results are applicable to the design of practical high-performance feedback controllers for plants whose dynamics vary significanly during operation. Many aerospace systems fall into this category.

Moerder, Daniel D.↗

NASA NTRS1985

The Feinheit method: A phase-independent formulation of the period-luminosity relation for Cephids

The tight correlation between color and luminosity during the cyclical variation of a Cepheid is calibrated and shown to be a direct means by which most of the phase-dependent temperature-induced amplitude of the variable can be transformed away. The resulting Feinheit function, F = V - alpha (B-V) gives rise to a random-phase period-luminosity relation with an rms scatter of less than 0.20 mag. A comparison shows that the blue Feinheit method for distance determinations method from single-phase observations is as accurate as near-infrared photometry but has the added advantage of being able to use panoramic detectors for the data acquisition. The Feinheit function is identified as the surface area variation of the Cepheid during its cycle.

Madore, B. F.↗

NASA NTRS1995

Reduction technique for tire contact problems

A reduction technique and a computational procedure are presented for predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of the reduction technique, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface.

Noor, Ahmed K.↗

NASA NTRS1995

Advances in reduction techniques for tire contact problems

Some recent developments in reduction techniques, as applied to predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities, are reviewed. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of reduction techniques, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface. Also, the research topics which have high potential for enhancing the effectiveness of reduction techniques are outlined.

Noor, Ahmed K.↗

NASA NTRS1992

An assumed-stress hybrid 4-node shell element with drilling degrees of freedom

An assumed-stress hybrid/mixed 4-node quadrilateral shell element is introduced that alleviates most of the deficiencies associated with such elements. The formulation of the element is based on the assumed-stress hybrid/mixed method using the Hellinger-Reissner variational principle. The membrane part of the element has 12 degrees of freedom including rotational or 'drilling' degrees of freedom at the nodes. The bending part of the element also has 12 degrees of freedom. The bending part of the element uses the Reissner-Mindlin plate theory which takes into account the transverse shear contributions. The element formulation is derived from an 8-node isoparametric element by expressing the midside displacement degrees of freedom in terms of displacement and rotational degrees of freedom at corner nodes. The element passes the patch test, is nearly insensitive to mesh distortion, does not 'lock', possesses the desirable invariance properties, has no hidden spurious modes, and for the majority of test cases used in this paper produces more accurate results than the other elements employed herein for comparison.

Aminpour, M. A.↗