Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal form methods”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1935

Spinning characteristics of wings II : rectangular Clark Y biplane cellule: 25 percent stagger; 0 degree decalage; gap/chord 1.0

General methods of theoretical analysis of airplane spinning characteristics have been available for some time. Some of these methods of analysis might be used by designers to predict the spinning characteristics of proposed airplane designs if the necessary aerodynamic data were known. The present investigation, to determine the spinning characteristics of wings, is planned to include variations in airfoil sections, plan forms, and tip shapes of monoplane wings and variations in stagger, gap, and decalage for biplane cellules. The first series of tests, made on a rectangular Clark Y monoplane wing, are reported in reference 1. That report also gives an analysis of the data for predicting the probable effects of various important parameters on the spin for normal airplanes using such a wing. The present report is the second of the series. It gives the aerodynamic characteristics of a rectangular Clark Y biplane cellule in spinning attitudes and includes a discussion of the data, using the method of analysis given in reference 1.

Bamber, M J↗
NASA NTRS1990

The simultaneous integration of many trajectories using nilpotent normal forms

Taylor's formula shows how to approximate a certain class of functions by polynomials. The approximations are arbitrarily good in some neighborhood whenever the function is analytic and they are easy to compute. The main goal is to give an efficient algorithm to approximate a neighborhood of the configuration space of a dynamical system by a nilpotent, explicitly integrable dynamical system. The major areas covered include: an approximating map; the generalized Baker-Campbell-Hausdorff formula; the Picard-Taylor method; the main theorem; simultaneous integration of trajectories; and examples.

Grayson, Matthew A.↗
NASA NTRS1943

The exact solution of shear-lag problems in flat panels and box beams assumed rigid in the transverse direction

A mathematical procedure is herein developed for obtaining exact solutions of shear-lag problems in flat panels and box beams: the method is based on the assumption that the amount of stretching of the sheets in the direction perpendicular to the direction of essential normal stresses is negligible. Explicit solutions, including the treatment of cut-outs, are given for several cases and numerical results are presented in graphic and tabular form. The general theory is presented in a from which further solutions can be readily obtained. The extension of the theory to cover certain cases of non-uniform cross section is indicated. Although the solutions are obtained in terms of infinite series, the present developments differ from those previously given in that, in practical cases, the series usually converge so rapidly that sufficient accuracy is afforded by a small number of terms. Comparisons are made in several cases between the present results and the corresponding solutions obtained by approximate procedures devised by Reissner and by Kuhn and Chiarito.

STRESSED SKIN-PLATES, FLAT↗
NASA NTRS1993

Uniform convergence of multigrid V-cycle iterations for indefinite and nonsymmetric problems

In this paper, we present an analysis of a multigrid method for nonsymmetric and/or indefinite elliptic problems. In this multigrid method various types of smoothers may be used. One type of smoother which we consider is defined in terms of an associated symmetric problem and includes point and line, Jacobi, and Gauss-Seidel iterations. We also study smoothers based entirely on the original operator. One is based on the normal form, that is, the product of the operator and its transpose. Other smoothers studied include point and line, Jacobi, and Gauss-Seidel. We show that the uniform estimates for symmetric positive definite problems carry over to these algorithms. More precisely, the multigrid iteration for the nonsymmetric and/or indefinite problem is shown to converge at a uniform rate provided that the coarsest grid in the multilevel iteration is sufficiently fine (but not depending on the number of multigrid levels).

Bramble, James H.↗
DOE OSTI2021

Evaluation of residual gas fraction estimation methods for cycle-to-cycle combustion variability analysis and modeling

Cycle-to-cycle combustion variability in spark-ignition engines during normal operation is mainly caused by random perturbations of the in-cylinder conditions such as the flow velocity field, homogeneity of the air-fuel distribution, spark energy discharge, and turbulence intensity of the flame front. Such perturbations translate into the variability of the energy released observed at the end of the combustion process. During normal operating conditions, the cycle-to-cycle variability (CCV) of the energy release behaves as random uncorrelated noise. However, during diluted combustion, in either the form of exhaust gas recirculation (EGR) or excess air (lean operation), the CCV tends to increase as dilution increases. Moreover, when the ignition limit is reached at high dilution levels, the combustion CCV is exacerbated by sporadic occurrences of incomplete combustion events, and the uncorrelation assumption no longer holds. The low or null energy released by partial burns and misfires has an impact on the following combustion event due to the residual gas that carries burned and unburned gases, which contributes to the deterministic coupling between engine cycles. Many residual gas fraction estimation methods, however, only address the nominal case where complete combustion occurs and combustion events are uncorrelated. Here we evaluate the efficacy of such methods on capturing the effects of partial burns and misfires on the residual gas estimate for high-EGR operation. The advantages and disadvantages of each method are discussed based on their ability to generate cycle-to-cycle estimates. Finally, a comparison between the different estimation techniques is presented based on their usefulness for control-oriented modeling.

42 ENGINEERING↗
DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improving nanoparticle superlattice stability with deformable polymer gels

The self-assembly of colloidal nanoparticles into ordered superlattices typically uses dynamic interactions to govern particle crystallization, as these non-permanent bonds prevent the formation of kinetically trapped, disordered aggregates. However, while the use of reversible bonding is critical in the formation of highly ordered particle arrangements, dynamic interactions also inherently make the structures more prone to disassembly or disruption when subjected to different environmental stimuli. Thus, there is typically a trade-off between the ability to initially form an ordered colloidal material and the ability of that material to retain its order under different conditions. Here, we present a method for embedding colloidal nanoparticle superlattices into a polymer gel matrix. This encapsulation strategy physically prevents the nanoparticles from dissociating upon heating, drying, or the introduction of chemicals that would normally disrupt the lattice. However, the use of a gel as the embedding medium still permits further modification of the colloidal nanoparticle lattice by introducing stimuli that deform the gel network (as this deformation in turn alters the nanoparticle lattice structure in a predictable manner). Moreover, encapsulation of the lattice within a gel permits further stabilization into fully solid materials by removing the solvent from the gel or by replacing the solvent with a liquid monomer that can be photopolymerized. This embedding method therefore makes it possible to incorporate ordered colloidal arrays into a polymer matrix as either dynamic or static structures, expanding their potential for use in responsive materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Confidence regions of planar cardiac vectors

A method for plotting the confidence regions of vectorial data obtained in electrocardiology is presented. The 90%, 95% and 99% confidence regions of cardiac vectors represented in a plane are obtained in the form of an ellipse centered at coordinates corresponding to the means of a sample selected at random from a bivariate normal distribution. An example of such a plot for the frontal plane QRS mean electrical axis for 80 horses is also presented.

Dubin, S.↗
NASA NTRS1993

Fluorescent-Antibody Measurement Of Cancer-Cell Urokinase

Combination of laboratory techniques provides measurements of amounts of urokinase in and between normal and cancer cells. Includes use of fluorescent antibodies specific against different forms of urokinase-type plasminogen activator, (uPA), fluorescence microscopy, quantitative analysis of images of sections of tumor tissue, and flow cytometry of different uPA's and deoxyribonucleic acid (DNA) found in suspended-tumor-cell preparations. Measurements provide statistical method for indicating or predicting metastatic potentials of some invasive tumors. Assessments of metastatic potentials based on such measurements used in determining appropriate follow-up procedures after surgical removal of tumors.

Morrison, Dennis R.↗
DOE OSTI2024

Unifying Combinatorial and Graphical Methods in Artificial Intelligence

Recently, a new graph Laplacian, called the inner product Laplacian, was introduced which generalizes many existing Laplacians, including the normalized and combinatorial Laplacian and their weighted variants. The key observation behind the inner product Laplacian is that by defining appropriate inner product spaces on the vertices and edges, the standard Laplacians can be recovered as Hodge Laplacians over the simplicial complex formed by the edges and vertices. These inner product spaces form a natural way to incorporate non-combinatorial information into the definition of a domain-specific Laplacian. In particular, in contrast to current domain-specific weighting schemes which rely solely on edge weights, information regarding the similarity of non-adjacent vertices and arbitrary pairs of edges can be effectively incorporated into the Laplacian. In order to illustrate this approach we consider the problem of calculating the potential energy of an atomistic configuration using Graph Neural Networks. In comparison with start-of-the-art approaches, such as SchNet, our approach replaces a learned (via auto-encoder) representation of the atom types with an inner product space on atoms based on scientific knowledge (e.g., electronegativity). We will illustrate how this approach captures key chemical properties of the molecules and compare the energy calculations with state-of-the-art neural network approaches. However, to compute the resulting Laplacian involves a mixture of sparse and dense matrix computation and yields a dense matrix as the basis for the graph convolution. This dense convolutional kernel necessitates moving away from the standard message passing framework for graph neural networks and increases the computational cost of applying the kernel. In order to mitigate these costs we investigate means of leveraging the mixed sparse and dense computations to reduce the overall computational cost and how these approaches can be automatically transferred to energy efficient hardware (e.g., field programmable gate arrays (FPGAs)).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Tidal estimation in the Atlantic and Indian Oceans, 3 deg x 3 deg solution

An estimation technique was developed to extrapolate tidal amplitudes and phases over entire ocean basins using existing gauge data and the altimetric measurements provided by satellite oceanography. The technique was previously tested. Some results obtained by using a 3 deg by 3 deg grid are presented. The functions used in the interpolation are the eigenfunctions of the velocity (Proudman functions) which are computed numerically from a knowledge of the basin's bottom topography, the horizontal plan form and the necessary boundary conditions. These functions are characteristic of the particular basin. The gravitational normal modes of the basin are computed as part of the investigation; they are used to obtain the theoretical forced solutions for the tidal constituents. The latter can provide the simulated data for the testing of the method and serve as a guide in choosing the most energetic functions for the interpolation.

Sanchez, Braulio V.↗
DOE OSTI2023

Birefringent waveplate and method for forming a waveplate having a birefringent metasurface

The present disclosure relates to a waveplate having a substrate forming an optic. The substrate may have an integral portion forming a plurality of angled columnar features on an exposed surface thereof. The plurality of angled columnar features may further be aligned parallel with a directional plane formed non-parallel to a reference plane, with the reference plane being normal to a surface of the substrate. The metasurface forms a birefringent metasurface.

Feigenbaum, Eyal↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗
DOE OSTI2024

Criteria for Grid-Forming Inverters Transitioning Between Current Limiting Mode and Normal Operation

Here, this letter focuses on the fault recovery analysis of inverters using multi-loop droop based grid-forming control, taking into account different current limiting strategies. A criterion is proposed to determine when the inverter can exit the current limiting mode after fault clearance. The theoretical findings are verified through electromagnetic transient simulations using the direct (d)-axis priority-based current limiting method, circular current limiting method, and virtual impedance method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1972

The acoustic response of rooms with open windows to airborne sounds.

The objective of the work described in this and the companion paper was to establish a theory for predicting the sound field generated in a room by a sonic boom incident on an open window. In this paper, some basic theoretical results are presented. First, the case of a normally incident harmonic wave was considered. Expressions for the pressure field were obtained by viewing the room as a terminated duct and by using a Green function method. The concept of mode excitation distribution functions was formulated and used to match the boundary conditions. This concept has been extended for oblique incidence. A modified form of Laplace transform technique was used to obtain expressions in the time domain for transient signals.

Vaidya, P. G.↗
NASA NTRS2001

Acoustical and Flowfield Characterization of a Scaled Tabletop Rocket

An analysis of the acoustical and flowfield environment for the scaled 1-pound-force (lbf) thrust tabletop motor was performed. This tabletop motor from NASA Stennis Space Center Is composed of Plexiglas burning In gaseous oxygen with a graphite insert for the nozzle portion. The nozzle has a throat diameter of 0.2 inch and an exit diameter of 0.38 Inch. With a chamber pressure at 55 pounds per square Inch absolute (psia), a normal shock is formed immediately downstream of the nozzle exit plane as the combustion products exhaust into the ambient at atmospheric pressure. The jet characterization Is based on computational fluid dynamics (CFD) in conjunction with Kirchhoff surface integral formulation and compared with correlations developed for measured rocket noise and a pressure fluctuation scaling (PFS) method. Predictions and comparisons are made for the overall sound pressure levels (OASPL's) and spectral dependence of sound pressure level (SPL). The overall objective of this effort is to develop methods for scaling the acoustic and flowfield environment of rockets with a wide range of thrust (1 lbf to 1 million lbf).

Kandula, Max↗
DOE OSTI2025

Vexcel Imaging's Suitability for Automatic Verification

Accurate, independently verified geospatial data is essential for automated calibration, validation, and operational decision-making. This study evaluated the positional accuracy of Vexcel Imaging™’s UltraCam® Osprey imagery (7.5 cm GSD) using globally distributed Continuously Operating Reference Stations (CORS) as independent control. Despite manufacturer claims of 15 cm horizontal accuracy, residual errors were consistently one to two orders of magnitude larger, with no subset of imagery meeting precision thresholds. These discrepancies cannot be explained by normal photogrammetric or environmental factors and raise concerns about the reliability of the imagery for high-precision tasks. The results demonstrate that Vexcel imagery, in its current form, is unsuitable for workflows requiring rigorous spatial accuracy or automated verification. At the same time, the reproducible validation framework developed in this study establishes a scalable method for assessing commercial imagery, ensuring that future products can be independently and objectively verified before operational adoption.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Development and experimental qualification of novel disruption prevention techniques on DIII-D

Novel disruption prevention solutions spanning a range of control regimes are being developed and tested on DIII-D to enable ITER success. First, a new real-time control algorithm has been developed and tested for regulating nearness to stability limits and maintaining safety-margins. Its first application has been for reliable prevention of vertical displacement events (VDEs) by adjusting plasma elongation (κ) and the inner-gap between the plasma and inner-wall in response to real- time open-loop VDE growth rate (γ) estimators. VDEs were robustly prevented up to average open-loop growth rates of 800 rad/s with initial tunings, with only applying shape modification when near safety limits. Second, the disruption risk during fast, emergency shutdown after large tearing and locked modes can be significantly improved by transitioning to a limited topology during shutdown. More than 50% of emergency limited shutdowns after locked modes reach a final normalized current I N < 0.3 before terminating, scaling to the 3 MA ITER requirement. This is in contrast to diverted shutdowns, the majority of which disrupt at I N > 0.8. Despite improvements, these results highlight the critical importance of early prevention. Third, a novel emergency shut down method has been developed which excites instabilities to form a warm, helical core post-thermal quench. The current quench extends to ~100ms and avoids VDEs and runaway electron generation. Novel real-time machine learning disruption prediction has been integrated with the DIII-D proximity controller, and a real- time compatible multi-mode MHD spectroscopy technique has been developed. Results presented here were enabled by a focused effort, the Disruption Free Protocol, in DIII-D’s 2019-20 campaign to complement disruption prevention experiments with a large piggy-back program. In addition to testing novel techniques, it is estimated to have helped avoid 32 potential disruptions in piggyback operations with rapid, early shutdowns after large rotating n=1 or locked modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗