Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elliptic problems”

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 829 records · Page 46

Integrated Aerodynamic and Control System Design of Oblique Wing Aircraft

An efficient high speed aircraft design must achieve a high lift to drag ratio at transonic and supersonic speeds. In 1952 Dr. R. T. Jones proved that for any flight Mach number minimum drag at a fixed lift is achieved by an elliptic wing planform with an appropriate oblique sweep angle. Since then, wind tunnel tests and numerical flow models have confirmed that the compressibility drag of oblique wing aircraft is lower than similar symmetrical sweep designs. At oblique sweep angles above thirty degrees the highly asymmetric planform gives rise to aerodynamic and inertia couplings which affect stability and degrade the aircraft's handling qualities. In the case of the NASA-Rockwell Oblique Wing Research Aircraft, attempts to improve the handling qualities by implementing a stability augmentation system have produced unsatisfactory results because of an inherent lack of controllability in the proposed design. The present work focuses on improving the handling qualities of oblique wing aircraft by including aerodynamic configuration parameters as variables in the control system synthesis to provide additional degrees of freedom with which to further decouple the aircraft's response. Handling qualities are measured using a quadratic cost function identical to that considered in optimal control problems, but the controller architecture is not restricted to full state feedback. An optimization procedure is used to simultaneously solve for the aircraft configuration and control gains which maximize a handling qualities measure, while meeting imposed constraints on trim. In some designs wing flexibility is also modeled and reduced order controllers are implemented. Oblique wing aircraft synthesized by this integrated design method show significant improvement in handling qualities when compared to the originally proposed closed loop aircraft. The integrated design synthesis method is then extended to show how handling qualities may be traded for other types of mission performance (drag, weight, etc.). Examples are presented which show how performance can be maximized while maintaining a desired level of handling quality.

Morris, Stephen James↗

Predicting Chandra CCD Degradation with the Chandra Radiation Model

Not long after launch of the Chandra X-Ray Observatory, it was discovered that the Advanced CCD Imaging Spectrometer (ACIS) detector was rapidly degrading due to radiation. Analysis by Chandra personnel showed that this degradation was due to 10w energy protons (100 - 200 keV) that scattered down the optical path onto the focal plane. In response to this unexpected problem, the Chandra Team developed a radiation-protection program that has been used to manage the radiation damage to the CCDs. This program consists of multiple approaches - scheduled sating of the ACIS detector from the radiation environment during passage through radiation belts, real-time monitoring of space weather conditions, on-board monitoring of radiation environment levels, and the creation of a radiation environment model for use in computing proton flux and fluence at energies that damage the ACIS detector. This radiation mitigation program has been very successful. The initial precipitous increase in the CCDs' charge transfer inefficiency (CTI) resulting from proton damage has been slowed dramatically, with the front-illuminated CCDS having an increase in CTI of only 2.3% per year, allowing the ASIS detector's expected lifetime to exceed requirements. This paper concentrates on one aspect of the Chandra radiation mitigation program, the creation of the Chandra Radiation Model (CRM). Because of Chandra's highly elliptical orbit, the spacecraft spends most of its time outside of the trapped radiation belts that present the severest risks to the ACIS detector. However, there is still a proton flux environment that must be accounted for in all parts of Chandra's orbit. At the time of Chandra's launch there was no engineering model of the radiation environment that could be used in the outer regions of the spacecraft's orbit, so the CRM was developed to provide the flux environment of 100 - 200 keV protons in the outer magnetosphere, magnetosheath, and solar wind regions of geospace. This presentation describes CRM, its role in Chandra operations, and its prediction of the ACIS CTI increase.

Minow, Joseph I.↗

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

CONTRIBUTIONS OF THE ELASTICITY TO THE PRECESSION OF A TWO-LAYER EARTH MODEL

We focus on the updating of a specific contribution to the precession of the equator in longitude, usually named as “second order.” It stems from the crossing of certain terms of the lunisolar gravitational potential. The IAU2006 precession theory assigns it the value of −46.8 mas/cy that was derived for a rigid Earth model. Instead of that model, we consider a two-layer Earth composed of an elastic mantle and a liquid core, working out the problem within the Hamiltonian framework developed by Getino and Ferrándiz. The targeted effect is obtained without further simplifying assumptions through Hori’s canonical perturbation method applied up to the second order of perturbation. On account of using a more realistic Earth model, the revised value of the second-order contribution is significantly changed and reaches −55.29 mas/cy. That variation of the second-order contribution is larger than other contributions included in IAU2006. It must be compensated with an increase of −8.51 mas/cy in the value of the lunisolar first-order component p{sub A}{sup ′} of the precession of the equator rate, which is derived from the total rate by subtracting the remaining contributions accounted for in IAU2006 precession. The updating of the second-order contribution implies that the p{sub A}{sup ′} parameter has to be changed, from 5040684.593 to 5040693.104 mas/cy in absence of potential revisions of other contributions. It entails a proportional variation of Earth’s dynamical ellipticity H{sub d}, for which the estimation associated with IAU2006, 0.00327379448, should be updated to 0.00327380001, about 1.7 ppm larger.

79 ASTRONOMY AND ASTROPHYSICS↗

Electromagnetically induced absorption and electromagnetically induced transparency for optical transitions F{sub g} → F{sub e} in the field of elliptically polarized waves

Nonlinear laser spectroscopy is considered in the scheme of two collinear waves with arbitrary elliptical polarizations. Emphasis is placed on investigating the nonlinear corrections in the absorption spectrum of one of the waves. The spontaneous transfer of low-frequency Zeeman coherence is shown to affect the sign of the subnatural-width resonance. For a closed transition, the direction of the resonance profile has been found to depend only on the angular momenta F{sub e} and F{sub g}. On this basis, a classification has been developed for various transitions by the direction of the subnatural-width resonance profile.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient Bayesian inference with latent Hamiltonian neural networks in No-U-Turn Sampling

When sampling for Bayesian inference, one popular approach in the computational field is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS), which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice when relying on computationally expensive forward models. We propose Latent Hamiltonian neural networks (L-HNNs) with HMC and NUTS for solving Bayesian inference problems. Once trained, L-HNNs do not require numerical gradients of the target density during sampling, and hence numerous evaluations of the forward computational model. Moreover, L-HNNs satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well-suited for use within HMC and NUTS because stationarity can be shown. We also propose the integration of L-HNNs in an online error monitoring scheme, in which numerical gradients of the target density are used for a few samples whenever the L-HNNs prediction errors are large. This online error monitor scheme prevents sample degeneracy in regions of low probability density and ensures robust uncertainty quantification. We demonstrate L-HNNs in NUTS with online error monitoring on several analytical examples involving complex, heavy-tailed, and high-local-curvature probability densities. We then demonstrate the applicability of L-HNNs in NUTS to two computational case studies, namely the Allen-Cahn stochastic partial differential equation and an elliptic partial differential equation with 25 and 50 inference parameters, respectively. Overall, the L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. In conclusion, compared to traditional NUTS, L-HNNs in NUTS with online error monitoring required 1–2 orders of magnitude fewer numerical gradients of the target density and improved the effective sample size (ESS) per gradient (which is a measure of both the sampling quality and the computational expense) by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

ABOUT THE LINEARITY OF THE COLOR-MAGNITUDE RELATION OF EARLY-TYPE GALAXIES IN THE VIRGO CLUSTER

We revisit the color-magnitude relation of Virgo Cluster early-type galaxies in order to explore its alleged nonlinearity. To this aim, we reanalyze the relation already published from data obtained within the ACS Virgo Cluster Survey of the Hubble Space Telescope and perform our own photometry and analysis of the images of 100 early-type galaxies observed as part of this survey. In addition, we compare our results with those reported in the literature from data of the Sloan Digital Sky Survey. We have found that when the brightest galaxies and untypical systems are excluded from the sample, a linear relation arises in agreement with what is observed in other groups and clusters. The central regions of the brightest galaxies also follow this relation. In addition, we notice that Virgo contains at least four compact elliptical galaxies besides the well-known object VCC 1297 (NGC 4486B). Their locations in the ({mu}{sub eff})-luminosity diagram define a trend different from that followed by normal early-type dwarf galaxies, setting an upper limit in effective surface brightness and a lower limit in the effective radius for their luminosities. Based on the distribution of different galaxy sub-samples in the color-magnitude and ({mu}{sub eff})-luminosity diagrams, we draw some conclusions on their formation and the history of their evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Vortex equations governing the fractional quantum Hall effect

An existence theory is established for a coupled non-linear elliptic system, known as “vortex equations,” describing the fractional quantum Hall effect in 2-dimensional double-layered electron systems. Via variational methods, we prove the existence and uniqueness of multiple vortices over a doubly periodic domain and the full plane. In the doubly periodic situation, explicit sufficient and necessary conditions are obtained that relate the size of the domain and the vortex numbers. For the full plane case, existence is established for all finite-energy solutions and exponential decay estimates are proved. Quantization phenomena of the magnetic flux are found in both cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Probabilistic Approach to Interior Regularity of Fully Nonlinear Degenerate Elliptic Equations in Smooth Domains

We consider the value function of a stochastic optimal control of degenerate diffusion processes in a domain D. We study the smoothness of the value function, under the assumption of the non-degeneracy of the diffusion term along the normal to the boundary and an interior condition weaker than the non-degeneracy of the diffusion term. When the diffusion term, drift term, discount factor, running payoff and terminal payoff are all in the class of C{sup 1,1}( D-bar ) , the value function turns out to be the unique solution in the class of C{sub loc}{sup 1,1}(D) Intersection C{sup 0,1}( D-bar ) to the associated degenerate Bellman equation with Dirichlet boundary data. Our approach is probabilistic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Post-Keplerian Corrections to the Orbital Periods of a Two-body System and Their Application to the Galactic Center

We perform detailed numerical analyses of the orbital motion of a test particle around a spinning primary, with the aim of investigating the possibility of using the post-Keplerian (pK) corrections to the orbiter’s periods (draconitic, anomalistic, and sidereal) as a further opportunity to perform new tests of post-Newtonian gravity. As a specific scenario, the S-stars orbiting the massive black hole (MBH) supposedly lurking in Sgr A* at the center of the Galaxy are adopted. We first study the effects of the pK Schwarzchild, Lense–Thirring, and quadrupole moment accelerations experienced by a target star for various possible initial orbital configurations. It turns out that the results of the numerical simulations are consistent with the analytical ones in the small eccentricity approximation for which almost all the latter ones were derived. For highly elliptical orbits, the sizes of the three pK corrections considered turn out to increase remarkably. The periods of the observed S2 and S0-102 stars as functions of the MBH’s spin axis orientation are considered as well. The pK accelerations lead to corrections of the orbital periods of the order of 1–100 days (Schwarzschild), 0.1–10 hr (Lense–Thirring), and 1–10{sup 3} s (quadrupole) for a target star with a = 300–800 au and e ≈ 0.8, which could be measurable with future facilities.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark halos and elliptical galaxies as marginally stable dynamical systems

The origin of equilibrium gravitational configurations is sought in terms of the stability of their trajectories, as described by the curvature of their Lagrangian configuration manifold of particle positions—a context in which subtle spurious effects originating from the singularity in the two-body potential become particularly clear. We focus on the case of spherical systems, which support only regular orbits in the collisionless limit, despite the persistence of local exponential instability of N-body trajectories in the anomalous case of discrete point particle representation even as N → ∞. When the singularity in the potential is removed, this apparent contradiction disappears. In the absence of fluctuations, equilibrium configurations generally correspond to positive scalar curvature and thus support stable trajectories. A null scalar curvature is associated with an effective, averaged equation of state describing dynamically relaxed equilibria with marginally stable trajectories. The associated configurations are quite similar to those of observed elliptical galaxies and simulated cosmological halos and are necessarily different from the systems dominated by isothermal cores, expected from entropy maximization in the context of the standard theory of violent relaxation. It is suggested that this is the case because a system starting far from equilibrium does not reach a 'most probable state' via violent relaxation, but that this process comes to an end as the system finds and (settles in) a configuration where it can most efficiently wash out perturbations. We explicitly test this interpretation by means of direct simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗

Formation Flying In Highly Elliptical Orbits Initializing the Formation

In this paper several methods are examined for initializing formations in which all spacecraft start in a common elliptical orbit subsequent to separation from the launch vehicle. The tetrahedron formation used on missions such as the Magnetospheric Multiscale (MMS), Auroral Multiscale Midex (AMM), and Cluster is used as a test bed Such a formation provides full three degrees-of-freedom in the relative motion about the reference orbit and is germane to several missions. The type of maneuver strategy that can be employed depends on the specific initial conditions of each member of the formation. Single-impulse maneuvers based on a Gaussian variation-of-parameters (VOP) approach, while operationally simple and intuitively-based, work only in a limited sense for a special class of initial conditions. These 'tailored' initial conditions are characterized as having only a few of the Keplerian elements different from the reference orbit. Attempts to achieve more generic initial conditions exceed the capabilities of the single impulse VOP. For these cases, multiple-impulse implementations are always possible but are generally less intuitive than the single-impulse case. The four-impulse VOP formalism discussed by Schaub is examined but smaller delta-V costs are achieved in our test problem by optimizing a Lambert solution.

Mailhe, Laurie↗

Mixing of an Airblast-atomized Fuel Spray Injected into a Crossflow of Air

The injection of a spray of fuel droplets into a crossflow of air provides a means of rapidly mixing liquid fuel and air for combustion applications. Injecting the liquid as a spray reduces the mixing length needed to accommodate liquid breakup, while the transverse injection of the spray into the air stream takes advantage of the dynamic mixing induced by the jet-crossflow interaction. The structure of the spray, formed from a model plain-jet airblast atomizer, is investigated in order to determine and understand the factors leading to its dispersion. To attain this goal, the problem is divided into the following tasks which involve: (1) developing planar imaging techniques that visualize fuel and air distributions in the spray, (2) characterizing the airblast spray without a crossflow, and (3) characterizing the airblast spray upon injection into a crossflow. Geometric and operating conditions are varied in order to affect the atomization, penetration, and dispersion of the spray into the crossflow. The airblast spray is first characterized, using imaging techniques, as it issues into a quiescent environment. The spray breakup modes are classified in a liquid Reynolds number versus airblast Weber number regime chart. This work focuses on sprays formed by the "prompt" atomization mode, which induces a well-atomized and well-dispersed spray, and which also produces a two-lobed liquid distribution corresponding to the atomizing air passageways in the injector. The characterization of the spray jet injected into the crossflow reveals the different processes that control its dispersion. Correlations that describe the inner and outer boundaries of the spray jet are developed, using the definition of a two-phase momentum-flux ratio. Cross-sections of the liquid spray depict elliptically-shaped distributions, with the exception of the finely-atomized sprays which show kidney-shaped distributions reminiscent of those obtained in gaseous jet in crossflow systems. A droplet trajectory analysis overpredicts the liquid mass penetration, and indicates a need for a more rigorous model to account for the three-dimensional mixing field induced by the jet-crossflow interaction. Nonetheless, the general procedures and criteria that are outlined can be used to efficiently assess and compare the quality of sprays formed under different conditions.

Leong, May Y.↗

Nonlinear birefringence in plasmas: Polarization dynamics, vector modulational instability, and vector solitons

The propagation of an elliptically polarized intense laser pulse in an unmagnetized collisionless uniform plasma is considered. A multiple scale perturbation theory is employed to show that in a weakly relativistic regime, evolution of the components of vector potential is governed by two coherently coupled nonlinear Schrödinger (NLS) equations. A set of equations describing the evolution of Stokes parameters is derived within the continuous wave approximation and the dynamics of components of field amplitude is studied. The polarization dynamics of a pulse is investigated by means of numerical solution of the coupled NLS equations. A detailed analysis of vector modulational instability is presented and the possibility for occurrence of various kinds of vector solitary waves is addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Blended-Wing-Body (BWB) Fuselage Structural Design for Weight Reduction

Structural analysis and design of efficient pressurized fuselage configurations for the advanced Blended-Wing-Body (BWB) flight vehicle is a challenging problem. Unlike a conventional cylindrical pressurized fuselage, stress level in a box type BWB fuselage is an order of magnitude higher, because internal pressure primarily results in bending stress instead of skin-membrane stress. In addition, resulting deformation of aerodynamic surface could significantly affect performance advantages provided by lifting body. The pressurized composite conformal multi-lobe tanks of X-33 type space vehicle also suffered from similar problem. In the earlier BWB design studies, Vaulted Ribbed Shell (VLRS), Flat Ribbed Shell (FRS); Vaulted shell Honeycomb Core (VLHC) and Flat sandwich shell Honeycomb Core (FLHC) concepts were studied. The flat and vaulted ribbed shell concepts were found most efficient. In a recent study, a set of composite sandwich panel and cross-ribbed panel were analyzed. Optimal values of rib and skin thickness, rib spacing, and panel depth were obtained for minimal weight under stress and buckling constraints. In addition, a set of efficient multi-bubble fuselage (MBF) configuration concept was developed. The special geometric configuration of this concept allows for balancing internal cabin pressure load efficiently, through membrane stress in inner-stiffened shell and inter-cabin walls, while the outer-ribbed shell prevents buckling due to external resultant compressive loads. The initial results from these approximate finite element analyses indicate progressively lower maximum stresses and deflections compared to the earlier study. However, a relative comparison of the FEM weight per unit floor area of the segment unit indicates that the unit weights are still relatively higher that the conventional B777 type cylindrical or A380 type elliptic fuselage design. Due to the manufacturing concern associated with multi-bubble fuselage, a Y braced box-type fuselage alternative with special resin-film injected (RFI) stitched carbon composite with foam-core was designed by Boeing under a NASA research contract for the 480 passenger version. It is shown that this configuration can be improved to a modified multi-bubble fuselage which has better stress distribution, for same material and dimension.

Mukhopadhyay, V.↗

Confocal Annular Josephson Tunnel Junctions with Large Eccentricity

Confocal annular Josephson tunnel junctions (CAJTJs), which are the natural generalization of the circular annular Josephson tunnel junctions, have a rich nonlinear phenomenology due to the intrinsic non-uniformity of their planar tunnel barrier delimited by two closely spaced confocal ellipses. In the presence of a uniform magnetic field in the barrier plane, the periodically changing width of the elliptical annulus generates a asymmetric double-well for a Josephson vortex trapped in a long and narrow CAJTJ. The preparation and readout of the vortex pinned in one of the two potential minima, which are important for the possible realization of a vortex qubit, have been numerically and experimentally investigated for CAJTJs with the moderate aspect ratio 2:1. In this work, we focus on the impact of the annulus eccentricity on the properties of the vortex potential profile and study the depinning mechanism of a fluxon in more eccentric samples with aspect ratio 4:1. We also discuss the effects of the temperature-dependent losses as well as the influence of the current and magnetic noise.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

On quasiconformal maps and semilinear equations in the plane

Assume that Ω is a domain in the complex plane ℂ and A(z) is a symmetric 2×2 matrix function with measurable entries, detA = 1; and such that 1/K|ξ|{sup 2} ≤ 〈A(z)ξ, ξ〉 ≤ K|ξ|{sup 2}, ξ ∈ ℝ{sup 2}, 1 ≤ K < ∞ . In particular, for semilinear elliptic equations of the form div (A(z)∇u(z)) = f(u(z)) in Ω; we prove a factorization theorem that asserts that every weak solution u to the above equation can be expressed as the composition u = To𝜔; where 𝜔 : Ω → G stands for a K−quasiconformal homeomorphism generated by the matrix function A(z); and T(w) is a weak solution of the semilinear equation ∇T(w) = J(w)f(T(w)) in G: Here, the weight J(w) is the Jacobian of the inverse mapping 𝜔{sup −1}: Similar results hold for the corresponding nonlinear parabolic and hyperbolic equations. Some applications of these results to anisotropic media are given.

97 MATHEMATICS AND COMPUTING↗