Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational 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 757 records · Page 42

Methods of determining complete sensor requirements for autonomous mobility

A method of determining complete sensor requirements for autonomous mobility of an autonomous system includes computing a time variation of each behavior of a set of behaviors of the autonomous system, determining mobility sensitivity to each behavior of the autonomous system, and computing a change in mobility based upon the mobility sensitivity to each behavior and the time variation of each behavior. The method further includes determining the complete sensor requirements of the autonomous system through analysis of the relative magnitude of the change in mobility, the mobility sensitivity to each behavior, and the time variation of each behavior, wherein the relative magnitude of the change in mobility, the mobility sensitivity to each behavior, and the time variation of each behavior are characteristic of the stability of the autonomous system.

Curtis, Steven A.↗

Three-dimensional self-adaptive grid method for complex flows

A self-adaptive grid procedure for efficient computation of three-dimensional complex flow fields is described. The method is based on variational principles to minimize the energy of a spring system analogy which redistributes the grid points. Grid control parameters are determined by specifying maximum and minimum grid spacing. Multidirectional adaptation is achieved by splitting the procedure into a sequence of successive applications of a unidirectional adaptation. One-sided, two-directional constraints for orthogonality and smoothness are used to enhance the efficiency of the method. Feasibility of the scheme is demonstrated by application to a multinozzle, afterbody, plume flow field. Application of the algorithm for initial grid generation is illustrated by constructing a three-dimensional grid about a bump-like geometry.

Djomehri, M. Jahed↗

Three-dimensional self-adaptive grid method for complex flows

A self-adaptive grid procedure for efficient computation of three-dimensional complex flow fields is described. The method is based on variational principles to minimize the energy of a spring system analogy which redistributes the grid points. Grid control parameters are determined by specifying maximum and minimum grid spacing. Multidirectional adaptation is achieved by splitting the procedure into a sequence of successive applications of a unidirectional adaptation. One-sided, two-directional constraints for orthogonality and smoothness are used to enhance the efficiency of the method. Feasibility of the scheme is demonstrated by application to a multinozzle, afterbody, plume flow field. Application of the algorithm for initial grid generation is illustrated by constructing a three-dimensional grid about a bump-like geometry.

Djomehri, M. Jahed↗

A study of the threshold method utilizing raingage data

The threshold method for estimation of area-average rain rate relies on determination of the fractional area where rain rate exceeds a preset level of intensity. Previous studies have shown that the optimal threshold level depends on the climatological rain-rate distribution (RRD). It has also been noted, however, that the climatological RRD may be composed of an aggregate of distributions, one for each of several distinctly different synoptic conditions, each having its own optimal threshold. In this study, the impact of RRD variations on the threshold method is shown in an analysis of 1-min rainrate data from a network of tipping-bucket gauges in Darwin, Australia. Data are analyzed for two distinct regimes: the premonsoon environment, having isolated intense thunderstorms, and the active monsoon rains, having organized convective cell clusters that generate large areas of stratiform rain. It is found that a threshold of 10 mm/h results in the same threshold coefficient for both regimes, suggesting an alternative definition of optimal threshold as that which is least sensitive to distribution variations. The observed behavior of the threshold coefficient is well simulated by assumption of lognormal distributions with different scale parameters and same shape parameters.

Short, David A.↗

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗

More Genetic Engineering With Cloned Hemoglobin Genes

Cells modified to enhance growth and production of proteins. Method for enhancing both growth of micro-organisms in vitro and production of various proteins or metalbolites in these micro-organisms provides for incorporation of selected chromosomal or extrachormosomal deoxyribonucleic acid (DNA) sequences into micro-organisms from other cells or from artificial sources. Incorporated DNA includes parts encoding desired product(s) or characteristic(s) of cells and parts that control expression of productor characteristic-encoding parts in response to variations in environment. Extended method enables increased research into growth of organisms in oxygen-poor environments. Industrial applications found in enhancement of processing steps requiring oxygen in fermentation, enzymatic degradation, treatment of wastes containing toxic chemicals, brewing, and some oxidative chemical reactions.

Bailey, James E.↗

Application and development of the Schwinger multichannel scattering theory and the partial differential equation theory of electron-molecule scattering

One version of the multichannel theory for electron-target scattering based on the Schwinger variational principle, the SMC method, requires the introduction of a projection parameter. The role of the projection parameter a is investigated and it is shown that the principal-value operator in the SMC equation is Hermitian regardless of the value of a as long as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix representation of this operator is also Hermitian. The use of such basis is consistent with the Schwinger variational principle because the Lippmann-Schwinger equation automatically builds in the correct boundary conditions. Otherwise, an auxiliary condition needs to be introduced, and Takatsuka and McKoy's original value of a is one of the three possible ways to achieve Hermiticity. In all cases but one, a can be uncoupled from the Hermiticity condition and becomes a free parameter. An equation for a based on the variational stability of the scattering amplitude is derived; its solution has an interesting property that the scattering amplitude from a converged SMC calculation is independent of the choice of a even though the SMC operator itself is a-dependent. This property provides a sensitive test of the convergence of the calculation. For a static-exchange calculation, the convergence requirement only depends on the completeness of the one-electron basis, but for a general multichannel case, the a-invariance in the scattering amplitude requires both the one-electron basis and the N plus 1-electron basis to be complete. The role of a in the SMC equation and the convergence property are illustrated using two examples: e-CO elastic scattering in the static-exchange approximation, and a two-state treatment of the e-H2 Chi(sup 1)Sigma(sub g)(+) yields b(sup 3)Sigma(sub u)(+) excitation.

Weatherford, Charles A.↗

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon↗

Tip vortex geometry of a hovering helicopter rotor in ground effect

The wide-field shadowgraph method has been used to photograph the tip vortices of a hovering helicopter rotor in ground effect. The shadowgraphs were used to obtain quantitative measurements of the rotor tip vortex geometry both in and out of ground effect. Many important phenomena are visible in the rotor wake using this method. These include the variation in descent and contraction rates of the tip vortices in ground effect, and the interaction between tip vortices in the far wake. The tip vortex geometry from the shadowgraphs is compared with the tip vortex geometry predicted using a free wake hover performance analysis. The free wake analysis accurately predicts the tip vortex geometry both in and out of ground effect. Performance data from the test is compared with the performance predicted using several methods, including the free wake analysis. All methods provided reasonable predictions for the helicopter performance in ground effect.

Light, Jeffrey S.↗

Estimating Transmitted-Signal Phase Variations for Uplink Array Antennas

A method of estimating phase drifts of microwave signals distributed to, and transmitted by, antennas in an array involves the use of the signals themselves as phase references. The method was conceived as part of the solution of the problem of maintaining precise phase calibration required for proper operation of an array of Deep Space Network (DSN) antennas on Earth used for communicating with distant spacecraft at frequencies between 7 and 8 GHz. The method could also be applied to purely terrestrial phased-array radar and other radio antenna array systems. In the DSN application, the electrical lengths (effective signal-propagation path lengths) of the various branches of the system for distributing the transmitted signals to the antennas are not precisely known, and they vary with time. The variations are attributable mostly to thermal expansion and contraction of fiber-optic and electrical signal cables and to a variety of causes associated with aging of signal-handling components. The variations are large enough to introduce large phase drifts at the signal frequency. It is necessary to measure and correct for these phase drifts in order to maintain phase calibration of the antennas. A prior method of measuring phase drifts involves the use of reference-frequency signals separate from the transmitted signals. A major impediment to accurate measurement of phase drifts over time by the prior method is the fact that although DSN reference-frequency sources separate from the transmitting signal sources are stable and accurate enough for most DSN purposes, they are not stable enough for use in maintaining phase calibrations, as required, to within a few degrees over times as long as days or possibly even weeks. By eliminating reliance on the reference-frequency subsystem, the present method overcomes this impediment. In a DSN array to which the present method applies (see figure), the microwave signals to be transmitted are generated by exciters in a signal-processing center, then distributed to the antennas via optical fibers. At each antenna, the signals are used to drive a microwave power-amplifier train, the output of which is coupled to the antenna for transmission. A small fraction of the power-amplifier-train output is sent back to the signal-processing center along another optical fiber that is part of the same fiber-optic cable used to distribute the transmitted signal to the antenna. In the signal-processing center, the signal thus returned from each antenna is detected and its phase is compared with the phase of the signal sampled directly from the corresponding exciter. It is known, from other measurements, that the signal-propagation path length from the power-amplifier-train output port to the phase center of each antenna is sufficiently stable and, hence, that sampling the signal at the power-amplifier-train output port suffices for the purpose of characterizing the phase drift of the transmitted signal at the phase center of the antenna

Paal, Leslie↗

A spectral element method for the simulation of unsteady incompressible flows with heat transfer

The spectral element method is a high-order finite element technique for solution of the Navier-Stokes and energy equations. In the isoparametric spectral element discretization, the domain is broken up into general brick elements, and the dependent and independent variables represented as high-order tensor-product Lagrangian interpolants through Chebyshev collocation points. The nonlinear and convective terms in the governing equations are treated with explicit collocation, while the pressure and diffusive contributions are handled implicitly using variational projection operators. The method is applied to flow past a cylinder, flow in grooved channels, and natural convection in an enclosure.

Karniadakis, George E.↗

Discussion summary: Fictitious domain methods

Fictitious Domain methods are constructed in the following manner: Suppose a partial differential equation is to be solved on an open bounded set, Omega, in 2-D or 3-D. Let R be a rectangle domain containing the closure of Omega. The partial differential equation is first solved on R. Using the solution on R, the solution of the equation on Omega is then recovered by some procedure. The advantage of the fictitious domain method is that in many cases the solution of a partial differential equation on a rectangular region is easier to compute than on a nonrectangular region. Fictitious domain methods for solving elliptic PDEs on general regions are also very efficient when used on a parallel computer. The reason is that one can use the many domain decomposition methods that are available for solving the PDE on the fictitious rectangular region. The discussion on fictitious domain methods began with a talk by R. Glowinski in which he gave some examples of a variational approach to ficititious domain methods for solving the Helmholtz and Navier-Stokes equations.

Glowinski, Rowland↗

Combined forced and free convection in a curved duct

The purpose of this study is to investigate the flow and heat transfer characteristics of a combined forced and free convection flow in a curved duct. Solutions are obtained by solving the low Mach number model of the Navier-Stokes equation using a control volume method. The finite-volume method was developed with the use of a predictor-corrector numerical scheme and some new variations of the classical projection method. Solutions indicated that the existence of a buoyancy force has changed the entire flow structure inside a curved duct. Reversed flow at both inner and outer bend is observed. For moderate Reynolds number, the upstream section of the duct was significantly influenced by the free convection processes. In general, heat transfer is strong at the inner bend of the beginning of the heated section and at the outer bend on the last half of the heated section. The maximum velocity location is strongly influenced by the combined effects of buoyancy and centrifugal forces. A strong buoyancy force can reduce the strength of the secondary flow where it plays an important role in mixing.

Yam, Clement G.↗

Porosity and Variations in Microgravity Aerogel Nano-Structures: New Laser Speckle Characterization Methods - 2

The transition from sol to gel is a process that is critical to the properties of engineered nanomaterials, but one with few available techniques for observing the dynamic processes occurring during the evolution of the gel network. Specifically, the observation of various cluster aggregation models, such as diffusion-limited and reaction-limited cluster growth can be quite difficult. This can be rather important as the actual aggregation model can dramatically influence the mechanical properties of gels, and is significantly affected by the presence of convective flows, or their absence in microgravity. We have developed two new non-intrusive optical methods for observing the aggregation processes within gels in real time. These make use of the dynamic behavior of laser speckle patterns produced when an intense laser source is passed through a gelling sol. The first method is a simplified time-correlation measurement, where the speckle pattern is observed using a CCD camera and information on the movement of the scattering objects is readily apparent. This approach is extremely sensitive to minute variations in the flow field as the observed speckle pattern is a diffraction-based image, and is therefore sensitive to motions within the sol on the order of the wavelength of the probing light. Additionally, this method has proven useful in determining a precise time for the gel-point, an event often difficult to measure. Monitoring the evolution of contrast within the speckle field is another method that has proven useful for studying aeration. In this case, speckle contrast is dependent upon the size (correlation length) and number of scattering centers, increasing with increasing size, and decreasing with increasing numbers. The dynamic behavior of cluster growth in gels causes both of these to change simultaneously with time, the exact rate of which is determined by the specific aggregation model involved. Actual growth processes can now be observed, and the effects of varying gravity fields on the growth processes qualitatively described. Results on preliminary ground-based measurements have been obtained.

Hunt, A. J.↗

Vertical profiles of wind and temperature by remote acoustical sounding

An acoustical method was investigated for obtaining meteorological soundings based on the refraction due to the vertical variation of wind and temperature. The method has the potential of yielding horizontally averaged measurements of the vertical variation of wind and temperature up to heights of a few kilometers; the averaging takes place over a radius of 10 to 15 km. An outline of the basic concepts and some of the results obtained with the method are presented.

Fox, H. L.↗

Mechanically fastened composite laminates subjected to combined bearing-bypass and shear loading

Bolts and rivets provide a means of load transfer in the construction of aircraft. However, they give rise to stress concentrations and are often the source and location of static and fatigue failures. Furthermore, fastener holes are prone to cracks during take-off and landing. These cracks present the most common origin of structural failures in aircraft. Therefore, accurate determination of the contact stresses associated with such loaded holes in mechanically fastened joints is essential to reliable strength evaluation and failure prediction. As the laminate is subjected to loading, the contact region, whose extent is not known, develops between the fastener and the hole boundary through this contact region, which consists of slip and no-slip zones due to friction. The presence of the unknown contact stress distribution over the contact region between the pin and the composite laminate, material anisotropy, friction between the pin and the laminate, pin-hole clearance, combined bearing-bypass and shear loading, and finite geometry of the laminate result in a complex non-linear problem. In the case of bearing-bypass loading in compression, this non-linear problem is further complicated by the presence of dual contact regions. Previous research concerning the analysis of mechanical joints subjected to combined bearing-bypass and shear loading is non-existent. In the case of bearing-bypass loading only, except for the study conducted by Naik and Crews (1991), others employed the concept of superposition which is not valid for this non-linear problem. Naik and Crews applied a linear finite element analysis with conditions along the pin-hole contact region specified as displacement constraint equations. The major shortcoming of this method is that the variation of the contract region as a function of the applied load should be known a priori. Also, their analysis is limited to symmetric geometry and material systems, and frictionless boundary conditions. Since the contact stress distribution and the contact region are not known a priori, they did not directly impose the boundary conditions appropriate for modelling the contact and on-contact regions between the fastener and the hole. Furthermore, finite element analysis is not suitable for iterative design calculations for optimizing laminate construction in the presence of fasteners under complex loading conditions. In this study, the solution method developed by Madenci and Ileri (1992a,b) has been extended to determine the contact stresses in mechanical joints under combined bearing-bypass and shear loading, and bearing-bypass loading in compression resulting in dual contact regions.

Madenci, Erdogan↗

Analytic marginalization of N(z) uncertainties in tomographic galaxy surveys

In this paper, we present a new method to marginalize over uncertainties in redshift distributions, N(z), within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess N(z) governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions. In principle, this is marginalization over hundreds or thousands of new parameters describing potential deviations as a function of redshift and tomographic bin. However, by linearly expanding the theory predictions around a fiducial model, this marginalization can be performed analytically, resulting in a modified data covariance matrix that effectively downweights the modes of the data vector that are more sensitive to redshift distribution variations. We showcase this method by applying it to the galaxy clustering measurements from the Hyper Suprime-Cam first data release. We illustrate how to marginalize over sample variance of the calibration sample and a large general systematic uncertainty in photometric estimation methods, and explore the impact of priors imposing smoothness in the redshift distributions.

79 ASTRONOMY AND ASTROPHYSICS↗