Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,423 records · Page 79

Numerical methods for the design and analysis of wings at supersonic speeds

Numerical methods for the design and analysis of arbitrary-planform wings at supersonic speeds are reviewed. Certain deficiencies are revealed, particularly in application to wings with slightly subsonic leading edges. Recently devised numerical techniques which overcome the major part of these deficiencies are presented. The original development as well as the more recent revisions are subjected to a thorough review.

Carlson, H. W.↗

Determination of mineral dissolution regimes using flow-through time-resolved analysis (FT-TRA) and numerical simulation

Flow-through time resolved analysis (FT-TRA) involves subjecting small mineral samples (< 10 mg) inserted in a miniature flow-through cell (50 μL) to controlled flows of eluent analyzed on-line by ICP-MS. In this study, FT-TRA is used to empirically determine the dissolution regimes for the two well-studied minerals forsterite and calcite, representing minerals with relatively slow and fast dissolution kinetics. A proportional increase in steady-state effluent [Mg, Si] concentrations with increasing flow-through cell eluent residence times confirms a dominantly surface-controlled dissolution regime for a powdered forsterite sample at pH 2.3, implying that transport limitations are negligible. In contrast, the relationship between flow rates and dissolution rates for single grain calcite samples at pH 2.3-4 reveals that transport limitations affect the rate of calcite dissolution. In order to provide a quantitative and process-based assessment of the effect of diffusive transport limitations, simulations of the calcite experiments were performed with a high resolution, pore-scale model that considers the geometry of the calcite grain and the FT-TRA flow-through reactor. The pore-scale model reproduces the observed effluent [Ca] concentrations for all experimental conditions using a single set of surface kinetic parameters, by accounting for the formation of a diffusive boundary layer (DBL) that varies in thickness as a function of flow rates. These results demonstrate that combining FT-TRA with pore-scale modeling makes it possible to obtain unprecedented insights not achievable by either method separately, including quantification of DBL thicknesses and the determination of transport controls as a function of pH, flow velocity and residence times.

58 GEOSCIENCES↗

Numerical approach for the aerodynamic analysis if airfoils with laminar separation

A numerical method for simultaneously and efficiently coupling an external subsonic potential flow and an interior viscous flow such that the two flows match at an interfacing boundary is discussed. Both a panel method and a simple point compressible vortex model are used for the outer potential field. The interior flow solvers which were used are the Navier-Stokes and Euler codes of T. J. Coakley and the Euler code of A. Verhoff. In order to test compatibility, the panel method is coupled to the less expensive Euler codes since the coupling procedure is identical with the Navier-Stokes code. The results show significant efficiency improvements can be obtained over the uncoupled approach. Results also indicate the outer potential flow is best represented by the simple point compressible vortex model. The panel method couples smoothly to Coakley's implicit code but is numerically incompatible as coupled with the explicit Euler code. An improved Navier-Stokes code is under initial development which extends the Euler code to include the necessary viscous terms. Results are shown for all infinite length channel with one wavy periodic wall with and without laminar separation.

Halt, D. W.↗

A numerical solution of the Navier-Stokes equations for supercritical fluid thermodynamic analysis

An explicit numerical solution of the compressible Navier-Stokes equations is applied to the thermodynamic analysis of supercritical oxygen in the Apollo cryogenic storage system. The wave character is retained in the conservation equations which are written in the basic fluid variables for a two-dimensional Cartesian coordinate system. Control-volume cells are employed to simplify imposition of boundary conditions and to ensure strict observance of local and global conservation principles. Non-linear real-gas thermodynamic properties responsible for the pressure collapse phenomonon in supercritical fluids are represented by tabular and empirical functions relating pressure and temperature to density and internal energy. Wall boundary conditions are adjusted at one cell face to emit a prescribed mass flowrate. Scaling principles are invoked to achieve acceptable computer execution times for very low Mach number convection problems. Detailed simulations of thermal stratification and fluid mixing occurring under low acceleration in the Apollo 12 supercritical oxygen tank are presented which model the pressure decay associated with de-stratification induced by an ordinary vehicle maneuver and heater cycle operation.

Heinmiller, P. J.↗

Numerical methods for the analysis of sampled-data systems and for the computation of system zeros

MARSYAS is a computer aided control system design package for the simulation and analysis of open loop and closed loop dynamic systems. Outlined here is the numerical and theoretical basis of the MARSYAS functions developed during the summer of 1991. In particular, the numerical computation of the matrix exponential e(exp A) = I + A + A(exp 2)/2! + A(exp 3)/3! + ... and the numerical computation of the finite system zeros of a dynamic system (continuous time or discrete time) are presented.

Hodel, A. Scottedward↗

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗

FPCAS2D user's guide, version 1.0

The FPCAS2D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady two-dimensional full potential equation which is solved for a cascade of blades. The structural analysis is based on a two degree-of-freedom rigid typical section model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS2D code. A complete description of the input data is provided in this report. In addition, four test cases, including inputs and outputs, are provided.

Bakhle, Milind A.↗

FPCAS3D User's guide: A three dimensional full potential aeroelastic program, version 1

The FPCAS3D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady three-dimensional full potential equation which is solved for a blade row. The structural analysis is based on a finite-element model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS3D code. A complete description of the input data is provided in this report. In addition, six examples, including inputs and outputs, are provided.

Bakhle, Milind A.↗

Multi-variate Factorisation of Numerical Simulations

Factorisation is widely used in the analysis of numerical simulations. It allows changes in properties of a system to be attributed to changes in multiple variables associated with that system. There are many possible factorisation methods; here we discuss three previously-proposed factorisations that have been applied in the field of climate modelling: the linear factorisation, the Stein and Alpert (1993) factorisation, and the Lunt et al (2012) factorisation. We show that, when more than two variables are being considered, none of these three methods possess all three properties of "uniqueness", "symmetry", and "completeness". Here, we extend each of these factorisations so that they do possess these properties for any number of variables, resulting in three factorisations – the "linear-sum" factorisation, the "shared-interaction" factorisation, and the "scaled-total" factorisation. We show that the linear-sum factorisation and the shared-interaction factorisation reduce to be identical. We present the results of the factorisations in the context of studies that used the previously-proposed factorisations. This reveals that only the linear-sum/shared-interaction factorisation possesses a fourth property – "boundedness", and as such we recommend the use of this factorisation in applications for which these properties are desirable.

Numerican simulations↗

Hybrid Symbolic-Numerical Modeling and Parametric Stability Analysis of DC–AC Power Systems

Hybrid DC-AC power systems integrating diverse inverter-based resources (IBRs) and multi-terminal high-voltage direct current (MTDC) networks represent a promising architecture for future power grids, while introducing challenges for modeling, stability analysis, and control design. This paper develops a hybrid symbolic-numerical modeling framework and tool to characterize the parametric small-signal stability of DC-AC coupled power systems. The proposed approach constructs parametric state-space models to enable efficient representation of system dynamics under varying control parameters and network configurations, with target parameters retained as symbolic variables and the remainder treated numerically. The stability analysis framework covers eigenvalue, sensitivity, and stability boundary and region characterization. Enhanced linear matrix inequality (LMI) techniques are proposed to directly certify small-signal stability over regions of parameter space while also reducing the conservativeness and computational burden. The resulting tools and frameworks enable rapid parametric model construction across diverse grid conditions, thereby facilitating stability-informed control and operation in complex DC–AC power systems.

DC–AC power systems↗

Stability of mixing layers

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

Tam, Christopher↗

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

FDTD modeling of induced polarization phenomena in transient electromagnetics

The finite-difference time-domain scheme is augmented in order to treat the modelling of transient electromagnetic signals containing induced polarization effects from 3-D distributions of polarizable media. Compared to the non-dispersive problem, the discrete dispersive Maxwell system contains costly convolution operators. Key components to our solution for highly digitized model meshes are Debye decomposition and composite memory variables. We revert to the popular Cole-Cole model of dispersion to describe the frequency-dependent behaviour of electrical conductivity. Its inversely Laplace-transformed Debye decomposition results in a series of time convolutions between electric field and exponential decay functions, with the latter reflecting each Debye constituents' individual relaxation time. Finally, these function types in the discrete-time convolution allow for their substitution by memory variables, annihilating the otherwise prohibitive computing demands. Numerical examples demonstrate the efficiency and practicality of our algorithm.

58 GEOSCIENCES↗

Sonic-point capturing

A prototype scheme that produces perfectly smooth transonic solutions to nozzle-flow problems is derived and tested. The basic upwind scheme is described as well as satisfying the entropy condition, treatment of the source term, and numerical verification. The analysis yielded a numerical flux function for use near a sonic point, which is based on a full model of a transonic expansion wave, and a matched treatment for the source term.

Van Leer, Bram↗

A numerical method for the stress analysis of stiffened-shell structures under nonuniform temperature distributions

A numerical method is presented for the stress analysis of stiffened-shell structures of arbitrary cross section under nonuniform temperature distributions. The method is based on a previously published procedure that is extended to include temperature effects and multicell construction. The application of the method to practical problems is discussed and an illustrative analysis is presented of a two-cell box beam under the combined action of vertical loads and a nonuniform temperature distribution.

Heldenfels, Richard R↗

Asymptotic analysis of dissipative waves with applications to their numerical simulation

Various problems involving the interplay of asymptotics and numerics in the analysis of wave propagation in dissipative systems are studied. A general approach to the asymptotic analysis of linear, dissipative waves is developed. It was applied to the derivation of asymptotic boundary conditions for numerical solutions on unbounded domains. Applications include the Navier-Stokes equations. Multidimensional traveling wave solutions to reaction-diffusion equations are also considered. A preliminary numerical investigation of a thermo-diffusive model of flame propagation in a channel with heat loss at the walls is presented.

Hagstrom, Thomas↗