Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 631 records · Page 35

An experimental evaluation of metallic diaphragms for positive fuel expulsion in the atmosphere explorer hydrazine propulsion subsystem

Four Arde conospheroid metallic diaphragms were tested to evaluate their capability for use in the orbit adjust propulsion subsystem (OAPS) of the Explorer spacecraft. The diaphragms will be used for positive propellant expulsion and spacecraft center of mass control. A leak-free cycle life capability of nine reversals was demonstrated. The diaphragms rolled smoothly from ring to ring in a predictable manner on the first reversal. Varying amounts of diaphragm cocking and ring skipping were observed on subsequent reversals. The diaphragm pressure differential did not exceed 7 N/sq cm during any reversal. Cycle life capability, reversal mode, and pressure differential were not affected by sudden reversals, environmental tests, or 18,000 partial reversals. An expulsion efficiency of approximately 97 percent was demonstrated. The results of these tests show that metallic diaphragms can be used as an effective means of positive fuel expulsion; however, to achieve spacecraft center of mass control, the diaphragm must not be reversed prior to flight.

Woodruff, W. L.↗

Microwave inversion of leaf area and inclination angle distributions from backscattered data

The backscattering coefficient from a slab of thin randomly oriented dielectric disks over a flat lossy ground is used to reconstruct the inclination angle and area distributions of the disks. The disks are employed to model a leafy agricultural crop, such as soybeans, in the L-band microwave region of the spectrum. The distorted Born approximation, along with a thin disk approximation, is used to obtain a relationship between the horizontal-like polarized backscattering coefficient and the joint probability density of disk inclination angle and disk radius. Assuming large skin depth reduces the relationship to a linear Fredholm integral equation of the first kind. Due to the ill-posed nature of this equation, a Phillips-Twomey regularization method with a second difference smoothing condition is used to find the inversion. Results are obtained in the presence of 1 and 10 percent noise for both leaf inclination angle and leaf radius densities.

Lang, R. H.↗

Optical and structural characterization of InAs/GaAs quantum wells

Three InAs/GaAs single quantum wells of two-, three-, and four-monolayer thickness were characterized using optical and structural techniques. The results of high-resolution transmission electron (HRTEM) microscopy and optical studies which combine absorption, photoluminescence (PL), photoreflectance, and cathodoluminescence are presented. Using the polarization modulated absorptance technique, we observed two absorption features in our samples at 77 K. On the basis of their polarization properties and comparison with an envelope function calculation, these structures are assigned to transitions between the confined heavy-hole and confined and unconfined electron levels. Photoreflectance spectra of the three-monolayer sample in 77-300 K range show only the fundamental quantum well transition. The temperature dependence of this transition is approximately linear with a slope of 2.2 x 10 exp -4 eV/K, which is significantly lower than in both constituent materials. Comparison to the absorption data reveals that the PL spectra are affected by the carrier diffusion and therefore do not provide direct measure of the exciton density of states. The HRTEM images indicate that, while the interfaces of the two-monolayer sample are smooth and the well thickness is uniform, the four-monolayer sample has uneven interfaces and contains domains of two, three, and four monolayers.

Ksendzov, A.↗

Summation by parts, projections, and stability

We have derived stability results for high-order finite difference approximations of mixed hyperbolic-parabolic initial-boundary value problems (IBVP). The results are obtained using summation by parts and a new way of representing general linear boundary conditions as an orthogonal projection. By slightly rearranging the analytic equations, we can prove strict stability for hyperbolic-parabolic IBVP. Furthermore, we generalize our technique so as to yield strict stability on curvilinear non-smooth domains in two space dimensions. Finally, we show how to incorporate inhomogeneous boundary data while retaining strict stability. Using the same procedure one can prove strict stability in higher dimensions as well.

Olsson, Pelle↗

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗

Optimal Polynomial Smoothers and One‐Sided V‐Cycles for Poisson Problems

The solution to the Poisson equation arising from the spectral element discretization of the incompressible Navier‐Stokes equations needs robust preconditioning strategies. One such strategy is multigrid. To realize the potential of multigrid methods, effective smoothing strategies are needed. Chebyshev polynomial smoothers, in conjunction with pointwise Jacobi or additive Schwarz methods (ASMs), prove to be an effective smoother. Other polynomial smoothers, however, may provide superior convergence to the multigrid preconditioner. The authors compare the standard Chebyshev polynomial smoothers to both the novel fourth‐kind Chebyshev polynomial smoothers proposed by Lottes as well as smoothers based on the polynomial of best uniform approximation to as proposed by Kraus, Vassilevski, and Zikatanov. At the cost of symmetry, further improvements may be made. For example, a order polynomial smoother on both sides of the V‐cycle may be substituted with an order polynomial smoother on one side at no additional cost. The choice of omitting the postsmoother in favor of higher‐order polynomial presmoothing is advantageous in cases where the multigrid approximation property constant is large. The authors consider a 2D model problem based on finite differences to motivate the choice of polynomial smoother, order, and whether to apply postsmoothing for the target application of high‐order ‐geometric multigrid methods for GPU architectures. Results from both domains demonstrate the substantial improvement of these approaches over the standard Chebyshev polynomial smoother with a symmetric V‐cycle.

97 MATHEMATICS AND COMPUTING↗

An efficient second-order adaptive procedure for inserting CAD geometries into hexahedral meshes using volume fractions

Here, this paper is concerned with inserting three-dimensional computer-aided design (CAD) geometries into meshes composed of hexahedral elements using a volume fraction representation. An adaptive procedure for doing so is presented. The procedure consists of two steps. The first step performs spatial acceleration using a k-d tree. The second step involves subdividing individual hexahedra in an adaptive mesh refinement (AMR)-like fashion and approximating the CAD geometry linearly (as a plane) at the finest subdivision. The procedure requires only two geometric queries from a CAD kernel: determining whether or not a queried spatial coordinate is inside or outside the CAD geometry and determining the closest point on the CAD geometry’s surface from a given spatial coordinate. We prove that the procedure is second-order accurate for sufficiently smooth geometries and sufficiently refined background meshes. We demonstrate the expected order of accuracy is achieved with several verification tests and illustrate the procedure’s effectiveness for several exemplar CAD geometries.

Adaptive↗

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗

Development of silicon nitride of improved toughness

The application of reaction sintered Si2N4 energy absorbing surface layers to hot-pressed Si3N4 was investigated. The surface layer was formed by in-place nitridation of silicon powder. It was found that reaction sintered Si3N4 layers of 1 mm thickness, fabricated from either -100, +200, -200, or -325 mesh Si powder and nitrided in 96% N2/4% H2 so that approximately 20-25 vol % unnitrided Si remained in the layer, resulted in a sevenfold increase in ballistic impact resistance of a 0.64 cm thick hot-pressed SI3N4 substrate from RT 1370 C. Both NC-132 SI3N4, with MgO additive, and NCX-34 Si3N4, with Y2O3 additive, were evaluated as substrate material. The finer grain size -200 and -325 mesh nitrided Si layers were for their smoothness and relatively high density. It was found that nitriding in N2/H2 mixtures, rather than pure N2, resulted in a microstructure that did not substantially degrade the strength of the hot-pressed Si3N4 substrate. Thermal cycling tests on the RSSN/HPSN combinations from 200 C to 1370 C for 75 cycles in air did not degrade the impact resistance nor the interfacial bonding, although a large amount of internal silica formation occurred within the RSSN layer. Mach 0.8, 5 hr, hot gas erosion tests showed no surface recession of RSSN layers at 1200 C and slight surface recession at 1370 C.

Brennan, J. J.↗

Digital control of highly augmented combat rotorcraft

Proposed concepts for the next generation of combat helicopters are to be embodied in a complex, highly maneuverable, multiroled vehicle with avionics systems. Single pilot and nap-of-the-Earth operations require handling qualities which minimize the involvement of the pilot in basic stabilization tasks. To meet these requirements will demand a full authority, high-gain, multimode, multiply-redundant, digital flight-control system. The gap between these requirements and current low-authority, low-bandwidth operational rotorcraft flight-control technology is considerable. This research aims at smoothing the transition between current technology and advanced concept requirements. The state of the art of high-bandwidth digital flight-control systems are reviewed; areas of specific concern for flight-control systems of modern combat are exposed; and the important concepts are illustrated in design and analysis of high-gain, digital systems with a detailed case study involving a current rotorcraft system. Approximate and exact methods are explained and illustrated for treating the important concerns which are unique to digital systems.

Tischler, Mark B.↗

Optimum filters and smoothers design for carrier phase and frequency tracking

The report presents the application of fixed lag smoothing algorithms to the problem of estimation of the phase and frequency of a sinusoidal carrier received in the presence of process noise and additive observation noise. A suboptimal structure consists of a phase-locked loop (PLL) followed by post-loop correction to the phase and frequency estimates. When the PLL is operating under high signal-to-noise ratio, the phase detector is approximately linear, and the smoother equations then correspond to the optimal linear equations for an equivalent linear signal model. The performance of such a smoother can be predicted by linear filtering theory. However, if the PLL is operating near the threshold region of the signal-to-noise ratio, the phase detector cannot be assumed to be linear. Then the actual performance of the smoother can significantly differ from that predicted by linear theory. In this report we present both the theoretical and simulated performance of such smoothers derived on the basis of various models for the phase and frequency processes.

Kumar, Rajendra↗

Velocity distributions among colliding asteroids

The probability distribution for impact velocities between two given asteroids is wide, non-Gaussian, and often contains spikes according to our new method of analysis in which each possible orbital geometry for collision is weighted according to its probability. An average value would give a good representation only if the distribution were smooth and narrow. Therefore, the complete velocity distribution we obtain for various asteroid populations differs significantly from published histograms of average velocities. For all pairs among the 682 asteroids in the main-belt with D greater than 50 km, we find that our computed velocity distribution is much wider than previously computed histograms of average velocities. In this case, the most probable impact velocity is approximately 4.4 km/sec, compared with the mean impact velocity of 5.3 km/sec. For cases of a single asteroid (e.g., Gaspra or Ida) relative to an impacting population, the distribution we find yields lower velocities than previously reported by others. The width of these velocity distributions implies that mean impact velocities must be used with caution when calculating asteroid collisional lifetimes or crater-size distributions. Since the most probable impact velocities are lower than the mean, disruption events may occur less frequently than previously estimated. However, this disruption rate may be balanced somewhat by an apparent increase in the frequency of high-velocity impacts between asteroids. These results have implications for issues such as asteroidal disruption rates, the amount/type of impact ejecta available for meteoritical delivery to the Earth, and the geology and evolution of specific asteroids like Gaspra.

Bottke, William F., Jr.↗

Impact of comet Shoemaker-Levy 9 on Jupiter

We have employed three-dimensional numerical simulations of the impact of Comet Shoemaker-Levy 9 (SL9) on Jupiter and the resulting vapor plume expansion using the smoothed particle hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disc. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km; however, no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Takata, Toshiko↗

srMO-BO-3GP: A sequential regularized multi-objective Bayesian optimization for constrained design applications using an uncertain Pareto classifier

Bayesian optimization (BO) is an efficient and flexible global optimization framework that is applicable to a very wide range of engineering applications. To leverage the capability of the classical BO, many extensions, including multi-objective, multi-fidelity, parallelization, and latent-variable modeling, have been proposed to address the limitations of the classical BO framework. In this work, we propose a novel multi-objective BO formalism, called srMO-BO-3GP, to solve multi-objective optimization problems in a sequential setting. Three different Gaussian processes (GPs) are stacked together, where each of the GPs is assigned with a different task. The first GP is used to approximate a single-objective computed from the multi-objective definition, the second GP is used to learn the unknown constraints, and the third one is used to learn the uncertain Pareto frontier. At each iteration, a multi-objective augmented Tchebycheff function is adopted to convert multi-objective to single-objective, where the regularization with a regularized ridge term is also introduced to smooth the single-objective function. Finally, we couple the third GP along with the classical BO framework to explore the convergence and diversity of the Pareto frontier by the acquisition function for exploitation and exploration. The proposed framework is demonstrated using several numerical benchmark functions, as well as a thermomechanical finite element model for flip-chip package design optimization.

42 ENGINEERING↗

The time development of a blast wave with shock heated electrons

Accurate approximations are presented for the time development of both edge conditions and internal structures of a blast wave with shock heated electrons, and equal ion and electron temperatures at the shock. The cases considered evolve in cavities with power law ambient densities (including the uniform ambient density case) and have negligible external pressure. Account is taken of possible saturation of the thermal conduction flux. The structures evolve smoothly to the adiabatic structures.

Edgar, R. J.↗

The time development of a blast wave with shock-heated electrons

Accurate approximations are presented for the time development of both edge conditions and internal structures of a blast wave with shock heated electrons, and equal ion and electron temperatures at the shock. The cases considered evolve in cavities with power law ambient densities (including the uniform ambient density case) and have negligible external pressure. Account is taken of possible saturation of the thermal conduction flux. The structures evolve smoothly to the adiabatic structures.

Edgar, R. J.↗

Formation, structure, and stability of MHD intermediate shocks

It was recently shown by Wu (1987) that intermediate shocks are admissible and can be formed through nonlinear wave steepening from continuous waves. In this paper, the formation, structure, and stability of intermediate shocks in dissipative MHD are considered in detail. The differences between the conventional theory and the present one are pointed out and clarified. It is shown that all four types of intermediate shocks can be formed from smooth waves. It is also shown that there are free parameters in the structure of the intermediate shocks, and that these parameters are related to the shock stability. In addition, the paper shows that a rotational discontinuity can not exist with finite width, indicates how this is related to the existence of time-dependent intermediate shocks, and shows why the conventional theory is not a good approximation to dissipative MHD solutions whenever there is rotation in magnetic field.

Wu, C. C.↗

Impact on comet Shoemaker-Levy 9 on Jupiter

Three-dimensional numerical simulations of the impact of Comet Shoemaker - Levy 9 on Jupiter and the resulting vapor plume expansion were conducted using the Smoothed Particle Hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 km to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disk. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km, but no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Ahrens, Thomas J.↗