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 235 records · Page 13

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↗

TCC in the interior of moduli space and its implications for the string landscape and cosmology

We consider the classical Friedmann-Robertson-Walker solutions that describe a universe undergoing a transition from an accelerating expansion phase in the past to an eternal decelerating expansion phase in the future, driven by a scalar field evolving in a potential energy landscape. We show that any solution for which the accelerating phase violates the Trans-Planckian Censorship Conjecture (TCC), even in the interior of moduli space, never approaches the asymptotic vacuum with zero particles. Based on the assumption that the effective field theory must be valid for the vacuum on the asymptotic boundary, as motivated by holography and string theory, we argue that (multi-field) scalar potentials with such solutions are disallowed, thus strengthening the case for TCC. In particular, assuming the regularity of the future vacuum state in the string landscape, we derive results that imply a new set of highly-nonlinear constraints across the string landscape which in the absence of certain meta-stable vacua make realizing inflation impossible.

Cosmological models↗

SmallTAL: Real-Time Egocentric Online Temporal Action Localization for the Data-Impoverished

Abstract We propose a real-time, online temporal action localization system that requires a small amount of annotated data. The main challenges we address are high intra-class variability and a large and diverse background class. We address these using a flexible frame descriptor, dynamic time warping, and a novel approach to database construction. Our solution receives egocentric RGB-D streams as input and makes predictions at regular temporal intervals. We validate our approach by localizing actions in a digital twin of an electrical substation, in which certain objects have been replaced by functional virtual replicas.

Computer Science↗

The eXtended virtual element method for elliptic problems with weakly singular solutions

This paper introduces a novel eXtended virtual element method, an extension of the conforming virtual element method. The X-VEM is formulated by incorporating appropriate enrichment functions in the local spaces. The method is designed to handle highly generic enrichment functions, including singularities arising from fractured domains. By achieving consistency on the enrichment space, the method is proven to achieve arbitrary approximation orders even in the presence of singular solutions. The paper includes a complete convergence analysis under general assumptions on mesh regularity, and numerical experiments validating the method’s accuracy on various mesh families, demonstrating optimal convergence rates in the L 2 - and H 1 - norms on fractured or L-shaped domains.

97 MATHEMATICS AND COMPUTING↗

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↗

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.↗

Anisotropic physics-regularized interpretable machine learning of microstructure evolution

Anisotropic Physics-Regularized Interpretable Machine Learning Microstructure Evolution (APRIMME) is a general-purpose machine learning solution for grain growth simulations. In prior work, PRIMME employed a deep neural network to predict site-specific migration as a function of its neighboring sites to model normal, isotropic, grain growth behavior. This work aims to extend this method by incorporating grain boundary misorientation-based grain growth behavior. APRIMME is trained on anisotropic simulations created using the Monte Carlo-Potts (MCP) model. Furthermore, the results of this work are compared statistically using grain radius, number of sides per grain, mean neighborhood misorientations, and the standard deviation of triple junction dihedral angles, and are found to match in most cases. The exceptions are small and seem to be related to two causes: (1) the deterministic model of APRIMME is learning from the stochastic simulations of MCP, which seems to accentuate triple junction behaviors; and, (2) a bias against very small grains is made evident in a quicker decrease in grains than expected at the beginning of an APRIMME simulation. APRIMME is also evaluated for its general ability to capture anisotropic grain growth behavior by first investigating different test case initial conditions, including a circle grain, three grain, and hexagonal grain microstructures.

36 MATERIALS SCIENCE↗

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↗

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.↗

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↗

Using Apptainer in a Pilot-based Distributed Workload

GlideinWMS is a pilot and pressure-based workload manager for distributed scientific computing. Many experiments like CMS and Fermilab’s Neutrino experiments use it to provision elastic clusters for their analysis and simulations, split into close to a million concurrent jobs. Most user jobs require containers, and the pilots use Apptainer to set up the desired platform. For the pilots that run as regular batch jobs, Apptainer is safer, lighter, and easier to use than other containerization solutions. Many images used by the pilots are expanded SIF images distributed via the CernVM-FS: this combination is very efficient. At Fermilab, for example, we store on GitHub Dockerfiles that mimic the platform in the worker nodes of local clusters. GitHub workflows build and push the images to Docker Hub, and a service periodically pulls and converts them to the expanded SIF images in the CernVM-FS, so the scientists can find a familiar environment everywhere. Apptainer has also been used to run services inside the pilot jobs, like benchmarks that characterize the worker node being used, or a Triton Inference Server that allows sharing a GPU with all the jobs that run in parallel on a node.

Mambelli, Marco [Fermilab] (ORCID:0000000294892681↗

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.↗

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.↗

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.↗

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↗

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↗

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↗

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↗