# 4 Mathematics for Engineering

# Introduction

“*Omnia de mathematica agitur!*” Indeed, many will say that mathematics is the engineer’s language, so the first thing people need to do before they can study engineering is to learn the essential mathematics. While not all types of engineering problem-solving require mathematics, the structure and rigor that mathematics brings to engineering analyses usually shed much light on the eventual solution to any given situation.

Most people reading this book will already be familiar with most of the basic mathematics needed to progress into most engineering subjects. However, they will still need to review fundamental mathematical concepts to progress successfully. Geometry, algebra, trigonometry, calculus, and vectors give engineers the essential mathematical tools that keep track of processes and can be used to solve problems. Differential equations, for example, also arise in engineering problem-solving. To this end, all engineers must become well-versed in solving various differential equations. In particular, when doing arithmetical and algebraic problems, engineers must do them cold without making any mistakes. *Errores in engineeringo plerumque sunt catastrophes.*

In engineering, quantities with both magnitude and direction are often required, such as forces, velocities, accelerations, and their manipulations involving scalar and vector products. Vector quantities are most conveniently expressed using vector notation, a shorthand version of the corresponding scalar equations. For example, the equations of motion for a fluid are often written as vector equations. Therefore, engineers must become comfortable using vector quantities, including shorthand versions of other vector operators such as the gradient operator, the Laplace operator, and the substantial derivative. *Exercitatio perfectos efficit*.

Learning Objectives

- Remember the essential mathematics that may have already been forgotten.
- Review algebraic manipulations using some simple examples.
- Revise the rules of calculus (differentiation and integration) with examples.
- Comfortably use the concepts associated with vectors.
- Understand the meaning of the primary vector operations such as scalar product, cross-product, gradient, divergence, and curl.
- Know how to interpret line integrals, surface integrals, and volume integrals.
- Appreciate the meaning of the substantial derivative and the Laplace operator.
- Know about the need to use ordinary and partial differential equations in engineering.
- Review the numerical rules of accuracy, significant digits, and rounding.

# Algebra

There is a lot of algebra in engineering. Algebra is a branch of mathematics that deals with symbols and the arithmetic operations of these symbols. Like most things, doing algebra is an acquired skill that takes practice. A review session for students taking their first classes in engineering might start with them writing down several algebraic equations and going about rearranging or simplifying them. To this end, consider the following examples.

Check Your Understanding #1 – Review problems in algebra

- Solve for in the equation .
- Simplify the expression given by .
- Solve the equation given by .
- Simplify the expression .
- Solve for in the expression .

## Show solution/hide solution.

- Distributing the 4, then . Adding 8 to both sides gives . Finally, dividing both sides by 4, then .
- Removing parentheses and combining like terms gives .
- Distributing the 2 gives . Combining like terms gives . Subtracting 6 from both sides gives . Finally, dividing both sides by 3 gives .
- Distributing the coefficients gives . Combining like terms leads to .
- First, notice that the expression is undefined when because it would result in a division by zero. So, is not a valid solution. Factor the numerator to get
Now, each factor can be set equal to zero, so

Solving each equation separately gives the solution to the original expression as and .

# Calculus

Engineering requires calculus, and in the use of calculus, there is a need to know how to do differentiation and integration. Both processes have rules, and success with calculus means first learning the rules, i.e., the mechanical processes usually taught to engineers by the mathematics department. But to succeed with calculus in engineering means much more, in that the process of doing the mathematics must also be physically interpreted. For example, differentiation is akin to finding a rate of change, a slope, or a gradient, and integration is about finding a sum, such as an area, a volume, or a contribution of something that acts or is distributed along a particular path or area. Tangible interpretation of mathematics within the framework of engineering is an acquired skill.

## Ordinary Differentiation

Finding the change in one thing with respect to another is the essence of differentiation. The fundamental rules of differentiation are:

1. Power Rule:

(1)

2. Constant Rule:

(2)

3. Sum/Difference Rule:

(3)

4. Product Rule:

(4)

5. Quotient Rule:

(5)

6. Chain Rule:

(6)

Having laid down these basic rules, which are the first things to learn, they are necessary but insufficient for engineering. Interpreting and understanding the meaning of these processes is more challenging, but it becomes the only way to be successful in engineering problem-solving. Experienced engineers take these “rules” for granted and just do the mathematics* valor facialis*.

Check Your Understanding #2 – Review problems using differentiation

- Find the derivative of the function .
- Find the derivative of the function .
- Find the derivative of the function .
- Find the derivative of the function .
- Find the derivative of the function .

## Show solution/hide solution.

- By taking the derivative term by term, then

2. Applying the chain rule leads to

3. Using the power rule gives

4. Applying the chain rule gives

5. Using the power rule and the chain rule leads to

## Partial Differentiation

Partial differentiation is a mathematical operation used to find the rate of change of a multivariable function with respect to one specific variable while keeping all other variables constant. In the context of a function of two variables, , the partial derivative of with respect to is denoted as , and this represents the rate of change of with respect to while treating as a constant. Similarly, the partial derivative of with respect to is denoted as , which represents the rate of change of with respect to while treating as a constant.

In general, to compute a partial derivative, the function is differentiated with respect to the variable of interest, treating all other variables as constants. Partial differentiation is a fundamental tool in calculus. It is extensively used in engineering to analyze how a function changes with respect to specific variables in complex systems with multiple independent variables or degrees of freedom.

Check Your Understanding #3 – Determining partial derivatives

Find the partial derivatives of the function .

## Show solution/hide solution.

First, differentiate the function with respect to while treating as a constant, giving

Next, differentiate the function with respect to while treating as a constant, giving

# Integration

Integration is a summation often explained as the opposite of differentiation, i.e., anti-differentiation. Integration is more complex than differentiation, but it comes up in some form in most engineering problem-solving. Integration techniques involve using the power rule, trigonometric and logarithmic integrals, and substitution. Most engineers only remember some of the rules of integration. Still, they can look them up as standard integrals in books, or now they can be solved online using AI (artificial intelligence) such as WolfgramAlpha. The ability to successfully use calculus should never be considered a quiz. *Labora callide, non fortiter.*

The fundamental rules of integration are:

1. Power Rule:

(7)

2. Constant Multiple Rule:

(8)

3. Sum/Difference Rule:

(9)

4. Integration by Parts:

(10)

5. Substitution Rule:

(11)

Again, the rules are the first things to learn, but the meaning of the integration processes must also be understood in engineering.

Check Your Understanding #4 – Review problems using integration

- Evaluate the integral .
- Evaluate the integral .
- Evaluate the integral .
- Evaluate the integral .
- Evaluate the integral .

## Show solution/hide solution.

1. Integrating term by term gives

where is the constant of integration.

2. Integrating each term separately leads to

3. Applying the power rule of integration gives

4. Using the substitution , the integral becomes

Integrating then gives

By substituting back , the final result is

and simplifying the expression gives

5. Applying appropriate substitutions and integration techniques gives

# Working with Vectors

A vector is like a direction, force, or something else with a magnitude and direction. Vectors come up everywhere in engineering. *Si vectores nescis, tunc errabis.* Vectors are represented by arrows with a length proportional to the magnitude, and the orientation denotes the direction. Vectors are also used to point to something in space, called a *position vector*. The location of an arbitrary point in space, say , can be defined by specifying the values of three coordinates, i.e., in terms of in a standard Cartesian coordinate system. As shown in the figure below, the point can be located by the position vector, , where

(12)

In general, if is a given vector in a Cartesian coordinate system and , and are the components of in the and directions, then

(13)

as shown in the figure below.

Vectors can be velocities (a speed in a given direction), a force, an acceleration, or something else. Vectors can be added to find a total effect (the total velocity, force, acceleration, etc.), done component-by-component, nose to tail. If is directed from point to point and a second vector is defined as

(14)

that points from to , then the resultant vector that points from to is

(15)

Alternative coordinate systems may describe any problem, e.g., in a cylindrical or spherical coordinate system instead of a Cartesian coordinate system. Usually, deciding on a suitable coordinate system is a matter of convenience, and it is also important to determine whether or not the mathematics can be appropriately simplified in an alternative coordinate system. For example, in introductory engineering problems, Cartesian coordinate systems are used primarily, but sometimes polar coordinates are used when convenient.

# Scalar and Vector Fields

A scalar quantity given as a function of coordinate space (and perhaps time ) is called a scalar field. For example, fluid properties such as pressure , density , and temperature are all scalar quantities, i.e.,

(16)

Similarly, a vector quantity given as a function of coordinate space and time is called a vector field. e.g., velocity is a vector quantity. A velocity can be written in terms of its scalar components, i.e.,

(17)

where the components are

(18)

More often than not, the velocity vector is written as

(19)

where , and are the components in the , and direction, respectively.

Similar expressions can be written for cylindrical and spherical coordinates. In many theoretical aerodynamic and other engineering problems, various scalar and vector fields will be unknowns that must subsequently be obtained in a solution with prescribed initial and other boundary conditions.

# Scalar Products

Let the vectors and be given by

(20)

Then the scalar or “dot” product or “ *dot* ” is given by

(21)

Notice that the “centered dot” or “” represents the scalar product operator.

As shown in the figure below, physically, the scalar product is the projection of one vector onto another, i.e., the component of one vector in the direction of the other. For example, the scalar projection (or scalar component) of a vector in the direction of a vector is given by

(22)

where is the angle between the vectors and .

Notice that because is the projection of onto , then is equal to the projection of onto times the magnitude of . Similarly, because is the projection of onto , then also equals the projection of onto times the magnitude of . However, the easiest way to remember what the dot product is just .

Need some help with scalar products? Here is a short video lesson on how to perform scalar products and an interpretation of the physical meaning of this vector operation from Dr. Leishman’s “Math & Physics Hints and Tips” series.

# Vector Products

The vector product (cross product) is given by

(23)

where the vector product or “cross product” is denoted by the operator, is the angle between the two vectors, and is a unit vector perpendicular to the plane containing the two vectors, as shown in the figure below. Notice that the cross-product is also the area of the shaded parallelogram. Like the dot or scalar product, the concept of a cross-product has many uses and applications in engineering.

In general, the vector product can also be written as

(24)

Expanding this operation out in terms of the third-order determinant gives

(25)

Similar expressions for and can be written in cylindrical and spherical coordinate systems.

Need some help with vector products? Here is a short video lesson on how to perform vector products as well as the interpretation of the physical meaning of this vector operation from Dr. Leishman’s “Math & Physics Hints and Tips” series.

Check Your Understanding #5 – Review problems with vector operations

Given two vectors and , then evaluate:

## Show solution/hide solution.

# Differentiating a Vector

To differentiate a vector, the vector must depend on time, i.e., its magnitude changes. Therefore, the vector’s derivative will be its time rate of change. For example, let the vector be , a position in space. The time derivative of will be a velocity, i.e.,

(26)

Therefore, the derivative of a vector is a vector.

In other possible engineering examples, if were the momentum, then the time rate of change of momentum would be a force. If the parameter were to be a velocity, then the time rate of change of velocity is acceleration. If the parameter were to be work, then the time derivative would be power. Therefore, differentiating a vector can lead to valuable outcomes for engineering applications.

Check Your Understanding #6 – Differentiating a vector

Differentiate the vector function

## Show solution/hide solution.

To differentiate this vector function, it is necessary to take the derivative of each component with respect to , i.e.,

By simplifying the derivatives of each component, then

This example demonstrates how to differentiate a vector function component-wise; the result is another vector.

# Line Integrals

Consider a vector field . Also, consider a curve in space connecting two points and , as shown in the figure below. Let be an elemental length of the curve.

The line integral along the curve from to is

(27)

which is simply a statement that the integral is the sum of the components of the vector field along the direction of the curve to .

If the curve is closed, i.e., points and are coincident, then the line integral is

(28)

where the counterclockwise direction is considered positive according to the convention used in the field of mathematics, as shown in the figure below.

Check Your Understanding #7 – Line integrals

Evaluate the integral along the curve defined by the straight line segment from point (0, 0) to (2, 4).

## Show solution/hide solution.

To calculate the line integral, the curve can be parameterized by

where varies from 0 to 1, corresponding to the endpoints of the line segment. The line integral is given by

Substituting the parameterizations and simplifying them gives

And then evaluating the integral gives

# Surface Integrals

Consider now a surface bounded by curve , as shown in the figure below. Let be an elemental surface area and let be a unit normal vector, i.e., one perpendicular to the surface. Let be a scalar field in this space and a vector field. A surface integral over the surface can be defined as

(29)

where is called the elemental vector area or the *unit normal vector area*. If the surface is closed, then the integral is written as

(30)

Line integrals can also be based on scalar functions, e.g., the scalar integrated along a curve given by

where represents an infinitesimal length along the curve .

# Unit Normal Vector

The meaning of a *unit normal vector*, also known as a normal unit vector or simply a *normal vector*, should be understood because it comes up in many engineering problems. The unit normal is a vector perpendicular (or orthogonal) to a surface at a particular point. It has a magnitude of 1, hence the term “unit.” The normal vector points always away from the surface in a perpendicular direction, as shown in the figure below.

The unit normal vector is often denoted by the symbol . For a surface defined by a function , where are the coordinates of a point on the surface, then the unit normal vector at that point can be obtained by taking the gradient of the function and normalizing it, as given by

(31)

where denotes the Euclidean *L*^{2} norm or magnitude of the vector. By normalizing the gradient vector, the resulting unit normal vector will have a magnitude of 1. This normalization is achieved by dividing the gradient vector by its magnitude.

Different methods exist to find the unit normal from discrete surface points, such as the cross-product approach or the least-squares fitting procedure. The cross-product approach is the easiest. The steps are:

- Select three neighboring points on the surface defined by the vectors , and .
- Calculate the vectors between these points, i.e., and .
- Compute the cross product of these vectors giving .
- Normalize the resulting vector to obtain the unit normal vector, i.e., .

Remember that the unit normal vector is a crucial concept in vector calculus and is used in various mathematical and physical applications, such as determining the direction of a force on a surface and solving problems involving surface integrals or differential equations.

Check Your Understanding #8 – Calculating unit normal vectors

- In a three-dimensional Cartesian coordinate system, what is the unit normal vector for the direction of the x-axis?
- Consider a two-dimensional plane defined by . What is the unit normal to this plane?
- Consider a three-dimensional plane defined by . What is the unit normal to this plane?

## Show solution/hide solution.

- This vector is parallel to the -axis and points in the positive -direction, so .
- In this case, . Therefore, at any point on this plane, the unit normal vector will be
- In this case, . Therefore, at any point on this plane, the unit normal vector will be

# Volume Integrals

Consider a volume in space, as shown in the figure below, which contains a small elemental volume .

Let be a scalar field in this space. The volume integral over the volume of the quantity is written as

(32)

The above result is a scalar. If is a vector field in space, then the volume integral is written as

(33)

and the result in this latter case will be a vector.

Need to brush up on your understanding of integration? Here is a short video lesson on the processes of integration from Dr. Leishman’s “Math & Physics Hints and Tips” series.

# Integral Relations

Some important relations between line, surface, and volume integrals occasionally arise. Consider an area bounded by a closed curve . Let be a vector field. The line integral of over is related to the surface integral of over by using Stokes’ Theorem, i.e.,

(34)

Again, consider the volume enclosed by the closed surface . Then, the surface and volume integrals of the vector field are related by using Gauss’s divergence theorem, i.e.,

(35)

which connects a volume integral to a surface integral. Finally, if represents a scalar field, then the gradient theorem states that

(36)

# Gradient of a Scalar Field

If , then the gradient (or grad ) is given by

(37)

where the operator (which is called nabla or del) is defined as

(38)

Physically, is the normal to the surface defined by = constant, as illustrated in the figure below.

Related to the gradient is the concept of a directional derivative. Consider some scalar quantity defined at some point over a surface. Choose some arbitrary direction away from the point and let be a unit vector in that direction. The rate of change of per unit length in the direction is

(39)

The term is called the *directional derivative*. Therefore, it can be appreciated that the rate of change of in any arbitrary direction is simply the component of in that particular direction.

Check Your Understanding #9 – Gradient of a scalar field

A two-dimensional pressure field is given by

If an isobar is a constant pressure line, plot the isobars for = 0, 2, 4, 6, and 8.

## Show solution/hide solution.

The equations of the isobars are

which are all straight lines, as shown in the figure below.

The gradient is given by

so that

and the magnitude of the pressure gradient is

The pressure gradient vector represents the magnitude and direction of the pressure change at a specific point in space. It points toward the steepest change in pressure, and its magnitude indicates the rate of pressure change in that direction. In this case, the pressure gradient vector is given by , which means that the pressure changes most rapidly in the direction , and the rate of decrease is 2 units in the direction and 1 unit in the direction.

To visualize the pressure gradient vectors, they are typically plotted as arrows; the tail of each arrow is placed at a specific point in space, and the direction and length of the arrows indicate the direction and magnitude of the pressure gradient at that point. Some arbitrary points can be selected, and then the corresponding pressure gradient vectors are plotted, e.g.,

Notice that the gradient vector runs perpendicular to the isobars, and at each of the points, A, B, and C, the pressure gradient vector is the same.

# Divergence of a Vector Field

Consider a vector field , then

(40)

which is a scalar quantity.

In aerodynamics and fluid dynamics, the divergence can be interpreted as a flux density or the amount of flux entering or leaving a point in space. Mass fluxes such as , , etc. are often used. Divergence is the rate of flux expansion (positive divergence) or flux contraction (negative divergence). By convention, the flux is defined as positive when it leaves a closed surface. So, the divergence is just the net flux per unit volume or flux density of the contents of a given region of space.

# Curl of a Vector Field

Again, consider the vector field defined by

(41)

Then the curl of is defined as

(42)

and expanding this determinant gives

(43)

A physical interpretation of the curl is that it represents the angular velocity of the rotation or “spin” of the contents of a given region of space. In aerodynamics, the curl of a velocity field is called *vorticity*, a physical concept that can help interpret the behavior of a flow. If the vorticity is zero, then the flow is said to be *irrotational*, and if the vorticity is non-zero, then the flow is *rotational*.

Check Your Understanding #10 – Finding the curl of a vector field

In a solution to a specific engineering problem, the solution of the vector field is

Find the curl of this field.

## Show solution/hide solution.

The curl of a vector requires the evaluation of a determinant, i.e.,

Evaluating this determinant gives

After differentiation then

or

# Physical Laws in Engineering

Taking a temporary break from the mathematics review allows one to recall some of the physical laws used in engineering. Aerospace engineers must be well-versed in understanding the laws of physics and how these are applied in the form of principles and the corresponding equations. In the following chapters of this ebook, the laws and principles of physics will be well used to explain fluid dynamics and other physical behavior. If you have yet to read or listen to the Feynman Lectures on physics, you should do so.

## Conservation Principles

Conservation principles in physics refer to fundamental laws that describe the preservation or constancy of specific quantities in physical systems. These laws state that specific physical quantities remain constant in isolated systems or undergo specific transformations under certain conditions.

1. Conservation of Energy (*Conservatio Energiae*):

“Energy cannot be created or destroyed; it can only be converted from one form to another or transferred between different objects or systems.” This principle states that the total energy of an isolated system remains constant over time, i.e.,

where denotes energy.

2. Conservation of Mass (*Conservatio Massae*):

“Mass is neither created nor destroyed in a closed system; it remains constant.” This principle states the total mass of a closed system remains unchanged in any physical or chemical process, i.e.,

where denotes mass.

3. Conservation of Momentum (*Conservatio Momentum*):

“The total momentum of an isolated system remains constant if no external forces act upon it.” Momentum, defined as the product of an object’s mass and velocity, is conserved in a closed system where no external forces are present, i.e.,

where denotes linear momentum.

4. Conservation of Angular Momentum (*Conservatio Momenti Angulares*):

“The total angular momentum of an isolated system remains constant without external torques.” Angular momentum, which is related to an object’s rotational motion, is conserved when no external torques are applied to the system, i.e.,

where denotes angular momentum.

## NEWTON’S LAWS

Isaac Newton formulated three fundamental laws of motion in his *Philosophiæ Naturalis Principia Mathematica*. These laws laid the foundation for explaining classical mechanics and have been essential in understanding various fields, including engineering.

1. Newton’s First Law of Motion (Law of Inertia). *Lex prima motus Newtonii (Lex inertiae)*:

“An object at rest tends to stay at rest, and an object in motion tends to stay in motion with the same speed and in the same direction unless acted upon by an external force.” Inertia is the property of an object that resists changes in its motion. According to the first law, an object will maintain its state of rest or uniform motion unless a net external force acts on it.

2. Newton’s Second Law of Motion (Law of Acceleration). *Lex secunda motus Newtonii (Lex accelerationis)*:

“The acceleration of an object is directly proportional to the net force acting on it and inversely proportional to its mass.” This law can be mathematically expressed as

where is a force, is a mass, and is an acceleration, or more specifically

where is the momentum of “something.”

3. Newton’s Third Law of Motion (Law of Action-Reaction). *Lex tertia motus Newtonii (Lex actionis et reactionis)*:

“For every action, there is an equal and opposite reaction.” According to the third law, whenever an object exerts a force on another object, the second object exerts an equal and opposite force on the first object. Forces always occur in pairs, and they act on different objects.

# Laplace Operator

The Laplace operator appears in several of the equations of aerodynamics and is a second-order differential operator defined as the divergence (i.e., ) of the gradient (i.e., ). The Laplace operator can be applied to both scalar and vector fields.

In the case of a scalar field, say , the Laplace operator would be written in Cartesian coordinates as

(44)

Of course, the Laplace operator appears in the familiar Laplace equation, i.e.,

(45)

In aerodynamics, the Laplace equation is, in fact, the governing equation for an incompressible, irrotational flow. To this end, the linearity of the Laplace equation is helpful in that if and are both flow solutions to the Laplace equation, then is also a solution, i.e., the idea that incompressible, irrotational flows can be combined using the principle of *linear superposition*, which is a significant advantage when dealing with complex flows.

The Laplace operator can also be applied to a vector field, i.e., . In this case, then

(46)

where

(47)

Check Your Understanding #11 – Finding a Laplacian

A two-dimensional function , where and are the Cartesian coordinates. Determine the Laplacian of this function.

## Show solution/hide solution.

1. To find the Laplacian, the second partial derivatives must be obtained, followed by taking their sum. Calculate the first partial derivatives:

2. Calculate the second partial derivatives:

3. Sum the second partial derivatives to obtain the Laplacian:

So, in this example, the Laplacian of the function is 4.

# Substantial Derivative Operator

The substantial derivative operator in Cartesian coordinates is written as

(48)

or because the velocity components are given by

(49)

then

(50)

Recall that vector gradient operator is defined as

(51)

Hence, the substantial derivative operator can be written as

(52)

The substantial derivative can be applied to any flow field variable, scalar, or vector. For example, if is the component of the velocity field then can be written as

(53)

The term is called the *local derivative*, which can be interpreted physically as the time rate of change of a given quantity at a fixed point. The term is called the c*onvective derivative*, which can be interpreted as the time rate of change from the movement of the given quantity from one point to another where its properties are spatially different. It will be apparent that terms like have units of time.

Notice that for a vector then is

(54)

which, in terms of the scalar components, is

(55)

Check Your Understanding #12 – Calculating a substantial derivative

A fluid flows through a pipe, and the temperature distribution is given by . What is the change in temperature at the point where the fluid velocity is ?

## Show solution/hide solution.

The temporal change and the convective terms involving the gradients and velocity components must be obtained to calculate the substantial temperature derivative at that point. The substantial derivative is

The contributing terms are

(56)

Therefore,

So, the substantial temperature derivative at the point in this example is 40.

Notice that the substantial derivative incorporates both the local temporal change in temperature and the convective transport from fluid motion. In this case, the local temporal change in temperature alone contributes a rate of 3 units per unit of time. The remaining terms involving convective transport, which consider the temperature gradients and the fluid velocity, contribute to a total rate of 37 units per unit of time.

# Ordinary Differential Equations (ODEs)

An ordinary differential equation (ODE) is an equation that involves ordinary derivatives of a function. A common engineering problem is solving an ODE, i.e., determining what function (or functions) will satisfy a given differential equation or a set of such equations. For example, given the ODE

(57)

then what will be the function ?

The antiderivative of is , so

(58)

where is some unknown constant. Therefore, solving an ODE is more complicated than just anti-differentiation because the value of the constant must also be determined. To this end, additional information, such as an initial or boundary condition, is required. Commonly, the value of will be known at some specific value of , e.g., . In this case, it will be apparent that

(59)

for all values of .

Check Your Understanding #13 – Solving an ODE by the method of separation of variables

Use the method of separation of variables to solve the ODE given by

to determine the function .

## Show solution/hide solution.

The separation of variables is a method frequently encountered in engineering practice. In this case, the separation of variables gives

Integrating both sides gives

If this ODE has an initial condition, say , then must satisfy

So, the constant of integration must be

Therefore, the particular solution to this ODE is

# Partial Differential Equations (PDEs)

Partial differential equations (PDEs) deal with equations that have multiple independent variables and their partial derivatives. PDEs model various engineering phenomena, including heat conduction, fluid dynamics, structural dynamics, elasticity, aeroelasticity, and acoustics.

PDEs can be classified into several types based on their order and linearity. The highest order of the partial derivatives involved determines the order of a PDE. For example, a second-order PDE involves second-order partial derivatives. Linearity refers to whether the equation is linear or nonlinear regarding the unknown function and its derivatives.

Examples of PDEs include:

Navier-Stokes Equations (incompressible, convective form):

Solving PDEs can be challenging, and different techniques are employed depending on the type of equation and its properties. PDEs are used in incompressible flow theory to help find velocity potentials and stream functions. Analytical methods involve finding exact solutions using techniques such as the separation of variables, Fourier series, or Laplace transforms.

However, analytical solutions are not always possible, especially for nonlinear PDEs or complex boundary conditions. To this end, numerical methods are often used to approximate solutions for PDEs. These methods involve discretizing the domain and approximating the derivatives using finite difference, finite element, or finite volume techniques.

Check Your Understanding #14 – PDEs in incompressible flow theory

In a specific incompressible potential flow, the velocity field in the plane is given by the pair of partial differential equations

Find the velocity potential .

## Show solution/hide solution.

Integrating the first equation with respect to gives

For the second equation, integrating with respect to gives

So, there are two equations

By inspection then and , so the velocity potential is

which are the equations of a set of parabolas. The final result can be verified by partial differentiation of .

# Working with Numbers

In engineering problems, information is often presented as numbers that must be manipulated to find a final answer. Much engineering work is done with numbers, and it is frequently said that anything else can only be an opinion, i.e., the truth can only be expressed in numbers or data. In this process, questions may be raised about whether it is necessary to round the numbers to a certain number of decimal places and/or significant digits. Common sense rules generally prevail, but some guidelines should be followed.

## Rounding Decimals

Sometimes, a number must be rounded to a certain number of decimal places or even a whole number. To this end, some numbers must be removed from the decimals at the end of the number. The method is to look at the value of the following number after the one to stop at, then apply the rules for rounding:

- Round up if the first digit to be discarded exceeds 5. e.g., 3.14159 is 3.1416 to 4 decimal places.
- Round down if the first digit to be discarded is below 5. e.g., 3.141592 is 3.14159 to 5 decimal places.
- If the first discarded digit is 5, round up if a non-zero digit follows it, but round down if it is followed by a zero. For example, 3.14159 is 3.142 to 3 decimal places.

In multi-step calculations, rounding up or down of numbers is generally permitted only at the end of the calculation; the appropriate number of decimal places (or otherwise the accepted accuracy of the known quantities) must be carried forth and then finally round off the number(s), as needed. However, rounding off intermediate numbers may be avoided altogether in many cases. Instead, the known numbers with an acceptable level of accuracy to a certain number of decimal places are used. However, using numbers to 8 decimal places is usually unnecessary.

So, remember that if rounding off numbers is needed, which is often the case, this is usually done only at the end of the calculation. To this end, common-sense rules with rounding numbers generally prevail. Rounding 9.81 to 9.8 may be acceptable (1 decimal place or two significant digits), but rounding to 10.0 is unacceptable. Rounding 0.002378 to 0.00238 may be fine, but rounding to 0.02 is generally unacceptable. Rounding 3.14156 to 3.142 may be acceptable, but rounding to 3.0 is unacceptable.

## Significant Digits

Significant digits are the total number used to express a measured or calculated quantity. Like decimal places, significant digits measure how “precise” a number is given. The value with the fewest significant digits sets the standard for the degree of precision.

Generally, a goal is to carry through numerical calculations containing values that have (more or less) the same number of significant digits to arrive at an answer with the needed (or maximum possible) accuracy. There are also some simple rules for significant digits:

- All non-zero digits in a number are always considered significant.
- Any zeros appearing in a number between two non-zero digits are considered significant.
- A final zero or trailing zeros in the decimal portion are considered significant.

Check Your Understanding #15 – Significant digits

- How many significant digits are in 14.638?
- How many significant digits are in 0.002378?
- How many significant digits are in 273.70?

## Show solution/hide solution.

- Answer = 5; 2. Answer = 4; 3. Answer = 5

## Manipulating Numbers

There are a few basic rules that should be followed when manipulating numbers:

- Addition and subtraction. Generally, the answer should have the same number of decimal places as the term with the fewest decimal places, e.g., 4.123 + 2.56 = 6.68.
- Multiplication and division. Generally, the answer should have the same number of significant digits as the term with the fewest number of significant digits, e.g., 4.121 x 6.62 = 27.3.
- Integers are always treated as if they have infinite significant digits; e.g., 5 is considered for operations as 5.000…

Check Your Understanding #16 – Rounding & accuracy

## Show solution/hide solution.

A ruler can be read easily to about 0.1 of a millimeter, so if the disk’s diameter were reported as 50.1 mm, that would be a suitable measurement. However, if the diameter is reported as 50.13 mm, this value implies that you can read the ruler to 0.01 mm, which is perhaps more unlikely. Therefore, the area is = 1971.3568, but if this area value is reported, it implies that the diameter was known to 8 significant figures. So, the area using the ruler must only be reported as 1971 mm. If the diameter could be measured to 2 decimal places, its area could be reported as 1971.4 mm.

Finally, common-sense rules should also be applied when working with numbers in engineering calculations. For example, it would be inappropriate to round up “g” (acceleration under gravity) from 9.81 m/s to 10.0 when simultaneously using the value of to 4 decimal places in the same calculation. Sometimes, when doing homework, it is necessary to put down the calculator and pencil, apply common sense, and say: “Based on what I just did, does this numerical answer seem right?” Ultimately, getting the “right” numerical answers in problem-solving means understanding the basic concepts of significant digits and rounding and learning how to spot likely numerical errors. Spotting errors requires experience, and experience is gained by practice.

# Summary & Closure

The ability of engineers to bring to bear the necessary mathematical concepts is an essential aspect of problem-solving. Geometry, algebra, trigonometry, calculus, and differential equations give engineers the crucial mathematical tools to help engineers solve complex problems, quantify processes, and make informed decisions. For example, many engineering problems involve vector quantities (e.g., forces, velocities, accelerations, etc.), which may require various vector operations such as scalar and vector products, line and surface integrals, volume integrals, and gradients. In addition, calculus brings other mathematical tools to the table that are often needed for engineering problem-solving. For example, besides pure mathematics, engineers nearly always need numerical results, which must be evaluated accurately and consistently, and a consistent number of significant digits or decimal places based on the available data.

5-Question Self-Assessment Quickquiz

For Further Thought or Discussion

- Explain the fundamentals of a scalar quantity. Make a list of some familiar scalar quantities used in engineering.
- Explain the fundamentals of a vector quantity. Make a list of some examples of vector quantities used in engineering.
- Explain the physical meaning of the curl of a vector field. Is the curl of a uniform flow positive, negative, or zero?
- What happens when the simple rules of significant digits are not used in basic arithmetic calculations?
- What happens when numbers are prematurely rounded up or down, and there is still a need to carry forward the needed number of significant digits in complex arithmetic calculations?
- If the function A(x,y) = K(x^2 + y^2) where K is a constant, then what is A(r, theta)? Your answer may depend on whether you are biased toward pure mathematics, physics, or engineering.

Other Useful Online Resources

There is more to explore with these mathematical concepts:

- This video goes into more depth on vector fields, divergence, and curl.
- If you want to learn more about the Laplace equation, this website has a lot of information on that and other mathematical operations.
- Here are some more short video lessons on differentiation and partial differentiation from Dr. Leishman’s “Math & Physics Hints and Tips” series:

- Here is a short, fun video on the history of bookbag calculators from Dr. Leishman’s “Math & Physics Hints and Tips” series.