Search NASA⌕ Search

SEARCH · Search NASA

Results for “second moment 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 181 records · Page 10

Finite-field Calculations of Molecular Polarizabilities Using Field-Induced Polarization Functions: Second- and Third-order Perturbation Correlation Corrections to the Coupled Hartree-Fock Polarizability of H2O

Ordinary Rayleigh-Schroudinger perturbation theory with Moller-Plesset (RSMP) partitioning is used to calculate second- and third-order correlation corrections to the CHF polarizability and dipole moment of the water molecule by a finite-field procedure. Pade approximants are found to be useful in accelerating the convergence of the property perturbation expansions. Field-induced polarization functions suitable for polarizability calculations are determined. The average polarizability calculated, neglecting vibrational averaging, with Dunning's (9s5p/4s-4s2p/2s) contracted GTO basis set augmented by field-induced lslp2d/lp polarization functions is within 3 per cent of the experimental result. Correlation corrections to the dipole moment and polarizability of the water molecule calculated by the finite-field RSMP and single + double excitation CI(SDCI) methods for the same basis set are found to be in close agreement. The RSMP approach has the advantages of being size-consistent and of being capable of greater efficiency than the SCDI method. Comparative calculations carried out using Epstein-Nesbet partitioning show that through third order RSEN correlation perturbation expansions for the dipole moment and polarizability are less rapidly convergent than RSMP expansions. However, reasonable accord with RSMP results can be achieved by using Pade approximants to accelerate the convergence of RSEN energy perturbation expansions. The convergence of RSMP property correlation expansions based on the zeroth-order uncoupled-Hartree-Fock (UCHF) and coupled-Hartree-Fock (CHF) approximations are compared through third order. Whereas the CHF + RSMP expansions are for practical purposes fully converged, the UCHF + RSMP expansions are not adequately converged.

Langhoff, S. R.↗

Variational Coupled Loads Analysis using the Hybrid Parametric Variation Method

Time-domain coupled loads analysis (CLA)is used to determine the response of a launch vehicle and payload system to transient forces, such as liftoff, engine ignitions and shutdowns, jettison events, and atmospheric flight loads, such as buffet. CLA, using Hurty/Craig-Bampton (HCB)component models, is the accepted method for the establishment of design-level loads for launch systems. However, uncertainty in the component models flows into uncertainty in predicted system results. Uncertainty in the structural responses during launch is a significant concern because small variations in launch vehicle and payload mode shapes and their interactions can result in significant variations in system loads. Uncertainty quantification (UQ)is used to determine statistical bounds on prediction accuracy based on model uncertainty. In this paper uncertainty is treated at the HCB component-model level. In an effort to account for model uncertainties and statistically bound their effect on CLA predictions, this work combines CLA with UQ in a process termed variational coupled loads analysis (VCLA). The modeling of uncertainty using a parametric approach, in which input parameters are represented by random variables, is common, but its major drawback is the resulting uncertainty is limited to the form of the nominal model. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. Model-form uncertainty can be represented using a nonparametric approach based on random matrix theory (RMT). In this work, UQ is performed using the hybrid parametric variation (HPV)method, which combines parametric with nonparametric uncertainty at the HCB component model level. The HPV method requires the selection of dispersion values for the HCB fixed-interface (FI)eigenvalues, and the HCB mass and stiffness matrices. The dispersions are based upon component test-analysis modal correlation results. During VCLA, random component models are assembled into an ensemble of random systems using a Monte Carlo (MC)approach. CLA is applied to each of the ensemble members to produce an ensemble of system-level responses for statistical analysis. The proposed methodology is demonstrated through its application to a buffet loads analysis of NASA’s Space Launch System (SLS)during the transonic regime fifty seconds after liftoff. Core stage (CS)section shears and moments are recovered, and statistics are computed.

Uncertainty Quantification↗

An investigation on the effect of second-order additional thickness distributions to the upper surface of an NACA 64-206 airfoil

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of an NACA 64-206 airfoil. Additional thickness distributions employed were in the form of two second-order polynomial arcs which have a specified thickness at a given chordwise location. The forward arc disappears at the airfoil leading edge, the aft arc disappears at the airfoil trailing edge. At the juncture of the two arcs, x = x, continuity of slope is maintained. The effect of varying the maximum additional thickness and its chordwise location on airfoil lift coefficient, pitching moment, and pressure distribution was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic NACA 64-206 airfoil, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.

Merz, A. W.↗

Continental and oceanic crustal magnetization modelling

Inversion of magnetic data from the MAGSAT satellite, to arrive at intensities of magnetization of the Earth's crust, was performed by two different methods. The first method uses a spherical harmonic model of the magnetic field. The coefficients believed to represent sources in the Earth's crust can then be inverted to arrive at vertical dipole moments per unit area at the Earth's surface. The spherical harmonic models contain coefficients of degrees of harmonics up to 23. The dipole moment per unit area for a surface element can then be determined by summing the contribution for each individual degree of harmonic. The magnetic moments were calculated for continental and oceanic areas separately as well as over certain latitudinal segments. Of primary concern was to determine whether there are any differences between continental and oceanic areas. The second analysis with magnetization intensities was made using narrower ranges of degrees of harmonics, assuming that higher degrees are present in the core field signal.

Harrison, C. G. A.↗

An investigation on the effect of second-order additional thickness distributions to the upper surface of an NACA 64 sub 1-212 airfoil

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of an NACA 64 sub 1 - 212 airfoil. Additional thickness distributions employed were in the form of two second-order polynomial arcs which have a specified thickness at a given chordwise location. The forward arc disappears at the airfoil leading edge, the aft arc disappears at the airfoil trailing edge. At the juncture of the two arcs, x = x, continuity of slope is maintained. The effect of varying the maximum additional thickness and its chordwise location on airfoil lift coefficient, pitching moment, and pressure distribution was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic NACA 64 sub 1 - 212 airfoil, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.

Hague, D. S.↗

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

Growth of thin films of organic nonlinear optical materials by vapor growth processes - An overview and examination of shortfalls

Research on the growth of second- and third-order nonlinear optical (NLO) organic thin film by vapor deposition is reviewed. Particular attention is given to the experimental methods for growing thin films of p-chlorophenylurea, diacetylenes, and phthalocyanines; characteristics of the resulting films; and approaches for advancing thin film technology. It is concluded that the growth of NLO thin films by vapor processes is a promising method for the fabrication of planar waveguides for nonlinear optical devices. Two innovative approaches are proposed including a method of controlling the input beam frequency to maximize nonlinear effects in thin films and single crystals, and the alternate approach to the molecular design of organic NLO materials by increasing the transition dipole moment between ground and excited states of the molecule.

Frazier, D. O.↗

Hit rates and radiation doses to nuclei of bone lining cells from alpha-particle-emitting radionuclides

Factors relating the local concentration of a bone-seeking alpha-particle emitter to the mean hit rate have been determined for nuclei of bone lining cells using a Monte Carlo procedure. Cell nuclei were approximated by oblate spheroids with dimensions and location taken from a previous histomorphometric study. The Monte Carlo simulation is applicable for planar and diffuse labels at plane or cylindrical bone surfaces. Additionally, the mean nuclear dose per hit, the dose mean per hit, the mean track segment length and its second moment, the percentage of stoppers, and the frequency distribution of the dose have been determined. Some basic features of the hit statistics for bone lining cells have been outlined, and the consequences of existing standards of radiation protection with regard to the hit frequency to cell nuclei are discussed.

NASA Discipline Musculoskeletal↗

A viscous/potential flow interaction analysis method for multi-element infinite swept wings, volume 1

An analysis method and computer program have been developed for the calculation of the viscosity dependent aerodynamic characteristics of multi-element infinite swept wings in incompressible flow. The wing configuration consisting at the most of a slat, a main element and double slotted flap is represented in the method by a large number of panels. The inviscid pressure distribution about a given configuration in the normal chord direction is determined using a two dimensional potential flow program employing a vortex lattice technique. The boundary layer development over each individual element of the high lift configuration is determined using either integral or finite difference boundary layer techniques. A source distribution is then determined as a function of the calculated boundary layer displacement thickness and pressure distributions. This source distribution is included in the second calculation of the potential flow about the configuration. Once the solution has converged (usually after 2-5 iterations between the potential flow and boundary layer calculations) lift, drag, and pitching moments can be determined as functions of Reynolds number.

Dvorak, F. A.↗

Backscatter RCS for TE and TM excitations of dielectric-filled cavity-backed apertures in two-dimensional bodies

Transverse electric (TE) and transverse magnetic (TM) scattering from dielectric-filled, cavity-backed apertures in two-dimensional bodies are treated using the method of moments technique to solve a set of combined-field integral equations for the equivalent induced electric and magnetic currents on the exterior of the scattering body and on the associated aperture. Results are presented for the backscatter radar cross section (RCS) versus the electrical size of the scatterer for two different dielectric-filled cavity-backed geometries. The first geometry is a circular cylinder of infinite length which has an infinite length slot aperture along one side. The cavity inside the cylinder is dielectric filled and is also of circular cross section. The two cylinders (external and internal) are of different radii and their respective longitudinal axes are parallel but not collocated. The second is a square cylinder of infinite length which has an infinite length slot aperture along one side. The cavity inside the square cylinder is dielectric-filled and is also of square cross section.

Goggans, Paul M.↗

Kinematics of velocity and vorticity correlations in turbulent flow

The kinematic problem of calculating second-order velocity moments from given values of the vorticity covariance is examined. Integral representation formulas for second-order velocity moments in terms of the two-point vorticity correlation tensor are derived. The special relationships existing between velocity moments in isotropic turbulence are expressed in terms of the integral formulas yielding several kinematic constraints on the two-point vorticity correlation tensor in isotropic turbulence. Numerical evaluation of these constraints suggests that a Gaussian curve may be the only form of the longitudinal velocity correlation coefficient which is consistent with the requirement of isotropy. It is shown that if this is the case, then a family of exact solutions to the decay of isotropic turbulence may be obtained which contains Batchelor's final period solution as a special case. In addition, the computed results suggest a method of approximating the integral representation formulas in general turbulent shear flows.

Bernard, P. S.↗

NASA Experience with Pogo in Human Spaceflight Vehicles

An overview of more than 45 years of NASA human spaceflight experience is presented with respect to the thrust axis vibration response of liquid fueled rockets known as pogo. A coupled structure and propulsion system instability, pogo can result in the impairment of the astronaut crew, an unplanned engine shutdown, loss of mission, or structural failure. The NASA history begins with the Gemini Program and adaptation of the USAF Titan II ballistic missile as a spacecraft launch vehicle. It continues with the pogo experienced on several Apollo-Saturn flights in both the first and second stages of flight. The defining moment for NASA s subsequent treatment of pogo occurred with the near failure of the second stage on the ascent of the Apollo 13 mission. Since that time NASA has had a strict "no pogo" philosophy that was applied to the development of the Space Shuttle. The "no pogo" philosophy lead to the first vehicle designed to be pogo-free from the beginning and the first development of an engine with an integral pogo suppression system. Now, more than 30 years later, NASA is developing two new launch vehicles, the Ares I crew launch vehicle propelling the Orion crew excursion vehicle, and the Ares V cargo launch vehicle. A new generation of engineers must again exercise NASA s system engineering method for pogo mitigation during design, development and verification.

Larsen, Curtis E.↗

Ozone and nitrogen dioxide ground based monitoring by zenith sky visible spectrometry in Arctic and Antarctic

Unattended diode array spectrometers have been designed for ground based stratospheric trace species monitoring by zenith sky visible spectrometry. Measurements are performed with a 1.0 nm resolution between 290 nm and 590 nm in order to allow simultaneous evaluations of column densities of ozone, nitrogen dioxide. Field tests have shown that the species can be monitored with a precision of + or - 2 Dobson for the first and + or - 2.10 to the 15th mol/sq cm for the second, although the absolute accuracy of the method is limited by the error of the estimation of the atmospheric optical path of the scattered light. Two identical instruments were set up in January 1988, one in Antarctica at Dumont d'Urville (66 S, 140 E) to be operated all year and another one in the Arctic at ESRANGE at Kiruna (68 N; 22 E) which will operate to the final warming of spring 1988. The data are processed in real time at both stations. O3 and NO2 columns are transmitted together with surface and stratospheric temperature and winds. They are also recorded for further treatment and search for OClO and BrO. Only one month of data from Antarctica is available at the moment. Obtained during polar summer, they cannot show more than stable columns of O3 and NO2 and for the last species, the buildup of its diurnal variation.

Pommereau, J. P.↗

Exploratory Studies in Generalized Predictive Control for Active Gust Load Alleviation

The results of numerical simulations aimed at assessing the efficacy of Generalized Predictive Control (GPC) for active gust load alleviation using trailing- and leading-edge control surfaces are presented. The equations underlying the method are presented and discussed, including system identification, calculation of control law matrices, and calculation of commands applied to the control effectors. Both embedded and explicit feedforward paths for inclusion of disturbance effects are addressed. Results from two types of simulations are shown. The first used a 3-DOF math model of a mass-spring-dashpot system subject to user-defined external disturbances. The second used open-loop data from a wind-tunnel test in which a wing model was excited by sinusoidal vertical gusts; closed-loop behavior was simulated in post-test calculations. Results obtained from these simulations have been decidedly positive. In particular, results of closed-loop simulations for the wing model showed reductions in root moments by factors as high as 1000, depending on whether the excitation is from a constant- or variable-frequency gust and on the direction of the response.

Kvaternik, Raymond G.↗

Distribution of SLS Integrated Load Uncertainty to Surface Pressures and Sectional Loads

Aerodynamic loads that are important to launch vehicle programs such as NASA’s Space Launch System (SLS) include both integrated loads such as the force & moment on the entire vehicle and distributed loads. In this work two examples of distributed loads are considered: the pressure field on the surface of the vehicle and sectional loads, which are one-dimensional distributions along the axis of the launch vehicle. In some modern flight programs, the integrated loads, such as lift and drag, used to design the guidance and control laws for the vehicle come from wind tunnel testing, while distributed loads are produced using Computational Fluid Dynamics(CFD). The first task that is addressed in this paper, then, is to provide a formal method to adjust the distributed loads so that integrating them matches the prescribed integrated load. In addition, the integrated loads in a launch vehicle typically include an uncertainty estimate. The second task is to distribute this prescribed integrated uncertainty to each point in a distributed load. Both tasks are addressed using the same technique, which is to create distributed load profiles that isolate adjustments to one integrated load while leaving the others unaffected. These adjustments are informed by Proper Orthogonal Decomposition (POD) of the entire CFD-based distributed load database. Once applied, the adjusted distributed loads can be used to evaluate any scalar quantity of interest that might be needed by downstream users such as structural analysis or trajectory modelers.

SLS↗

Nearfield Summary and Statistical Analysis of the Second AIAA Sonic Boom Prediction Workshop

A summary is provided for the Second AIAA Sonic Boom Workshop held 8-9 January 2017 in conjunction with AIAA SciTech 2017. The workshop used three required models of increasing complexity: an axisymmetric body, a wing body, and a complete configuration with flow-through nacelle. An optional complete configuration with propulsion boundary conditions is also provided. These models are designed with similar nearfield signatures to isolate geometry and shock/expansion interaction effects. Eleven international participant groups submitted nearfield signatures with forces, pitching moment, and iterative convergence norms. Statistics and grid convergence of these nearfield signatures are presented. These submissions are propagated to the ground, and noise levels are computed. This allows the grid convergence and the statistical distribution of a noise level to be computed. While progress is documented since the first workshop, improvement to the analysis methods for a possible subsequent workshop are provided. The complete configuration with flow-through nacelle showed the most dramatic improvement between the two workshops. The current workshop cases are more relevant to vehicles with lower loudness and have the potential for lower annoyance than the first workshop cases. The models for this workshop with quieter ground noise levels than the first workshop exposed weaknesses in analysis, particularly in convective discretization.

Park, Michael A.↗

Magnetopause shape as a bivariate function of interplanetary magnetic field B(sub z) and solar wind dynamic pressure

We present a new method for determining the shape of the magnetopause as a bivariate function of the hourly averaged solar wind dynamic pressure (p) and the north-south component of the interplanetary magnetic field (IMF) B(sub z). We represent the magnetopause (for X(sub GSE) greater than -40 R(sub E)) as an ellipsoid of revolution in solar-wind-aberrated coordinates and express the (p, B(sub z)) dependence of each of the three ellipsoid parameters as a second-order (6-term) bivariate expansion in Inp and B(sub z). We define 12 overlapping bins in a normalized dimensionless (p, B(sub z)) `control space' and fit an ellipsoid to those magnetopause crossings having (p, B(sub z)) values within each bin. We also calculate the bivariate (Inp, B(sub z)) moments to second order over each bin in control space. We can then calculate the six control-space expansion coefficients for each of the three ellipsoid parameters in configuration space. From these coefficients we can derive useful diagnosis of the magnetopause shape as joint functions of p and B(sub z): the aspect ratio of the ellipsoid's minor-to-major axes; the flank distance, radius of curvature, and flaring angle (at X(sub GSE) = 0); and the subsolar distance and radius of curvature. We confirm and quantify previous results that during periods of southward B(sub z) the subsolar magnetopause moves inward, while at X(sub GSE) = 0 the flank magnetopause moves outward and the flaring angle increases.

Roelof, Edmond C.↗

Nearfield Summary and Statistical Analysis of the Second AIAA Sonic Boom PredictionWorkshop

A summary is provided for the Second AIAA Sonic Boom Workshop held 8–9 January 2017 in conjunction with AIAA SciTech 2017. The workshop used three required models of increasing complexity: an axisymmetric body, a wing body, and a complete configuration with flow-through nacelle. An optional complete configuration with propulsion boundary conditions is also provided. These models are designed with similar nearfield signatures to isolate geometry and shock/expansion interaction e ects. Eleven international participant groups submitted nearfield signatures with forces, pitching moment, and iterative convergence norms. Statistics and grid convergence of these nearfield signatures are presented. These submissions are propagated to the ground, and noise levels are computed. This allows the grid convergence and the statistical distribution of a noise level to be computed. While progress is documented since the first workshop, improvement to the analysis methods for a possible subsequent workshop are provided. The complete configuration with flow-through nacelle showed the most dramatic improvement between the two workshops. The current workshop cases are more relevant to vehicles with lower loudness and have the potential for lower annoyance than the first workshop cases. The models for this workshop with quieter ground noise levels than the first workshop exposed weaknesses in analysis, particularly in convective discretization.

Park, Michael A.↗