@misc{MTMT:37444899, title = {Exact analytical solution for the non-selfadjoint problem of Maxwell–Cattaneo–Vernotte heat conduction with heat-transfer boundary condition}, url = {https://m2.mtmt.hu/api/publication/37444899}, author = {Szücs, Mátyás and Fülöp, Tamás}, unique-id = {37444899}, year = {2026}, orcid-numbers = {Szücs, Mátyás/0000-0002-2492-0392; Fülöp, Tamás/0000-0003-2708-7065} } @article{MTMT:36074073, title = {GENERIC-motivated extended symplectic numerical methods for dissipative mechanical systems: A concrete example and general messages}, url = {https://m2.mtmt.hu/api/publication/36074073}, author = {Fülöp, Tamás and Szücs, Mátyás and Takács, Donát M.}, doi = {10.32973/jcam.2025.001}, journal-iso = {JCAM}, journal = {JOURNAL OF COMPUTATIONAL AND APPLIED MECHANICS}, volume = {20}, unique-id = {36074073}, issn = {1586-2070}, abstract = {The GENERIC (General Equation for the Non-Equilibrium Reversible-Irreversible Coupling) framework is used, first, to present the irreversible thermodynamical formulation of the Rayleigh potential, and, subsequently, to realize the dynamical equations of motion of a supercritical van der Waals fluid separated from the environment by a movable rigid piston. Motivated by the latter, a quasi-symplectic generalization of the symplectic Euler finite-difference numerical scheme is introduced for this setup. A remarkable advantage over the explicit Euler scheme regarding artificial numerical antidissipation is illustrated. The possibility of controlled artificial numerical damping is addressed. The quasi-symplectic scheme is proven to enable efficient simulation of the considerably nonlinear and sensitive processes near the liquid-vapour critical point. A numerically feasible quantitative measure of nonlinearity of time-dependent processes is introduced and applied.}, year = {2025}, eissn = {2732-0189}, pages = {3-33}, orcid-numbers = {Fülöp, Tamás/0000-0003-2708-7065; Szücs, Mátyás/0000-0002-2492-0392; Takács, Donát M./0000-0002-8463-745X} } @article{MTMT:37071195, title = {The Kluitenberg–Verhás model of viscoelasticity and its irreversible thermodynamical consistency via rheological energy}, url = {https://m2.mtmt.hu/api/publication/37071195}, author = {Fülöp, Tamás}, doi = {10.32973/jcam.2025.004}, journal-iso = {JCAM}, journal = {JOURNAL OF COMPUTATIONAL AND APPLIED MECHANICS}, volume = {20}, unique-id = {37071195}, issn = {1586-2070}, abstract = {Rocks and plastics are two examples of classes of materials where theKluitenberg–Verh´as model of viscoelasticity has proved adequate for describing the observedmechanical behaviour. This model cannot be interpreted as a generalized Maxwell model,and requires in fact a novel rheological circuit element beyond the spring and the dashpot.Moreover, its thermodynamical consistency leads to nontrivial inequality conditions on themodel coefficients, ignoring which would lead to violation of the second law of thermodynamics and would thus be a principal source of instability of numerical solutions. The presentpaper explores the role of the rheological energy contribution embedded in the Kluitenberg–Verh´as model, in a didactic and application-oriented way. The explicit meaning of theinvolved irreversible thermodynamical additional state variable – internal degree of freedom– is identified. Neither thermal expansion nor heat conduction is neglected, enhancing theapplicability of the constructed framework.}, year = {2025}, eissn = {2732-0189}, pages = {93-117}, orcid-numbers = {Fülöp, Tamás/0000-0003-2708-7065} } @article{MTMT:34726196, title = {Improving the accuracy of the Newmark method through backward error analysis}, url = {https://m2.mtmt.hu/api/publication/34726196}, author = {Takács, Donát M. and Fülöp, Tamás}, doi = {10.1007/s00466-024-02580-3}, journal-iso = {COMPUT MECH}, journal = {COMPUTATIONAL MECHANICS}, volume = {75}, unique-id = {34726196}, issn = {0178-7675}, abstract = {We use backward error analysis for differential equations to obtain modified or distorted equations describing the behaviour of the Newmark scheme applied to the transient structural dynamics equation. Based on the newly derived distorted equations, we give expressions for the numerically or algorithmically distorted stiffness and damping matrices of a system simulated using the Newmark scheme. Using these results, we show how to construct compensation terms from the original parameters of the system, which improve the performance of Newmark simulations. The required compensation terms turn out to be slight modifications to the original system parameters (e.g. the damping or stiffness matrices), and can be applied without changing the time step or modifying the scheme itself. Two such compensations are given: one eliminates numerical damping, while the other achieves fourth-order accurate calculations using the traditionally second-order Newmark method. The performance of both compensation methods is evaluated numerically to demonstrate their validity, and they are compared to the uncompensated Newmark method, the generalized-α method and the 4th-order Runge–Kutta scheme. © The Author(s) 2024.}, year = {2025}, eissn = {1432-0924}, pages = {1585-1606}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065} } @article{MTMT:36161803, title = {The piston effect in supercritical fluids investigated via a reversible–irreversible vector field splitting-based explicit time integration scheme}, url = {https://m2.mtmt.hu/api/publication/36161803}, author = {Takács, Donát M. and Fülöp, Tamás and Kovács, Róbert Sándor and Szücs, Mátyás}, doi = {10.1063/5.0271333}, journal-iso = {PHYS FLUIDS}, journal = {PHYSICS OF FLUIDS}, volume = {37}, unique-id = {36161803}, issn = {1070-6631}, abstract = {In the vicinity of the liquid–vapor critical point, supercritical fluids behave strongly compressible and, in parallel, thermophysical properties have strong state dependence. These lead to various peculiar phenomena, one of which is the piston effect where a sudden heating induces a mechanical pulse. The coupling between thermal and mechanical processes, in the linear approximation, yields a non-trivially rich thermoacoustics. The numerous applications of supercritical fluids raise the need for a reliable yet fast and efficient numerical solution for thermoacoustic time and space dependence in this sensitive domain. Here, we present a second-order accurate, fully explicit staggered space–time grid finite difference method for such coupled linear thermoacoustic problems. Time integration is based on the splitting of the state space vector field representing the interactions that affect the dynamics into reversible and irreversible parts, which splitting procedure leads to decoupled wave and heat equations. The former is a hyperbolic partial differential equation, while the latter is a parabolic one; therefore, different time integration algorithms must be amalgamated to obtain a reliable, dispersion error-free, and dissipation error-free numerical solution. Finally, the thermoacoustic approximation of the supercritical piston effect is investigated via the developed method.}, year = {2025}, eissn = {1089-7666}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065; Kovács, Róbert Sándor/0000-0001-5822-6035; Szücs, Mátyás/0000-0002-2492-0392} } @article{MTMT:36178057, title = {Improving Discrete Numerical Methods for Dynamics Using Continuous Mathematical Tools}, url = {https://m2.mtmt.hu/api/publication/36178057}, author = {Takács, Donát M. and Fülöp, Tamás}, doi = {10.4028/p-11TtO6}, journal-iso = {ADV SCI TECHNOL}, journal = {ADVANCES IN SCIENCE AND TECHNOLOGY}, volume = {165}, unique-id = {36178057}, issn = {1662-8969}, abstract = {We briefly review backward error analysis as a useful mathematical technique for improving numerical methods used for solving ordinary differential equations describing dynamical systems. Then, we show how backward error analysis-based compensation, an approach recently introduced by the authors, can be applied to the second-order Newmark method for eliminating numerical damping and achieving fourth-order convergence. The presented improvements only require modifying the physical parameters of the system and the excitation, while the Newmark method is left intact. We compare the performance of the resulting improvements with that of several other numerical methods, including a novel partitioned method based on local extrapolation.}, year = {2025}, eissn = {1662-0356}, pages = {21-30}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065} } @article{MTMT:36214178, title = {On the coordinate system-dependence of the accuracy of symplectic numerical methods}, url = {https://m2.mtmt.hu/api/publication/36214178}, author = {Takács, Donát M. and Fülöp, Tamás}, doi = {10.33993/jnaat542-1577}, journal-iso = {J NUMER ANA. APPROX THEORY}, journal = {JOURNAL OF NUMERICAL ANALYSIS AND APPROXIMATION THEORY}, volume = {54}, unique-id = {36214178}, issn = {2457-6794}, abstract = {Symplectic numerical methods have become a widely-used choice for the accurate simulation of Hamiltonian systems in various fields, including celestial mechanics, molecular dynamics and robotics. Even though their characteristics are well-understood mathematically, relatively little attention has been paid in general to the practical aspect of how the choice of coordinates affects the accuracy of the numerical results, even though the consequences can be computationally significant. The present article aims to fill this gap by giving a systematic overview of how coordinate transformations can influence the results of simulations performed using symplectic methods. We give a derivation for the non-invariance of the modified Hamiltonian of symplectic methods under coordinate transformations, as well as a sufficient condition for the non-preservation of a first integral corresponding to a cyclic coordinate for the symplectic Euler method. We also consider the possibility of finding order-compensating coordinate transformations that improve the order of accuracy of a numerical method. Various numerical examples are presented throughout.}, year = {2025}, eissn = {2501-059X}, pages = {315-344}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065} } @{MTMT:36227399, title = {Improving Discrete Numerical Methods for Dynamics Using Continuous Mathematical Tools}, url = {https://m2.mtmt.hu/api/publication/36227399}, author = {Takács, Donát M. and Fülöp, Tamás}, booktitle = {The 10th International Scientific Conference on Advances in Mechanical Engineering (ISCAME)}, unique-id = {36227399}, abstract = {We briefly review backward error analysis as a useful mathematical technique for improving numerical methods used for solving ordinary differential equations describing dynamical systems. Then, we show how backward error analysis-based compensation, an approach recently introduced by the authors, can be applied to the second-order Newmark method for eliminating numerical damping and achieving fourth-order convergence. The presented improvements only require modifying the physical parameters of the system and the excitation, while the Newmark method is left intact. We compare the performance of the resulting improvements with that of several other numerical methods, including a novel partitioned method based on local extrapolation.}, year = {2025}, pages = {235-244}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065} } @article{MTMT:34259085, title = {Szimplektikus numerikus módszerek és műszaki alkalmazásuk néhány gyakorlati szempontja [Symplectic numerical methods and some practical aspects of their application for engineering problems]}, url = {https://m2.mtmt.hu/api/publication/34259085}, author = {Takács, Donát M. and Fülöp, Tamás}, doi = {10.35925/j.multi.2024.4.1}, journal-iso = {MULTIDISZCIPLINÁRIS TUDOMÁNYOK}, journal = {MULTIDISZCIPLINÁRIS TUDOMÁNYOK: A MISKOLCI EGYETEM KÖZLEMÉNYE}, volume = {14}, unique-id = {34259085}, issn = {2062-9737}, year = {2024}, eissn = {2786-1465}, pages = {3-23}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065} } @article{MTMT:34575282, title = {Thermodynamically extended symplectic numerical simulation of viscoelastic, thermal expansion and heat conduction phenomena in solids}, url = {https://m2.mtmt.hu/api/publication/34575282}, author = {Takács, Donát M. and Pozsár, Áron and Fülöp, Tamás}, doi = {10.1007/s00161-024-01280-w}, journal-iso = {CONTINUUM MECH THERM}, journal = {CONTINUUM MECHANICS AND THERMODYNAMICS}, volume = {36}, unique-id = {34575282}, issn = {0935-1175}, abstract = {Symplectic numerical schemes for reversible dynamical systems predict the solution reliably over large times as well, and are a good starting point for extension to schemes for simulating irreversible situations like viscoelastic wave propagation and heat conduction coupled via thermal expansion occuring in rocks, plastics, biological samples etc. Dissipation error (artificial nonpreservation of energies and amplitudes) of the numerical solution should be as small as possible since it should not be confused with the real dissipation occurring in the irreversible system. In addition, the other well-known numerical artefact, dispersion error (artificial oscillations emerging at sharp changes), should also be minimal to avoid confusion with the true wavy behavior. The continuum thermodynamical aspects (respect for balances with fluxes, systematic constitutive relationships between intensive quantities and fluxes, the second law of thermodynamics with positive definite entropy production, and the spacetime-based kinematic viewpoint) prove valuable for obtaining such extended schemes and for monitoring the solutions. Generalizing earlier works in this direction, here, we establish and investigate such a numerical scheme for one-dimensional viscoelastic wave propagation in the presence of heat conduction coupled via thermal expansion, demonstrating long-term reliability and the applicability of thermodynamics-based quantities in supervising the quality of the solution.}, year = {2024}, eissn = {1432-0959}, pages = {525-538}, orcid-numbers = {Takács, Donát M./0000-0002-8463-745X; Fülöp, Tamás/0000-0003-2708-7065} }