Articles (NTIS)
Permanent URI for this collection
Browse
Recent Submissions
Item On the Accurate Determination of the Orthometric Correction to Levelled Height Differences - A Case Study in Hong Kong(2025) Tenzer, Robert; Ababio, Albertini Nsiah; Foroughi, Ismael; Pitoňák, Martin; Novák, Pavel; Chen, Wenjin; Ghomsi, Franck Eitel KemgangOrthometric heights are practically determined from levelling and gravity measurements by applying orthometric corrections to levelled height differences. Currently, Helmert’s definition of orthometric heights is mostly used, with the mean gravity computed only approximately from observed surface gravity by applying the Poincaré-Prey gravity reduction. In this study, we compute orthometric corrections by applying the state-of-the-art method that utilizes advanced numerical procedures to account for the topographic relief and 3D mass density variations, while adopting the Earth’s spherical approximation. The non-topographic contribution of masses inside the geoid is evaluated by solving geodetic boundary-values problems. We apply this numerical procedure for the first time to practically determine orthometric heights of levelling benchmarks from levelling and gravity measurements and digital terrain and rock density models. Results obtained after readjustment of newly determined orthometric heights at the levelling network covering the Hong Kong territories are compared with Helmert’s orthometric heights. Comparison reveals that errors in Helmert’s orthometric heights vary between -3.13 and 0.95 cm. Such errors are very significant when compared to accurate values of the cumulative orthometric correction between -1.88 and 0.84 cm. Moreover, large errors (up to 1 cm) occur already at levelling benchmarks at very low elevations (< 100 m). These findings demonstrate that an accurate determination of orthometric heights is crucial even for regions with moderately elevated topography.Item Geoid determination using airborne gravity vectors(2025) Goli, Mehdi; Foroughi, Ismael; Pagiatakis, Spiros; Ferguson, Stephen; Novák, PavelIn traditional airborne gravimetry, the vertical component of the gravity vector is used as an approximation of the measured magnitude of the gravity vector, which enters the determination of the local geoid. In this study, a comprehensive computational scheme for determining the local geoid using three components of the airborne gravity vector is presented. Our approach extends the existing one-step method for local geoid modelling by incorporating the full gravity vector measured by airborne sensors as boundary values in the gravimetric boundary-value problem. We derive integral kernel functions along with far-zone contributions for the three components of the airborne gravity vector and apply deterministic modifications to them. To validate our derivations, we use GGM-based airborne gravity vectors burdened with realistic coloured noise at one of the most challenging test sites for geoid determination, the 1-cm geoid test area in Colorado (USA). Results of closed-loop tests confirm that applying all three components of the GGM-based airborne gravity vector improves the internal accuracy of the geoid by 50 % compared to using only the vertical component. We further use real airborne gravity vectors observed at a test site in the same region and show that the RMS of the estimated geoid heights evaluated against the reference geoidal heights along the GSVS17 line is 2.3 cm using the “traditional approach” and 1.3 cm including the horizontal components. This indicates a significant improvement in the external accuracy (~46 %) of the geoid when the full gravity vector is used, without using other heterogeneous observations.Item Three Years of VoiceMOS Challenges: Lessons Learned by the UWB-NTIS-TTS Team(2025) Kunešová, Marie; Matoušek, Jindřich; Lehečka, Jan; Švec, Jan; Tihelka, Daniel; Hanzlíček, ZdeněkAutomatic prediction of mean-opinion scores (MOS) promises a faster, cheaper alternative to listening tests, yet robust generalization across speakers, languages, and domains remains a significant challenge. This article presents our system designs and experimental results from three years of participation in the VoiceMOS Challenges (2022–2024), covering MOS prediction for synthesized or voice-converted speech and singing voice, including out-of-domain and cross-language conditions. We evaluate six neural architectures – wav2vec 2.0, QuartzNet, CNN-RNN, LDNet, RawNet3, and HiFi-GAN – and their ensembles. Across all tasks, we find that 1) self-supervised acoustic encoders are the most consistently reliable foundation, 2) ensembling yields rapidly diminishing returns once complementary representations are covered, and 3) the diversity and balance of training data outweigh architectural complexity. Notably, the indiscriminate fusion strategy that performed well in 2022 degrades under the mismatched French TTS conditions of 2023, emphasizing the importance of out-of-domain validation. Further experiments show that carefully pruned ensembles can modestly outperform the best single model while remaining within real-time constraints. We conclude with several observations to guide the development of computationally efficient, domain-robust MOS prediction systems.Item Modeling of groundwater flow in porous medium layered over inclined impermeable beds(2025) Girg, Petr; Kotrla, LukášWe propose a new mathematical model of groundwater flow in porous medium layered over inclined impermeable beds. In its full generality, this is a free-surface problem. To obtain analytically tractable model, we use generalized Dupuit-Forchheimer assumption for inclined impermeable bed. In this way, we arrive at parabolic partial differential equation which is a generalization of the classical Boussinesq equation. The novelty of our approach consists in considering nonlinear constitutive law of the power type. Thus introducing p-Laplacian-like differential operator into the Boussinesq equation. Unlike in the classical case of the Boussinesq equation, the convective term cannot be set aside from the main part of the diffusive term and remains incorporated within it. In this article, we analyze qualitative properties of the stationary solutions of our model. In the case of p > 2, we study validity of Weak and Strong Maximum Principles as well. We use methods based on the linearization of the p-Laplacian-type problems in the vicinity of known solution, error estimates, and analysis of Green's function of the linearized problem.Item Functionally-Equivalent Formalization and Automated Model Checking of Function Block Diagrams(2025) Ausberger, Tomáš; Kubíček, Karel; Medvecová, PavlaIn the development and verification of safety-critical and safety-related Instrumentation and Control (I&C) systems, it is essential to ensure there is no deviation from the requirements of the assignment during development. Model checking is a method of formal verification which can be used to prove whether a formal model satisfies its formal requirement. Since algorithms of I&C systems are generally informal, they can not be verified by model checking directly, but they must be carefully translated. This article presents a new method based on functionally-equivalent formalisation and model checking. This method can be used for automatic verification of I&C algorithms by model checking while preserving obtained proofs from a formalised model in the original algorithm. There are several problems associated with the verification of PLCs by model checking: 1) State space explosion, 2) Model consistency, 3) Specifying Properties to be Checked, 4) Representing PLC execution cycle, 5) TONs timers representation. This aims of this article is to address points 1), 2), 4) and 5). The article also presents conditions for implementing these algorithms in a target I&C system under which the obtained proofs can also be expected in the physical I&C system. This article primarily focuses on formalisation and model checking of Function Block Diagram (FBD) algorithms. However, the presented methods can also be extended to other programming languages.Item First-principles investigation of gas adsorption on bilayer transition metal dichalcogenides for sensing toxic gases(2025) Damte, Jemal Yimer; Ataalite, HassanTransition metal dichalcogenides (TMDs) have shown significant promise in gas sensing applications due to their high catalytic activity and unique electronic properties, which facilitate effective interactions with various gas molecules. This makes them ideal candidates for high-performance gas sensors. In this study, we investigated the sensing properties of nitrogen-containing gases (NCGs) on several heterostructures-namely, MoS2/WTe2, MoTe2/WS2, MoS2/TiO2 and MoS2/IrO2-using density functional theory calculations. The results indicate that NH3 and NOX exhibit weak electronic interactions with MoS2/WTe2 and MoTe2/WS2 heterostructures, while strong electronic interactions are observed with MoS2/TiO2 and MoS2/IrO2 heterostructures. Electron transport properties were further assessed using Non-Equilibrium Green’s Function calculations, revealing promising gas sensing characteristics for NH3 detection across all heterostructures and particularly effective NOX detection with MoS2/WTe2 and MoTe2/WS2 heterostructures. These findings highlight the potential of MoS2/WTe2 and MoTe2/WS2 as sensitive and selective gas sensors for both NH3 and NOX, providing valuable insights for developing advanced gas-sensing technologies with diverse practical applications.Item Challenges in Automatic Differentiation and Numerical Integration in Physics-Informed Neural Networks Modelling(2025) Daněk, Josef; Pospíšil, JanIn this paper, we numerically examine the precision challenges that emerge in automatic differentiation and numerical integration in various tasks now tackled by physics-informed neural networks (PINNs). Specifically, we illustrate how ill-posed problems or inaccurately computed functions can cause serious precision issues in differentiation and integration. A major difficulty lies in detecting these problems. A simple large-scale view of the function or good-looking loss functions or convergence results may not reveal any potential errors, and the resulting outcomes are often mistakenly considered correct. To address this, it is often critical to determine whether standard double-precision arithmetic suffices or if higher precision is necessary, but using higher precision arithmetic with neural networks does not have to bring an improvement at all. Three problematic use-cases for solving differential equations using PINNs are analyzed in detail. For the case requiring numerical integration, we also evaluate several numerical quadrature methods and suggest particular numerical analysis steps to choose the most suitable method.Item Dihedral tilings of the sphere by regular polygons and quadrilaterals I: squares and rhombi(2025) Luk, Hoi PingWe classify the edge-to-edge dihedral tilings of the sphere by squares and rhombi with a systematic procedure applicable to similar problems.Item Some Bounds on the Threshold Dimension of Graphs(2025) Francis, Mathew C.; Majumder, Atrayee; Mathew, RogersThe threshold dimension of a graph G is the minimum number of threshold graphs whose intersection yields G. We give tight or nearly tight upper bounds for the threshold dimension of a graph in terms of various graph parameters including treewidth, maximum degree, degeneracy, number of vertices, and vertex cover number. We also study threshold dimension of random graphs and graphs with high girth.Item WO3/CuWO4 nanocomposite thin films for humidity resilient acetone gas sensing(2025) Kumar, Nirmal; Kumar, Akash; Čapek, Jiří; Comini, Elisabetta; Haviar, StanislavWe demonstrate a high-performing and selective acetone gas sensor based on WO3/CuWO4 nanocomposites produced by sequentially sputter-deposited WO3 thin films and CuO nanoparticles, engineered to reduce the humidity interference. By optimizing deposition order and layer thicknesses, we harnessed synergistic p-n/n-n heterojunctions and the formation of a catalytic CuWO4 ternary phase. The best performing configuration (20 nm of tungsten oxide film on top of nanoparticles) exhibits a high response (S = 23) to 10 ppm acetone at 300 °C, fast response in dry/humid conditions (38 s/58 s), and low detection limit (0.6 ppm). More importantly, the sensor exhibited > 95 % retention of its response in 90 % relative humidity compared to a loss of >50 % for pristine WO3. The reduced humidity interference is assigned to heterojunction formation at the WO3/CuWO4 interface, Lewis acid sites that allow for acetone selective adsorption, and bulk-dominated conduction. This noble-metal-free acetone sensor overcomes a known shortcoming of metal oxide-based sensors, enabling accurate acetone detection in humid environments for breath-based disease diagnosis (e.g., diabetes) and industrial safety monitoring.Item Ultra-low-resistivity nitrogen-doped p-type Cu2O thin films fabricated by reactive HiPIMS(2025) Rezek, Jiří; Koloros, Jan; Houška, Jiří; Čerstvý, Radomír; Haviar, Stanislav; Kolenatý, David; Damte, Jemal Yimer; Baroch, PavelWe have successfully fabricated the nitrogen-doped cuprous oxide thin films on the amorphous standard soda-lime glass by reactive high-power impulse magnetron sputtering. The energy of film-forming particles was controlled by the value of pulse-averaged target power density, which has a significant impact on the elemental composition, structure and optoelectrical properties of the films. We have shown that the high-energy regime is more suitable for preserving Cu2O structure and leads to continuous substitution of oxygen by nitrogen compared with the low-energy regime. Moreover, in the high-energy regime, it is possible, to some extent, to independently control the electrical resistivity and optical properties. The electrical resistivity decreases down to ≈ 0.05 ohm.cm (upper bound of the hole mobility of 0.08 ± 0.05 cm²/Vs) at the optical band gap 2.0–2.3 eV. Special attention is paid to the formation of nitrogen molecules and their ability to form shallow acceptor states. Experimental results supported by our Density Functional Theory calculations indicate that N2 replacing Cu in the Cu2O lattice is one possible (but not the only possible) acceptor. We have also found that the formation of nitrogen molecules is preferred in a high-energy regime.Item Strongly thermochromic W-doped VO2 films with a large temperature coefficient of electrical resistance near room temperature(2025) Farrukh, Sadoon; Vlček, Jaroslav; Rezek, Jiří; Houška, Jiří; Čerstvý, Radomír; Kozák, TomášWe report the crystal structure, surface morphology, electronic band structure, optical and electrical properties, and semiconductor-metal transition characteristics of strongly thermochromic W-doped VO2 films with a large (up to −16 % 1/K) temperature coefficient of electrical resistance at a small hysteresis width of electrical resistivity (down to 3 °C) near room temperature, and with a wide temperature operation range at a high detection sensitivity (≥ 8 % 1/K) and low values of the electrical resistivity. They were deposited at a reduced substrate temperature of 350 °C onto soda-lime glass (SLG) with two versions of yttria-stabilized zirconia (YSZ) interlayers possessing different cubic crystal orientations, and onto bare SLG and monocrystalline YSZ and Al2O3 substrates for comparison. The W-doped VO2 depositions were performed using reactive deep oscillation magnetron sputtering with feedback pulsed O2 flow control, allowing us to increase the deposition rate of films up to 20–30 nm/min for a target-substrate distance of 100 mm. The results are important for further improvement of the thermochromic performance of VO2-based coatings for energy-saving smart windows and for a new design of high-performance infrared detectors and temperature sensors prepared by a fast, low-temperature, scalable synthesis.Item The radial integral of the geopotential(2025) Tenzer, Robert; Novák, Pavel; Eshagh, MehdiIn Newtonian theory of gravitation, used in Earth’s and planetary sciences, gravitational acceleration is standardly regarded as the most fundamental parameter that describes any vectorial gravitational field. Considering only conservative gravitational field, the vectorial field can be described by a scalar function of 3D position called the gravitational potential from which other parameters (particularly the gravitational attraction and the gravitational gradient) are derived by applying the gradient operators. Gradients of the Earth’s gravity potential are nowadays measured with high accuracy and applied in various geodetic and geophysical applications. In geodesy, the gravity and gravity gradient measurements are used to determine the Earth’s gravity potential (i.e., the geopotential) that is related to geometry of equipotential surfaces, most notably the geoid approximating globally the mean sea surface. Reversely to the application of gradient operator, the application of radial integral to gravity yields the gravity potential differences and the same application to gravity gradient yields the gravity differences. This procedure was implemented in definitions of rigorous orthometric heights and differences between normal and orthometric heights (i.e., the geoid-to-quasigeoid separation). Following this concept, we introduce the radially integrated gravity potential (i.e., the geopotential), and provide mathematical definitions of this functional in spatial and spectral domains. We also define its relationship with other parameters of the Earth’s gravity field via Poisson, Hotine, and Stokes integrals. We then discuss prospects of using this functional in gravimetric geophysics in the context of interpreting the Earth’s inner structure. In numerical examples, we demonstrate that the indefinite radial integral of the disturbing potential (i.e., difference between actual and normal gravity potentials) has a spatial pattern that better exhibits a long-wavelength signature of deep mantle than the global geoidal geometry. This finding is explained by the fact that a more detailed spatial pattern attributed mainly to a lithospheric structure is filtered out proportionally with increasing degree of spherical harmonics in this functional. The global geoidal geometry, on the other hand, comprises not only a deep mantle signature but eventually also a gravitational signature of lithosphere, most notably across large orogens, even after applying spectral decompensation or filtering.Item Nonlinear dynamics of ball vibration absorber considering stability, stationarity and rolling-condition boundaries(2025) Dyk, Štěpán; Bulín, Radek; Rendl, JanThe paper presents a detailed nonlinear analysis of a ball vibration absorber (BVA), which consists of a harmonic oscillator with a spherical cavity and a rolling ball as an absorber. Frequency response curves are calculated using the harmonic balance method and pseudo arc length continuation, and stability is assessed using stability analysis applied to modulation equations. This is particularly important at higher excitation amplitudes where the modulation equations provide information on the presence of strongly modulated response regimes. Codimension-2 continuation is used to identify the onset of instability and non-stationary regions with respect to all key design parameters. The study highlights the critical role of the rolling and contact conditions in maintaining the validity of the solution and provides conditions for their satisfaction. The results provide valuable insights into the non-linear dynamic behaviour of the BVA, revealing its effectiveness in vibration reduction and its limitations due to parameter selection and design constraints.Item Convolution neural network for fluid flow simulations in cascade with oscillating blades(2025) Bublík, Ondřej; Heidler, Václav; Vimmr, JanThis paper aims to design a computational model for simulating the unsteady flow field in a cascade of oscillating blades. The core of the new model is a convolutional neural network, which is trained on a simplified cascade consisting of three blades. The primary advantage lies in significantly reducing the computational cost, as the new model is several orders of magnitude faster than traditional CFD methods for evaluations, though training the model remains computationally intensive. The convolutional neural network can accurately predict the unsteady flow field, as demonstrated in validation examples. In the next step, a composition algorithm is proposed to combine several simplified cases, enabling the solution of a cascade with any number of blades.Item Thermally-induced microstructural evolution in nanoparticle-based CuO, WO3 and CuO–WO3 thin films for hydrogen gas sensing(2025) Shaji, Kalyani; Haviar, Stanislav; Zeman, Petr; Procházka, Michal; Čerstvý, Radomír; Kumar, Nirmal; Čapek, JiříThis study systematically investigates the microstructural evolution of nanoparticle-based CuO, WO3, and composite CuO–WO3 thin films induced by their post-deposition annealing. The films were reactively deposited using a magnetron-based gas aggregation technique, with the composite films consisting of alternating monolayers of CuO and WO3 nanoparticles. After deposition, the films were annealed in synthetic air at temperatures ranging from 200 to 400 °C and characterized using scanning electron microscopy, X-ray diffraction, Raman spectroscopy, and X-ray photoelectron spectroscopy. Annealing of the CuO films led to the most pronounced changes associated with a gradual enhancement of crystallinity accompanied by significant particle growth with increasing annealing temperature, while the WO3 and CuO–WO3 films were more thermally stable to crystallization and particle growth. Notably, at 400 °C, the CuO–WO3 films crystallized into a novel 𝛾-CuWO4 phase. The annealed films were further evaluated for their gas-sensing performance upon H2 exposure, and the obtained results were analysed in relation to film properties and the microstructural evolution induced by annealing.Item Homogenization of the acoustic streaming in periodic rigid porous structures(2025) Rohan, Eduard; Moravcová, FannyThe paper presents a new model of the acoustic streaming (AS) in rigid porous media. The modelling is based on the classical perturbation approach combined with the periodic homogenization. The first one enables to linearize the Navier–Stokes equations for a barotropic fluid using the decomposition into the first and the second order subproblems governing the fluid dynamics in the rigid period scaffolds. The acoustic wave captured by the first order problem provides the Reynolds stress which appears in the second order problem as the streaming source term. Both the subproblems are treated by the homogenization resulting in the dynamic Darcy flow equations. Using spectral analysis of the characteristic microscopic dynamic Stokes flow and the associated spectral decomposition of the responses, the dynamic permeability is derived and also the driving force for the time-averaged permanent flow is evaluated. The AS can be observed at both the macroscopic and the microscopic levels. While the acoustics-driven microflows are observed for any microstructure, the macroscopic AS depends on the porous microstructure geometry, its nonsymmetry and boundary conditions. For porous particulate structures, the forces and force moments acting on the suspended particles are computed and the influence of the wave frequency and geometrical features is examined. All these effects are illustrated using 2D examples of periodic porous microstructures.Item Consistent determination of the gravimetric geoid and orthometric height(2025) Tenzer, Robert; Novák, PavelVarious computational methods have been developed and applied to determine regional gravimetric geoid models with high accuracy using surface gravity and terrain data, while also often taking into consideration topographic mass density information. Helmert’s orthometric height is, on the other hand, until now solely used for practical realization of vertical geodetic controls in countries where the orthometric height is adopted for the definition of official height systems. Whereas small errors (at the level of a few centimeters) are reported for accurately determined regional gravimetric geoid models, errors in Helmert’s orthometric height reach several centimeters and decimeters already at levelling networks realized in lowlands and regions with moderately elevated topography. In mountainous regions with extremely elevated topography, these errors reach several meters. In Helmert’s definition of the orthometric height, the mean value of gravity within topographic masses is computed approximately from observed surface gravity by applying the Poincaré-Prey gravity gradient reduction, without applying complex computational methods that are used in the gravimetric geoid modelling. This approximation introduces errors due to assuming a constant topographic mass density and disregarding terrain geometry and mass density heterogeneities inside the geoid. Consequently, values of Helmert’s orthometric heights are not consistent with accurately determined regional gravimetric geoid models and should not be fitted or combined with GNSS/levelling data. To address this theoretical inconsistency, we propose a computational scheme based on applying developed methods for consistent determination of the regional gravimetric geoid and orthometric height to achieve their full compatibility by means of improving the accuracy of the orthometric height. We demonstrate that computational methods applied in the regional gravimetric geoid modelling can be modified to determine also the accurate orthometric height, so that both quantities are computed consistently and simultaneously. We also show that the proposed computational scheme can be also used for an accurate conversion of normal to orthometric heights by means of applying the geoid-to-quasigeoid separation. This allows an independent validation of regional gravimetric geoid models.Item Unveiling effects of Zr alloying on structure and properties of nanocrystalline Cu–Zr films(2025) Zhadko, Mariia; Benediktová, Anna; Čerstvý, Radomír; Houška, Jiří; Čapek, Jiří; Kolenatý, David; Minár, Jan; Baroch, Pavel; Zeman, PetrNanocrystalline Cu–Zr films with Zr content in the range of 0.3–2.7 at.% were deposited by direct current magnetron sputter deposition. Effects of Zr alloying on the structure, surface, mechanical, and electrical properties were systematically investigated using X-ray diffraction, electron microscopy, atomic force microscopy, indentation, and the four-point probe method. The experimental results revealed that the Zr content significantly affects the structural and functional characteristics of the films, with the most notable changes observed between 0.3 and ≈1.5 at.% Zr. Beyond this range, further increase in the Zr content results in only minor changes in the microstructure and mechanical properties, while the solubility, electrical resistivity, and surface roughness continue to rise. The alloyed Cu–Zr films exhibit hardness values between 3.2 and 4.2 GPa, exceeding 2.5 GPa measured for the unalloyed Cu film, which is attributed to the combined effect of grain boundary strengthening due to structural refinement and Zr segregation, along with solid solution strengthening.Item Multilayer design of sustainable multifunctional Zr–Cu–N coatings: A route for enhanced mechanical and antibacterial performance(2025) Daniel, Rostislav; Ziegelwanger, Tobias; Zítek, Michal; Červená, Michaela; Haviar, Stanislav; Meindlhumer, Michael; Baroch, Pavel; Keckes, Jozef; Zeman, PetrWear-resistant protective coatings with antimicrobial activity are essential for durability and hygiene in healthcare, public spaces, food industry, consumer products, and industrial environments. This study developed sustainable multifunctional Zr–Cu–N coatings with exceptional damage tolerance, and antibacterial properties using non-reactive and reactive sputtering of only two elemental Zr and Cu targets without external heating. The coatings’ superior performance stems from a sophisticated multilayer architecture combining elastic ZrCu metallic glass, hard and stiff ZrN ceramic, and hard and tough ZrN–Cu nanocomposite coatings. Each constituent was optimized for composition and mechanical properties before integration into multilayer structures to provide high damage tolerance and antibacterial functionality. Antibacterial efficacy was tested in a high-traffic environment over 60 days, showing consistent antimicrobial performance. Fracture stress and toughness were assessed through in situ bending experiments on microcantilever beams fabricated by focused ion beam milling. Results revealed that optimizing the thicknesses of ductile and stiff sublayers significantly enhances damage tolerance while maintaining high hardness and wear resistance. The incorporation of Cu in an unbonded state within the ZrN–Cu nanocomposite facilitates sustainable and scalable production of these multifunctional coatings with antibacterial properties, making them ideal for large surface applications in high-traffic environments like hospitals, office buildings, and public transport.