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 703 records · Page 39

Micropatterning tractional forces in living cells

Here we describe a method for quantifying traction in cells that are physically constrained within micron-sized adhesive islands of defined shape and size on the surface of flexible polyacrylamide gels that contain fluorescent microbeads (0.2-microm diameter). Smooth muscle cells were plated onto square (50 x 50 microm) or circular (25- or 50-microm diameter) adhesive islands that were created on the surface of the gels by applying a collagen coating through microengineered holes in an elastomeric membrane that was later removed. Adherent cells spread to take on the size and shape of the islands and cell tractions were quantitated by mapping displacement fields of the fluorescent microbeads within the gel. Cells on round islands did not exhibit any preferential direction of force application, but they exerted their strongest traction at sites where they formed protrusions. When cells were confined to squares, traction was highest in the corners both in the absence and presence of the contractile agonist, histamine, and cell protrusions were also observed in these regions. Quantitation of the mean traction exerted by cells cultured on the different islands revealed that cell tension increased as cell spreading was promoted. These results provide a mechanical basis for past studies that demonstrated a similar correlation between spreading and growth within various anchorage-dependent cells. This new approach for analyzing the spatial distribution of mechanical forces beneath individual cells that are experimentally constrained to defined sizes and shapes may provide additional insight into the biophysical basis of cell regulation. Copyright 2002 Wiley-Liss, Inc.

NASA Discipline Cell Biology↗

Microwave Treatment for Cardiac Arrhythmias

NASA seeks to transfer the NASA developed microwave ablation technology, designed for the treatment of ventricular tachycardia (irregular heart beat), to industry. After a heart attack, many cells surrounding the resulting scar continue to live but are abnormal electrically; they may conduct impulses unusually slowly or fire when they would typically be silent. These diseased areas might disturb smooth signaling by forming a reentrant circuit in the muscle. The objective of microwave ablation is to heat and kill these diseased cells to restore appropriate electrical activity in the heart. This technology is a method and apparatus that provides for propagating microwave energy into heart tissues to produce a desired temperature profile therein at tissue depths sufficient for thermally ablating arrhythmogenic cardiac tissue while preventing excessive heating of surrounding tissues, organs, and blood. A wide bandwidth double-disk antenna is effective for this purpose over a bandwidth of about six gigahertz. A computer simulation provides initial screening capabilities for an antenna such as antenna, frequency, power level, and power application duration. The simulation also allows optimization of techniques for specific patients or conditions. In comparison with other methods that involve direct-current pulses or radio frequencies below 1 GHz, this method may prove more effective in treating ventricular tachycardia. This is because the present method provides for greater control of the location, cross-sectional area, and depth of a lesion via selection of the location and design of the antenna and the choice of microwave power and frequency.

Hernandez-Moya, Sonia↗

Modeling Optical Properties of Mineral Aerosol Particles by Using Nonsymmetric Hexahedra

We explore the use of nonsymmetric geometries to simulate the single-scattering properties of airborne dust particles with complicated morphologies. Specifically, the shapes of irregular dust particles are assumed to be nonsymmetric hexahedra defined by using the Monte Carlo method. A combination of the discrete dipole approximation method and an improved geometric optics method is employed to compute the single-scattering properties of dust particles for size parameters ranging from 0.5 to 3000. The primary optical effect of eliminating the geometric symmetry of regular hexahedra is to smooth the scattering features in the phase function and to decrease the backscatter. The optical properties of the nonsymmetric hexahedra are used to mimic the laboratory measurements. It is demonstrated that a relatively close agreement can be achieved by using only one shape of nonsymmetric hexahedra. The agreement between the theoretical results and their measurement counterparts can be further improved by using a mixture of nonsymmetric hexahedra. It is also shown that the hexahedron model is much more appropriate than the "equivalent sphere" model for simulating the optical properties of dust particles, particularly, in the case of the elements of the phase matrix that associated with the polarization state of scattered light.

Bi, Lei↗

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗

Adjoint methods for aerodynamic wing design

A model inverse design problem is used to investigate the effect of flow discontinuities on the optimization process. The optimization involves finding the cross-sectional area distribution of a duct that produces velocities that closely match a targeted velocity distribution. Quasi-one-dimensional flow theory is used, and the target is chosen to have a shock wave in its distribution. The objective function which quantifies the difference between the targeted and calculated velocity distributions may become non-smooth due to the interaction between the shock and the discretization of the flowfield. This paper offers two techniques to resolve the resulting problems for the optimization algorithms. The first, shock-fitting, involves careful integration of the objective function through the shock wave. The second, coordinate straining with shock penalty, uses a coordinate transformation to align the calculated shock with the target and then adds a penalty proportional to the square of the distance between the shocks. The techniques are tested using several popular sensitivity and optimization methods, including finite-differences, and direct and adjoint discrete sensitivity methods. Two optimization strategies, Gauss-Newton and sequential quadratic programming (SQP), are used to drive the objective function to a minimum.

Grossman, Bernard↗

Effects of turbulence in a planetary atmosphere on radio occultation

A formulation is given for the correlation functions for the log-amplitude and phase fluctuations of both spherical and plane waves propagating in a turbulent medium whose correlation function for refractive index fluctuations is described by the product of a function of the average coordinate and a function of the difference coordinate. The results are applied to radio occultation of a flyby space probe by the atmosphere of Venus, assuming that the turbulence in the atmosphere exists as a layer, that it is localized, isotropic, and smoothly varying, and that the localized turbulence is described by the Kolmogorov spectrum. Closed-form solutions for both variances and frequency spectra of the log-amplitude and phase fluctuations are obtained using Rytov's method, and it is seen that the shape of the frequency spectra depends greatly on the characteristics and extent of the turbulence.

Woo, R.↗

Long-term differential energy spectrum for solar-flare iron-group particles

A long-term solar-flare differential energy spectrum for iron-group nuclei from approximately 0.1 to approximately 600 MeV/amu is derived from track density profile measurements in sample 64455 and sample 68815. Measurements from uneroded surfaces were obtained from quench crystals of plagioclase in 64455, and a Kr-81/Kr method indicates that the exposure age of this sample is 2,010,000 yrs. The power laws which best fit the normalized track density data are reported; the energy spectrum consists of two power law curves smoothly joined together which in turn are smoothly connected to a modulated galactic cosmic-ray spectrum. Standard track production versus depth profiles can be used to determine solar-flare track exposure ages and erosion rates for lunar samples.

Blanford, G. E.↗

Two-stream Maxwellian kinetic theory of cloud droplet growth by condensation

A new growth rate formula (NGRF) is developed for the rate of growth of cloud droplets by condensation. The theory used is a modification of the Lees-Shankar theory in which the two-stream Maxwellian distribution function of Lees is used in Maxwell's method of moments to determine the transport of water vapor to and heat away from the droplet. Boundary conditions at the droplet are the usual conditions set in terms of accommodation coefficients, and the solution passes smoothly into diffusion flow in the far region. Comparisons are given between NGRF and the conventional formula showing close agreement (approximately 0.1%) for large radii with significant difference (approximately 5%) for small radii (not greater than 1 micron). Growth times for haze droplets in a Laktionov chamber are computed.

Robinson, N. F.↗

Multigrid for hypersonic inviscid flows

The use of multigrid methods to solve the Euler equations for hypersonic flow is discussed. The steady state equations are considered with a Runge-Kutta smoother based on the time accurate equations together with local time stepping and residual smoothing. The effect of the Runge-Kutta coefficients on the convergence rate was examined considering both damping characteristics and convection properties. The importance of boundary conditions on the convergence rate for hypersonic flow is discussed. Also of importance are the switch between the second and fourth difference viscosity. Solutions are given for flow around the bump in a channel and flow around a biconic section.

Decker, Naomi H.↗

Design of transonic compressor cascades using hodograph method

The use of the Hodograph Method in the design of a transonic compressor cascade is discussed. The design of the flow mode in the transonic compressor cascade must be as follows: the flow in the nozzle part should be uniform and smooth; the location of the sonic line should be reasonable; and the aerodynamic character of the flow canal in the subsonic region should be met. The rate through cascade may be determined by the velocity distribution in the subsonic region (i.e., by the numerical solution of the Chaplygin equation). The supersonic sections A'C' and AD are determined by the analytical solution of the Mixed-Type Hodograph equation.

Chen, Zuoyi↗

Grid generation research at OSU

In the last two years, effort was concentrated on: (1) surface modeling; (2) surface grid generation; and (3) 3-D flow space grid generation. The surface modeling shares the same objectives as the surface modeling in computer aided design (CAD), so software available in CAD can in principle be used for solid modeling. Unfortunately, however, the CAD software cannot be easily used in practice for grid generation purposes, because they are not designed to provide appropriate data base for grid generation. Therefore, we started developing a generalized surface modeling software from scratch, that provides the data base for the surface grid generation. Generating surface grid is an important step in generating a 3-D space for flow space. To generate a surface grid on a given surface representation, we developed a unique algorithm that works on any non-smooth surfaces. Once the surface grid is generated, a 3-D space can be generated. For this purpose, we also developed a new algorithm, which is a hybrid of the hyperbolic and the elliptic grid generation methods. With this hybrid method, orthogonality of the grid near the solid boundary can be easily achieved without introducing empirical fudge factors. Work to develop 2-D and 3-D grids for turbomachinery blade geometries was performed, and as an extension of this research we are planning to develop an adaptive grid procedure with an interactive grid environment.

Nakamura, S.↗

New optimality criteria methods - Forcing uniqueness of the adjoint strains by corner-rounding at constraint intersections

In new, iterative continuum-based optimality criteria (COC) methods, the strain in the adjoint structure becomes non-unique if the number of active local constraints is greater than the number of design variables for an element. This brief note discusses the use of smooth envelope functions (SEFs) in overcoming economically computational problems caused by the above non-uniqueness.

Rozvany, G. I. N.↗

Task Space Angular Velocity Blending for Real-Time Trajectory Generation

The invention is embodied in a method of controlling a robot manipulator moving toward a target frame F(sub 0) with a target velocity v(sub 0) including a linear target velocity v and an angular target velocity omega(sub 0) to smoothly and continuously divert the robot manipulator to a subsequent frame F(sub 1) by determining a global transition velocity v(sub 1), the global transition velocity including a linear transition velocity v(sub 1) and an angular transition velocity omega(sub 1), defining a blend time interval 2(tau)(sub 0) within which the global velocity of the robot manipulator is to be changed from a global target velocity v(sub 0) to the global transition velocity v(sub 1) and dividing the blend time interval 2(tau)(sub 0) into discrete time segments (delta)t. During each one of the discrete time segments delta t of the blend interval 2(tau)(sub 0), a blended global velocity v of the manipulator is computed as a blend of the global target velocity v(sub 0) and the global transition velocity v(sub 1), the blended global velocity v including a blended angular velocity omega and a blended linear velocity v, and then, the manipulator is rotated by an incremental rotation corresponding to an integration of the blended angular velocity omega over one discrete time segment (delta)t.

Volpe, Richard A.↗

Causal Correlation Functions and Fourier Transforms: Application in Calculating Pressure Induced Shifts

By adopting a concept from signal processing, instead of starting from the correlation functions which are even, one considers the causal correlation functions whose Fourier transforms become complex. Their real and imaginary parts multiplied by 2 are the Fourier transforms of the original correlations and the subsequent Hilbert transforms, respectively. Thus, by taking this step one can complete the two previously needed transforms. However, to obviate performing the Cauchy principal integrations required in the Hilbert transforms is the greatest advantage. Meanwhile, because the causal correlations are well-bounded within the time domain and band limited in the frequency domain, one can replace their Fourier transforms by the discrete Fourier transforms and the latter can be carried out with the FFT algorithm. This replacement is justified by sampling theory because the Fourier transforms can be derived from the discrete Fourier transforms with the Nyquis rate without any distortions. We apply this method in calculating pressure induced shifts of H2O lines and obtain more reliable values. By comparing the calculated shifts with those in HITRAN 2008 and by screening both of them with the pair identity and the smooth variation rules, one can conclude many of shift values in HITRAN are not correct.

correlation function↗

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗

A design method for entrance sections of transonic wind tunnels with rectangular cross sections

A mathematical technique developed to design entrance sections for transonic or high-speed subsonic wind tunnels with rectangular cross sections is discribed. The transition from a circular cross-section setting chamber to a rectangular test section is accomplished smoothly so as not to introduce secondary flows (vortices or boundary-layer separation) into a uniform test stream. The results of static-pressure measurements in the transition region and of static and total-pressure surveys in the test section of a pilot model for a new facility at the Ames Research Center are presented.

Lionel, L.↗

Optimization methods applied to hybrid vehicle design

The use of optimization methods as an effective design tool in the design of hybrid vehicle propulsion systems is demonstrated. Optimization techniques were used to select values for three design parameters (battery weight, heat engine power rating and power split between the two on-board energy sources) such that various measures of vehicle performance (acquisition cost, life cycle cost and petroleum consumption) were optimized. The apporach produced designs which were often significant improvements over hybrid designs already reported on in the literature. The principal conclusions are as follows. First, it was found that the strategy used to split the required power between the two on-board energy sources can have a significant effect on life cycle cost and petroleum consumption. Second, the optimization program should be constructed so that performance measures and design variables can be easily changed. Third, the vehicle simulation program has a significant effect on the computer run time of the overall optimization program; run time can be significantly reduced by proper design of the types of trips the vehicle takes in a one year period. Fourth, care must be taken in designing the cost and constraint expressions which are used in the optimization so that they are relatively smooth functions of the design variables. Fifth, proper handling of constraints on battery weight and heat engine rating, variables which must be large enough to meet power demands, is particularly important for the success of an optimization study. Finally, the principal conclusion is that optimization methods provide a practical tool for carrying out the design of a hybrid vehicle propulsion system.

Donoghue, J. F.↗