2025-05-15 | | Total: 15
This work is concerned with constructing a second-order, invariant-domain preserving approximation of the compressible multi-species Euler equations where each species is modeled by an ideal gas equation of state. We give the full solution to the Riemann problem and derive its maximum wave speed. The maximum wave speed is used in constructing a first-order invariant-domain preserving approximation. We then extend the methodology to second-order accuracy and detail a convex limiting technique which is used for preserving the invariant domain. Finally, the numerical method is verified with analytical solutions and then validated with several benchmarks and laboratory experiments.
Surfactants have important effects on the dynamics of droplets on solid surfaces, which has inspired many industrial applications. Phase-field surfactant model with moving contact lines (PFS-MCL) has been employed to investigate the complex droplet dynamics with surfactants, while its numerical simulation remains challenging due to the coupling of gradient flows with respect to transport distances involving nonlinear and degenerate mobilities. We propose a novel structure-preserving variational scheme for PFS-MCL model with the dynamic boundary condition based on the minimizing movement scheme and optimal transport theory for Wasserstein gradient flows. The proposed scheme consists of a series of convex minimization problems and can be efficiently solved by our proposed primal-dual splitting method and its accelerated versions. By respecting the underlying PDE's variational structure with respect to the transport distance, the proposed scheme is proved to inherits the desirable properties including original energy dissipation, bound-preserving, and mass conservation. Through a suite of numerical simulations, we validate the performance of the proposed scheme and investigate the effects of surfactants on the droplet dynamics.
This paper presents a novel quasi-3D theory and the corresponding mixed beam element model to achieve accurate solutions for functionally graded beams. The key innovations include the development of equilibrium-based stress expressions, the modified cross-sectional stiffness matrix, and the mixed beam element model based on semi-analytical definition of internal force fields. In contrast to the conventional quasi-3D theory where stress expressions are derived from constitutive equations and geometric relations, the stress expressions in this study are derived from the differential equilibrium equations among stresses, ensuring strict adherence of stress solutions to equilibrium conditions. To incorporate the influence of equilibrium-derived stress distributions, the modified cross-sectional stiffness matrix is derived, enhancing the theoretical and practical feasibility of the beam model. For beam element construction, the mixed variational principle of two-field variables is employed, with generalized internal forces and generalized displacements regarded as two independent fields. Especially, semi-analytical internal force fields, which partially satisfy the differential equilibrium equations, are introduced to improve the element performance. Numerical examples are conducted to verify the accuracy and effectiveness of the proposed theory and beam element.
In this paper, a vertex-based auxiliary space multigrid(V-ASMG) method as a preconditioner of the PCG method is proposed for solving the large sparse linear equations derived from the linear elasticity equations. The main key of such V-ASMG method lies in an auxiliary region-tree structure based on the geometrically regular subdivision. The computational complexity of building such a region-tree is O(qNlog2N), where N is the number of the given original grid vertices and q is the power of the ratio of the maximum distance dmax to minimum distance dmin between the given original grid vertices. The process of constructing the auxiliary region-tree is similar to the method in [17], but the selection of the representative points is changed. To be more specific, instead of choosing the barycenters, the correspondence between each grid layer is constructed based on the position relationship of the grid vertices. There are two advantages for this approach: the first is its simplicity, there is no need to deal with hanging points when building the auxiliary region-tree, and it is possible to construct the restriction/prolongation operator directly by using the bilinear interpolation function, and it is easy to be generalized to other problems as well, due to all the information we need is only the grid vertices; the second is its strong convergence, the corresponding relative residual can quickly converge to the given tolerance(It is taken to be 10−6 in this paper), thus obtaining the desired numerical solution. Two- and three-dimensional numerical experiments are given to verify the strong convergence of the proposed V-ASMG method as a preconditioner of the PCG method.
In this paper, we propose a federated algorithm for solving large linear systems that is inspired by the classic randomized Kaczmarz algorithm. We provide convergence guarantees of the proposed method, and as a corollary of our analysis, we provide a new proof for the convergence of the classic randomized Kaczmarz method. We demonstrate experimentally the behavior of our method when applied to related problems. For underdetermined systems, we demonstrate that our algorithm can be used for sparse approximation. For inconsistent systems, we demonstrate that our algorithm converges to a horizon of the least squares solution. Finally, we apply our algorithm to real data and show that it is consistent with the selection of Lasso, while still offering the computational advantages of the Kaczmarz framework and thresholding-based algorithms in the federated setting.
The aim of this work is to investigate semi-Lagrangian approximation schemes on unstructured grids for viscous transport and conservative equations with measurable coefficients that satisfy a one-sided Lipschitz condition. To establish the convergence of the schemes, we exploit the characterization of the solution for these equations expressed in terms of measurable time-dependent viscosity solution and, respectively, duality solution. We supplement our theoretical analysis with various numerical examples to illustrate the features of the schemes.
One of the fundamental mathematical models for studying random genetic drift is the Kimura equation, derived as the large-population limit of the discrete Wright-Fisher model. However, due to the degeneracy of the diffusion coefficient, it is impossible to impose a suitable boundary condition that ensures the Kimura equation admits a classical solution while preserving biological significance. In this work, we propose a modified model for random genetic drift that admits classical solutions by modifying the domain of the Kimura equation from (0,1) to (δ,1−δ) with δ being a small parameter, which allows us to impose a Robin-type boundary condition. By introducing two additional variables for the probabilities in the boundary region, we effectively capture the conservation of mass and the fixation dynamics in the original model. To numerically investigate the modified model, we develop a hybrid Eulerian-Lagrangian operator splitting scheme. The scheme first solves the flow map equation in the bulk region using a Lagrangian approach with a no-flux boundary condition, followed by handling the boundary dynamics in Eulerian coordinates. This hybrid scheme ensures mass conservation, maintains positivity, and preserves the first moment. Various numerical tests demonstrate the efficiency, accuracy, and structure-preserving properties of the proposed scheme. Numerical results demonstrate the key qualitative features of the original Kimura equation, including the fixation behavior and the correct stationary distribution in the small-δ limit.
We prove convergence and stability of the discrete exterior calculus (DEC) solutions for the Hodge-Laplace problems in two dimensions for families of meshes that are non-degenerate Delaunay and shape regular. We do this by relating the DEC solutions to the lowest order finite element exterior calculus (FEEC) solutions. A Poincaré inequality and a discrete inf-sup condition for DEC are part of this proof. We also prove that under appropriate geometric conditions on the mesh the DEC and FEEC norms are equivalent. Only one side of the norm equivalence is needed for proving stability and convergence and this allows us to relax the conditions on the meshes.
Waveform inversion seeks to estimate an inaccessible heterogeneous medium from data gathered by sensors that emit probing signals and measure the generated waves. It is an inverse problem for a second order wave equation or a first order hyperbolic system, with the sensor excitation modeled as a forcing term and the heterogeneous medium described by unknown, spatially variable coefficients. The traditional ``full waveform inversion" (FWI) formulation estimates the unknown coefficients via minimization of the nonlinear, least squares data fitting objective function. For typical band-limited and high frequency data, this objective function has spurious local minima near and far from the true coefficients. Thus, FWI implemented with gradient based optimization algorithms may fail, even for good initial guesses. Recently, it was shown that it is possible to obtain a better behaved objective function for wave speed estimation, using data driven reduced order models (ROMs) that capture the propagation of pressure waves, governed by the classic second order wave equation. Here we introduce ROMs for vectorial waves, satisfying a general first order hyperbolic system. They are defined via Galerkin projection on the space spanned by the wave snapshots, evaluated on a uniform time grid with appropriately chosen time step. Our ROMs are data driven: They are computed in an efficient and non-iterative manner, from the sensor measurements, without knowledge of the medium and the snapshots. The ROM computation applies to any linear waves in lossless and non-dispersive media. For the inverse problem we focus attention on acoustic waves in a medium with unknown variable wave speed and density. We show that these can be determined via minimization of an objective function that uses a ROM based approximation of the vectorial wave field inside the inaccessible medium.
We provide a framework for interpreting Discrete Exterior Calculus (DEC) numerical schemes in terms of Finite Element Exterior Calculus (FEEC). We demonstrate the equivalence of cochains on primal and dual meshes with Whitney and generalized Whitney forms which allows us to analyze DEC approximations using tools from FEEC. We demonstrate the applicability of our framework by rigorously proving convergence with rates for the Hodge-Laplacian problem in full k-form generality on well-centered meshes on contractible domains. We also provide numerical results illustrating optimality of our derived convergence rates. Moreover, we demonstrate how superconvergence phenomena can be explained in our framework with corresponding numerical results.
A Jacobian free Newton Krylov (JFNK) method with a globalization scheme is introduced to solve large and complex nonlinear systems of equations that arise in groundwater flow models of multi-layer aquifer systems. We explore the advantages of the JFNK method relative to the Newton-Krylov (NK) method and identify the circumstances in which the JFNK method demonstrates computing efficiency. We perform the validation and efficiency of the JFNK method on various test cases involving an unconfined single-layer aquifer and a two-layer aquifer with both confined and unconfined conditions. The results are validated by the NK method. The JFNK method is incorporated in Integrated Water Flow Model (IWFM), an integrated hydrologic model developed and maintained by California Department of Water Resources. We examine the determinacy of the JFNK's adaptability on practical models such as the California Central Valley Groundwater-Surface Water Simulation Model (C2VSim).
The recently proposed affine frequency division multiplexing (AFDM) is a new transmission waveform that has shown excellent performance in high-mobility environments, making it a sensible option for the next-generation wireless networks. In this paper, we investigate an energy-efficient index modulation scheme for AFDM by leveraging spread spectrum, referred to as IM-AFDM-SS, to combat the interference caused by the doubly dispersive channels. Specifically, the information bits are conveyed by the transmitted symbols as well as the indices of the selected spreading codes in our proposed IM-AFDM-SS scheme. To avoid extensive computations, we also develop a lowcomplexity maximal ratio combining (MRC) detector algorithm, which recovers the spreading codes first and demodulates the symbols afterwards. Moreover, an upper bound on the bit error rate (BER) of the proposed IM-AFDM-SS system with maximumlikelihood (ML) detection is derived. Numerical results demonstrate the superiority of the proposed IM-AFDM-SS system over the classical AFDM spread spectrum (AFDM-SS) and the existing index modulated AFDM (IM-AFDM) systems.
This study proves that the injectivity condition for the integral operator of the direct-indirect mixed Burton-Miller (BM) boundary integral equation (BIE) for Helmholtz transmission problems is identical to that for the ordinary BM BIE for Helmholtz transmission problems with a transmissive circular inclusion. Although some numerical methods based on the direct-indirect mixed BM BIE can be computed faster than the ordinary BM BIE, its well-posedness has been unclear. This study resolves a part of the well-posedness, namely the injectivity of the integral operator with a transmissive circular inclusion.
Aerodynamic design optimization is an important problem in aircraft design that depends on the interplay between a numerical optimizer and a high-fidelity flow physics solver. Derivative-based, first and (quasi) second order, optimization techniques are the de facto choice, particularly given the availability of the adjoint method and its ability to efficiently compute gradients at the cost of just one solution of the forward problem. However, implementation of the adjoint method requires careful mathematical treatment, and its sensitivity to changes in mesh quality limits widespread applicability. Derivative-free approaches are often overlooked for large scale optimization, citing their lack of scalability in higher dimensions and/or the lack of practical interest in globally optimal solutions that they often target. However, breaking free from an adjoint solver can be paradigm-shifting in broadening the applicability of aerodynamic design optimization. We provide a systematic benchmarking of a select sample of widely used derivative-based and derivative-free optimization algorithms on the design optimization of three canonical aerodynamic bodies, namely, the NACA0012 and RAE2822 airfoils, and the ONERAM6 wing. Our results demonstrate that derivative-free methods are competitive with derivative-based methods, while outperforming them consistently in the high-dimensional setting. These findings highlight the practical competitiveness of modern derivative-free strategies, offering a scalable and robust alternative for aerodynamic design optimization when adjoint-based gradients are unavailable or unreliable.
We introduce Flexi-VAE, a data-driven framework for efficient single-shot forecasting of nonlinear parametric partial differential equations (PDEs), eliminating the need for iterative time-stepping while maintaining high accuracy and stability. Flexi-VAE incorporates a neural propagator that advances latent representations forward in time, aligning latent evolution with physical state reconstruction in a variational autoencoder setting. We evaluate two propagation strategies, the Direct Concatenation Propagator (DCP) and the Positional Encoding Propagator (PEP), and demonstrate, through representation-theoretic analysis, that DCP offers superior long-term generalization by fostering disentangled and physically meaningful latent spaces. Geometric diagnostics, including Jacobian spectral analysis, reveal that propagated latent states reside in regions of lower decoder sensitivity and more stable local geometry than those derived via direct encoding, enhancing robustness for long-horizon predictions. We validate Flexi-VAE on canonical PDE benchmarks, the 1D viscous Burgers equation and the 2D advection-diffusion equation, achieving accurate forecasts across wide parametric ranges. The model delivers over 50x CPU and 90x GPU speedups compared to autoencoder-LSTM baselines for large temporal shifts. These results position Flexi-VAE as a scalable and interpretable surrogate modeling tool for accelerating high-fidelity simulations in computational fluid dynamics (CFD) and other parametric PDE-driven applications, with extensibility to higher-dimensional and more complex systems.