Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,279 records · Page 71

An installed nacelle design code using a multiblock Euler solver. Volume 1: Theory document

An efficient multiblock Euler design code was developed for designing a nacelle installed on geometrically complex airplane configurations. This approach employed a design driver based on a direct iterative surface curvature method developed at LaRC. A general multiblock Euler flow solver was used for computing flow around complex geometries. The flow solver used a finite-volume formulation with explicit time-stepping to solve the Euler Equations. It used a multiblock version of the multigrid method to accelerate the convergence of the calculations. The design driver successively updated the surface geometry to reduce the difference between the computed and target pressure distributions. In the flow solver, the change in surface geometry was simulated by applying surface transpiration boundary conditions to avoid repeated grid generation during design iterations. Smoothness of the designed surface was ensured by alternate application of streamwise and circumferential smoothings. The capability and efficiency of the code was demonstrated through the design of both an isolated nacelle and an installed nacelle at various flow conditions. Information on the execution of the computer program is provided in volume 2.

Chen, H. C.↗

An installed nacelle design method using multiblock Euler solver

An efficient multiblock Euler design method was developed for designing a nacelle installed on geometrically complex airplane configurations. This approach employed a design module based on a direct iterative surface curvature method developed at NASA-Langley. A multiblock Euler flow solver was used for computing flow around complex geometries. The flow solver used a finite-volume formulation with explicit time-stepping to solve the Euler equations. It used a multiblock version of the multigrid method to accelerate the convergence of the calculations. The design module successively updated the surface geometry to reduce the difference between the computed and target pressure distributions. In the flow solver, the change in surface geometry was simulated by applying surface transpiration boundary conditions to avoid repeated grid generation during design iterations. Smoothness of the designed surface was ensured by alternate application of streamwise and circumferential smoothings.

Chen, H. C.↗

A Real-Time Method for Estimating Viscous Forebody Drag Coefficients

This paper develops a real-time method based on the law of the wake for estimating forebody skin-friction coefficients. The incompressible law-of-the-wake equations are numerically integrated across the boundary layer depth to develop an engineering model that relates longitudinally averaged skin-friction coefficients to local boundary layer thickness. Solutions applicable to smooth surfaces with pressure gradients and rough surfaces with negligible pressure gradients are presented. Model accuracy is evaluated by comparing model predictions with previously measured flight data. This integral law procedure is beneficial in that skin-friction coefficients can be indirectly evaluated in real-time using a single boundary layer height measurement. In this concept a reference pitot probe is inserted into the flow, well above the anticipated maximum thickness of the local boundary layer. Another probe is servomechanism-driven and floats within the boundary layer. A controller regulates the position of the floating probe. The measured servomechanism position of this second probe provides an indirect measurement of both local and longitudinally averaged skin friction. Simulation results showing the performance of the control law for a noisy boundary layer are then presented.

Whitmore, Stephen A.↗

Rough-Wall Channel Analysis Using Suboptimal Control Theory

The original aim of this work was to shed some light on the physics of turbulence over rough walls using large-eddy simulations and the suboptimal-control wall boundary conditions introduced by Nicoud et al. It was hoped that, if that algorithm was used to fit the mean velocity profile of the simulations to that of a rough-walled channel, instead of to a smooth one, the wall stresses introduced by the control algorithm would give some indication of what aspects of rough walls are most responsible for the modification of the flow in real turbulence. It was similarly expected that the structure of the resulting velocity fluctuations would share some of the characteristics of rough-walled flows, thus again suggesting what is intrinsic and what is accidental in the effect of geometric wall roughness. A secondary goal was to study the effect of 'unphysical' boundary conditions on the outside flow by observing how a relatively major change of the target velocity profile, and therefore presumably of the applied wall stresses, modifies properties such as the dominant length scales of the velocity fluctuations away from the wall. As will be seen below, this secondary goal grew more important during the course of the study, which was carried out during a short summer visit of the first two authors to the CTR. It became clear that there are open questions about the way in which the control algorithm models the boundary conditions, even for smooth walls, and that these questions make the physical interpretation of the results difficult. Considerable more work in that area seems to be needed before even relatively advanced large-eddy simulations, such as these, can be used to draw conclusions about the physics of wall-bounded turbulent flows. The numerical method is the same as in Nicoud et al. The modifications introduced in the original code are briefly described in section 2, but the original paper should be consulted for a full description of the algorithm. The results are presented in section 3 and summarized in section 4. The elementary properties of turbulence over rough walls which are used in the text have been taken from recent reviews such as Raupach et al. or Jimenez.

Flores, O.↗

Combining Disparate Measures of Metabolic Rate During Simulated Spacewalks

Scientists from NASA's Extravehicular Activities (EVA) Physiology Systems and Performance Project help design space suits for future missions, during which astronauts are expected to perform EVA activities on the Lunar or Martian surface. During an EVA, an astronaut's integrated metabolic rate is used to predict how much longer the activity can continue and still provide a safe margin of remaining consumables. For EVAs in the Apollo era, NASA physicians monitored live data feeds of heart rate, O 2 consumption, and liquid cooled garment (LCG) temperatures, which were subjectively combined or compared to produce an estimate of metabolic rate. But these multiple data feeds sometimes provided conflicting estimates of metabolic rate, making real-time calculations of remaining time difficult for physician/monitors. Currently, designs planned for the Constellation Program EVAs utilize an automated, but largely heuristic methodology for incorporating the above three measurements, plus an additional one - CO 2 production, ignoring data that appears in conflict; however a more rigorous model-based approach is desirable. In this study, we show how principal axis factor analysis, in combination with OLS regression and LOWESS smoothing can be used to estimate metabolic rate as a data-driven weighted average of heart rate, O 2 consumption, LCG temperature data, and CO 2 production. Preliminary results suggest less sensitivity to occasional spikes in observed data feeds, and reasonable within-subject reproducibility when applied to subsequent tasks. These methods do not require physician monitoring and as such can be automated in the electronic components of future space suits. With additional validation, our models show promise for increasing astronaut safety, while reducing the need for and potential errors associated with human monitoring of multiple systems.

Alan H Feiveson↗

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗

Real Surface Effects in Elastohydrodynamic Lubrication

Following a general discussion of real effects in EHL, consideration is given in detail to the role played by stochastic roughness superimposed on the nominal solid boundaries. Particular attention is given to the full-film EHL regime where incipient asperity contact bears a negligible fraction of the load. Flow, from which an averaged Reynolds equation can be formed, is nonetheless modified by the amplitude and texture of the roughness patterns. In describing these effects by means of certain well-defined flow factors, it is found that only two extra parameters are needed - the rms surface height and the ratio of the correlation lengths in the two principal roughness directions. In a lowest order perturbation expansion of these flow factors in powers of the ratio of rms roughness to nominal film thickness, no other properties of the roughness appear. In most cases of practical interest, the factors describing Poiseuille flow are separable into the sum of two single-surface flow factors which means that a combination of a single equivalent rough surface versus an ideal smooth surface can always be found. For flow entrained by slip velocity the factors separate instead into a difference of single-surface factors and it becomes significant which of the two surfaces carries the equivalent roughness. Results are discussed of some applications of the averaged Reynolds equation based on the flow factor method to the EHL line contact problem. Finally, the partial EHL regime is considered where comparable load fractions are carried by the hydrodynamic film and by incipient mechanical contact. An extension of the method into this regime by combining it with asperity contact models appears most encouraging.

John H Tripp↗

Photoelectrochemical fabrication of spectroscopic diffraction gratings

Photoelectrochemical etching was demonstrated as a means of fabricating a variety of periodic structures in semiconductors. The semiconductor is used as an electrode in an electrochemical cell, and is in contact with a liquid electrolyte. When the crystal is held at a positive voltage and illuminated, etching occurs in only the illuminated regions to a depth proportional to the illumination intensity and exposure time. In Phase 1, it was determined that diffraction gratings could be produced in gallium arsenide crystals by this method, using either a scanned focused laser beam or by uniform illumination of a ruling mask defined in metal or photoresist on the crystal surface. The latter approach was determined to produce V-grooves if the mask is oriented along certain crystallographic directions. These V-grooves were produced with an exceedingly smooth crystal morphology due to the highly controllable nature of the process and the mild electrolytes involved. The results form the basis for photoelectrochemical fabrication of deep, low pitch Eschelle gratings for use in high orders in NASA spectrographic instrumentation such as the Space Telescope Imaging Spectrograph.

Rauh, R. David↗

Modal Correction Method For Dynamically Induced Errors In Wind-Tunnel Model Attitude Measurements

This paper describes a method for correcting the dynamically induced bias errors in wind tunnel model attitude measurements using measured modal properties of the model system. At NASA Langley Research Center, the predominant instrumentation used to measure model attitude is a servo-accelerometer device that senses the model attitude with respect to the local vertical. Under smooth wind tunnel operating conditions, this inertial device can measure the model attitude with an accuracy of 0.01 degree. During wind tunnel tests when the model is responding at high dynamic amplitudes, the inertial device also senses the centrifugal acceleration associated with model vibration. This centrifugal acceleration results in a bias error in the model attitude measurement. A study of the response of a cantilevered model system to a simulated dynamic environment shows significant bias error in the model attitude measurement can occur and is vibration mode and amplitude dependent. For each vibration mode contributing to the bias error, the error is estimated from the measured modal properties and tangential accelerations at the model attitude device. Linear superposition is used to combine the bias estimates for individual modes to determine the overall bias error as a function of time. The modal correction model predicts the bias error to a high degree of accuracy for the vibration modes characterized in the simulated dynamic environment.

Buehrle, R. D.↗

Method For Identifying Sedimentary Bodies From Images And Its Application To Mineral Exploration

A method is disclosed for identifying a sediment accumulation from an image of a part of the earth's surface. The method includes identifying a topographic discontinuity from the image. A river which crosses the discontinuity is identified from the image. From the image, paleocourses of the river are identified which diverge from a point where the river crosses the discontinuity. The paleocourses are disposed on a topographically low side of the discontinuity. A smooth surface which emanates from the point is identified. The smooth surface is also disposed on the topographically low side of the point.

Wilkinson, Murray Justin↗

Images of the Extended Outer Regions of the Debris Ring around HR 4796 A

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (approx 9 M(sub Jup) at 0".5, approx 3 M(sub Jup) at 1").

Thalmann, C.↗

A method of calibrating wind velocity sensors with a modified gas flow calibrator

A procedure was described for calibrating air velocity sensors in the exhaust flow of a gas flow calibrator. The average velocity in the test section located at the calibrator exhaust was verified from the mass flow rate accurately measured by the calibrator's precision sonic nozzles. Air at elevated pressures flowed through a series of screens, diameter changes, and flow straighteners, resulting in a smooth flow through the open test section. The modified system generated air velocities of 2 to 90 meters per second with an uncertainty of about two percent for speeds below 15 meters per second and four percent for the higher speeds. Wind tunnel data correlated well with that taken in the flow calibrator.

Stump, H. P.↗

Marching iterative methods for the parabolized and thin layer Navier-Stokes equations

Downstream marching iterative schemes for the solution of the Parabolized or Thin Layer (PNS or TL) Navier-Stokes equations are described. Modifications of the primitive equation global relaxation sweep procedure result in efficient second-order marching schemes. These schemes take full account of the reduced order of the approximate equations as they behave like the SLOR for a single elliptic equation. The improved smoothing properties permit the introduction of Multi-Grid acceleration. The proposed algorithm is essentially Reynolds number independent and therefore can be applied to the solution of the subsonic Euler equations. The convergence rates are similar to those obtained by the Multi-Grid solution of a single elliptic equation; the storage is also comparable as only the pressure has to be stored on all levels. Extensions to three-dimensional and compressible subsonic flows are discussed. Numerical results are presented.

Israeli, M.↗

Marching iterative methods for the Parabolized and Thin Layer Navier-Stokes equations

Downstream marching iterative schemes for the solution of the Parabolized or Thin Layer (PNS or TL) Navier-Stokes equations are described. Modifications of the primitive equation global relaxation sweep procedure result in efficient second-order marching schemes. These schemes take full account of the reduced order of the approximate equations as they behave like the SLOR for a single elliptic equation. The improved smoothing properties permit the introduction of Multi-Grid acceleration. The proposed algorithm is essentially Reynolds number independent and therefore can be applied to the solution of the subsonic Euler equations. The convergence rates are similar to those obtained by the Multi-Grid solution of a single elliptic equation; the storage is also comparable as only the pressure has to be stored on all levels. Extensions to three-dimensional and compressible subsonic flows are discussed. Numerical results are presented.

Israeli, M.↗

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

Polymer thin films as universal substrates for extreme ultraviolet absorption spectroscopy of molecular transition metal complexes

Polystyrene and polyvinyl chloride thin films are explored as sample supports for extreme ultraviolet (XUV) spectroscopy of molecular transition metal complexes. Thin polymer films prepared by slip-coating are flat and smooth, and transmit much more XUV light than silicon nitride windows. Analytes can be directly cast onto the polymer surface or co-deposited within it. The M-edge XANES spectra (40–90 eV) of eight archetypal transition metal complexes (M = Mn, Fe, Co, Ni) are presented to demonstrate the versatility of this method. The films are suitable for pump/probe transient absorption spectroscopy, as shown by the excited-state spectra of Fe(bpy) 3 2+ in two different polymer supports.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of thermodynamic fluctuations and pattern size on the nucleation behavior during area selective deposition

In this work, the authors explore the impact that thermodynamic fluctuations have on the spontaneous appearance of defects during area selective deposition on patterned surfaces that are fully covered by adsorbates under local thermodynamic equilibrium, such as self-assembled monolayers. By using a simple lattice gas model for the adsorbed monolayer, the authors were able to track the spontaneous formation of defects in the monolayer as a function of the pattern width. The results indicate that, for pattern widths of the order of tens of nanometers, roughening effects at the pattern edge can be the leading source for the spontaneous appearance of nucleation defects. This leads to an enhancement of the density of defects that can be up to three orders of magnitude higher than those expected in uniform (not patterned) surfaces. The model also predicts a density of defects that is inversely proportional to the pattern width. Finally, if the dynamic nucleation of defects during area selective deposition is driven by thermodynamic fluctuations, the model predicts that the nucleation rate should be proportional to the total precursor fluency and independent of purge times. Moreover, a tight confinement of the monolayer through a high quality smooth interface in the patterned substrate and strong cohesive interactions between adsorbates should each contribute to a reduction of the overall defect density.

36 MATERIALS SCIENCE↗

Finite element analysis of notch behavior using a state variable constitutive equation

The state variable constitutive equation of Bodner and Partom was used to calculate the load-strain response of Inconel 718 at 649 C in the root of a notch. The constitutive equation was used with the Bodner-Partom evolution equation and with a second evolution equation that was derived from a potential function of the stress and state variable. Data used in determining constants for the constitutive models was from one-dimensional smooth bar tests. The response was calculated for a plane stress condition at the root of the notch with a finite element code using constant strain triangular elements. Results from both evolution equations compared favorably with the observed experimental response. The accuracy and efficiency of the finite element calculations also compared favorably to existing methods.

Dame, L. T.↗