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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1995

A far-field non-reflecting boundary condition for two-dimensional wake flows

Far-field boundary conditions for external flow problems have been developed based upon long-wave perturbations of linearized flow equations about a steady state far field solution. The boundary improves convergence to steady state in single-grid temporal integration schemes using both regular-time-stepping and local-time-stepping. The far-field boundary may be near the trailing edge of the body which significantly reduces the number of grid points, and therefore the computational time, in the numerical calculation. In addition the solution produced is smoother in the far-field than when using extrapolation conditions. The boundary condition maintains the convergence rate to steady state in schemes utilizing multigrid acceleration.

Danowitz, Jeffrey S.↗

NASA NTRS2003

Prospect of Using Numerical Dynamo Model for Prediction of Geomagnetic Secular Variation

Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Kuang, Weijia↗

NASA NTRS1989

Evaluation of VSAERO in prediction of aerodynamic characteristics of helicopter hub fairings

A low-order panel code, VSAERO, was used to predict the aerodynamic characteristics of helicopter hub fairings. Since the simulation of this kind of bluff body by VSAERO was not documented before, the VSAERO solutions were correlated with experimental data to establish their validity. The validation process revealed that simulation of the aerodynamic environment around a hub fairing was sensitive to several modeling parameters. Some of these parameters are body and wake panels arrangement, streamwise and spanwise separation location, and the most prominent one-the wake modeling. Three wake models were used: regular wake, separated wake, and jet model. The regular wake is a wake with negligible thickness (thin wake). It is represented by a single vortex sheet. The separated wake and the jet model in the present application are wakes with finite thickness (thick wake). They consist of a vortex sheet enclosing a region of low-energy flow. The results obtained with the reqular wake were marginally acceptable for sharp-edged hub fairings. For all other cases under consideration, the jet model results correlated slightly better. The separated wake, which seemed to be the most appropriate model, caused the solution to diverge. While the regular wake was straight-forward to apply in simulations, the jet model was not. It requires the user to provide information about the doublet strength gradient on wake panels by guessing the efflux velocities at the wake shedding location. In summary, VSAERO neither predicts accurately the aerodynamic characteristics of helicopter hub fairings nor was cost effective.

Louie, Alexander↗

NASA NTRS2011

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms.

Knapmeyer, Martin↗

NASA NTRS1972

An optimum settling problem for time lag systems.

A solution is presented to an optimization problem for time lag systems by the classical method of Lagrange multipliers in a Banach space. Following terminology and assumption definitions, the regularity and controllability of the Lagrange multipliers problem is discussed, and a set of necessary conditions for an optimal control is derived. In conclusion, the solution existence, uniqueness, and sufficiency are established.

Jacobs, M. Q.↗

NASA NTRS2000

Observational Model for Precision Astrometry with the Space Interferometry Mission

The Space Interferometry Mission (SIM) is a space-based 10-m baseline Michelson optical interferometer operating in the visible waveband that is designed to achieve astrometric accuracy in the single digits of the microarcsecond domain. Over a narrow field of view SIM is expected to achieve a mission accuracy of 1 microarcsecond. In this mode SIM will search for planetary companions to nearby stars by detecting the astrometric "wobble" relative to a nearby reference star. In its wide-angle mode, SIM will provide 4 microarcsecond precision absolute position measurements of stars, with parallaxes to comparable accuracy, at the end of its 5-year mission. The expected proper motion accuracy is around 3 microarcsecond/year, corresponding to a transverse velocity of 10 m/ s at a distance of 1 kpc. The basic astrometric observable of the SIM instrument is the pathlength delay. This measurement is made by a combination of internal metrology measurements that determine the distance the starlight travels through the two arms of the interferometer, and a measurement of the white light stellar fringe to find the point of equal pathlength. Because this operation requires a non-negligible integration time, the interferometer baseline vector is not stationary over this time period, as its absolute length and orientation are time varying. This paper addresses how the time varying baseline can be "regularized" so that it may act as a single baseline vector for multiple stars, as required for the solution of the astrometric equations.

TILE OBSERVATIONS↗

NASA NTRS2005

Measuring, Enabling and Comparing Modularity, Regularity and Hierarchy in Evolutionary Design

For computer-automated design systems to scale to complex designs they must be able to produce designs that exhibit the characteristics of modularity, regularity and hierarchy - characteristics that are found both in man-made and natural designs. Here we claim that these characteristics are enabled by implementing the attributes of combination, control-flow and abstraction in the representation. To support this claim we use an evolutionary algorithm to evolve solutions to different sizes of a table design problem using five different representations, each with different combinations of modularity, regularity and hierarchy enabled and show that the best performance happens when all three of these attributes are enabled. We also define metrics for modularity, regularity and hierarchy in design encodings and demonstrate that high fitness values are achieved with high values of modularity, regularity and hierarchy and that there is a positive correlation between increases in fitness and increases in modularity. regularity and hierarchy.

Hornby, Gregory S.↗

NASA NTRS1994

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗

NASA NTRS1977

The variety of variables in mathematical expressions

Various distinctions are discussed which can be made regarding evaluation of mathematical expressions: regular evaluation vs. infinite evaluation vs. finite evaluation, regular variables vs. mathematical variables vs. shadow variables vs. labels, simplification vs. evaluation vs. solution of equations. The unsatisfactory state of evaluation strategies in symbolic systems is due to insufficient use of such distinctions in the past.

Moses, J.↗

NASA NTRS1990

Implementation of two projection methods on a shared memory multiprocessor - DEC VAX 6240

The relative performance of two iterative schemes, based on projection techniques, is compared on a shared memory multiprocessor - VAX 6240. The CG accelerated Block-SSOR method and the CG accelerated Symmetric-Kaczmarz method are considered for the solution of large sparse nonsymmetric systems of linear equations. It is shown that the regular structure of many matrices can be exploited by the CG-accelerated Block-SSOR method to provide good speedup in a multiprocessing environment. However, the CG accelerated Symmetric-Kaczmarz method, while being a viable alternative on a scalar machine, is unable to benefit from multiprocessing.

Kamath, C.↗

NASA NTRS1991

Energy functions for regularization algorithms

Regularization techniques are widely used for inverse problem solving in computer vision such as surface reconstruction, edge detection, or optical flow estimation. Energy functions used for regularization algorithms measure how smooth a curve or surface is, and to render acceptable solutions these energies must verify certain properties such as invariance with Euclidean transformations or invariance with parameterization. The notion of smoothness energy is extended here to the notion of a differential stabilizer, and it is shown that to void the systematic underestimation of undercurvature for planar curve fitting, it is necessary that circles be the curves of maximum smoothness. A set of stabilizers is proposed that meet this condition as well as invariance with rotation and parameterization.

Delingette, H.↗

NASA NTRS2004

Statistical analysis of geodetic networks for detecting regional events

We present an application of hidden Markov models (HMMs) to analysis of geodetic time series in Southern California. Our model fitting method uses a regularized version of the deterministic annealing expectation-maximization algorithm to ensure that model solutions are both robust and of high quality.

geodetic networks↗

NASA NTRS2013

The BKG/IGGB VLBI Analysis Center

In 2012, the activities of the BKG/IGGB VLBI Analysis Center, as in previous years, consisted of routine computations of Earth orientation parameter (EOP) time series and of a number of research topics in geodetic VLBI. The VLBI group at BKG continued its regular submissions of time series of tropospheric parameters and the generation of daily SINEX (Solution INdependent EXchange format) files. Quarterly updated solutions have been computed to produce terrestrial reference frame (TRF) and celestial reference frame (CRF) realizations. Routine computations of the UT1-UTC Intensive observations include all sessions of the Kokee-Wettzell and Tsukuba-Wettzell baselines and the networks Kokee-Svetloe-Wettzell and Ny-degAlesund-Tsukuba-Wettzell. The VLBI group at BKG developed a procedure to get the most probable station positions of Tsukuba after the earthquake on March 11, 2011 for the epochs of the Intensive sessions. The analysis of the Intensive sessions with station Tsukuba could be resumed in February 2012. At IGGB, the emphasis has been placed on individual research topics.

Thorandt, Volkmar↗

NASA NTRS2023

DPOD2020: A DORIS Extension of the ITRF2020 for Precise Orbit Determination

As one of the tracking systems used to determine orbits of the altimeter mission satellites (such as TOPEX/Poseidon, Envisat, Jason-1/-2/-3, CryoSat-2, Saral/Altika, Sentinel-3A/-3B, HY-2A/C/D, Jason-CS/Sentinel-6A, SWOT), DORIS (Doppler Orbitography Radiopositionning Integrated by Satellite) allows to determine positions and velocities of tracking stations that define a stable reference for the estimation of the precise orbits and thus are fundamental for the quality of the altimeter data and derived mean sea level products. Due to the time evolution of the DORIS ground network, some stations included in the 2020 realization of the International Terrestrial Reference Frame (ITRF2020) have been decommissioned and since 2021.0 a few new stations were added to the tracking network. Therefore, to satisfy operational requirements for POD (Precise Orbit Determination) and routine delivery of geodetic products, the International DORIS Service (IDS) regularly updates the DPOD (DORIS terrestrial reference frame for Precise Orbit Determination). The DPOD solutions include mean positions and velocities of all the DORIS stations since 1993.0 derived from the stacking of the latest IDS weekly combined series aligned to the current ITRF. In this paper, we first present the stacking process of the DPOD2020 version 1.0. Then, we address the validation procedure of the DPO2020 including comparison with ITRF2020 and POD tests. For eighty percent of all the time segments of all the DORIS stations, the station position differences between DPOD2020 version 1.0 and ITRF2020 are smaller than ten millimeters. The major position differences between these two solutions are associated with the DORIS sites either localized in the South Atlantic Anomaly region or with time spans smaller than one year. Compared to DPOD2014, the DPOD2020 shows reduction of the main statistics of the DORIS-to-DORIS tie residuals (differences between the estimated and measured ties). The POD tests showed similar results for DPOD2020 and DPOD2014 for most of the altimetric satellites (TOPEX/Poseidon, Jason-1, CryoSat-2, Jason-3). In addition, we observed better POD results with DPOD2020 for the latest altimetric satellite Sentinel-3A as well as a slight degradation for Jason-2. That degradation was fully explained by slightly worst results for the stations localized in the South Atlantic Anomaly region.

DORIS↗

NASA NTRSDate not supplied

Evaluating SAR Radiometric Terrain Correction Solutions: Optimal products for applied users

Operational applications of Synthetic Aperture Radar (SAR) are under development around the world, driven by the regularly-acquired, free-and-open source C-band SAR observations provided by ESA’s Sentinel-1 sensor constellation since 2014. Groups like SERVIR, a joint NASA and USAID initiative, are at the forefront of remote sensing applications for societal benefit. A takeaway from SERVIR’s experience is the need for appropriately geocoded and fully calibrated SAR data that is ready to use for a range of ecosystems-related applications. Radiometric Terrain Corrected (RTC) data are key entry-level products for multiple applications that range from ecosystems to hazards. This work fills a gap in current research by evaluating several RTCs produced by open-source software solutions (SNAP-7 and ISCE-2), the gold standard commercial software (GAMMA), a Google Earth Engine (GEE) based workflow, and the uncorrected GRD products currently available in GEE. RTCs were analyzed for geolocation quality, absolute radiometric calibration, and fidelity of the radiometric terrain flattening over ten sites representing varied terrains. In addition, a time series analysis was conducted over two locations. Overall, no significant differences for radiometric calibration were found across RTC products. However, all RTCs performed better than uncorrected GRD products. The main differences between products were found in geolocation quality. These results not only demonstrate the need for the uptake and distribution of RTC products for ecosystems applications, but demonstrate the ability to do so with open source methods, adding value to developing affordable operational applications.

Helen Blue Parache↗

NASA NTRS2007

Process for High-Rate Fabrication of Alumina Nanotemplates

An anodizing process, at an early stage of development at the time of reporting the information for this article, has shown promise as a means of fabricating alumina nanotemplates integrated with silicon wafers. Alumina nanotemplates are basically layers of alumina, typically several microns thick, in which are formed approximately regular hexagonal arrays of holes having typical diameters of the order of 10 to 100 nm. Interest in alumina nanotemplates has grown in recent years because they have been found to be useful as templates in the fabrication of nanoscale magnetic, electronic, optoelectronic, and other devices. The present anodizing process is attractive for the fabrication of alumina nanotemplates integrated with silicon wafers in two respects: (1) the process involves self-ordering of the holes; that is, the holes as formed by the process are spontaneously arranged in approximately regular hexagonal arrays; and (2) the rates of growth (that is, elongation) of the holes are high enough to make the process compatible with other processes used in the mass production of integrated circuits. In preparation for fabrication of alumina nanotemplates in this process, one first uses electron-beam evaporation to deposit thin films of titanium, followed by thin films of aluminum, on silicon wafers. Then the alumina nanotemplates are formed by anodizing the aluminum layers, as described below. In experiments in which the process was partially developed, the titanium films were 200 A thick and the aluminum films were 5 m thick. The aluminum films were oxidized to alumina, and the arrays of holes were formed by anodizing the aluminum in aqueous solutions of sulfuric and/or oxalic acid at room temperature (see figure). The diameters, spacings, and rates of growth of the holes were found to depend, variously, on the composition of the anodizing solution, the applied current, or the applied potential, as follows: In galvanostatically controlled anodizing, regardless of the chemical composition of the solution, relatively high current densities (50 to 100 mA/cm2) resulted in arrays of holes that were more nearly regular than were those formed at lower current densities. . The rates of elongation of the holes were found to depend linearly on the applied current density: the observed factor of proportionality was 1.2 (m/h)/(mA/cm2). For a given fixed current density and room temperature, the hole diameters were found to depend mainly on the chemical compositions of the anodizing solutions. The holes produced in sulfuric acid solutions were smaller than those produced in oxalic acid solutions. The arrays of holes produced in sulfuric acid were more ordered than were those produced in oxalic acid. . The breakdown voltage was found to decrease logarithmically with increasing concentration of sulfuric acid. The breakdown voltage was also found to decrease with temperature and to be accompanied by a decrease in hole diameter. The hole diameter was found to vary linearly with applied potential, with a slope of 2.1 nm/V. This slope differs from slopes (2.2 and 2.77 nm/V) reported for similar prior measurements on nanotemplates made from bulk aluminum. The differences among these slopes may be attributable to differences among impurities and defects in bulk and electron-beam-evaporated aluminum specimens.

Myung, Nosang↗

NASA NTRS1994

Echo mapping of active galactic nuclei broad-line regions: Fundamental algorithms

We formulate and test a series of algorithms for echo mapping the emission-line regions near active galactic nuclei from measurements of correlated variability in their line and continuum light curves. The linear regularization method (LRM) employs a direct inversion of evenly spaced light-curve data, with a regularization parameter that can be used to control the trade-off between noise and resolution. Matrix formulas express the formal solution as well as its variance and covariance in terms of uncertainties in the measurements. Unlike the maximum-entropy method (MEM), LRM applies to kernels with both positive and negative values, but the results are somewhat limited by ringing effects. A positivity constraint proves effective in controlling the ringing. MEM combines regularization and positivity in a natural way, but similar results are also found using positivity constraints with nonentropic regularization functions. Direct inversions of unevenly sampled light curves require interpolating the noisy data. In this case better results are found by solving for both the continuum light curve and kernel function in a simultaneous fit to the data. Our conclusion is that while echo mapping currently gives ambiguous results, the algorithms are not the limiting factor. Progress depends on efforts to increase the accuracy and completeness of sampling of the observed light curves.

Vio, Roberto↗

NASA NTRS1972

Mariner 9 celestial mechanics experiment - Gravity field and pole direction of Mars.

Analysis of the Mariner 9 radio-tracking data shows that the Martian gravity field is rougher than that of earth or the moon, and that the accepted direction of the Mars rotation axis is in error by about 0.5 deg. Contours of equivalent surface heights deduced from a sixth-degree solution for the Martian gravity field are presented. These contours represent the deviations from sphericity of a uniformly dense body with an external potential which is given by the first sixth-degree solution. In addition to Doppler observations, ranging or group-delay measurements have been made regularly since orbit insertion.

Lorell, J.↗