Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical analysis”

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

Sensitivity Analysis of Numerical Modeling Input Parameters on Wind Turbine Loads in Deterministic Transient Load Cases

Aero-hydro-elastic-servo numerical models used to design and analyze wind turbines are based on thousands of variable input parameters that dictate the inflow, aerodynamic, structural, and control characteristics of the system as well as sea state, hydrodynamic, and mooring characteristics for fixed-bottom and floating offshore wind turbines. Each of these parameters has some level of uncertainty, which can significantly impact the predicted loads. Understanding the uncertainty in the inputs is critical to understanding the uncertainty in the outputs. This work demonstrates a screening technique to identify which parameters ultimate loads are most sensitive to so that more focus can be given to quantifying the possible range of those parameters. This technique has been demonstrated previously for different turbine and load case types and is extended here for a floating offshore wind turbine in design load cases with transient events both in the inflow and operations. Each load case features a deterministic gust, including variations in wind speed, direction, and shear. Load cases are considered with an operating turbine as well as with prescribed fault, startup, and shutdown procedures. The study found that key input parameters with a large impact on loads include the length of the gust, the magnitude of direction change and speed in the gust, the initial wind speed, and the shape of the gust profile.

17 WIND ENERGY↗

On numerical nonlinear analysis of highly flexible spinning cantilevers

The general nonlinear discretized equations of motion of spinning elastic solids and structures are derived as a set of nonlinear ordinary differential equations for the case when the strain-displacement and velocity-displacement relations are nonlinear up to the second order. It is shown that the cost of generation of such equations is proportional to the fourth power of the number of degrees of freedom. A computer program is written to automatically generate the equations for the case of spinning cantilevers with initial imperfections. The types and the number of the coordinate functions used in the trial solution are parameters of the program.

Utku, S.↗

Equilibrium gas flow computations. II - An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, equilibrium gas laws. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer flux-vector splittings, and Roe's approximate Riemann solver are presented for three-dimensional, time-varying grids. The approximations inherent in previous generalizations are discussed.

Vinokur, Marcel↗

The use of available potential energy to evaluate the impact of satellite data on numerical model analysis during FGGE

To evaluate the effect of the FGGE satellite observing system, the following two data sets were compared by examining the available potential energy (APE) and extratropical cyclone activity within the entire global domain during the first Special Observing Period: (1) the complete FGGE IIIb set, which incorporates satellite soundings, and (2) a NOSAT set which incorporates only conventional data. The time series of the daily total APEs indicate that NOSAT values are larger than the FGGE values, although in the Northern Hemisphere the differences are negligible. Analyses of cyclone scale features revealed only minor differences between the Northern Hemisphere FGGE and NOSAT analyses. On the other hand, substantial differences were revealed in the two Southern Hemisphere analyses, where the satellite soundings apparently add detail to the FGGE set.

Horn, Lyle H.↗

Application of the p-version of the finite-element method to global-local problems

A brief survey is given of some recent developments in finite-element analysis technology which bear upon the three main research areas under consideration in this workshop: (1) analysis methods; (2) software testing and quality assurance; and (3) parallel processing. The variational principle incorporated in a finite-element computer program, together with a particular set of input data, determines the exact solution corresponding to that input data. Most finite-element analysis computer programs are based on the principle of virtual work. In the following, researchers consider only programs based on the principle of virtual work and denote the exact displacement vector field corresponding to some specific set of input data by vector u(EX). The exact solution vector u(EX) is independent of the design of the mesh or the choice of elements. Except for very simple problems, or specially constructed test problems, vector u(EX) is not known. Researchers perform a finite-element analysis (or any other numerical analysis) because they wish to make conclusions concerning the response of a physical system to certain imposed conditions, as if vector u(EX) were known.

Szabo, Barna A.↗

The systems improved numerical fluids analysis code (SINFAC)

A computer program has been developed to model the steady-state performance of single and two-phase heat transport systems. The development effort had the intent of making the program general, relatively easy to implement and easy to integrate into currently available software. The primary use of the program is for, but not limited to, heat transport system for spacecraft. In order to show the program has achieved its intended goals, the program was first applied to a small loop and a quick parametric study done. Next a verification study was done in which the temperature data from a thermal loop flown on the space shuttle was compared to modeler predictions.

Schweickart, Russell B.↗

A homotopy approach for combined control-structure optimization - Constructive analysis and numerical examples

This paper outlines the development of methods for the combined control-structure optimization of physical systems encountered in the technology of large space structures. The objectives of the approach taken in this paper is not to produce the 'best' optimized design, but rather to efficiently produce a family of design options so as to assist in early trade studies, typically before hard design constraints are imposed. The philosophy is that these are candidate designs to be passed on for further considerations, and their function is more to guide the development of the system design rather than to represent the ultimate product. A homotopy approach involving multi-objective functions is developed for this purpose. Analytical and numerical examples are also presented.

Scheid, R. E.↗

Fourier analysis of numerical algorithms for the Maxwell equations

The Fourier method is used to analyze the dispersive, dissipative, and isotropy errors of various spatial and time discretizations applied to the Maxwell equations on multi-dimensional grids. Both Cartesian grids and non-Cartesian grids based on hexagons and tetradecahedra are studied and compared. The numerical errors are quantitatively determined in terms of phase speed, wave number, propagation direction, gridspacings, and CFL number. The study shows that centered schemes are more efficient than upwind schemes. The non-Cartesian grids yield superior isotropy and higher accuracy than the Cartesian ones. For the centered schemes, the staggered grids produce less errors than the unstaggered ones. A new unstaggered scheme which has all the best properties is introduced. The study also demonstrates that a proper choice of time discretization can reduce the overall numerical errors due to the spatial discretization.

Liu, Yen↗

Radiative properties of a plasma moving across a magnetic field. I - Theoretical analysis. II - Numerical results

The early-time evolution of plasmas moving across a background magnetic field is addressed with a 2D model in which a plasma cloud is assumed to have formed instantaneously with a velocity across a uniform background magnetic field and with a Gaussian density profile in the two dimensions perpendicular to the direction of motion. This model treats both the dynamics associated with the formation of a polarization field and the generation and propagation of electromagnetic waves. In general, the results indicate that, to zeroth order, the plasma cloud behaves like a large dipole antenna oriented in the direction of the polarization field which oscillates at frequencies defined by the normal mode of the system. Radiation damping is shown to play an important role in defining the plasma cloud evolution, causing a rapid decay of the polarizaiton field and a loss of plasma kinetic energy and momentum on time scales comprable to several ion gyroperiods. Scaling laws are derived for the plasma momentum and energy loss rates, and predictions for the braking time, the amplitude and spectrum of the radiation field, and the total radiated power are presented for conditions relevant to the recent Combined Release and Radiation Effects Satellite experiments.

Roussel-Dupre, Robert↗

Effect of the thickness variation and initial imperfection on buckling of composite cylindrical shells: Asymptotic analysis and numerical results by BOSOR4 and PANDA2

This study is an extension of a previous investigation of the combined effect of axisymmetric thickness variation and axisymmetric initial geometric imperfection on buckling of isotropic shells under uniform axial compression. Here the anisotropic cylindrical shells are investigated by means of Koiter's energy criterion. An asymptotic formula is derived which can be used to determine the critical buckling load for composite shells with combined initial geometric imperfection and thickness variation. Results are compared with those obtained by the software packages BOSOR4 and PANDA2.

Li, Yi-Wei↗

Interpolation Method Needed for Numerical Uncertainty Analysis of Computational Fluid Dynamics

Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors in an unstructured grid, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors. Nomenclature

CFD↗

Analysis of Numerical Simulation Database for Pressure Fluctuations Induced by High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 and Reynolds number of Re(sub T) approximately 460 are conducted at two wall temperatures (Tw/Tr = 0.25, 0.76) to investigate the generated pressure fluctuations and their dependence on wall temperature. Simulations indicate that the influence of wall temperature on pressure fluctuations is largely limited to the near-wall region, with the characteristics of wall-pressure fluctuations showing a strong temperature dependence. Wall temperature has little influence on the propagation speed of the freestream pressure signal. The freestream radiation intensity compares well between wall-temperature cases when normalized by the local wall shear; the propagation speed of the freestream pressure signal and the orientation of the radiation wave front show little dependence on the wall temperature.

Duan, Lian↗

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Hu, Fang↗