Homogeneous Boundary Conditions Research Articles

Construction of Short-Time Heat Conduction Solutions in One-Dimensional Finite Rectangular Bodies

Construction of Short-Time Heat Conduction Solutions in One-Dimensional Finite Rectangular Bodies

Optimal distributed control for a Cahn-Hilliard-Darcy system with mass sources, unmatched viscosities and singular potential

We study a Cahn-Hilliard-Darcy system with mass sources, which can be considered as a basic, though simplified, diffuse interface model for the evolution of tumor growth. This system is equipped with an impermeability condition for the averaged velocity u as well as homogeneous Neumann boundary conditions for the phase function ϕ and the chemical potential µ. The source term in the convective Cahn-Hilliard equation contains a control R that can be thought as a drug or a nutrient in applications. Our goal is to study a distributed optimal control problem in the two-dimensional setting with a cost functional of tracking-type. In the physically relevant case with unmatched viscosities for the binary fluid mixtures and a singular potential, we first prove the existence and uniqueness of a global strong solution with ϕ being strictly separated from the pure phases ±1. This well-posedness result enables us to characterize the control-to-state mapping S: R→ϕ. Then we obtain the existence of an optimal control, the Frechet differentiability of S and first-order necessary optimality conditions expressed through a suitable variational inequality for the adjoint variables. Finally, we prove the differentiability of the control-to-costate operator and establish a second-order sufficient condition for the strict local optimality.

In-flight anti-icing simulation of electrothermal ice protection systems with inhomogeneous thermal boundary condition

Conventional anti-icing computational solvers calculate the convective heat transfer coefficient using a homogeneous thermal boundary condition, assuming a constant temperature across the surface. However, this approach can lead to inaccuracies in regions with significant temperature variations. To address this limitation, the present study proposes a novel approach that utilizes an inhomogeneous thermal boundary condition, which updates the convective heat transfer coefficient based on the transient wall temperature distribution. We developed a unified finite volume framework that efficiently integrates various in-house solvers, including a compressible Navier-Stokes-Fourier (NSF) airflow solver, a Eulerian droplet impingement solver, a PDE-based ice solver, and a multilayer heat conduction solver. Thermal interaction between the different solvers was modeled using the conjugate heat transfer method. The effects of updating airflow on anti-icing results were investigated by comparing the results obtained by coupled and decoupled solvers. While the decoupled solver computes airflow once and remains unchanged, the coupled solver updates the airflow during the anti-icing simulation. Our findings show that the decoupled solver has the highest deviations from the coupled solver in evaporative anti-icing regimes, where dry regions form within the protection limits. Using coupled solvers in designing ice protection systems in evaporative regimes improves temperature prediction accuracy, enabling designers to reduce safety factors and save energy.

Asymptotic Behavior of Global Solution of an Anomalous Gray Scott Model

The object of this paper is to prove that existence global and asymptotic behavior of solutions for anomalous coupled reaction diffusion system (Gray Scott model) with homogeneous Neumann boundary conditions. The existence and uniqueness of the local solution are given by the Banach fixed point theorem. However, the asymptotic behavior is investiged by technique semi group estimates and the Sobolev embedding theorem

Existence of solutions for anisotropic parabolic Ni-Serrin type equations originated from a capillary phenomena with nonstandard growth nonlinearity

We consider an initial boundary value problem for a class of anisotropic parabolic Ni-Serrin type equations with nonstandard nonlinearity in a bounded smooth domain with homogeneous Dirichlet boundary condition. Because the nonlinear perturbation leads to difficulties (it does not have a definite sign) in obtaining a priori estimates in the energy method, we had to modify the Tartar method significantly. Under suitable assumptions, we obtain the global existence, decay and extinction of solutions.

Analytic and Gevrey class regularity for parametric semilinear reaction-diffusion problems and applications in uncertainty quantification

We investigate a class of parametric elliptic semilinear partial differential equations of second order with homogeneous essential boundary conditions, where the coefficients and the right-hand side (and hence the solution) may depend on a parameter. This model can be seen as a reaction-diffusion problem with a polynomial nonlinearity in the reaction term. The efficiency of various numerical approximations across the entire parameter space is closely related to the regularity of the solution with respect to the parameter. We show that if the coefficients and the right-hand side are analytic or Gevrey class regular with respect to the parameter, the same type of parametric regularity is valid for the solution. The key ingredient of the proof is the combination of the alternative-to-factorial technique from our previous work [1] with a novel argument for the treatment of the power-type nonlinearity in the reaction term. As an application of this abstract result, we obtain rigorous convergence estimates for numerical integration of semilinear reaction-diffusion problems with random coefficients using Gaussian and Quasi-Monte Carlo quadrature. Our theoretical findings are confirmed in numerical experiments.

Solutions to a coupled Cahn–Hilliard-phase-field-crystal model for grain boundary motion

In this paper, we consider an initial–boundary value problem (IBVP) of a coupled Cahn–Hilliard-phase-field-crystal (CH-PFC) system subject to homogeneous Neumann boundary conditions in one spatial dimension. This CH-PFC model couples the composition field of a diffusing species with the crystallographic and can be used to model the diffusion-induced grain boundary migration in crystalline materials. Under suitable assumptions on the coefficients and initial data, we prove that the IBVP possesses a global weak solution. Our existence proof, which contributes to the verification of the model, is only valid in one space dimension.

Global weak solvability in a self-consistent chemotaxis-Navier-Stokes system involving Dirichlet boundary conditions for the signal

Global weak solvability in a self-consistent chemotaxis-Navier-Stokes system involving Dirichlet boundary conditions for the signal

The Conservative and Efficient Numerical Method of 2-D and 3-D Fractional Nonlinear Schrödinger Equation Using Fast Cosine Transform

This paper introduces a novel approach employing the fast cosine transform to tackle the 2-D and 3-D fractional nonlinear Schrödinger equation (fNLSE). The fractional Laplace operator under homogeneous Neumann boundary conditions is first defined through spectral decomposition. The difference matrix Laplace operator is developed by the second-order central finite difference method. Then, we diagonalize the difference matrix based on the properties of Kronecker products. The time discretization employs the Crank–Nicolson method. The conservation of mass and energy is proved for the fully discrete scheme. The advantage of this method is the implementation of the Fast Discrete Cosine Transform (FDCT), which significantly improves computational efficiency. Finally, the accuracy and effectiveness of the method are verified through two-dimensional and three-dimensional numerical experiments, solitons in different dimensions are simulated, and the influence of fractional order on soliton evolution is obtained; that is, the smaller the alpha, the lower the soliton evolution.

Existence and stability of steady-state bifurcations in a prey–predator system with a nonmonotonic functional response

A cross-diffusion prey–predator system exhibiting the prey group defense under homogeneous Neumann boundary conditions is studied. By considering the diffusive rate of the prey as a bifurcation parameter, we investigate sudden changes in the population dynamics of the prey and predator which can have a substantial effect on population size of the species. First a priori estimate for positive steady states is obtained. Next we prove the existence of a pitchfork bifurcation of positive steady states at a simple eigenvalue. The structure of the global steady-state bifurcation is discussed. We also investigate the stability of the trivial solution line and nontrivial steady-state solutions via the eigenvalue perturbation theory. To illustrate our theoretical results some numerical simulations are given. Numerical examples contain a supercritical and a subcritical pitchfork bifurcation.

Integral Transform Solution for Natural Convection Within Horizontal Concentric Annular Cavities

Abstract Natural convection inside horizontal concentric annular cavities is dealt with through the generalized integral transform technique (GITT), offering a hybrid numerical-analytical solution of the continuity, Navier–Stokes, and energy equations in cylindrical coordinates. The flow is in steady-state, laminar regime, two-dimensional, buoyancy-induced, and the governing equations are written in the streamfunction-only formulation. Two strategies of integral transformation are adopted to verify the best computational performance, namely, the usual one with eigenvalue problems for both streamfunction and temperature defined in the radial variable, and a novel alternative with eigenvalue problems defined in the azimuthal angular coordinate. First, the eigenfunction expansions convergence behavior is analyzed to critically compare the two integral transform solution strategies. Then, test cases for different aspect ratios and Rayleigh numbers are validated with experimental data from the classical work of Kuehn and Goldstein. A maximum relative deviation of 5% is found comparing the GITT results for the average Nusselt number against the experimental data, while an 8% maximum relative deviation is found comparing against an empirical correlation by the same authors. It is concluded that the GITT solution with the eigenvalue problem in the angular coordinate yields better convergence rates than the more usual eigenfunction expansion in the radial variable. This is due to the originally homogeneous boundary conditions in the angular direction, which do not require filtering for convergence enhancement, as opposed to the required filter in the radial direction that introduces a source term in the filtered equation for the streamfunction.

Solving Boundary Value Problems by Sinc Method and Geometric Sinc Method

This paper introduces an efficient numerical method for approximating solutions to geometric boundary value problems. We propose the multiplicative sinc–Galerkin method, tailored specifically for solving multiplicative differential equations. The method utilizes the geometric Whittaker cardinal function to approximate functions and their geometric derivatives. By reducing the geometric differential equation to a system of algebraic equations, we achieve computational efficiency. The method not only proves to be computationally efficient but also showcases a valuable symmetric property, aligning with inherent patterns in geometric structures. This symmetry enhances the method’s compatibility with the often-present symmetries in geometric boundary value problems, offering both computational advantages and a deeper understanding of geometric calculus. To demonstrate the reliability and efficiency of the proposed method, we present several examples with both homogeneous and non-homogeneous boundary conditions. These examples serve to validate the method’s performance in practice.

Steady States of a Diffusive Population-Toxicant Model with Negative Toxicant-Taxis

This paper is dedicated to studying the steady state problem of a population-toxicant model with negative toxicant-taxis, subject to homogeneous Neumann boundary conditions. The model captures the phenomenon in which the population migrates away from regions with high toxicant density towards areas with lower toxicant concentration. This paper establishes sufficient conditions for the non-existence and existence of non-constant positive steady state solutions. The results indicate that in the case of a small toxicant input rate, a strong toxicant-taxis mechanism promotes population persistence and engenders spatially heterogeneous coexistence (see, Theorem 2.3). Moreover, when the toxicant input rate is relatively high, the results unequivocally demonstrate that the combination of a strong toxicant-taxis mechanism and a high natural growth rate of the population fosters population persistence, which is also characterized by spatial heterogeneity (see, Theorem 2.4).

Solvability of an Inverse Problem for an Elliptic-Type Equation

In this study, we consider an inverse problem of determining an unknown source function in the right-hand side of an elliptic equation which is ill-posed in the Hadamard sense. To investigate the solvability of the problem, we reduce it to a Dirichlet problem for a third-order partial differential equation with homogeneous boundary condition. Since the problem is linear, the proof of the uniqueness theorem is based on the Fredholm Alternative Theorem. We prove the existence of the solution to the problem by using the Galerkin method.

The engineering elastic constants of bio-inspired sandwich cores with wavy cylinders

Beetle Elytron Plates (BEPs) represent a new class of biomimetic sandwich cores with excellent mechanical properties and inspired by the internal structure of the beetle elytra. The cores have a hexagonal centre-symmetric configuration with through-thickness cylinders in the unit cell. The present work describes the behaviour of the engineering elastic constants of an evolution of these BEP topologies that feature cylinders with waviness. These BEP core configurations are simulated using Finite Element models, with periodic boundary conditions for volumetric homogenization. The models are then used to perform a parametric analysis against the geometry characteristics of the new cellular configurations. This work further explores the potential design space of BEP cellular metamaterials by altering the geometry of the trabeculae and facilitates further investigation into the mechanical properties of these new BEP topologies with wave-shaped units. Although – and as expected - the specific nondimensional out-of-plane tensile modulus has values here close to unity, BEP topologies with out-of-plane wave-shaped cylinders show enhanced in-plane moduli values than classical straight cylinder configurations, except for the Young's modulus of MBEP cores with only one semi-wave. The same trend is observed for the nondimensional Young's modulus of EBEP configurations, whereas the situation is opposite for the MBEP topologies. The nondimensional out-of-plane shear modulus decreases for both EBEP and MBEP configurations compared to the straight cylinder cores.

On a class of doubly nonlinear evolution equations in Musielak–Orlicz spaces

AbstractThis paper is concerned with a parabolic evolution equation of the form , settled in a smooth bounded domain of , , and complemented with the initial conditions and with (for simplicity) homogeneous Dirichlet boundary conditions. Here, stands for a diffusion operator, possibly nonlinear, which may range in a very wide class, including the Laplacian, the ‐Laplacian for suitable , the “variable‐exponent” ‐Laplacian, or even some fractional order operators. The operator is assumed to be in the form with being measurable in and maximal monotone in . The main results are devoted to proving existence of weak solutions for a wide class of functions that extends the setting considered in previous results related to the variable exponent case where . To this end, a theory of subdifferential operators will be established in Musielak–Orlicz spaces satisfying structure conditions of the so‐called ‐type, and a framework for approximating maximal monotone operators acting in that class of spaces will also be developed. Such a theory is then applied to provide an existence result for a specific equation, but it may have an independent interest in itself. Finally, the existence result is illustrated by presenting a number of specific equations (and, correspondingly, of operators , ) to which the result can be applied.

On the space-time analyticity of the inhomogeneous heat equation on the half space with Neumann boundary conditions

We consider the inhomogeneous heat equation on the half-space R + d with Neumann boundary conditions. We prove a space-time Gevrey regularity of the solution, with a radius of analyticity uniform up to the boundary of the half-space. We also address the case of homogeneous Robin boundary conditions. Our results generalize the case of homogeneous Dirichlet boundary conditions established by Kukavica and Vicol [A direct approach to Gevrey regularity on the half-space. In: Partial differential equations in fluid mechanics. Cambridge: Cambridge Univ. Press; 2018. p. 268–288. (London math. soc. lecture note ser.; vol. 452).].

SH waves in a weakly inhomogeneous half space with a nonlinear thin layer coating

We investigate the self-modulation of Love waves propagating in a nonlinear half-space covered by a nonlinear layer. We assume that the constituent material of the layer is nonlinear, homogeneous, isotropic, compressible, and hyperelastic, whereas for the half-space, it is nonlinear, heterogeneous, compressible and a different hyperelastic material. By employing the nonlinear thin layer approximation, the problem of wave propagation in a layered half-space is reduced to the one for a nonlinear heterogeneous half-space with a modified nonlinear homogeneous boundary condition on the top surface. This new problem is analyzed by a relevant perturbation method, and a nonlinear Schrödinger (NLS) equation defining the self-modulation of waves asymptotically is obtained. The dispersion relation is derived for different heterogeneous properties of the half-space and the thin layer. Then the results of the thin layer approximation are compared with the ones for the finite layer obtained in Teymur et al. (Int J Eng Sci 85:150–162, 2014). The solitary solutions of the derived NLS equation are obtained for selected real material models. It has been discussed how these solutions are influenced by the heterogeneity of the semi-infinite space.

Stability of Solutions to a Caginalp Phase-Field Type Equations

This paper is concerned with the study of the asymptotic behavior of a generalization of the Caginalp phase-field model subject to homogeneous Neumann boundary conditions and regular potentials involving two temperatures. This work follows on from a paper in which the well-posedness of the problem, the dissipativity of the system, and the existence of global and exponential attractors were demonstrated. In addition, a study on the semi-infinite cylinder was also carried out. Indeed, if it is true that the existence of a global attractor makes it possible to predict the asymptotic behavior of solutions on a bounded domain, it does not say that these solutions converge. After having shown the existence of the global attractor, it is therefore important to look at the convergence of the solutions over time. There are several methods for determining the asymptotic behavior of the solutions of a differential system. We can mention the one that consists of transforming the given differential equations into integral equations and then applying the classical Picard successive approximation procedure to them. This work is devoted to the study of the convergence of solutions to steady states, adapting a well-known result concerning Lojasiewicz-Simon’s inequality.

Global boundedness and asymptotic stabilization in a chemotaxis system with density-suppressed motility and nonlinear signal production

In this paper, we study the following chemotaxis model with density-suppressed motility and nonlinear production{ut=Δ(φ(v)u)+ru−μuα,x∈Ω,t>0,vt=Δv−v+wβ,x∈Ω,t>0,wt=Δw−w+uγ,x∈Ω,t>0 under homogeneous Neumann boundary conditions in a bounded domain Ω⊂Rn(n≥2) with smooth boundary, where r∈R,μ,β,γ>0 and α>1. The positive motility function satisfies φ(s)∈C3([0,∞)) and φ′(s)≤0 for all s≥0. It is showed that the system admits a globally bounded and classical solution under some conditions on α,β and γ. Then, under stricter constraints on φ, we obtained that the parameter β has a wider range than before, which is enough to ensure the global boundedness of the solution. Furthermore, if μ is sufficiently large, we proved that the solution converges to ((r+μ)1α−1,(r+μ)βγα−1,(r+μ)γα−1) in L∞(Ω) as t→∞ for all r∈R,μ,β,γ>0 and α>1, where r+:=max⁡{r,0}. Finally, we showed that the convergence is exponential in cases r>0 and r<0.