Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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,657 records · Page 92

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Bershader, D.↗

Wires, Solitons and the Big Bang. Final report

The overarching aim of the project was to improve our understanding of inflation and its end. The following were the main goals of the project: 1. Develop a novel statistical framework to understand particle production during inflation and reheating in scenarios with many poorly constrained components and calculate potential observational signatures. The goal was to be achieved in the following steps (i) Develop the theoretical framework (ii) Apply the framework to inflation and reheating 2. Develop analytical and numerical tools to explore the end of inflation and the energy transfer from inflation to daughter fields during reheating. The goal was to be achieved in two steps: (i) Characterize the eq. of state at the end of inflation using non-perturbative techniques and lattice simulations. (ii) Develop a novel numerical algorithm for gauge fields in a cosmological setting.

42 ENGINEERING↗

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers ($\mathscr{M}$ son ≲ 10 –2 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant–Friedrichs–Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume SEVEN-LEAGUE HYDRO (SLH) code, and it makes use of a low-Mach version of the five-wave Harten–Lax–van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit–explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at $\mathscr{M}$ son ~ 10 –3 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.

79 ASTRONOMY AND ASTROPHYSICS↗

Trust-Region Approximation of Extreme Trajectories in Power System Dynamics

In this work we present a novel technique, based on a trust-region optimization algorithm and second-order trajectory sensitivities, to compute the extreme trajectories of power system dynamic simulations given a bounded set that represents parametric uncertainty. Furthermore, we show how this method, while remaining computationally efficient compared with sampling-based techniques, overcomes the limitations of previous sensitivity-based techniques to approximate the bounds of the trajectories when the local approximation loses validity because of the nonlinearity. We present several numerical experiments that showcase the accuracy and scalability of the technique, including a demonstration on the IEEE New England test system.

42 ENGINEERING↗

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗

Development of a numerical procedure to map a general 3-d body onto a near-circle

Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.

Hommel, M. J.↗

Numerical evidence against advantage with quantum fidelity kernels on classical data

Quantum machine learning techniques are commonly considered one of the most promising candidates for demonstrating practical quantum advantage. In particular, quantum kernel methods have been demonstrated to be able to learn certain classically intractable functions efficiently if the kernel is well aligned with the target function. In the more general case, quantum kernels are known to suffer from exponential “flattening” of the spectrum as the number of qubits grows, preventing generalization and necessitating the control of the inductive bias by hyperparameters. We show that the general-purpose hyperparameter-tuning techniques proposed to improve the generalization of quantum kernels lead to the kernel becoming well approximated by a classical kernel, removing the possibility of quantum advantage. We provide extensive numerical evidence for this phenomenon utilizing multiple previously studied quantum feature maps and both synthetic and real data. Our results show that unless novel techniques are developed to control the inductive bias of quantum kernels, they are unlikely to provide a quantum advantage on classical data that lacks special structure.

Quantum Information, Science & Technology↗

Dynamic loads analyses of flexible airplanes - New and existing techniques

This paper reviews existing techniques for calculating dynamic loads for flexible airplanes and presents a new technique. The new technique involves the summation-of-forces method of writing dynamic loads equations. Until now this form of the dynamic loads equations has been formulated in the frequency domain. The new technique uses s-plane approximation methods (previously applied only to the equations of motion) to transform the dynamic loads equations from a second-order frequency-domain formulation with frequency-dependent coefficients into a linear-time-invariant state-space formulation. Several numerical examples demonstrate the usefulfness of the new technique and the high quality of the results. In addition, a convergence investigation establishes that the summation-of-forces method converges more quickly (that is, with fewer modes) than does the mode displacement method.

Pototzky, A. S.↗

Parallel pivoting combined with parallel reduction and fill-in control

Parallel algorithms for triangularization of large, sparse, and unsymmetric matrices are presented. The method combines the parallel reduction with a new parallel pivoting technique, control over generation of fill-ins and check for numerical stability, all done in parallel with the work being distributed over the active processes. The parallel pivoting technique uses the compatibility relation between pivots to identify parallel pivot candidates and uses the Markowitz number of pivots to minimize fill-in. This technique is not a preordering of the sparse matrix and is applied dynamically as the decomposition proceeds.

Alaghband, Gita↗

A Polar Scaling Technique for the Regularization of Strongly Singular and Strongly Near-Singular Helmholtz Surface Integrals Evaluated Over 2-D Domains

The numerical integration of expressions containing strong singularities or strong near-singularities has long been a challenging problem in the electromagnetics community. Much attention has been paid to this problem, as strong $1/R^{{2}}$ singularities routinely appear when implementing electromagnetic simulation techniques like the method of moments (MoM). To date, several techniques, from singularity extraction to singularity cancellation (SC), have been employed to deal with problems that require the evaluation of 2-D strongly singular integrals. However, no single technique has been proposed that can deal with both strong singularities and strong near-singularities in a fully numerical manner for arbitrary 2-D domains. Moreover, it has been claimed that the Helmholtz-type strongly singular integral found in the MoM is convergent in a principal value sense, but this convergence value has yet to be proven mathematically. In this work, we will conduct the convergence proof and introduce a “polar scaling” change of variables method that may be used to evaluate Helmholtz integrals with both strong and weak singularities/near-singularities. The technique is fully numerical and can in principle be applied to any planar or curved polygon and any nonsingular basis function. We will also provide numerical results showing useful convergence behavior for integrals involving both exact and near-singularities.

47 OTHER INSTRUMENTATION↗

Annual Research Briefs: 1995

This report contains the 1995 annual progress reports of the Research Fellows and students of the Center for Turbulence Research (CTR). In 1995 CTR continued its concentration on the development and application of large-eddy simulation to complex flows, development of novel modeling concepts for engineering computations in the Reynolds averaged framework, and turbulent combustion. In large-eddy simulation, a number of numerical and experimental issues have surfaced which are being addressed. The first group of reports in this volume are on large-eddy simulation. A key finding in this area was the revelation of possibly significant numerical errors that may overwhelm the effects of the subgrid-scale model. We also commissioned a new experiment to support the LES validation studies. The remaining articles in this report are concerned with Reynolds averaged modeling, studies of turbulence physics and flow generated sound, combustion, and simulation techniques. Fundamental studies of turbulent combustion using direct numerical simulations which started at CTR will continue to be emphasized. These studies and their counterparts carried out during the summer programs have had a noticeable impact on combustion research world wide.

Source record↗

Evolution of Local Microstructures: Spatial Instabilities of Coarsening Clusters

This work examines the diffusional growth of discrete phase particles dispersed within a matrix. Engineering materials are microstructurally heterogeneous, and the details of the microstructure determine how well that material performs in a given application. Critical to the development of designing multiphase microstructures with long-term stability is the process of Ostwald ripening. Ripening, or phase coarsening, is a diffusion-limited process which arises in polydisperse multiphase materials. Growth and dissolution occur because fluxes of solute, driven by chemical potential gradients at the interfaces of the dispersed phase material, depend on particle size. The kinetics of these processes are "competitive," dictating that larger particles grow at the expense of smaller ones, overall leading to an increase of the average particle size. The classical treatment of phase coarsening was done by Todes, Lifshitz, and Slyozov, (TLS) in the limit of zero volume fraction, V(sub v), of the dispersed phase. Since the publication of TLS theory there have been numerous investigations, many of which sought to describe the kinetic scaling behavior over a 0 range of volume fractions. Some studies in the literature report that the relative increase in coarsening rate at low (but not zero) volume fractions compared to that predicted by TLS is proportional to v(sub v)(exp 1/2), whereas others suggcest V(sub v)(exp 1/3). This issue has been resolved recently by simulation studies at low volume fractions in three dimensions by members of the Rensselaer/MSFC team. Our studies of ripening behavior using large-scale numerical simulations suggest that although there are different circumstances which can lead to either scaling law, the most important length scale at low volume fractions is the diffusional analog of the Debye screening length. The numerical simulations we employed exploit the use of a recently developed "snapshot" technique, and identifies the nature of the coarsening dynamics at various volume fractions. Preliminary results of numerical and experimental investigations, focused on the growth of finite particle clusters, provide important insight into the nature of the transition between the two scaling regimes. The companion microgravity experiment centers on the growth within finite particle clusters, and follows the temporal dynamics driving microstructural evolution, using holography.

Frazier, Donald O.↗

Composite laminates with spatially varying fiber orientations - 'Variable stiffness panel concept'

A solution has been obtained to the plane elasticity problem for a symmetrically laminated composite panel with spatially varying fiber orientations. Since variation of the fiber angles along the length of a composite laminate results in stiffness properties that change as a function of location, the laminates are called variable stiffness panels. An analysis of the stiffness variation and its effect on the elastic response of the panel is presented here. A numerical solution has been obtained using an iterative collocation technique. Corresponding closed-form solutions are given for three different sets of boundary conditions. Two of the cases considered have exact solutions and thus serve to validate the numerical model.

Guerdal, Zafer↗

Development of iterative techniques for the solution of unsteady compressible viscous flows

During the past two decades, there has been significant progress in the field of numerical simulation of unsteady compressible viscous flows. At present, a variety of solution techniques exist such as the transonic small disturbance analyses (TSD), transonic full potential equation-based methods, unsteady Euler solvers, and unsteady Navier-Stokes solvers. These advances have been made possible by developments in three areas: (1) improved numerical algorithms; (2) automation of body-fitted grid generation schemes; and (3) advanced computer architectures with vector processing and massively parallel processing features. In this work, the GMRES scheme has been considered as a candidate for acceleration of a Newton iteration time marching scheme for unsteady 2-D and 3-D compressible viscous flow calculation; from preliminary calculations, this will provide up to a 65 percent reduction in the computer time requirements over the existing class of explicit and implicit time marching schemes. The proposed method has ben tested on structured grids, but is flexible enough for extension to unstructured grids. The described scheme has been tested only on the current generation of vector processor architecture of the Cray Y/MP class, but should be suitable for adaptation to massively parallel machines.

Hixon, Duane↗

Diabatic forcing and initialization with assimilation of cloud and rain water in a forecast model: Methodology

The focus of this part of the investigation is to find one or more general modeling techniques that will help reduce the time taken by numerical forecast models to initiate or spin-up precipitation processes and enhance storm intensity. If the conventional data base could explain the atmospheric mesoscale flow in detail, then much of our problem would be eliminated. But the data base is primarily synoptic scale, requiring that a solution must be sought either in nonconventional data, in methods to initialize mesoscale circulations, or in ways of retaining between forecasts the model generated mesoscale dynamics and precipitation fields. All three methods are investigated. The initialization and assimilation of explicit cloud and rainwater quantities computed from conservation equations in a mesoscale regional model are examined. The physical processes include condensation, evaporation, autoconversion, accretion, and the removal of rainwater by fallout. The question of how to initialize the explicit liquid water calculations in numerical models and how to retain information about precipitation processes during the 4-D assimilation cycle are important issues that are addressed. The explicit cloud calculations were purposely kept simple so that different initialization techniques can be easily and economically tested. Precipitation spin-up processes associated with three different types of weather phenomena are examined. Our findings show that diabatic initialization, or diabatic initialization in combination with a new diabatic forcing procedure, work effectively to enhance the spin-up of precipitation in a mesoscale numerical weather prediction forecast. Also, the retention of cloud and rain water during the analysis phase of the 4-D data assimilation procedure is shown to be valuable. Without detailed observations, the vertical placement of the diabatic heating remains a critical problem.

Raymond, William H.↗

An optimal numerical filter for wide-field-of-view measurements of earth-emitted radiation

A technique is described in which all data points along an arc of the orbit may be used in an optimal numerical filter for wide-field-of-view measurements of earth emitted radiation. The statistical filter design is derived whereby the filter is required to give a minimum variance estimate of the radiative exitance at discrete points along the ground track of the satellite. An equation for the optimal numerical filter is given by minimizing the estimate error variance equation with respect to the filter weights, resulting in a discrete form of the Wiener-Hopf equation. Finally, variances of the errors in the radiant exitance can be computed along the ground track and in the cross track directions.

Smith, G. L.↗

Observation simulation experiments with regional prediction models

Research efforts in FY 1990 included studies employing regional scale numerical models as aids in evaluating potential contributions of specific satellite observing systems (current and future) to numerical prediction. One study involves Observing System Simulation Experiments (OSSEs) which mimic operational initialization/forecast cycles but incorporate simulated Advanced Microwave Sounding Unit (AMSU) radiances as input data. The objective of this and related studies is to anticipate the potential value of data from these satellite systems, and develop applications of remotely sensed data for the benefit of short range forecasts. Techniques are also being used that rely on numerical model-based synthetic satellite radiances to interpret the information content of various types of remotely sensed image and sounding products. With this approach, evolution of simulated channel radiance image features can be directly interpreted in terms of the atmospheric dynamical processes depicted by a model. Progress is being made in a study using the internal consistency of a regional prediction model to simplify the assessment of forced diabatic heating and moisture initialization in reducing model spinup times. Techniques for model initialization are being examined, with focus on implications for potential applications of remote microwave observations, including AMSU and Special Sensor Microwave Imager (SSM/I), in shortening model spinup time for regional prediction.

Diak, George↗

A nonlinear programming approach for optimizing two-stage lifting vehicle ascent to orbit

An optimal atmospheric flight branched trajectory-shaping capability is presented based on the Davidon-Fletcher-Powell variable metric parameter optimization technique. Gradient information is generated using finite difference methods. A typical atmospheric flight branched optimization problem is analyzed which requires the determination of 31 parameters. This parameter set includes the three-dimensional description of vehicle attitude control angles for three branches of flight: first-stage ascent, second-stage ascent, and first-stage flyback. The important inflight inequality contraints required to maintain the integrity of the vehicles are considered. Some of the numerical methods employed are discussed, along with several new auxiliary techniques developed to improve the compatibility of the numerical gradient and iterator.

Kamm, J. L.↗