Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “asymptotic 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 163 records · Page 9

A nonlinear structural concept for drag-reducing compliant walls

Two mechanisms of drag reduction for flow over flat plates were investigated. The first mechanism employs Bushnell's hypothesis that compliant walls produce drag reduction by interfering with the formation of the turbulent spots in a turbulent boundary layer. It is shown that the amplitudes and frequencies of compliant wall motions for drag reduction might be achieved by using slightly curved walls and the resulting large amplitude motions of snap buckling. A simple structural model of an arch is used in the analysis, and an asymptotic method is developed. The required wall motions can be obtained by using materials like mylar. In addition, the delay of transition from laminar to turbulent flow by driven walls was studied for Poiseuille channel flow. The walls are driven by a periodic traveling wave. A significant increase in the transitional Reynolds number is obtained by appropriately prescribing the wavelength and phase velocity of the wall motion. Previously developed asymptotic methods are used in the analysis. Previously announced in STAR as N83-11061

Reiss, E. L.↗

Optimal guidance of aero-assisted transfer vehicles based on matched asymptotic expansions

This paper addresses a number of issues related to the matched asymptotic expansions analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is rigorously followed in developing a complete algebraic solution to the problem of inclination change with minimum energy loss. Repeated solution of these algebraic equations along the trajectory, treating each current state as an initial state, constitute a feedback guidance algorithm.

Calise, A. J.↗

Computer analysis of shells of revolution using asymptotic results

It is suggested that asymptotic results for the behavior of thin shells can be incorporated in a general computer code for the analysis of a complex shell structure. The advantage when compared to existing finite difference or finite element codes is a substantial reduction in computational labor with the capability of working to a specified level of accuracy. A reduction in user preparation time and dependance on user judgment is also gained, since mesh spacing can be internally generated. The general theory is described in this paper, as well as the implementation in the computer code FAST 1 (Functional Algorithm for Shell Theory) for the analysis of the general axisymmetric shell structure with axisymmetric loading.

Steele, C. R.↗

A Coupled Nonlinear Spacecraft Attitude Controller and Observer with an Unknown Gyro Misalignment and Gyro Bias

A nonlinear control scheme for attitude control of a spacecraft is combined with a nonlinear gyro misalignment and bias observer for the case of constant gyro misalignment and bias. A persistency of excitation analysis shows the observer gyro bias estimates converge to the true bias values exponentially fast. The convergence of the misalignment estimates is also presented. Then; the resulting coupled, closed loop dynamics are proven by a Lyapunov analysis to be globally stable, with asymptotically perfect tracking. The analysis is extended to consider the effects of noise in addition to the gyro misalignment and bias. A simulation of the proposed observer-controller design is given for a rigid spacecraft tracking a specified, time-varying attitude sequence to illustrate the theoretical claims.

Thienel, Julie↗

Asymptotic boundary conditions for dissipative waves - General theory

An outstanding issue in computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗

Asymptotic boundary conditions for dissipative waves: General theory

An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗

Performing Bayesian Analyses With AZURE2 Using BRICK: An Application to the 7 Be System

Phenomenological R-matrix has been a standard framework for the evaluation of resolved resonance cross section data in nuclear physics for many years. It is a powerful method for comparing different types of experimental nuclear data and combining the results of many different experimental measurements in order to gain a better estimation of the true underlying cross sections. Yet a practical challenge has always been the estimation of the uncertainty on both the cross sections at the energies of interest and the fit parameters, which can take the form of standard level parameters. Frequentist (χ 2 -based) estimation has been the norm. In this work, a Markov Chain Monte Carlo sampler, emcee, has been implemented for the R-matrix code AZURE2, creating the Bayesian R-matrix Inference Code Kit (BRICK). Bayesian uncertainty estimation has then been carried out for a simultaneous R-matrix fit of the 3 He (α,γ) 7 Be and 3 He (α,α) 3 He reactions in order to gain further insight into the fitting of capture and scattering data. Both data sets constrain the values of the bound state α-particle asymptotic normalization coefficients in 7 Be. The analysis highlights the need for low-energy scattering data with well-documented uncertainty information and shows how misleading results can be obtained in its absence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ground-based infrared observations of variable IRAS sources as candidates for late asymptotic giant branch stars

Analysis of the color distribution of OH/IR stars and IRAS low-resolution spectra class 30 objects suggests the presence of a well-defined evolutionary sequence which is populated by late asymptotic giant branch (LAGB) stars. The paper reports ground-based identification and infrared photometry of 10 candidates of news LAGB stars. None of the selected sources are found to have optical counterparts, and eight of the 10 show a strong 10-micron silicate absorption feature. It is suggested that these stars represent an invisible extension of extreme Mira variables and are some of the most evolved stars observed to date.

Kwok, Sun↗

Stability analysis of nonlinear systems by multiple time scaling

The asymptotic solution for the transient analysis of a general nonlinear system in the neighborhood of the stability boundary was obtained by using the multiple-time-scaling asymptotic-expansion method. The nonlinearities are assumed to be of algebraic nature. Terms of order epsilon to the 3rd power (where epsilon is the order of amplitude of the unknown) are included in the solution. The solution indicates that there is always a limit cycle which is stable (unstable) and exists above (below) the stability boundary if the nonlinear terms are stabilizing (destabilizing). Extension of the solution to include fifth order nonlinear terms is also presented. Comparisons with harmonic balance and with multiple-time-scaling solution of panel flutter equations are also included.

Morino, L.↗

Theoretical characterization of the potential energy surface for the reversible reaction H + O2 yields HO2(asterisk) yields OH + O. III - Computed points to define a global potential energy surface

Computed energies and geometries are reported which, combined with previously published calculations, permit a global representation of the potential energy surface for the reaction H + O2 yields HO2(asterisk) yields OH + O. These new calculations characterize the potential energy surface (PES) for all H atom angles of approach to O2 and for the region of the inner repulsive wall. The region of the T-shaped H-O2 exchange saddle point is connected with the constrained energy minimum (CEM) path, and a new collinear H-O2 exchange saddle point is characterized which lies only 9 kcal/mol above the H + O2 asymptote. A vibrational analysis which utilizes local cubic and quartic polynomial representations of the PES along the CEM path has been carried out. Optimal geometries, energies, and harmonic frequencies are reported along with anharmonic analyses for the O2 and OH asymptotes and for the HO2 minimum region of the PES.

Walch, Stephen P.↗

Generation of Tollmien-Schlichting waves by free-stream disturbances at low Mach numbers

The method of matched asymptotic expansions is used to study the generation of Tollmien-Schlichting waves by free stream disturbances incident on a flat plate boundary layer. Near the leading edge, the motion is governed by the unsteady boundary layer equation, while farther downstream it is governed (to lowest order) by the Orr-Sommerfeld equation with slowly varying coefficients. It is shown that there is an overlap domain where the Tollmien-Schlichting wave solutions to the Orr-Sommerfeld equation and an appropriate asymptotic solution of the unsteady boundary layer equation match, in the matched asymptotic expansion sense. The analysis leads to a set of scaling laws for the asymptotic structure of the unsteady boundary layer.

Goldstein, M. E.↗

Random motion analysis of flexible satellite structures

A singular perturbation formulation is used to study the responses of a flexible satellite when random measurement errors can occur. The random variables, at different instants of time, are assumed to be uncorrelated. Procedures for obtaining maxima and minima are described, and a variation of the linear method is developed for the formal solution of the two-point boundary-value problems represented by the variational equations. Random and deterministic solutions for the structural position coordinates are studied, and an analytic algorithm for treating the force equation of motion is developed. Since the random system indicated by the variational equation will always be asymptotically unstable, any analysis of stability must be based on the deterministic system.

Huang, T. C.↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Guidance law development for aeroassisted transfer vehicles using matched asymptotic expansions

This report addresses and clarifies a number of issues related to the Matched Asymptotic Expansion (MAE) analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is developed and rigorously followed to obtain a set of algebraic equations for the problem of inclination change with minimum energy loss. A detailed evaluation of the zeroth order guidance algorithm for aeroassisted orbit transfer is performed. It is shown that by exploiting the structure of the MAE solution procedure, the original problem, which requires the solution of a set of 20 implicit algebraic equations, can be reduced to a problem of 6 implicit equations in 6 unknowns. A solution that is near optimal, requires a minimum of computation, and thus can be implemented in real time and on-board the vehicle, has been obtained. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the zeroth order MAE problem to obtain the feedback controls. Finally, a general procedure is developed for constructing a MAE solution up to first order, of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is valid for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is shown to be inappropriate since it is not valid over a narrow range of the independent variable. That is, it is not uniformly valid. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished-one where the left boundary condition coincides with, or lies to the right of, the singular region, and another one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure where the obtained solution is uniformly valid to O(Epsilon(exp 2)). The potential application of this procedure to aeroassisted plane change is also described and partially evaluated.

Calise, Anthony J.↗

Pseudospectral collocation methods for fourth order differential equations

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Malek, Alaeddin↗

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗