This paper presents a novel shooting algorithm for solving two-point boundary value problems (BVPs) for differential equations of one independent variable. This algorithm includes the following steps: First, a shooting step is performed, in which a guess for the initial condition is made and a forward numerical integration of the differential equation is performed so that an Initial Value Probl...

We consider a one-phase Stefan problem for the heat equation with a superlinear reaction term. It is known from a previous work (Ghidouche, Souplet, & Tarzia [5]) that all global solutions are bounded and decay uniformly to 0. Moreover, it was shown in Ghidouche, Souplet, & Tarzia [5] that either: (i) the free boundary converges to a finite limit and the solution decays at an exponential rate, ...

A method for reducing the spurious pressure oscillations observed when simulating moving boundary flow problems with sharp-interface immersed boundary methods (IBMs) is proposed. By first identifying the primary cause of these oscillations to be the violation of the geometric conservation law near the immersed boundary, we adopt a cut-cell based approach to strictly enforce geometric conservati...

A patch representation differing from the traditional treatments in the boundary element method (BEM) is presented, which we call the constant "node patch" method. Its application to solving the Poisson-Boltzmann equation (PBE) demonstrates considerable improvement in speed compared with the constant element and linear element methods. In addition, for the node-based BEMs, we propose an efficie...

In this paper, a new initial value method for solving a class of nonlinear singularly perturbed boundary value problems with a boundary layer at one end is proposed. The method is designed for the practicing engineer or applied mathematician who needs a practical tool for these problems (easy to use, modest problem preparation and ready computer implementation). Using singular perturbation anal...

Nonlinear convective roll cells that develop in thin layers of magnetized ferrofluids heated from above are examined in the limit as the wavenumber of the cells becomes large. Weakly nonlinear solutions of the governing equations are extended to solutions that are valid at larger distances above the curves of marginal stability. In this region, a vortex flow develops where the fundamental vorte...

Boundary impurities are known to dramatically alter certain bulk properties of (1+1)-dimensional strongly correlated systems. The entanglement entropy of a zero temperature Luttinger liquid bisected by a single impurity is computed using a novel finite size scaling or bosonization scheme. For a Luttinger liquid of length 2L and UV cutoff epsilon, the boundary impurity correction (deltaSimp) to ...

The Casimir force due to thermal fluctuations (or pseudo-Casimir force) was recently calculated for the perfect Bose gas in the slab geometry for various boundary conditions. The Casimir pressure due to quantum fluctuations in a weakly-interacting dilute Bose-Einstein condensate (BEC) confined to a parallel plate geometry was also recently calculated for Dirichlet boundary conditions. In this p...

A trigonometric plate theory is assessed for the static bending analysis of plates resting on Winkler elastic foundation. The theory considers the effects of transverse shear and normal strains. The theory accounts for realistic variation of the transverse shear stress through the thickness and satisfies the traction free conditions at the top and bottom surfaces of the plate without using shea...

In this paper, dynamic modeling of a double layer cylindrical functionally graded (FG) microshell is considered. Modeling is based on the first-order shear deformation theory (FSDT), and the equations of motion are derived using the Hamilton's principle. It assumes that functionally graded length scale parameter changes along the thickness. Generalized differential quadrature method (GDQM) is u...