CHAPTER ONE

INTRODUCTION

1.1 Background

The advection-diffusion equation is a fundamental partial differential equation (PDE) used to describe the

transport of a scalar quantity (e.g., concentration, temperature) through a medium due to advection and

diffusion. Advection refers to the bulk movement of the transported quantity with the flow of the fluid, while

diffusion represents the random spreading due to molecular motion (Leveque, 2007). The equation combines

these two processes and serves as the foundation for understanding various transport phenomena in fields

ranging from environmental science to engineering.

The general form of the advection-diffusion equation in one dimension is given by:

where u(x,t) is the transported scalar field (such as concentration or temperature), v is the advection velocity,

D is the diffusion coefficient, and f(x,t) is a source term. This equation models the behavior of pollutants in

water bodies, heat conduction in materials, or the spread of chemicals in industrial processes (Crank, 1975).

Its importance is evident in multiple domains, particularly in environmental modeling, chemical

engineering, and biological transport.

In environmental engineering, for instance, the advection-diffusion equation is used to simulate the transport

of pollutants in rivers and lakes. It plays a critical role in understanding the fate and dispersion of pollutants,

such as heavy metals, chemicals, or nutrients, within natural water systems (Patankar, 1980). Similarly, in

atmospheric sciences, the equation is applied to model the spread of pollutants or particulate matter in the

atmosphere, helping in the prediction of air quality and the evaluation of environmental risks (Yamamoto &

Ragsdale, 2017).

The solution of the advection-diffusion equation is often not straightforward, especially when dealing with

non-linearities, time-varying boundary conditions, or irregular domains. In many cases, analytical solutions

are either non-existent or difficult to obtain due to the complexity of the boundary conditions and the

geometry of the system (Wylie, 2017). This has led to the development of numerous numerical methods

designed to approximate solutions to the advection-diffusion equation.

Numerical methods such as Finite Difference Methods (FDM), Finite Element Methods (FEM), and Spectral

Methods have been extensively employed for the solution of the advection-diffusion equation (Fletcher,

2012). While these methods are highly effective for many standard problems, they often face challenges in

dealing with time-dependent boundary conditions, complex source terms, or irregular geometries. As such,

the development of alternative methods has been an area of active research.

One promising approach is the Laplace transform method, which is used to convert time-dependent

differential equations into algebraic equations in the complex frequency domain. The key advantage of this

method is its ability to handle time-dependent boundary conditions efficiently. By transforming the equation

from the time domain to the frequency (Laplace) domain, the time derivative is converted into an algebraic

term, simplifying the solution process (Abd-Alla, 1996). Once the solution is obtained in the Laplace

domain, the inverse Laplace transform is applied to return to the time domain, providing the solution to the

original PDE.

The use of the Laplace transform for solving PDEs is particularly advantageous when dealing with complex

initial and boundary conditions, as it provides a systematic way to incorporate these conditions into the

solution process. For instance, by transforming both the advection-diffusion equation and the boundary

conditions into the Laplace domain, we can simplify the process of solving the equation and ensure that the

solution respects the specified conditions (Wang & Yip, 2016). This is particularly useful in environmental

problems where boundary conditions often vary with time and space.

Although the Laplace transform simplifies many aspects of the solution process, the main challenge lies in

the numerical inversion of the Laplace transform to recover the solution in the time domain. Techniques

such as Euler’s method and more advanced methods like Padé approximants or Prony’s method are

commonly used for numerical inversion (Cohen & Spector, 2003). While these methods are effective, they

require careful numerical implementation to avoid inaccuracies, especially in systems with highly dynamic

boundary conditions or irregular geometries.

The numerical Laplace transform method has found applications in a variety of fields. In fluid dynamics, it

has been applied to solve the advection-diffusion equation for the dispersion of contaminants in rivers or

lakes (Cohen & Spector, 2021). In heat transfer problems, it is used to model transient heat conduction in

complex materials (Li & Xu, 2016). Additionally, in biological modeling, the Laplace transform approach is
used to simulate the diffusion of substances like nutrients, drugs, or pathogens through tissues, which is

important for understanding physiological processes and drug delivery (Bertalanffy, 1969).

One of the main benefits of the Laplace transform method is its ability to handle time-dependent boundary

conditions. In many real-world applications, the boundary conditions are not static but evolve with time. For

example, in pollution modeling, the rate of pollutant discharge into a water body may change with time,

making the boundary conditions dynamic. The Laplace transform provides a means of incorporating such

time-varying conditions into the solution process without introducing significant complexity (LeVeque,

2007).

Despite its advantages, the method does face limitations. The inversion of the Laplace transform can be

computationally expensive, especially for problems with non-linearities or very long time scales.

Furthermore, the need for appropriate discretization in both the spatial and Laplace domains adds an

additional layer of complexity to the numerical solution. Nonetheless, recent advances in numerical

techniques and the increased computational power available today make the Laplace transform method an

attractive choice for solving the advection-diffusion equation in many practical applications (Wang & Yip,

2016).

In summary, the advection-diffusion equation plays a crucial role in modeling various transport processes

across scientific and engineering domains. While exact solutions are often difficult to obtain, numerical

methods, particularly the Laplace transform method, offer an effective way to approximate solutions. The

ability of the Laplace transform method to handle time-dependent boundary conditions and complex initial

conditions makes it a valuable tool for solving the advection-diffusion equation in real-world scenarios, from

environmental science to industrial engineering.

1.2 Problem Statement

The advection-diffusion equation plays a crucial role in modeling the transport of scalar quantities, such as

heat, mass, or pollutants, across various scientific and engineering domains. This equation governs the

movement of these quantities under the combined effects of advection and diffusion, which are prevalent in

natural and industrial systems. However, obtaining analytical solutions to the advection-diffusion equation is

often difficult due to the complexity of boundary conditions, irregular geometries, and time-dependent

variables. In real-world applications, such as environmental modeling, fluid dynamics, and heat transfer,
these factors often prevent the use of exact solutions, making it necessary to rely on numerical methods. The

challenge is to develop efficient and accurate numerical approaches that can handle these complexities and

provide meaningful results.

Among the various numerical techniques available, the Laplace transform method offers a promising

approach for solving the advection-diffusion equation. By transforming the time-dependent problem into an

algebraic system in the complex frequency domain, the Laplace transform simplifies the solution process,

especially when dealing with time-varying boundary conditions or external forcing terms. The main

difficulty in applying this method lies in the need for numerical inversion of the Laplace transform to

recover the solution in the time domain. This inversion is computationally intensive, and the accuracy of the

resulting solution can be sensitive to the choice of inversion method. While several methods, such as Euler’s

method and Padé approximants, have been proposed for this task, they may not always provide the desired

precision or efficiency for complex advection-diffusion problems.

The problem addressed in this study is to explore and evaluate the numerical Laplace transform method for

solving the advection-diffusion equation in both simple and complex scenarios. Specifically, this work aims

to develop an efficient algorithm for discretizing the advection-diffusion equation in space, applying the

Laplace transform in the time domain, and then numerically inverting the Laplace transform to obtain a

solution in the time domain. The goal is to assess the accuracy, stability, and computational efficiency of the

method, comparing it to other traditional numerical methods, such as finite difference methods, in terms of

both performance and applicability to real-world transport phenomena. By providing a robust and

computationally feasible approach, this study seeks to contribute to the effective numerical modeling of

transport processes in various scientific and engineering fields.

1.3 Significance of the Study

This study holds significant value in the context of numerical modeling of transport phenomena, particularly

in the fields of environmental engineering, fluid dynamics, and heat transfer. The advection-diffusion

equation is pivotal in modeling the spread of pollutants, heat, and other scalar quantities through various

media, and accurate numerical solutions are essential for predicting and managing these processes.

Traditional methods, such as finite difference and finite element methods, often face limitations in handling

complex, time-dependent boundary conditions and irregular geometries. The use of the numerical Laplace
transform method can overcome many of these limitations by efficiently transforming the problem into an

algebraic form in the Laplace domain, where it becomes easier to solve. This study, by exploring the

applicability and accuracy of the Laplace transform approach, aims to offer a more reliable and

computationally efficient alternative to traditional methods, improving the modeling and prediction of

transport phenomena in real-world scenarios.

Furthermore, the results of this study could contribute significantly to the advancement of computational

techniques in scientific and engineering research. By developing a robust and efficient algorithm for solving

the advection-diffusion equation using the Laplace transform method, the research will provide a valuable

tool for predictive modeling in various domains such as environmental management, industrial process

optimization, and biological modeling. For example, in environmental science, the ability to accurately

model pollutant dispersion in rivers or the atmosphere is crucial for assessing risks and formulating effective

mitigation strategies. Similarly, in chemical engineering, optimizing the transport of reactants in industrial

reactors is critical for improving efficiency and safety. By offering a computationally feasible and accurate

solution, this study will support the development of better decision-making tools in these fields and

contribute to a deeper understanding of complex transport processes.

1.4 Justification

The need for efficient and accurate numerical methods to solve the advection-diffusion equation arises from

the complexities inherent in real-world transport processes. In many practical applications, analytical

solutions to the advection-diffusion equation are either difficult or impossible to obtain due to irregular

boundary conditions, complex geometries, and time-dependent factors. Traditional numerical methods, such

as finite difference methods (FDM) and finite element methods (FEM), are widely used to approximate

solutions, but they can be computationally expensive, especially in the case of high-dimensional or non-

linear problems. Furthermore, these methods often struggle with time-varying boundary conditions and

require significant computational resources for accurate results in large-scale simulations. The numerical

Laplace transform method, by transforming the time-dependent advection-diffusion equation into an

algebraic system in the frequency domain, offers a more efficient approach to tackle these challenges,

allowing for easier incorporation of time-dependent conditions and reducing computational complexity.
In addition, the Laplace transform method is particularly advantageous when dealing with time-dependent

boundary conditions—a common feature in many real-world transport problems. For example, in

environmental pollution modeling, pollutant sources may vary over time, or the velocity of fluid flow may

change with external factors such as wind patterns. These dynamic conditions are difficult to model directly

in the time domain but can be more effectively handled in the Laplace domain. By numerically inverting the

Laplace transform, it is possible to recover the solution in the time domain while accounting for these time-

dependent effects. This study justifies its approach by addressing the need for a robust and computationally

efficient numerical tool that can handle a wide range of transport problems, offering more accuracy and

flexibility than traditional methods. The potential for improved modeling of environmental, chemical, and

biological systems underscores the value of this research, as it will contribute to better predictions, decision-

making, and optimization in these critical fields.

1.5 AIMS AND OBJECTIVES

Aim:

The primary aim of solving the advection-diffusion equation using the numerical Laplace method is to find

an efficient and accurate solution to the partial differential equation that models physical phenomena

involving the transport and spread of substances (e.g., heat, pollutants, etc.) in a moving medium.

Objectives:

1. Derive the Laplace-Transformed Equation: Convert the advection-diffusion equation into the Laplace

domain, simplifying the complexity of solving a time-dependent PDE.

2. Solve the Transformed Equation: Develop and implement a numerical method to solve the

transformed advection-diffusion equation. This typically involves discretizing the spatial domain

using methods like finite differences or finite elements.

3. Apply Inversion of the Laplace Transform: Once the solution is obtained in the Laplace domain, use

a numerical Laplace inversion method (e.g., Stehfest’s algorithm) to retrieve the solution in the time

domain

1.6 Research Questions

1. How can the Laplace transform method be applied to the advection-diffusion equation to simplify its

solution process, particularly for time-dependent boundary conditions?

 2. What numerical methods can be developed and implemented to efficiently solve the Laplace-

transformed advection-diffusion equation in the frequency domain?

3. How accurate and computationally efficient is the numerical inversion of the Laplace transform in

recovering the solution to the advection-diffusion equation in the time domain?

4. How does the performance of the Laplace transform method compare to traditional numerical

methods (e.g., finite difference or finite element methods) in solving the advection-diffusion

equation, especially in scenarios with complex boundary conditions or varying transport

coefficients?

1.7 Scope and Limitations of the Study

Scope of the Study

This study focuses on the application of the numerical Laplace transform method to solve the advection-

diffusion equation, which models the transport and dispersion of scalar quantities (such as heat, mass, or

pollutants) in a moving medium. The scope of the research includes:

1. Mathematical Formulation: The study begins with deriving the Laplace-transformed version of the

advection-diffusion equation, transforming the time-dependent PDE into an algebraic equation in the

frequency domain. This approach aims to simplify the solution process by eliminating the time-

dependence and handling the complexity of boundary conditions in a more manageable form.

2. Numerical Solution: The study focuses on developing and implementing numerical methods for

solving the transformed equation. This involves discretizing the spatial domain using techniques

such as finite difference methods (FDM) or finite element methods (FEM), and solving the resulting

system of equations in the Laplace domain. The research will also explore the use of numerical

methods for the inversion of the Laplace transform, such as Stehfest’s algorithm, to recover the

solution in the time domain.

3. Performance Evaluation: The study will assess the accuracy, stability, and computational efficiency

of the Laplace transform method for different types of advection-diffusion problems, including those

with simple and complex boundary conditions. The method will be compared with traditional

numerical methods to evaluate its effectiveness in handling real-world transport scenarios.

 4. Applications: Although the primary focus is on the mathematical and computational aspects, the

study aims to demonstrate the applicability of the Laplace transform method in modeling various

transport phenomena, such as pollutant dispersion in water bodies, heat transfer in materials, and

other engineering and environmental processes.

Limitations of the Study

While the study presents the potential advantages of using the numerical Laplace transform method for

solving the advection-diffusion equation, several limitations must be considered:

1. Computational Complexity of Inversion: One of the primary limitations of the Laplace transform

method lies in the numerical inversion process. While algorithms like Stehfest’s algorithm are

effective, they can still be computationally expensive, particularly for long time periods or highly

dynamic systems. The quality of the inversion process also depends on the choice of numerical

method, which could introduce errors if not properly implemented.

2. Boundary Conditions and Initial Conditions: Although the Laplace transform method excels in

handling time-dependent boundary conditions, its effectiveness may diminish in cases involving non-

linear or highly irregular boundary conditions. The study assumes that the boundary conditions are

sufficiently regular and that the method will be most effective in scenarios where these conditions are

not overly complex or difficult to model.

3. Domain and Geometry Limitations: The study primarily considers 1D or 2D advection-diffusion

problems for simplicity, assuming relatively simple geometries such as rectangular or linear domains.

The method may face challenges when applied to more complex geometries or higher-dimensional

systems, where discretization and numerical inversion could become more computationally

demanding.

4. Linear Assumptions: The study focuses on linear advection-diffusion problems, where the transport

and diffusion coefficients are assumed to be constant or vary in a simple, predictable manner. While

this is appropriate for many practical scenarios, the approach may not be directly applicable to non-

linear advection-diffusion equations without modifications or approximations.

5. Validation and Comparisons: Although comparisons with traditional numerical methods are part of

the study, the research may be limited by the availability of benchmark problems or experimental
data for validation. This could impact the ability to fully assess the practical performance of the

Laplace transform method in real-world applications.

 CHAPTER TWO

LITETRATURE REVIEW

2.1 Introduction

This chapter provides an overview of the key concepts, methods, and previous research related to the

advection-diffusion equation, its solution techniques, and the use of the Laplace transform method in

solving partial differential equations (PDEs). The chapter also explores existing numerical methods,

including the finite difference method (FDM), finite element method (FEM), and the Laplace

transform, to highlight the evolution of numerical techniques used in solving transport equations. It

concludes by identifying the gaps in the literature that this study aims to address.

2.2 Conceptual Framework

The conceptual framework of this study revolves around the application of the Laplace transform

method to solve the advection-diffusion equation, which is central to modeling transport phenomena

in various scientific fields. The framework integrates key principles from mathematical modeling,

numerical methods, and physical processes involved in advection and diffusion. By transforming the

advection-diffusion equation into the Laplace domain, the study aims to simplify the complexities

associated with time-dependent solutions and improve computational efficiency.

At the core of the framework is the advection-diffusion equation, which describes the transport of a

substance (such as heat or pollutants) within a moving medium. This equation incorporates two main

processes: advection, where the substance moves with the velocity of the medium, and diffusion,

where the substance spreads due to random molecular motion. Both of these processes are often

coupled and need to be handled simultaneously, making the equation inherently complex. Traditional

numerical methods, such as finite difference and finite element methods, are often used to solve the

advection-diffusion equation. However, these methods can struggle with handling time-dependent

boundary conditions and may require significant computational resources, especially in high-

dimensional problems.

The Laplace transform method offers a solution to these challenges by converting the time-dependent

terms of the PDE into the frequency domain. The main advantage of this approach is its ability to

simplify time-dependent problems into algebraic equations, which are generally easier to solve. After
solving the transformed equation in the Laplace domain, the solution is converted back into the time

domain using numerical inversion techniques such as Stehfest's algorithm. The inversion process

plays a critical role in retrieving the solution in the time domain, which is essential for capturing the

time evolution of the transported substance.

One of the primary advantages of using the Laplace transform is its ability to handle complex

boundary conditions effectively. For instance, when external conditions change over time—such as

varying velocities or fluctuating external forces—the Laplace transform provides a systematic

approach to incorporate these variations into the solution. This is particularly valuable in

environmental engineering, where pollutant dispersion or heat transfer is often influenced by time-

varying factors. Furthermore, the Laplace transform allows for analytic continuation, which means

that solutions can be extended beyond the domain of known data, providing insights into the

behavior of the system at times or conditions that were not directly simulated.

The conceptual framework also acknowledges the limitations of the Laplace transform approach.

While it simplifies time-dependence, the numerical inversion of the Laplace transform introduces its

own set of challenges, including accuracy and computational efficiency. Numerical algorithms, such

as Euler's method or Padé approximants, are employed to approximate the inverse of the Laplace

transform, but their performance can degrade when the system is highly non-linear or when the

Laplace domain solution is difficult to represent accurately.

To address these challenges, the framework proposes a hybrid approach in which the Laplace

transform method is combined with conventional numerical techniques like finite differences or

finite elements. This hybrid approach aims to take advantage of the strengths of both methods—

using the Laplace transform for simplifying the time-dependence and traditional methods for

handling spatial discretization and boundary conditions. The study will also explore adaptations of

the Laplace transform method for non-linear advection-diffusion equations, where non-linearities in

the transport or diffusion terms can significantly impact the solution.

The validation of the proposed approach is a key component of the conceptual framework. The study

will compare the results obtained from the Laplace transform method with those derived from more

traditional numerical methods. Benchmark problems, such as pollutant dispersion in a river or heat
conduction in a rod, will be used to evaluate the accuracy, stability, and computational efficiency of

the proposed method. These comparisons will provide a robust assessment of the potential benefits

and drawbacks of using the Laplace transform method in practical applications.

Additionally, the framework considers the broader impact of this research on applied sciences. For

example, in environmental engineering, more accurate models of pollutant dispersion can improve

risk assessment and decision-making regarding environmental protection. In chemical engineering,

efficient modeling of reaction rates and transport phenomena can lead to better reactor design and

process optimization. In biological systems, the ability to accurately model the diffusion of chemicals

or drugs through tissues can improve therapeutic strategies.

2.2.1 The Advection-Diffusion Equation

The advection-diffusion equation is a fundamental equation in the modeling of transport phenomena,

such as the diffusion of heat or the advection of pollutants. The general form of the advection-

diffusion equation in one spatial dimension is given by:

In physical terms, the equation describes the change in the quantity u(x,t) over time due to the

combined effects of advection (the transport of the substance due to the movement of the fluid) and

diffusion (the spreading of the substance due to random motion of particles) (Crank, 1975). The

equation is widely used to model phenomena in various fields, such as environmental science

(pollutant dispersion), engineering (heat transfer), and biological systems (chemical diffusion

through tissues).

One of the main challenges in solving the advection-diffusion equation lies in the interaction

between advection and diffusion. Advection tends to move the quantity in a particular direction,

whereas diffusion causes it to spread out. When the advection velocity is much larger than the
 diffusion coefficient, the solution may exhibit sharp gradients, making it difficult to capture the

solution with traditional numerical methods. This issue is often referred to as the advection-dominant

problem (Leveque, 2007). On the other hand, when diffusion dominates, the solution tends to smooth

out, making the behavior easier to capture. These challenges emphasize the need for accurate and

efficient numerical methods for solving the advection-diffusion equation.

2.2.2 Numerical Methods for Solving the Advection-Diffusion Equation

Numerical methods are essential for solving the advection-diffusion equation, as analytical solutions

are often not available for most practical scenarios. Various numerical techniques have been

developed to address the complexities of these equations, each with its own strengths and limitations.

i. Finite Difference Method (FDM) is one of the most commonly used numerical methods for solving

PDEs, including the advection-diffusion equation. FDM approximates derivatives using finite

differences, discretizing the spatial and time domains. For example, in the explicit scheme, the time

derivative is approximated using forward differences, while the spatial derivative is approximated

using central differences (Patankar, 1980). However, FDM suffers from stability issues, especially in

the case of advection-dominated problems. The well-known Courant-Friedrichs-Lewy (CFL)

condition imposes a constraint on the time step to ensure stability. In the case of high advection

velocities, this condition may lead to very little time steps, making the method computationally

expensive and less efficient.

ii. Finite Element Method (FEM) is another popular method for solving PDEs, particularly in cases

with complex geometries and boundary conditions. FEM discretizes the spatial domain into smaller

elements and approximates the solution by constructing piecewise polynomial functions over each

element (Zienkiewicz & Taylor, 2005). While FEM is highly flexible and can handle irregular

geometries more easily than FDM, it is computationally more demanding. Additionally, it may also

suffer from instability in advection-dominated problems if not carefully designed.

iii. Spectral Methods offer a high-order approximation to the solution, achieving high accuracy with

fewer grid points compared to FDM or FEM. These methods approximate the solution as a sum of

basis functions, such as Fourier series or Chebyshev polynomials (Canuto et al., 2016). Spectral
methods are known for their accuracy and efficiency in smooth problems, but they face challenges in

handling sharp gradients or discontinuities, such as those often found in advection-dominated cases.

Despite the availability of these methods, each has limitations in terms of accuracy, stability, and

computational cost. As a result, researchers have sought alternative approaches to improve the

solution process, particularly for time-dependent and complex problems.

2.3 The Laplace Transform Method in Solving PDEs

The Laplace transform is a powerful mathematical tool for transforming time-dependent problems

into simpler algebraic equations. It is widely used in the solution of linear ordinary and partial

differential equations, particularly when dealing with complex boundary conditions or time-varying

external forces (Abd-Alla, 1996). The basic idea is to transform the time-dependent PDE into an

algebraic equation in the complex frequency domain, which is easier to solve. Once the equation is

solved in the Laplace domain, the inverse Laplace transform is applied to retrieve the solution in the

original time domain.

In the context of the advection-diffusion equation, applying the Laplace transform simplifies the time

derivative, turning it into an algebraic term. The resulting equation can then be solved using standard

numerical techniques in the Laplace domain, such as finite difference or finite element methods.

After obtaining the solution in the frequency domain, numerical inversion methods are used to

recover the time-dependent solution. This approach is particularly advantageous when dealing with

time-dependent boundary conditions, which are common in many real-world applications (Wang &

Yip, 2016).

However, the main challenge with the Laplace transform method is the numerical inversion process,

which is necessary to return to the time domain. Direct inversion of the Laplace transform is often

not feasible, and thus numerical inversion algorithms, such as Euler’s method, Padé approximants, or

Stehfest’s algorithm, are used (Cohen & Spector, 2022). These methods aim to approximate the

inverse transform with sufficient accuracy while keeping the computation feasible.

2.4 Previous Research on Numerical Laplace Transform in Advection-Diffusion Problem

The application of numerical methods to solve transport equations, including advection-diffusion

equations, has been widely researched in local academic settings, particularly in fields like
environmental engineering, fluid dynamics, and chemical engineering. Much of the research in these

areas has focused on refining traditional numerical methods, such as finite difference and finite

element methods, and comparing them with newer techniques, including the Laplace transform

method.

One significant study conducted at a university explored the application of the finite difference

method to model pollutant dispersion in rivers. The research aimed to simulate the movement and

spread of pollutants under varying environmental conditions. The study found that while the finite

difference method provided reasonable approximations for simple boundary conditions, it struggled

to handle time-varying sources and complex geometries commonly found in natural systems

(Nguyen & Tan, 2021). This limitation highlighted the need for more sophisticated numerical

methods, such as the Laplace transform, that could better accommodate dynamic and complex

boundary conditions.

In another study, researchers investigated the application of finite element methods to model the

transport of heat in non-homogeneous materials (Ibrahim & Hussein, 2017). While finite element

methods provided an efficient solution for problems involving complex geometries, the study noted

that these methods were computationally expensive, especially when dealing with multi-dimensional

problems. The research emphasized the need for more computationally efficient methods, which

could handle the time-dependence of the heat transfer process without compromising accuracy. This

gap in efficiency led to further exploration of alternative methods, including the Laplace transform

method.

A study on pollutant migration in soil systems found that traditional methods, such as finite

difference and finite element methods, were effective in steady-state simulations but often produced

inaccurate results in time-dependent cases. This limitation became especially pronounced when

external factors such as rainfall or soil permeability were included in the model. To overcome these

challenges, the study proposed combining finite difference methods with Laplace transforms for

time-dependent problems. While this combined method showed promise, the study acknowledged

the challenges of Laplace inversion and proposed further research into numerical algorithms for

more accurate and efficient inversion (Kamal & Ali, 2018).

In the realm of heat diffusion, a study focused on modeling thermal conductivity in composite

materials using the Laplace transform method. The researchers found that the Laplace transform was

able to handle complex boundary conditions, such as time-varying heat sources and insulated

surfaces, with ease. However, the study also pointed out that the method became computationally

expensive when dealing with highly non-linear materials, which often require extensive numerical

inversion techniques. This research highlighted the need for further improvements in Laplace

inversion algorithms and their efficiency in handling non-linear problems.

A significant portion of research has also explored the application of Laplace transforms in fluid

dynamics. For instance, a study on turbulent flow in pipes suggested that while the finite difference

method was effective for solving steady-state flow equations, it struggled with capturing the

transients involved in turbulent flow. The use of the Laplace transform method in conjunction with

finite difference methods provided an effective solution for time-dependent simulations, especially

when the system included time-varying boundary conditions such as fluctuating pressures and

velocities (Fakhruddin & Rahman, 2016).

In biological transport problems, local studies have applied the Laplace transform method to model

the diffusion of drugs through tissues. One such study used the method to simulate the time-

dependent diffusion of a drug in a tissue layer with varying permeability. The research found that

while traditional methods could provide approximate solutions, the Laplace transform method

offered better accuracy in capturing the time-dependent behavior of drug transport, particularly when

the permeability was not constant but varied with time (Sulaiman et al., 2020). This demonstrated the

effectiveness of the Laplace transform in handling time-varying biological processes.

Despite the promising results, most local research has noted the limitations of numerical inversion

methods when applying the Laplace transform to real-world problems. A study by Tan and Khoo

(2019) addressed this issue by proposing a modified inversion technique based on Padé

approximants to improve the accuracy and reduce computational costs. Their results indicated that

the modified approach significantly reduced the errors associated with numerical inversion while

maintaining the method’s computational efficiency.

 In conclusion, different researches have shown that while traditional numerical methods are widely

used in solving the advection-diffusion equation, they often struggle with handling time-varying

boundary conditions, complex geometries, and non-linear problems. The Laplace transform method

offers a promising alternative, but challenges related to numerical inversion and computational

efficiency remain. This study aims to build on the local body of research by improving the numerical

inversion techniques and extending the Laplace transform method to handle more complex and non-

linear advection-diffusion problems efficiently.

2.5 Gaps in the Literature

While the Laplace transform method has been shown to be effective in solving various transport

problems, there are several gaps in the literature that this study seeks to address:

1. Limited application to complex geometries: Most studies have focused on relatively simple

domains. This research will explore the application of the Laplace transform method in more

complex geometries and higher-dimensional problems.

2. Numerical inversion for non-linear cases: The inversion of the Laplace transform for non-linear

advection-diffusion equations is an area that requires further exploration. This study will investigate

the potential for adapting the Laplace transform method for non-linear systems and the challenges

involved.

3. Comparative studies with modern methods: While many studies have explored the Laplace

transform method in isolation, few have conducted comprehensive comparisons with more advanced

methods, such as modern finite element or spectral methods. This research will compare the Laplace

transform method to traditional methods, evaluating both accuracy and computational efficiency.
References

Abd-Alla, M. (1996). Laplace transforms and their applications to differential equations. Dover

Publications.

Bertalanffy, L. (1969). General systems theory: Foundations, development, applications. George

Braziller.

Canuto, C., Hussaini, M. Y., Quarteroni, A., & Zang, T. A. (2016). Spectral methods: Fundamentals in

single domains. Springer.

Cohen, H., & Spector, G. (2021). Numerical methods for Laplace transform inversion. Journal of

Computational Physics, 186(2), 435-457. [Link]

Cohen, R. E., & Spector, L. A. (2022). Numerical inversion of the Laplace transform: A comparison of

algorithms. Journal of Computational Physics, 190(1), 246-258.

Crank, J. (1975). The Mathematics of Diffusion. Oxford University Press.

Fakhruddin, M., & Rahman, M. M. (2016). A numerical study of transient turbulent flow in pipes

using Laplace transforms. Journal of Fluid Dynamics, 50(3), 89-101.

Fletcher, C. (2012). Computational techniques for fluid dynamics (3rd ed.). Springer.

Ibrahim, R., & Hussein, A. (2017). Numerical solutions of heat conduction in non-homogeneous

materials using finite element methods. International Journal of Heat and Mass Transfer,

114, 110-118.

Kamal, M. R., & Ali, R. (2018). The use of finite difference and Laplace transforms for time-

dependent pollutant migration in soils. Environmental Engineering Science, 35(4), 299-311.

Khamis, M., & Ali, M. (2020). Numerical inversion techniques for Laplace transforms in solving

transport equations: A review. Computational and Applied Mathematics, 39(1), 1-13.

LeVeque, R. J. (2007). Finite difference methods for ordinary and partial differential equations:

Steady-state and time-dependent problems. Society for Industrial and Applied Mathematics.

Li, L., & Xu, M. (2016). Laplace transform methods for heat conduction in multilayered structures.

Journal of Heat Transfer, 138(1), 110-118. [Link]

Nguyen, T., & Tan, H. (2021). Modeling pollutant dispersion in rivers using finite difference methods.

Journal of Environmental Engineering, 141(4), 04015012.

Patankar, S. V. (1980). Numerical heat transfer and fluid flow. Hemisphere Publishing Corporation.

Sulaiman, M. S., Yusof, R., & Tan, Y. S. (2020). Modeling drug diffusion through biological tissues: A

comparison of numerical methods. Journal of Biomedical Engineering, 43(5), 529-540.

Tan, S., & Khoo, J. (2019). Modified inversion technique for Laplace transforms in solving complex

transport equations. Applied Mathematical Modelling, 43(7), 2145-2157.

Wang, S., & Yip, R. (2016). The Laplace transform method in transport equations. Environmental

Modeling & Assessment, 11(2), 181-191.

Wylie, E. (2017). Advanced engineering mathematics (8th ed.). McGraw-Hill Education.

Wylie, R. C. (2017). Introduction to mass and heat transfer in biological systems. McGraw-Hill.

Yamamoto, S., & Ragsdale, R. (2017). Air pollution: Its origin and control. Cambridge University

Press.

Zienkiewicz, O. C., & Taylor, R. L. (2005). The Finite Element Method (6th ed.). Elsevier.