SEARCH · Engineering Papers

Results for “Perturbative 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Computation of turbulent boundary layers on curved surfaces, 1 June 1975 - 31 January 1976

An accurate method was developed for predicting effects of streamline curvature and coordinate system rotation on turbulent boundary layers. A new two-equation model of turbulence was developed which serves as the basis of the study. In developing the new model, physical reasoning is combined with singular perturbation methods to develop a rational, physically-based set of equations which are, on the one hand, as accurate as mixing-length theory for equilibrium boundary layers and, on the other hand, suitable for computing effects of curvature and rotation. The equations are solved numerically for several boundary layer flows over plane and curved surfaces. For incompressible boundary layers, results of the computations are generally within 10% of corresponding experimental data. Somewhat larger discrepancies are noted for compressible applications.

Wilcox, D. C.↗

NASA NTRS1988

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

NASA NTRS1986

A generalization of the method of harmonic balance

A procedure is provided for generalizing the method of harmonic balance to obtain higher-order approximations to the periodic solutions of differential equations for two systems, one a conservative system and the other nonconservative. The procedure is currently being applied to investigate the possible solution behaviors of singular, nonlinear oscillators, where the usual perturbation methods do not work.

Mickens, R. E.↗

NASA NTRS1973

Program for the solution of multipoint boundary value problems of quasilinear differential equations

Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.

Source record↗

NASA NTRS1972

A finite element procedure for nonlinear prebuckling and initial postbuckling analysis

A procedure cast in a form appropriate to the finite element method is presented for geometrically nonlinear prebuckling and postbuckling structural analysis, including the identification of snap-through type of buckling. The principal features of this procedure are the use of direct iteration for solution of the nonlinear algebraic equations in the prebuckling range, an interpolation scheme for determination of the initial bifurcation point, a perturbation method in definition of the load-displacement behavior through the postbuckling regime, and extrapolation in determination of the limit point for snap-through buckling. Three numerical examples are presented in illustration of the procedure and in comparison with alternative approaches.

Mau, S. T.↗

NASA NTRS1985

An Overview of Latest Model Reduction and Control Methods of Large Flexible Space Structures

The latest trends and theoretical developments involved with the modeling and control of Large Flexible Space Structures (LFSS) are described. The paper addresses first the basic problems, characteristics, and difficulties inherent in modeling and control of LFSS. Major sources of difficulties and errors are the stiffness and damping operators of the dynamic model. Extensions of Linear Quadratic Gaussian (LQG) theory as applied to LFSS are presented, including frequency-shaped cost functionals and perturbation methods. The minimum data/maximum entropy approach which uses a stochastic design model to overcome difficulties found in the LQG-based methods is described. Latest trends in system theory including balanced realization and singular-value analysis are used to determine reduced order controllers and models. Ad hoc methods such as component cost analysis and modal cost analysis are discussed in context with the closed-loop reduction problem of controller order versus performance. The minimum data/maximum entropy approach also addresses controller order versus performance. Those areas of control science and large scale systems that appear to have an important role in understanding and solving LFSS modeling and control are also identified.

Santiago, J. M.↗

NASA NTRS1989

On the method of matched asymptotic expansions

The present evaluation of the method of asymptotic expansions (MAE) indicates that the various terms of the common solution of MAE can be generated as polynomials in stretched variables, without actually solving them from the outer solution, as is currently the practice. It is also noted that the common solution of the MAE and the intermediate solution of the singular-perturbation method are the same; these methods therefore yield identical results for a certain class of problems. Two illustrative problems are treated.

Naidu, D. S.↗

DOE OSTI2022

Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman’s, Pearson’s and Kendall’s tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗

NASA NTRS1948

On similarity rules for transonic flows

A method used by Tsien to derive similarity rules for hypersonic flows is utilized to derive Von Karman's similarity rules for transonic flows. At the lower limit of the transonic region of flow the theory yields a formula for the critical stream Mach numbers of a given family of symmetrical profiles. It is further shown that this formula can also be obtained by means of the Prandtl-Glauert small-perturbation method. Investigation of the behavior of the similarity parameter in the region where the thickness coefficient approaches zero and the critical stream Mach number approaches unity shows that it possesses a limiting value characteristic of the prescribed family of shapes.

Kaplan, Carl↗

NASA NTRS1992

Backscattering from a randomly rough dielectric surface

A backscattering model for scattering from a randomly rough dielectric surface is developed based on an approximate solution of a pair of integral equations for the tangential surface fields. Both like and cross-polarized scattering coefficients are obtained. It is found that the like polarized scattering coefficients contain two types of terms: single scattering terms and multiple scattering terms. The single scattering terms in like polarized scattering are shown to reduce the first-order solutions derived from the small perturbation method when the roughness parameters satisfy the slightly rough conditions. When surface roughnesses are large but the surface slope is small, only a single scattering term corresponding to the standard Kirchhoff model is significant. If the surface slope is large, the multiple scattering term will also be significant. The cross-polarized backscattering coefficients satisfy reciprocity and contain only multiple scattering terms. The difference between vertical and horizontal scattering coefficients is found to increase with the dielectric constant and is generally smaller than that predicted by the first-order small perturbation model. Good agreements are obtained between this model and measurements from statistically known surfaces.

Fung, Adrian K.↗

NASA NTRS1996

A Numerical Simulation of Scattering from One-Dimensional Inhomogeneous Dielectric Random Surfaces

In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.

Sarabandi, Kamal↗

NASA NTRS2001

A Study of Green's Function Methods Applied to Space Radiation Protection

The purpose of this research was to study the propagation of galactic ions through various materials. Galactic light ions result from the break up of heavy ion particles and their propagation through materials is modeled using the one-dimensional Boltzmann equation. When ions enter materials there can occur (i) the interaction of ions with orbital electrons which causes ionization within the material and (ii) ions collide with atoms causing production of secondary particles which penetrate deeper within the material. These processes are modeled by a continuum model. The basic idea is to place a control volume within the material and examine the change in ion flux across this control volume. In this way on can derive the basic equations for the transport of light and heavy ions in matter. Green's function perturbation methods can then be employed to solve the resulting equations using energy dependent nuclear cross sections.

Heinbockel, John H.↗

NASA NTRS1984

Rational approach for assumed stress finite elements

A new method for the formulation of hybrid elements by the Hellinger-Reissner principle is established by expanding the essential terms of the assumed stresses as complete polynomials in the natural coordinates of the element. The equilibrium conditions are imposed in a variational sense through the internal displacements which are also expanded in the natural co-ordinates. The resulting element possesses all the ideal qualities, i.e. it is invariant, it is less sensitive to geometric distortion, it contains a minimum number of stress parameters and it provides accurate stress calculations. For the formulation of a 4-node plane stress element, a small perturbation method is used to determine the equilibrium constraint equations. The element has been proved to be always rank sufficient.

Pian, T. H. H.↗

NASA NTRS2010

Technique for Performing Dielectric Property Measurements at Microwave Frequencies

A paper discusses the need to perform accurate dielectric property measurements on larger sized samples, particularly liquids at microwave frequencies. These types of measurements cannot be obtained using conventional cavity perturbation methods, particularly for liquids or powdered or granulated solids that require a surrounding container. To solve this problem, a model has been developed for the resonant frequency and quality factor of a cylindrical microwave cavity containing concentric cylindrical samples. This model can then be inverted to obtain the real and imaginary dielectric constants of the material of interest. This approach is based on using exact solutions to Maxwell s equations for the resonant properties of a cylindrical microwave cavity and also using the effective electrical conductivity of the cavity walls that is estimated from the measured empty cavity quality factor. This new approach calculates the complex resonant frequency and associated electromagnetic fields for a cylindrical microwave cavity with lossy walls that is loaded with concentric, axially aligned, lossy dielectric cylindrical samples. In this approach, the calculated complex resonant frequency, consisting of real and imaginary parts, is related to the experimentally measured quantities. Because this approach uses Maxwell's equations to determine the perturbed electromagnetic fields in the cavity with the material(s) inserted, one can calculate the expected wall losses using the fields for the loaded cavity rather than just depending on the value of the fields obtained from the empty cavity quality factor. These additional calculations provide a more accurate determination of the complex dielectric constant of the material being studied. The improved approach will be particularly important when working with larger samples or samples with larger dielectric constants that will further perturb the cavity electromagnetic fields. Also, this approach enables the ability to have a larger sample of interest, such as a liquid or powdered or granulated solid, inside a cylindrical container.

Barmatz, Martin B.↗

NASA NTRS1972

Perturbation and harmonic balance methods for nonlinear panel flutter.

A systematic way of applying both perturbation methods and harmonic balance methods to nonlinear panel flutter problems is developed here. Results obtained by both these methods for two-dimensional simply supported and three-dimensional clamped-clamped plates with six modes agree well with those obtained by the straightforward direct integration method, yet require less computer time and provide better insight into the solutions. Effects of viscoelastic structural damping on the flutter stability boundary are generally found to be destabilizing and the postflutter behavior becomes more explosive. The methods developed here may be of interest in related vibration problems.

Kuo, C.-C.↗

DOE OSTI2022

Recombination between 13 C and 2 H to Form Acetylide ( 13 C 2 2 H – ) Probes Nanoscale Interactions in Lipid Bilayers via Dynamic Secondary Ion Mass Spectrometry: Cholesterol and GM 1 Clustering

Although it is thought that there is lateral heterogeneity of lipid and protein components within biological membranes, probing this heterogeneity has proven challenging. The difficulty in such experiments is due to both the small length scale over which such heterogeneity can occur, and the significant perturbation resulting from fluorescent or spin labeling on the delicate interactions within bilayers. Atomic recombination during dynamic nanoscale secondary ion imaging mass spectrometry (NanoSIMS) is a non-perturbative method for examining nanoscale bilayer interactions. Atomic recombination is a variation on conventional NanoSIMS imaging, whereby an isotope on one molecule combines with a different isotope on another molecule during the ionization process, forming an isotopically enriched polyatomic ion in a distance-dependent manner. In this work, we show that the recombinant ion, 13 C 2 2 H – , is formed in high yield from 13 C- and 2 H-labeled lipids. The low natural abundance of triply labeled acetylide also makes it an ideal ion to probe GM 1 clusters in model membranes and the effects of cholesterol on lipid–lipid interactions. We find evidence supporting the cholesterol condensation effect as well as the presence of nanoscale GM 1 clusters in model membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electromagnetic formalism

This is a companion paper to the previous work [P. H. Yoon, Phys. Plasmas 31, 032309 (2024)] in which the nonlinear susceptibilities of weakly turbulent magnetized plasma are derived under a simplifying assumption of electrostatic interaction. The present paper extends the analysis to a general situation of electromagnetic interaction. The main novelty of the previous and present papers is that by employing the Bessel function addition theorem, the mathematical definitions for the susceptibilities are substantially simplified, a procedure that has not been discussed in the existing literature. In the present paper, a full set of Maxwell’s equations are considered in conjunction with the nonlinear Vlasov equation, which is solved by a perturbative method. The result is a fully general nonlinear susceptibility, given in tensorial form, which is applicable for weakly turbulent magnetized plasmas.

Materials Science↗