SEARCH · Engineering Papers

Results for “numerical calculations: variational”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Implicit total variation diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme. Previously announced in STAR as N83-23085

Yee, H. C.↗

NASA NTRS1988

A scalar/vector potential formulation of an airfoil in nonuniform stream

A recent innovation in aerodynamic calculation techniques is the formulation based upon scalar and vector potentials (Stokes-Helmholtz decomposition) for the velocity field. This technique is presented: an implementation involving classical panel method is studied. The application of interest is to two-dimensional airfoils moving into a nonuniform approach flow. Comparisons with theoretical and numerical results are included. Furthermore the variations of lift and moment coefficients of quasi-steady simulated flight of an airfoil through microburst wind shear data are studied.

Wang, J. C. A.↗

NASA NTRS2006

NPSS Multidisciplinary Integration and Analysis

The objective of this task was to enhance the capability of the Numerical Propulsion System Simulation (NPSS) by expanding its reach into the high-fidelity multidisciplinary analysis area. This task investigated numerical techniques to convert between cold static to hot running geometry of compressor blades. Numerical calculations of blade deformations were iteratively done with high fidelity flow simulations together with high fidelity structural analysis of the compressor blade. The flow simulations were performed with the Advanced Ducted Propfan Analysis (ADPAC) code, while structural analyses were performed with the ANSYS code. High fidelity analyses were used to evaluate the effects on performance of: variations in tip clearance, uncertainty in manufacturing tolerance, variable inlet guide vane scheduling, and the effects of rotational speed on the hot running geometry of the compressor blades.

Hall, Edward J.↗

NASA NTRS1992

A time-accurate high-resolution TVD scheme for solving the Navier-Stokes equations

A total variation diminishing (TVD) scheme has been developed and incorporated into an existing time-accurate high-resolution Navier-Stokes code. The accuracy and the robustness of the resulting solution procedure have been assessed by performing many calculations in four different areas: shock tube flows, regular shock reflection, supersonic boundary layer, and shock boundary layer interactions. These numerical results compare well with corresponding exact solutions or experimental data.

Kim, Hyun Dae↗

NASA NTRS1993

A time-accurate high-resolution TVD scheme for solving the Navier-Stokes equations

Kim, Hyun D.↗

NASA NTRS2014

Solar Ion Processing of Itokawa Grains: Reconciling Model Predictions with Sample Observations

Analytical TEM observations of Itokawa grains reported to date show complex solar wind ion processing effects in the outer 30-100 nm of pyroxene and olivine grains. The effects include loss of long-range structural order, formation of isolated interval cavities or "bubbles", and other nanoscale compositional/microstructural variations. None of the effects so far described have, however, included complete ion-induced amorphization. To link the array of observed relationships to grain surface exposure times, we have adapted our previous numerical model for progressive solar ion processing effects in lunar regolith grains to the Itokawa samples. The model uses SRIM ion collision damage and implantation calculations within a framework of a constant-deposited-energy model for amorphization. Inputs include experimentally-measured amorphization fluences, a Pi steradian variable ion incidence geometry required for a rotating asteroid, and a numerical flux-versus-velocity solar wind spectrum.

Christoffersen, Roy↗

DOE OSTI2025

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1988

Hypersonic turbulent wall boundary layer computations

The Baldwin-Lomax algebraic turbulence model was modified for hypersonic flow conditions. Two coefficients in the outer layer eddy viscosity model were determined as functions of Mach number and temperature ratio. By matching the solutions from the Baldwin-Lomax model to those from the Cebeci-Smith model for a flat plate at hypersonic speed, the new values of the coefficient were obtained. The results show that the values of C sub cp and C sub kleb are functions of both Mach number and wall temperature ratio. The C sub cp and C sub kleb variations with Mach number and wall temperature were used for the calculations of both a 4 deg wedge flow at Mach 18 and an axisymmetric Mach 20 nozzle flow. The Navier-Stokes equations with thin layer approximation were solved for the above hypersonic flow conditions and the results were compared with existing experimental data. The agreement between the numerical solutions and the existing experimental data were good. The modified Baldwin-Lomax model thus is useful in the computations of hypersonic flows.

Kim, S. C.↗

NASA NTRS1988

Hypersonic turbulent wall boundary layer computations

The Baldwin-Lomax (1978) algebraic turbulence model was modified for hypersonic flow conditions. Two coefficients in the outer-layer eddy-viscosity model were determined as functions of Mach number and temperature ratio. By matching the solutions from the Baldwin-Lomax model to those from the Cebeci-Smith (1974) model for a flat plate at hypersonic speed, the new values of the coefficients were obtained. The results show that the values of C(cp) and C(kleb) are functions of both Mach number and wall temperature ratio. The C(cp) and C(kleb) variations with Mach number and wall temperature were used for the calculations of both a 4-deg wedge flow at Mach 18 and an axisymmetric Mach 20 nozzle flow. The Navier-Stokes equations with thin-layer approximation were solved for the above hypersonic flow conditions and the results were compared with existing experimental data. The agreement between the numerical solutions and the existing experimental data were good. The modified Baldwin-Lomax model thus is useful in the computations of hypersonic flows.

Kim, S. C.↗

NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

NASA NTRS1978

The relation between temporal and spatial stability in three-dimensional flows

An analysis is presented of the nonparallel spatial or temporal stability of three-dimensional incompressible, isothermal boundary-layer flows taking into account the transverse velocity component as well as the axial and crossflow variations of the mean flow. The method of multiple scales is used to derive partial differential equations that describe the axial and crossflow variations of the disturbance amplitude, phase and wavenumbers. This equation is used to derive the expressions that relate the temporal and spatial instabilities. These relations are functions of the complex group velocities. Moreover, this equation is used to derive the expression that relates the spatial amplification in any direction to a calculated amplification in any other direction. These relations are verified by numerical results obtained for two- and three-dimensional disturbances in two- and three-dimensional flows.

Nayfeh, A. H.↗

NASA NTRS1989

An improvement of convection fidelity in Euler calculations

A new solution procedure is presented to solve the Euler equations for steady, compressible, rotational, inviscid flows. The approach is aimed at achieving real inviscid solutions in Euler calculations by eliminating numerical diffusion. The variables in the Euler equations are divided into elliptic and convective quantities, using the Clebsch velocity decomposition. The convective quantities are then transported without numerical contamination using an efficient convection operator, while the elliptic quantities are integrated with a relaxation procedure. This approach provides a generalization of the full potential formulation to rotation Euler physics by allowing variations of convective quantities. Results are demonstrated for several transonic flows in two dimensions.

Lee, K. D.↗

NASA NTRS1990

Analysis of the single-blow transient testing technique for perforated plate heat exchangers

An analysis of the single-blow transient testing technique for heat exchanger surfaces is made for perforated plate and similar discontinuous surfaces. The model assumes that there is no temperature variation across each plate, and allows for axial conduction in spacers: the resulting axial fluid temperature profile is discontinuous, rather than the usual continuous profile. Numerical solutions are obtained to the resulting set of coupled first-order differential equations for a step change in inlet fluid temperature. Results are presented in tabular form, which allow the heat transfer coefficients to be calculated from test data using the maximum slope technique.

Rodriguez, Jose I.↗

DOE OSTI2025

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads in extreme idling conditions

Abstract. Floating offshore wind turbine (FOWT) systems are subject to complex environmental loads, with significant potential for damage in extreme storm conditions. Design simulations in these conditions are required to assess the survivability of the device with some level of confidence. Aero-hydro-servo-elastic engineering tools can be used with a reasonable balance of accuracy and computational efficiency. The models require many input parameters to describe the air and water conditions, the system properties, and the load calculations. Each of these parameters has some possible range, due to either statistical uncertainty or variations with time. Variation in the input parameters can have important effects on the uncertainty in the resulting loads, but it is not practical to perform detailed assessments of the impact of this uncertainty for every input parameter. This work demonstrates a method to identify the input parameters that have the most impact on the loads to focus further inspection. The process is done specifically for extreme storm load cases defined in the International Electrotechnical Commission design requirements for floating offshore wind turbines. The analysis was performed using the International Energy Agency Wind 15 MW offshore reference wind turbine atop the University of Maine VolturnUS-S reference platform in two US offshore wind regions, the Gulf of Maine and Humboldt Bay. It was found that the direction of incident waves and current, yaw misalignment, and the length of mooring line sections were among the primary sensitivities.

17 WIND ENERGY↗

DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗

NASA NTRS1985

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The novel implicit and unconditionally stable, high resolution Total Variation Diminishing (TVD) scheme whose application to steady state calculations is presently examined is a member of a one-parameter family of implicit, second-order accurate systems developed by Harten (1983) for the computation of weak solutions for one-dimensional hyperbolic conservation laws. The scheme will not generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments for a quasi-one-dimensional nozzle problem show that the experimentally determined stability limit correlates exactly with the theoretical stability limit for the nonlinear scalar hyberbolic conservation laws.

Yee, H. C.↗

NASA NTRS1985

Application of TVD schemes for the Euler equations of gas dynamics

First-order, second-order, and implicit total variation diminishing (TVD) schemes are reviewed using the modified flux approach. Some transient and steady-state calculations are then carried out to illustrate the applicability of these schemes to the Euler equations. It is shown that the second-order explicit TVD schemes generate good shock resolution for both transient and steady-state one-dimensional and two-dimensional problems. Numerical experiments for a quasi-one-dimensional nozzle problem show that the second-order implicit TVD scheme produces a fairly rapid convergence rate and remains stable even when running with a Courant number of 10 to the 6th.

Yee, H. C.↗

NASA NTRS2000

Magnetic Susceptibility Effects and Lorentz Damping in Diamagnetic Fluids

A great number of crystals (semi-conductor and protein) grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity and g-jitter. Both static and dynamic (rotating or travelling wave) magnetic fields can be used to reduce the effects of convection in materials processing. In semi-conductor melts, due to their relatively high electrical conductivity, the induced Lorentz force can be effectively used to curtail convective effects. In melts/solutions with reduced electrical conductivity, such as aqueous solutions used in solution crystal growth, protein crystal growth and/or model fluid experiments for simulating melt growth, however, the variation of the magnetic susceptibility with temperature and/or concentration can be utilized to better damp fluid convection than the Lorentz force method. This paper presents a comprehensive, comparative numerical study of the relative damping effects using static magnetic fields and gradients in a simple geometry subjected to a thermal gradient. The governing equations are formulated in general terms and then simplified for the numerical calculations. Operational regimes, based on the best damping technique for different melts/solutions are identified based on fluid properties. Comparisons are provided between the numerical results and available results from experiments in surveyed literature.

Ramachandran, Narayanan↗