Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution regularity”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Solar-wind structure determined by corotating coronal inhomogeneities. II - Arbitrary perturbations.

The problem of the solar-wind structure resulting from long-lived inhomogeneities in the solar corona has been extended within the framework of a linearized hydrodynamic approach to allow for arbitrary perturbations in the plasma parameters, namely, velocity, density, and temperature. By using a Parker model for the zero-order flow speed, solutions are given for perturbations imposed at the zero-order critical point. Using the restriction that the solutions must be regular at the critical point, it is possible to express all situations as a linear sum of the solutions for pure velocity, density, and temperature perturbations at the critical point. The results of the integrations for three inner-boundary conditions can be used as the elements of matrices that map arbitrary perturbations at one heliocentric distance to any other distance. One interesting feature of the solutions is that the amplitudes at earth for the pure initial temperature situation are larger than for the other two situations.

Siscoe, G. L.↗
NASA NTRS2006

Mascons, GRACE, and Time-variable Gravity

The GRACE mission has been in orbit now for three years and now regularly produces snapshots of the Earth s gravity field on a monthly basis. The convenient standard approach has been to perform global solutions in spherical harmonics. Alternative local representations of mass variations using mascons show great promise and offer advantages in terms of computational efficiency, minimization of problems due to aliasing, and increased temporal resolution. In this paper, we discuss the results of processing the GRACE KBRR data from March 2003 through August 2005 to produce solutions for GRACE mass variations over mid-latitude and equatorial regions, such as South America, India and the United States, and over the polar regions (Antarctica and Greenland), with a focus on the methodology. We describe in particular mascon solutions developed on regular 4 degree x 4 degree grids, and those tailored specifically to drainage basins over these regions.

Lemoine, F.↗
NASA NTRS2005

Computation of the Streamfunction and Velocity Potential for Limited and Irregular Domains

An algorithm is proposed for the computation of streamfunction and velocity potential from given horizontal velocity vectors based on solving a minimization problem. To guarantee the uniqueness of the solution and computational reliability of the algorithm, a Tikhonov regularization is applied. The solution implies that the obtained streamfunction and velocity potential have minimal magnitude, while the given velocity vectors can be accurately reconstructed from the computed streamfunction and velocity potential. Because the formulation of the minimization problem allows for circumventing the explicit specification of separate boundary conditions on the streamfunction and velocity potential, the algorithm is easily applicable to irregular domains. By using an advanced minimization algorithm with the use of adjoint techniques, the method is computationally efficient and suitable for problems with large dimensions. An example is presented for coastal oceans to illustrate the practical application of the algorithm.

velocity↗
NASA NTRS1992

Inverse problems: Fuzzy representation of uncertainty generates a regularization

In many applied problems (geophysics, medicine, and astronomy) we cannot directly measure the values x(t) of the desired physical quantity x in different moments of time, so we measure some related quantity y(t), and then we try to reconstruct the desired values x(t). This problem is often ill-posed in the sense that two essentially different functions x(t) are consistent with the same measurement results. So, in order to get a reasonable reconstruction, we must have some additional prior information about the desired function x(t). Methods that use this information to choose x(t) from the set of all possible solutions are called regularization methods. In some cases, we know the statistical characteristics both of x(t) and of the measurement errors, so we can apply statistical filtering methods (well-developed since the invention of a Wiener filter). In some situations, we know the properties of the desired process, e.g., we know that the derivative of x(t) is limited by some number delta, etc. In this case, we can apply standard regularization techniques (e.g., Tikhonov's regularization). In many cases, however, we have only uncertain knowledge about the values of x(t), about the rate with which the values of x(t) can change, and about the measurement errors. In these cases, usually one of the existing regularization methods is applied. There exist several heuristics that choose such a method. The problem with these heuristics is that they often lead to choosing different methods, and these methods lead to different functions x(t). Therefore, the results x(t) of applying these heuristic methods are often unreliable. We show that if we use fuzzy logic to describe this uncertainty, then we automatically arrive at a unique regularization method, whose parameters are uniquely determined by the experts knowledge. Although we start with the fuzzy description, but the resulting regularization turns out to be quite crisp.

Kreinovich, V.↗
NASA NTRS2003

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗
NASA NTRS1974

A finite element solution algorithm for the Navier-Stokes equations

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the steady-state kinematics and thermodynamics of a variable viscosity, compressible multiple-species fluid. For an incompressible fluid, the motion may be transient as well. The primitive dependent variables are replaced by a vorticity-streamfunction description valid in domains spanned by rectangular, cylindrical and spherical coordinate systems. Use of derived variables provides a uniformly elliptic partial differential equation description for the Navier-Stokes system, and for which the finite element algorithm is established. Explicit non-linearity is accepted by the theory, since no psuedo-variational principles are employed, and there is no requirement for either computational mesh or solution domain closure regularity. Boundary condition constraints on the normal flux and tangential distribution of all computational variables, as well as velocity, are routinely piecewise enforceable on domain closure segments arbitrarily oriented with respect to a global reference frame.

Baker, A. J.↗
NASA NTRS1994

Renormalization Group Theory of Bolgiano Scaling in Boussinesq Turbulence

Bolgiano scaling in Boussinesq turbulence is analyzed using the Yakhot-Orszag renormalization group. For this purpose, an isotropic model is introduced. Scaling exponents are calculated by forcing the temperature equation so that the temperature variance flux is constant in the inertial range. Universal amplitudes associated with the scaling laws are computed by expanding about a logarithmic theory. Connections between this formalism and the direct interaction approximation are discussed. It is suggested that the Yakhot-Orszag theory yields a lowest order approximate solution of a regularized direct interaction approximation which can be corrected by a simple iterative procedure.

Rubinstein, Robert↗
NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗
NASA NTRS2001

Calculation of Phase Equilibria in the Y2O3-Yb2O3-ZrO2 System

Rare earth oxide stabilized zirconias find a wide range of applications. An understanding of phase equilibria is essential to all applications. In this study, the available phase boundary data and thermodynamic data is collected and assessed. Calphad-type databases are developed to completely describe the Y2O3-ZrO2, Yb2O3-ZrO2, and Y2O3-Yb2O3 systems. The oxide units are treated as components and regular and subregular solution models are used. The resultant calculated phase diagrams show good agreement with the experimental data. Then the binaries are combined to form the database for the Y2O3-Yb2O3-ZrO2 psuedo-ternary.

Jacobson, Nathan S.↗
NASA NTRS2005

Dependence of Ru2O3 Activity on Composition of Silicate Melts: Using Statistical Correlations to Infer Thermodynamic Behavior in the Melt

Understanding variations in activity with composition is an essential step in improving prediction of partition coefficients during magma evolution. Variations in activity with composition are complex and do not generally exhibit ideal behavior relative to a traditional melt-component set. Although deviations from component ideality can be modeled numerically by simply fitting to compositional variables (such as in a regular or subregular solution model), such models have not been particularly successful for describing variations in trace component activities. A better approach might be to try to identify components that do a better job of describing the behavior of the species in the melt. Electrochemical Measurement of Ru2O3 activities: Electrodes were inserted into silicate melt beads of various compositions (Table 1) suspended on Ptwire loops in a 1-atm gas mixing furnace. An electrical potential was imposed between the electrodes, the imposed potential increasing along a step ramp with a pulse imposed on each step (Fig. 1). Current flows between electrodes when electroactive species in the melt are oxidized or reduced at the electrodes. The resulting current was measured at the top and bottom of the voltage pulse, and the difference (the differential current) was plotted against potential. The peak of the resulting curve is related to the activity coefficient for the particular electroactive species (Ru2O3) in the melt [1, 2, 3]. A significant part of the nonideal contribution to activity is due not to intrinsic properties of the component in the melt, but to our ignorance about the state and mixing properties of the component in the melt.

Colson, R. O.↗
NASA NTRS2021

GRACE Fast Mascons from Spherical Harmonics and a Regularization Design Trade Study

Mass concentration (mascon) solutions have become a prominent medium for investigating time-variable gravity recovered by GRACE and GRACE Follow-On. While GRACE Level-2 spherical harmonic products require various post-processing techniques to eliminate correlated noise, mascon formulations employ spatial regularization strategies during the estimation step to improve signal recovery. However, mascon estimation has traditionally required large computing resources and GRACE Level-1B processing capabilities. In this study, we show that a typical mascon estimation system can be reformulated to allow for the estimation of regularized mascons from Level-2 spherical harmonics. Provided that spherical harmonic solution covariances are available, the computed mascons will be mathematically equivalent to similar mascons estimated from Level-1B observations. This method is computationally efficient, better leverages GRACE spherical harmonics than past methods, and matches the performance of typical mascon solutions without locking scientists into pre-determined regularization designs. We develop a proof-of-concept solution using ITSG-Grace2018 and compare results with traditional mascons from JPL and NASA GSFC. We then assess the effects of spherical harmonic truncation and use of regularization correlations on basin signal recovery. We find that spherical expansions to degree and order 60 provide the minimum expansion necessary to study most basins, while larger expansions help further localize signals. We also find that diagonal regularizations (i.e., regularizations that do not contain inter-mascon correlations) are adversely affected by leakage, especially across boundaries such as coastlines where signals are not highly correlated, whereas including inter-mascon correlations and regional boundaries in the regularization greatly improves signal recovery.

M. J. Croteau↗
NASA NTRS2020

Development of a Daily GRACE Mascon Solution for Terrestrial Water Storage

The Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow On missions have provided a global history of terrestrial water storage changes since 2002. Traditional GRACE products resolve monthly time-variable gravity at spatial resolutions of 300-500 km, with many recent efforts focusing on regularized mass concentration (mascon) solutions to better resolve signals spatially. However, monthly resolution inhibits the applicability of GRACE to investigations of sub-monthly signals. This study presents a new daily mass change solution estimated as deviations from the NASA Goddard Space Flight Center (GSFC) converged monthly mascon product and quantifies the fundamental trade-off between temporal and spatial resolution in GRACE-only solutions. As an iteration of the GSFC product, this daily solution represents a mixing of monthly information at higher relative spatial resolution with daily information at lower relative spatial resolution, resolving each temporal scale at the best resolution achievable without requiring hydrological model dependencies in the estimation. The resolution of the daily mascons is 300-400 km at high latitudes and 600-1,000 km lower, depending on each daily orbital track and the proximity of each mascon to constraint region boundaries (e.g., coastlines). Through simulations and model comparisons, we demonstrate that daily signals over areas larger than 400,000 km(2) are recoverable, with basins larger than 800,000 km(2) exhibiting strong signal recovery relative to leakage. This analysis establishes baselines for daily signal recovery from GRACE in the context of longer time scale solutions, characterizes leakage inherent in daily GRACE information, and creates new opportunities for applying GRACE to investigations of sub-monthly signals.

M. J. Croteau↗
NASA NTRS2013

Advanced Signal Processing for High Temperatures Health Monitoring of Condensed Water Height in Steam Pipes

An advanced signal processing methodology is being developed to monitor the height of condensed water thru the wall of a steel pipe while operating at temperatures as high as 250deg. Using existing techniques, previous study indicated that, when the water height is low or there is disturbance in the environment, the predicted water height may not be accurate. In recent years, the use of the autocorrelation and envelope techniques in the signal processing has been demonstrated to be a very useful tool for practical applications. In this paper, various signal processing techniques including the auto correlation, Hilbert transform, and the Shannon Energy Envelope methods were studied and implemented to determine the water height in the steam pipe. The results have shown that the developed method provides a good capability for monitoring the height in the regular conditions. An alternative solution for shallow water or no water conditions based on a developed hybrid method based on Hilbert transform (HT) with a high pass filter and using the optimized windowing technique is suggested. Further development of the reported methods would provide a powerful tool for the identification of the disturbances of water height inside the pipe.

health monitoring↗
NASA NTRS1986

Regularity estimates up to the boundary for elliptic systems of difference equations

Regularity estimates up to the boundary for solutions of elliptic systems of finite difference equations were proved. The regularity estimates, obtained for boundary fitted coordinate systems on domains with smooth boundary, involve discrete Sobolev norms and are proved using pseudo-difference operators to treat systems with variable coefficients. The elliptic systems of difference equations and the boundary conditions which are considered are very general in form. The regularity of a regular elliptic system of difference equations was proved equivalent to the nonexistence of eigensolutions. The regularity estimates obtained are analogous to those in the theory of elliptic systems of partial differential equations, and to the results of Gustafsson, Kreiss, and Sundstrom (1972) and others for hyperbolic difference equations.

Strikwerda, J. C.↗
NASA NTRS1988

Supersonic combusting flow - A comparison of numerical and experimental results

The flow field of an unconfined coaxial supersonic burner (SSB) was obtained numerically, using a CFD code. The explicit MacCormack algorithm, with point implicit treatment of chemistry source terms, was employed in the CFD code. Combustion of hydrogen and air was simulated by a two-step finite-rate combustion model whereas turbulence was accounted for by a Prandtl mixing length scheme. Boundary conditions and some important features of the numerical scheme were discussed. The qualitative features of the numerical solution were compared with those obtained by regular as well as schlieren photographs. General agreements of the CFD solution and the experimental photographs were observed.

Chitsomboon, T.↗
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2019

Regularization and Error Characterization of GRACE Mascons

We present a new global time-variable gravity mascon solution derived from Gravity Recovery and Climate Experiment (GRACE) Level 1B data. The new product from the NASA Goddard Space Flight Center (GSFC) results from a novel approach that combines an iterative solution strategy with geographical binning of inter-satellite range-acceleration residuals in the construction of time-dependent regularization matrices applied in the inversion of mascon parameters. This estimation strategy is intentionally conservative as it seeks to maximize the role of the GRACE measurements on the final solution while minimizing the influence of the regularization design process. We fully reprocess the Level 1B data in the presence of the final mascon solution to generate true post-fit inter-satellite residuals, which are utilized to confirm solution convergence and to validate the mascon noise uncertainties. We also present the mathematical case that regularized mascon solutions are biased, and that this bias, or leakage, must be combined with the estimated noise variance to accurately assess total mascon uncertainties. The estimated leakage errors are determined from the monthly resolution operators. We present a simple approach to compute the total uncertainty for both individual mascon and regional analysis of the GSFC mascon product, and validate the results in comparison with independent mascon solutions and calibrated Stokes uncertainties. Lastly, we present the new solution and uncertainties with global analyses of the mass trends and annual amplitudes, and compute updated trends for the global ocean, and the respective contributions of the Greenland Ice Sheet, Antarctic Ice Sheet, Gulf of Alaska, and terrestrial water storage. This analysis highlights the successful closure of the global mean sea level budget, that is, the sum of global ocean mass from the GSFC mascons and the steric component from Argo floats agrees well with the total determined from sea surface altimetry.

B. D. Loomis↗