Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model selection”

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 685 records · Page 38

Opportunities for wave energy in bulk power system operations

Wave energy resources have high, yet largely untapped potential as candidate generation technology. In this paper, we perform a data-driven analysis to characterize the impact of wave energy integration on bulk-scale power systems and market operations. Through data-driven sensitivity studies centered on an optimization-based production cost modeling formulation, our work characterizes the inflection point beyond which wave integration starts impacting power system operations, considering present day transmission infrastructure. Furthermore, our analysis also considers the joint effects of wave energy integration and system-wide transmission expansion. Finally, potential resilience scenarios such as wildfire-driven transmission contingencies and heat wave events are investigated, whereby the contributions of grid-integrated wave energy in alleviating the effects of the resilience events are analyzed. As our demonstration test bed, we consider a reduced-order network topology for the U.S. Western Interconnection with wave energy generation integrated at carefully selected sites across the coastal areas of Washington, Oregon, and northern California. Our results indicate that over a representative year of operations, wave energy integration systematically reduces locational marginal prices (LMPs) of energy and price volatility, especially during periods of high wave resource availability (winter months for the U.S. west coast). Average, maximum, and minimum of hourly LMPs over a typical year of operation was reduced by 2.95, 51.28, and 1.13 $\$$/MWh respectively (over a baseline scenario with no wave energy integration), when the selected network model had a total of 5000 MW wave power installed capacity during the representative year of study. The effects of wave energy integration can remain localized with existing transmission infrastructure (identified to be most pronounced in the Pacific Northwest region in the example we studied). However, with concurrent transmission expansion, the impacts of wave energy integration are likely to have a higher geographical spread. Our results also indicate that wave energy may be able to assist power system operations during resilience events such as major transmission contingencies and heat wave events, although such benefits might be dependent on factors such as proximity of affected area to wave resources, availability of adequate resource potential and adequate transmission capacity.

16 TIDAL AND WAVE POWER↗

Effects of gamma-Radiation on Select Lipids and Antioxidants

Radiation encountered on an extended duration space mission (estimates of 3 Sieverts for a mission to Mars) poses a threat not only to human health, but also to the quality, nutritional value, and palatability of the food system. Free radicals generated by radiation interaction with foods may initiate many unwanted reactions including: 1) autoxidation in lipids that alters flavor, odor, and concentrations of essential fatty acids, and 2) depletion of antioxidants food products and dietary supplements. Studies have shown that antioxidants may provide long term health protection from oxidative stress caused by radiation exposure; therefore, consumption of antioxidants will be important. Stability of essential fatty acids is also important for astronauts long-term health status. The objectives of this study were to characterize the effects of low dose gamma-radiation on lipids and antioxidants by monitoring oxidation and reducing power, respectively, in model systems. Select oils and antioxidants were exposed to levels of gamma-radiation ranging from 0 to 1000 Gy (1 Gy = 1 Sv) using a Gammacell 220 and stored at ambient or elevated temperatures (65 C) for up to 3 months prior to analysis. A Fricke dosimeter was used to verify differences between the radiation doses administered. Primary and secondary products of lipid oxidation in soybean and peanut oils were monitored using conjugated diene and 2-thiobarbituric acid (TBARs) assays. Changes in fatty acid composition and formation and vitamin E levels were also measured. The reducing power of antioxidant compounds, including vitamins C and E and beta-carotene, was determined using the ferric reducing antioxidant power (FRAP) assay. Significant differences (alpha =0.05) were present between all radiation doses tested using the Fricke dosimeter. Increasing radiation doses above 3 Sv resulted in significantly (alpha =0.05) elevated levels of oxidation and free fatty acids in soybean and peanut oils. Decreases in concentrations of essential fatty acids upon increasing radiation doses were also observed. Increasing radiation doses caused significant (alpha =0.05) decreases in reducing power and hence the effectiveness of vitamins C and E as well as beta-carotene. This work establishes a need for quantifying the effects of space relevant radiation doses in the development of a food system for an extended duration mission and for identifying threshold radiation levels that will impact the useful shelf-life of the variety of foods that will be sent. Eventual rancidity of lipids and the loss of antioxidant bioprotective effects are major concerns for the acceptability and nutritional profile of a food system.

Gandolph, Jacob↗

Global Network of Slow Solar Wind

The streamer belt region surrounding the heliospheric current sheet (HCS) is generally treated as the primary or sole source of the slow solar wind. Synoptic maps of solar wind speed predicted by the Wang-Sheeley-Arge model during selected periods of solar cycle 23, however, show many areas of slow wind displaced from the streamer belt. These areas commonly have the form of an arc that is connected to the streamer belt at both ends. The arcs mark the boundaries between fields emanating from different coronal holes of the same polarity and thus trace the paths of belts of pseudostreamers, i.e., unipolar streamers that form over double arcades and lack current sheets. The arc pattern is consistent with the predicted topological mapping of the narrow open corridor or singular separator line that must connect the holes and, thus, consistent with the separatrix-web model of the slow solar wind. Near solar maximum, pseudostreamer belts stray far from the HCS-associated streamer belt and, together with it, form a global-wide web of slow wind. Recognition of pseudostreamer belts as prominent sources of slow wind provides a new template for understanding solar wind stream structure, especially near solar maximum.

solar wind↗

Comparison of linear regression, k-nearest neighbour and random forest methods in airborne laser-scanning-based prediction of growing stock

Abstract In this study, for five sites around the world, we look at the effects of different model types and variable selection approaches on forest yield modelling performances in an area-based approach (ABA). We compared ordinary least squares regression (OLS), k-nearest neighbours (kNN) and random forest (RF). Our objective was to test if there are systematic differences in accuracy between OLS, kNN and RF in ABA predictions of growing stock volume. The analyses are based on a 5-fold cross-validation at five study sites: an eucalyptus plantation, a temperate forest and three different boreal forests. Two completely independent validation datasets were also available for two of the boreal sites. For the kNN, we evaluated multiple measures of distance including Euclidean, Mahalanobis, most similar neighbour (MSN) and an RF-based distance metric. The variable selection approaches we examined included a heuristic approach (for OLS, kNN and RF), exhaustive search among all combinations (OLS only) and all variables together (RF only). Performances varied by model type and variable selection approaches among sites. OLS and RF had similar accuracies and were more efficient than any of the kNN variants. Variable selection did not affect RF performance. Heuristic and exhaustive variable selection performed similarly for OLS. kNN fared the poorest amongst model types, and kNN with RF distance was prone to overfitting when compared with a validation dataset. Additional caution is therefore required when building kNN models for volume prediction though ABA, being preferable instead to opt for models based on OLS with some variable selection, or RF with all variables together.

Cosenza, Diogo N.↗

Genomic selection in algae with biphasic lifecycles: A Saccharina latissima (sugar kelp) case study

Introduction Sugar kelp ( Saccharina latissima ) has a biphasic life cycle, allowing selection on both thediploid sporophytes (SPs) and haploid gametophytes (GPs). Methods We trained a genomic selection (GS) model from farm-tested SP phenotypic data and used a mixed-ploidy additive relationship matrix to predict GP breeding values. Topranked GPs were used to make crosses for further farm evaluation. The relationship matrix included 866 individuals: a) founder SPs sampled from the wild; b) progeny GPs from founders; c) Farm-tested SPs crossed from b); and d) progeny GPs from farm-tested SPs. The complete pedigree-based relationship matrix was estimated for all individuals. A subset of founder SPs ( n = 58) and GPs ( n = 276) were genotyped with Diversity Array Technology and whole genome sequencing, respectively. We evaluated GS prediction accuracy via cross validation for SPs tested on farm in 2019 and 2020 using a basic GBLUP model. We also estimated the general combining ability (GCA) and specific combining ability (SCA) variances of parental GPs. A total of 11 yield-related and morphology traits were evaluated. Results The cross validation accuracies for dry weight per meter ( r ranged from 0.16 to 0.35) and wet weight per meter ( r ranged 0.19 to 0.35) were comparable to GS accuracy for yield traits in terrestrial crops. For morphology traits, cross validation accuracy exceeded 0.18 in all scenarios except for blade thickness in the second year. Accuracy in a third validation year (2021) was 0.31 for dry weight per meter over a confirmation set of 87 individuals. Discussion Our findings indicate that progress can be made in sugar kelp breeding by using genomic selection.

59 BASIC BIOLOGICAL SCIENCES↗

Near Earth Object (NEO) Mitigation Options Using Exploration Technologies

This work documents the advancements in MSFC threat modeling and mitigation technology research completed since our last major publication in this field. Most of the work enclosed here are refinements of our work documented in NASA TP-2004-213089. Very long development times from start of funding (10-20 years) can be expected for any mitigation system which suggests that delaying consideration of mitigation technologies could leave the Earth in an unprotected state for a significant period of time. Fortunately there is the potential for strong synergy between architecture requirements for some threat mitigators and crewed deep space exploration. Thus planetary defense has the potential to be integrated into the current U.S. space exploration effort. The number of possible options available for protection against the NEO threat was too numerous for them to all be addressed within the study; instead, a representative selection were modeled and evaluated. A summary of the major lessons learned during this study is presented, as are recommendations for future work.

Arnold William↗

Shifts in hydroclimatology of US megaregions in response to climate change

Most of the population and economic growth in the United States occurs in megaregions as the clustered metropolitan areas, whereas climate change may amplify negative impacts on water and natural resources. This study assesses shifts in regional hydroclimatology of fourteen US megaregions in response to climate change over the 21st century. Hydroclimatic projections were simulated using the Variable Infiltration Capacity (VIC) model driven by three downscaled climate models from the Multivariate Adaptive Constructed Analogs (MACA) dataset to cover driest to wettest future conditions in the conterminous United States (CONUS). Shifts in the regional hydroclimatolgy and basin characteristics of US megaregions were represented as a combination of changes in the aridity and evaporative indices using the Budyko framework and Fu's equation. Changes in the climate types of US megaregions were estimated using the Fine Gaussian Support Vector Machine (SVM) method. The results indicate that Los Angeles, San Diego, and San Francisco are more likely to experience less arid conditions with some shifts from Continental to Temperate climate type while the hydroclimatology of Houston may become drier with some shifts from Temperate to Continental climate type. Additionally, water yield is likely to decrease in Seattle. Change in the hydroclimatology of Denver and Phoenix highly depends on the selected climate model. However, the basin characteristics of Phoenix have the highest sensitivity to climate change. Overall, the hydroclimatic conditions of Los Angeles, San Diego, Phoenix, Denver, and Houston have the highest sensitivity to climate change. Understanding of future shifts in hydroclimatology of megaregions can help decision-makers to attenuate negative consequences by implementing appropriate adaptation strategies, particularly in the water-scare megaregions.

54 ENVIRONMENTAL SCIENCES↗

Local Chain Dynamics in Sequence-Controlled Polymers as a Tunable Handle for Rare Earth Sequestration

Chain dynamics govern the intricate behaviors of proteins, underpinning functions such as catalysis, recognition, and stimulus response, and are an increasingly appreciated aspect of structure–function relationships. Analogously, manipulating chain dynamics and structure in abiotic polymers via sequence control is an exciting, yet underexplored, strategy for improving material functions. In this work, we report a systematic study relating the sequence of polymeric sequestrants to their structure and dynamics, as well as to their binding affinity and selectivity for model substrates, rare earth elements (REEs). A series of sequence-controlled polymers with metal chelating, solubilizing, and structure forming monomers was synthesized via multiblock polymerization, yielding compositionally identical polymers with spectroscopically resolved domains and distinct morphologies. Using a combination of small-angle X-ray scattering and 19 F NMR relaxometry measurements, we connected differences in polymer structure and dynamics to polymer sequence variables such as the patchiness (density) of the structure forming monomer and the location of the chelating monomer. Furthermore, we found that, relative to calcium, all polymers in the series collapse more and have slower dynamics when binding REEs (lanthanum and lutetium) , though the extent of these effects were sequence-dependent and localized to specific domains within the polymer. Notably, sequence-controlled polymers that exhibited the largest conformational and dynamic changes upon binding REEs also bound REEs with the greatest affinity and modest selectivity. Collectively, these results correlate monomer patterning with dynamics, morphology, and REE binding performance en route to the development of efficient and selective macromolecular chelators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Search for Long-Lived Particles with Muon Detector Shower Signature in the CMS Run-3 data

Many beyond standard model theories predict the existence of long-lived particles (LLPs). These LLPs can have sizable lifetimes and decay several meters from their production vertex. In this poster/talk, we present the analysis strategy for searching for LLPs using the Compact Muon Solenoid (CMS) Experiment. The proton-proton collision data used in this analysis were collected from 2022 to 2024 at a center-of-mass energy of 13.6 TeV, corresponding to an integrated luminosity of 170 fb^-1. The LLP decays are reconstructed as a high-multiplicity cluster of detector hits in the cathode strip chambers (CSC) of the muon system endcap. This signature is referred to as the Muon Detector Showers (MDS). This search requires events to contain at least one MDS cluster. The analysis focuses on LLP hadronic decays and LLP masses up to a few tens of GeV. We present the signal properties in MonteCarlo simulation, event selection, background modeling, and evaluation of the expected sensitivity. The results are interpreted under the Twin Higgs model benchmark.

Agyemang-Duah, Andrews [Grambling State U.]↗

A phase transition for finding needles in nonlinear haystacks with LASSO artificial neural networks

To fit sparse linear associations, a LASSO sparsity inducing penalty with a single hyperparameter provably allows to recover the important features (needles) with high probability in certain regimes even if the sample size is smaller than the dimension of the input vector (haystack). More recently learners known as artificial neural networks (ANN) have shown great successes in many machine learning tasks, in particular fitting nonlinear associations. Small learning rate, stochastic gradient descent algorithm and large training set help to cope with the explosion in the number of parameters present in deep neural networks. Yet few ANN learners have been developed and studied to find needles in nonlinear haystacks. Driven by a single hyperparameter, our ANN learner, like for sparse linear associations, exhibits a phase transition in the probability of retrieving the needles, which we do not observe with other ANN learners. To select our penalty parameter, we generalize the universal threshold of Donoho and Johnstone (Biometrika 81(3):425–455, 1994) which is a better rule than the conservative (too many false detections) and expensive cross-validation. In the spirit of simulated annealing, we propose a warm-start sparsity inducing algorithm to solve the high-dimensional, non-convex and non-differentiable optimization problem. We perform simulated and real data Monte Carlo experiments to quantify the effectiveness of our approach.

97 MATHEMATICS AND COMPUTING↗

Orographic enhancement of rainfall over the Congo Basin

Abstract The Congo rainforest located in central equatorial Africa is an important, yet understudied part of the globe surrounded by complex orographic features. A primitive understanding of precipitation processes such as mesoscale convective dynamics magnifies uncertainties in the future climate projections of the hydrological cycle over the Congo. Furthermore, the effects of orography, which is an important forcing for convection and precipitation, are poorly resolved by climate models, and ill‐conceptualized over the Congo. To address this knowledge gap, perturbed orographic forcing experiments are conducted using the Weather Research and Forecasting (WRF) mesoscale numerical model in a high‐resolution convection‐permitting model setup. The model simulated selected dates in November 2014. The thunderstorms and rainfall simulated in the control run for the case study analyzed in this article compared reasonably well to satellite‐derived brightness temperature and rainfall data. The results from this case study show that the dynamical impact of increasing the height of the East African Highlands is the blocking of the lower‐tropospheric tropical easterlies. This weakening of the lower‐tropospheric zonal winds increases the windshear over the Congo Basin resulting in slower propagating, more intense mesoscale convective systems with enhanced rainfall.

54 ENVIRONMENTAL SCIENCES↗

Image-based terrain modeling with thematic mapper applied to resolving the limit of Holocene Lake expansion in the Great Salt Lake Desert, Utah, part 1

The LANDSAT Thematic Mapper (TM) scanner records reflected solar energy from the earth's surface in six wavelength regions, or bands, and one band that records emitted energy in the thermal region, giving a total of seven bands. Useful research was extracted about terrain morphometry from remote sensing measurements and this information is used in an image-based terrain model for selected coastal geomorphic features in the Great Salt Lake Desert (GSLD). Technical developments include the incorporation of Aerial Profiling of Terrain System (APTS) data in satellite image analysis, and the production and use of 3-D surface plots of TM reflectance data. Also included in the technical developments is the analysis of the ground control point spatial distribution and its affects on geometric correction, and the terrain mapping procedure; using satellite data in a way that eliminates the need to degrade the data by resampling. The most common approach for terrain mapping with multispectral scanner data includes the techniques of pattern recognition and image classification, as opposed to direct measurement of radiance for identification of terrain features. The research approach in this investigation was based on an understanding of the characteristics of reflected light resulting from the variations in moisture and geometry related to terrain as described by the physical laws of radiative transfer. The image-based terrain model provides quantitative information about the terrain morphometry based on the physical relationship between TM data, the physical character of the GSLD, and the APTS measurements.

Merola, John A.↗

Evaluation of Spectral Band Adjustment Factors for Cross-Calibration of Visible Imagers

The CERES EBAF dataset provides TOA SW and LW fluxes for long-term monitoring of the Earth’s energy balance and to validate climate models. The EBAF products, based on the Terra and Aqua CERES instrument observed radiances, rely on coincident measurements from the onboard MODIS imagers to determine cloud properties that are used in the angular distribution model scene selection required to convert the CERES observed radiances into flux values. Once the Terra and Aqua orbits start drifting outside of their 15-minute window, the CERES project will rely solely on the NOAA-20 CERES observations. A seamless transition of fluxes and clouds can only occur if the analogous MODIS and VIIRS channels are properly inter-calibrated. The visible band (0.65 µm) in MODIS (band 1) and VIIRS (M5 or I1) is critical for retrieving the cloud mask and optical depth. The spectral response functions (SRF) of these bands differ noticeably and will require scene dependent spectral band adjustment factors (SBAF) for proper radiometric scaling between them. Fortunately, the VIIRS I1 and M5 bands are calibrated using the same solar diffuser as a reference, which implies that band reflectance differences must be due to the spectral disparity. The coincident I1 and M5 TOA reflectance measurements provide the optimal opportunity to validate SBAFs over many surface and cloud conditions. The CERES project maintains SCIAMACHY, GOME-2, Hyperion, and radiative transfer model-based scene-stratified hyper-spectral reflectance measurements that can be convolved with sensor pair SRFs to compute the corresponding SBAF. This study will evaluate the SCIAMACHY, GOME-2, Hyperion, and RTM based SBAFs over tropical ocean targets used in SNO inter-calibration of MODIS and VIIRS, including clear-sky ocean, deep convective clouds (DCC), and liquid cloud targets in terms of their applicability to absolute intercomparison of visible imagers. The large SCIAMACHY and GOME-2 footprints and the lack of Hyperion cloudy scenes will impact the resulting SBAFs and necessitates the need for future CLARREO hyper-spectral reflectances to decrease SBAF uncertainties Preliminary results indicate that the SCIAMACHY and GOME-2 based SBAFs for the VIIRS I1 and M5 band pairs may differ by 1.5% for some ATO scene types. Because most of the modern GEOs visible band SRFs encompass the VIIRS I1 band SRF, these evaluations are critical to ensure that the MODIS, VIIRS, and GEO clouds and fluxes are consistent. SBAFs for Libya-4, Dome-C and other Earth invariant target will also be evaluated.

D. R. Doelling↗

Correlating binding energies of adsorbed CO and H on model surfaces with CO/H 2 selectivity from co-electrolysis of CO 2 and H 2 O over copper–palladium bimetallic catalysts

Binding energies of adsorbed CO and H are key descriptors governing the activity and selectivity of the co-electrolysis of CO 2 and H 2 O to produce syngas with desired CO/H 2 ratios. Palladium hydride (PdH), which forms in situ at negative overpotentials, has been identified as the active Pd phase for CO 2 reduction to syngas. Herein, binding energies of CO and H are determined using temperature programmed desorption (TPD) of CO and H 2 from Pd(111), PdH/Pd(111), and Cu/PdH/Pd(111) under ultra-high vacuum (UHV) conditions. TPD results reveal that desorption of H 2 from subsurface PdH occurs at 460 K, while desorption from surface PdH is more facile at 320 K. CO desorption temperatures shift 20 K lower on PdH/Pd(111) compared to on Pd(111). The presence of 0.7 ML Cu further increases the desorption temperature of H 2 by 30 K while simultaneously reducing CO desorption temperatures by 70 K. Density functional theory (DFT) calculations show that CO adsorption onto Pd sites is hindered on the 0.7 ML Cu/PdH/Pd(111) surface while the kinetic barrier for H 2 desorption is increased. The trends in the binding energies of CO and H on model surfaces are consistent with electrochemical measurements of CuPd powder catalysts in a membrane electrode assembly (MEA), where H 2 evolution is reduced while CO production is enhanced compared to unmodified Pd catalysts. Overall, the results from model surface studies (TPD and DFT) provide a prediction and explanation for the activity and CO/H 2 ratios observed in electrochemical experiments. This study also demonstrates that CuPd is a promising catalyst with reduced Pd-loading to produce CO-rich syngas.

36 MATERIALS SCIENCE↗

A Generic and Multifunctional Electromagnetic Transient Model for Grid-Following Inverters

This paper presents a generic and multifunctional electromagnetic transient (EMT) dynamic model of grid-following (GFL) inverter-based resources (IBRs) using the PSCAD software platform. The features of the model include flexibility in selecting various types and combinations of DC sources covering photovoltaic modules, battery modules, and ideal DC-source modules as well as flexibility in selecting either switching or averaged models of the inverter. This model also covers exhaustive lists of controller algorithm, including open-loop/closed-loop PQ dispatch control, DC voltage and AC terminal voltage control, and conventional current control designed in the dq-domain, the ..alpha....beta.. -domain, and the positive-/negative-sequence domain. Moreover, this model is equipped with flexibility in selecting various types of current-limiting schemes, including saturation-based and latching-based current limiters, and anti-windup protection. Also, the EMT model is agnostic to the MVA rating and is suitable for interfacing transmission systems by complying with IEEE Std. 2800. The generality in the power circuits and the multifunctional options in the operation and control of the developed EMT model make it suitable for both academia and industry to study various power system aspects, including, but not limited to, the fault behavior of GFL IBRs, the impacts on the protection system, and the transient stability of a system interfaced with large numbers of GFL IBRs.

integrated circuit modeling↗

High-Resolution Wind Resource Data Set of the Greater Puerto Rico Region

In February 2022, the U.S. Department of Energy and six national laboratories launched the Puerto Rico Grid Resilience and Transitions to 100% Renewable Energy Study (PR100). PR100 aims to provide a comprehensive analysis of possible pathways for Puerto Rico's energy future, with a goal of 100% renewable energy by 2050. As a part of the renewable energy potential assessment in this project, we developed 20 years (2001-2020) of data using a numerical weather prediction (NWP) model for onshore and offshore wind resource assessment for the Puerto Rico. The research steps in developing the long-term wind resource data sets based on the NWP model were: 1. Model wind resource based on the Weather Research and Forecasting (WRF) model. 2. Develop WRF model configurations for Puerto Rico. 3. Test WRF with 11 different physics parameterizations for planetary boundary layer (PBL). 4. Assess WRF output from the different PBL schemes against observations. 5. Select a final model configuration which can produce the modeled wind speed with sufficient accuracy. 6. Produce 20 years wind resource data sets for Puerto Rico region. In the first stage of our framework for developing wind resource data, we developed a WRF model configuration using two nested domains (9 km and 3 km) to cover Puerto Rico and U.S. Virgin Islands and downscale the ERA5 reanalysis data (0.25 degrees x 0.25 degrees; hourly interval) to a 3-km domain. For the second stage, we implemented one-year simulations focused on using 11 different PBL physics parameterizations to find a combination of WRF physics parameterizations that could provide accurately modeled wind speed for Puerto Rico. We also analyzed the sensitivity of the modeled wind speed to PBL schemes for onshore and offshore locations. The WRF output resulting from the 11 WRF experiments using different PBL parameterizations were evaluated against observations obtained from the National Data Buoy Center (NDBC) as well as at hub height for a location for which measurements were available. A final model setup selected through the validation with observational data was used to produce 20 years of data with 3-km spatial and 5-minute temporal resolution. The WRF model output was post-processed to include wind profiles and basic atmospheric variables in a format that can be easily used for downstream modeling. The 20 years of wind resource data will be made available through NREL and support the estimation of wind energy development costs for the PR100 study.

17 WIND ENERGY↗

Identifying Differential Equations in Fourier Domain (FourierIdent)

We investigate identifying differential equations in the frequency domain. Fourier analysis is an important tool in theoretical analysis and numerical solvers of differential equations, yet there is limited work in exploring this connection in the identification of differential equations. This paper aims to identify the underlying differential equation in the frequency domain, from a given single realization of the differential equation perturbed by noise. Such setting imposes difficulties which are different from other identification methods where computation is carried out in the physical domain. We propose several ways to mitigate the challenges arising from noise in data and large differences in the magnitudes of frequency responses. The main takeaways are that identifying differential equations solely in the frequency domain is challenging, the method we propose is based on a form of domain partitions in the frequency domain, and this method shows benefits for complex data even with high level of noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain and to enhance the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms, and shows advantages on complex data with many frequency modes, even under high level of noise.

97 MATHEMATICS AND COMPUTING↗